{ "cells": [ { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import matplotlib.pyplot as plt\n", "import seaborn as sns" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Check out this page for various distributions to sample from:\n", "https://docs.scipy.org/doc/numpy-1.13.0/reference/routines.random.html" ] }, { "cell_type": "code", "execution_count": 11, "metadata": {}, "outputs": [], "source": [ "length=500\n", "randu=np.random.random_sample(length)\n", "records=np.zeros(length)\n", "records[0]=randu[0]\n", "k=1\n", "while (k]" ] }, "execution_count": 13, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXQAAAD4CAYAAAD8Zh1EAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nO19f5Bd1X3f57u7T9KTsLXCyIl5IJAdLGFFRRvWmETt2CI2IqaCHeEgMJkmHU+ZtHFqVLIZMfEgQd1BicaVk6mbhqbUtc0YQaAbYeisXUsZZ2jksMquDIu1icJPrewiGxbX0oLe7p7+8d5d3b3v/Lz3nHvPvXs+M6B999137/n5Pd/fX2KMISAgICCg/OgqugEBAQEBAXYQCHpAQEBARRAIekBAQEBFEAh6QEBAQEUQCHpAQEBARdBT1Isvuugidvnllxf1+oCAgIBS4ujRoz9mjK3mfVcYQb/88ssxMjJS1OsDAgICSgkiekX0XVC5BAQEBFQEgaAHBAQEVASBoAcEBARUBIGgBwQEBFQESoJORA8R0etE9LzgeyKiPyGiE0T0fSL6JfvNDAgICAhQQcfL5SsA/hOArwq+/zUAV7T/+wiAP23/G1AhDI1OYt/wBE5NTePi3joGt67DQF+j6GYFBATEoCTojLHvEtHlkltuBvBV1krbeISIeonofYyxH1pqY0DBGBqdxD1PPIfp5iwAYHJqGvc88RwAWCfq4eAIAMI6SAsbfugNAK/FPp9sX+sg6ER0J4A7AWDNmjUWXh2QB/YNT2C6OYvl56bxW0efxNKZcwCA1w4Bxzf8HNb//LutvOf4j36KUz94HbfOnU/pfOopwvEr32vtHQH+Y1Gsg23bgA9/2PpjbRB04lzjJllnjD0I4EEA6O/vt5KIPZzk7nFqahoAsO0H38Xvf7eleZuLpv3/AIz4i8AUHwTwQd6qeKa1oBja/2u/z8Y7A/yDbB1UBhdf7C1BPwng0tjnSwCcsvBcJfJUBSxmXNxbx+TUNK557XmcXt6LD3/2awCdJ6eN3jqe2XVd5vd8YNdTfE4AQL3WPT/P0ecHtm8M85wD8maaROuAALy090Zn760CbBD0gwA+S0SPoGUMfSsv/fm+4Qlc8OaPceTPfxsr3zlz/osv5PH2xYM4Y/TUus0LiDlwnoPPiujgSKKbaAExB4Dp5iz2DU8Egp4CJgS6CKZJtA4u7q2nfuZikeSVBJ2IvgHgYwAuIqKTAHYDqAEAY+y/AHgawCcBnABwFsC/dNXYJE5NTWPj/zuNle+cwV9e+VG8vOri+e8+9/Er8mrGosDxH/4Uwy+8joNX/rOO72xttJX1GmrdhObsef4syZnHYesgWUwwJdCR/SQO14fp4NZ1C9oItNbB4NZ1qZ63mCR5HS+X2xXfMwC/Y61FBri4t46ek61JevwXr8N33381gJYK4HMWVAA+oihOYz2A46OTOPXEc4CjjTY13USti7BqeQ1TZ5vz/ds3PMHl2FbWawueVTUOzEWfTAm06NB0eZhG7bDV9yIOpaJQWLZFGxjcug5PvNDK2NjsbnUlC4GxBVfEpWhOw8ZGi49NFxFmE0XKm3MMy5f0YPTe6xdcH3zsGJpzC+89c24GQ6OTAGBlXHw6FFzNNe9glF13of7QwUBfw9rYF3EoFYVSE/SBvgZW/8oa4GvAbFc3Gh5wZi6Jrg+cRpaNlhybJDGPkNxoA30N3PfkON4821xwvTnLsG94AgAyj4vuvNkk+rJn2Zzr+HtE6Ca+z5Bt9YcKLg7VPA4lX5iBUhN0ANh82UoAwKO/+1Hg2msLbo1bolt2ToM3NjzwNtpUgphHkPXdZFx05s3mYa16lq25Tr5HBNHhalv9YdJWW8yQ60OpaMk5jtITdDTbG71Wk99nCaqT2CXRLUr8tQWdMRBtNFXfs46LzrzZPKxVz9Kda9V61D1EG5Kxsqn+kMEVM2RyKKXhtH2QnCOUP9tijgQ9Ooknp6bBcP4kjvS4gJiI2CC6g1vXoV7rXnDNB5uBLkRj0E0EQouoiHzLZX23MS4682bzsFY9S6dPOusxyyGaN1wyQwN9DTyz6zq8tPdGPLPrOiExV40n7zci+0MRknMg6AaQncQRXBLdgb4GHti+EY3eupIAyjA0OonNew9h7a6nsHnvIemCtQnR2Hzx1qukGw2Q993GuOjMm83DWvUsnT6J1uNdB8bm5zXLISqDizWkO76u1q/O/k62Y/CxY8LnFSE5B5WLAUQnbvyEdi3eZRV/i9T3ZdXHyvqedVx02mZTF6vzLFWfZBxgNK+3XN3A40cnrUbZulhDQ6OTOPPOTMd1kVTiYv2aSgj7hic6PK8iEFCI1BMIugFEek1Ca6FFC0qHuJguTFtW9KL1fXnpY9NA1bYsBxJv/h7YvjHTnIrWY4Tp5iwOHz+d+T1J6K4h0ZpNXt+yfnXHoQMAq5bXsHvbBi2pxMb6NbVRyQ5UhmKClgJBN8Dg1nXYeWCsI88EA4wXlMnC/PzQc3j4yKvz783ClZTdU8YlZIdm8rv9OzZldol8YPvGTDlweFx+Eqempq0fojprSNTnkVfeWEC8J6emF6ztOJYv6elot8v1ayqByQ5UmZHZJYIO3QADfQ1h8ijTBaW7MIdGJ7kLXqbbk0GlpyxKv140ZAaxNMay+HPvfvSYkW5WF3E9uwgu9Lg6um4Rw/KN773WcV21p+JrskvgL2+jn6a2mMGt61Dr6mxPFwFnz80UsocCh26IRgrXQR7npyve7RuesHaIAHIuxCd/2ryhMoilEfOj8dQNoEqDiPvm+Zu78l7R4WRFfRONBQ8X99a1gtFs9tNEmonu23NwHFPTLTq0vNaF5hybD4LLew9VlkN3xWmaerGIuLst61drPUe26dNwJTIuxNTKXyXIJKa0Yr7KB9wm92zLA8rWu2TeNTwkr0Z7QTSGWb10bGGgr4Gx3dfj5b034uW9N2LViqULkssB+e6h0nPoL7zyE3wIwAf3fBur3/PueYLoitM0NYyJiCTPWLVl/WrsG57AzgNjSk4+ixVdxIUsZv26i8Al2bjxDu+shu88Dc6qd4m4eJHXzS1XN3D4+OmOvu88MMZ9/hxjXuZGL3oPlZqgD41O4uRzk/gQgHNdPZicmuYaLQG7nhwmG0c2wfHniNQdvA1AAO64do31zVv2SNQsUKkR0qgzZPndk1xl1dRdMsan/7ILtQ+usuVhKXoPlZqg7xuewO0zTcxQ13zRBZmGrghOUzTBK+s1bN57aH4RnT03o83Jy7wvsizGvBMx+QQdyct0nEXjyVMRFO1O6gIixseEISpbHpai91CpCfqpqWn0zM5gpluvG0VwmrwJrnURzpybmTekyHyJZW5nthdj1sCfMkN1MKZRZ5iMZ9Giuq/gjSFPNZl2jdo+SIveQ6Um6Bf31lGbm0Wzq1t5rwudpQ54E3z23ExHKlgRZIeQyWLU7avPgT+u4FLdoTueOqK6Lyla00A3yIjXJx3VZHSfKVwcpEXuoVIT9MGt63Dm8TnMdMm7wcuTnqfOMjnBa3c9pfU7lahm4steBf2sK4Lmg7pDJaqXeQ5Ngox2HhjDXQfGhLUNbM9V0Tpv2yi12+JAXwO/subdmOsRE/RVy2vcpE9FuuiJFktvvWbkdqabzMh2X4sIPsoS3KOCD+oOlStgmV1K0wQZiebX9lyVPYNpEqXm0AFg7colQO8KEPgGUdPCCHlsYhE3tuemDUZchq4BxmZfi+IUfcrh4QoyUd2HQ0cFkQSVNsiIN7+25ipZnHxZrWtBHVvfpR4RSk/Q0WwCS5YYT3SahcFbsIC5AcSW4UT3OTYJVlHqCZ9yeBQBXw4dEWQHvcx9U0XUk/NrY654xcnrtW6j/Dy+ohoEvVYznmjT+3kLdvCxYwBhPjIs0gGOvPIGvjCwUdpsW4YTnefYJFhFcYouCVrRnglJ8BgHXw4dERcuO+hNgoySSM6v6Vzx2uuDzcQVKkPQTSfa9P77nhzvWAS8XMgMwMNHXkX/ZRd6szhsEqyiOEXXBM0X7x5ZZsa84hFM2waoA+gAcZDR3Y8e43Lqomho3bkStVd0gPikvkoLYgbJcmyiv7+fjYyMZH/QTTcBJ08Cf/d32Z8lwNDoJO4ShCCL0OitZ0qN6itESaDyyKdRZrc9XWzee4h7YMrWU15zImsbwI+n0NkHa3c9JQwIfDlDeL+ovSJVT1n2LBEdZYz1876rDIduChPikMaToAqnPQ9Fqid84aJdIo1KKy8Vgqxt+3dsSi1BiaS+rDnFZcbYeq27cPWVCyxKgm7qqZGGOPtirHKBxUBYi0IalZZLu0ac8ekScLYX99YzHfSu1GmygyLSpVdN2islQY8vsidePI1L3rMCqw1+b8rRyCqT1GtdmG7OJa65P+1dqx8Wg3rDR6Qhbq7sGqa5yNMe9K6kPtlYJt8ZSeFlX+OlI+jxRUZsDuydd/D3P+nGM7GaniqYcjSyUl/TzTnUuggXLOvJzY/VZkk6HvL0NQ8Hx0KkIW6uOFxZLvI5xqzOlwupTzaWZY68laF0BD1aZP/mbx7F73/3qwCAb//CNUb6QlOOJr4weL9rzjEsX9KD0Xuv1+1GaqhK0pksxjQuaLalABubqmqHgilxc8XhihgcX3OR8yAay6q6LpaOoEeL7Iofv4o3l70L/73/JnznAx820hem4WiihSGyyKfVV5oSI1sl6dK6oNmEjU2V5lDQHfMyHRTR+ozavPPA2LwPeNo26zI+ZRqnCGWIvE2D0hH0aJH1zM3iJ8tX4k823w7AzCKehaMx4e5VCz0NMbJVkk5GTPPyNbexqUwPBd0xL6NIXkRu7zKOE1C+whm60ErORUQ3ENEEEZ0gol2c79cQ0WEiGiWi7xPRJ+03tYUomU733Cxm2mlz0+gLB/oaeGbXdXhp743c5F2q98chSs2rSiaVJuGSaMGZlqSTEVNeHwmtPthMxqWbXEwG00NBd8zLmAzLdptVCcNU7ywiiRsPvHaYJuUy7YvLZHIyKAk6EXUD+DKAXwPwIQC3E9GHErd9HsCjjLE+ALcB+M+2GxohWmQX9ACzXd25F4nVWeSA3uZKw6GKiK1pSToZMY33MXq+Kgtems1rI9Od6aGgO+ZlFMld5faWMT6iZ0frJG+CloSIsALQLqqdhjgXxRDoqFyuAXCCMfYiABDRIwBuBvBC7B4G4N3tv1cCOGWzkUkM9DWAtauA5ecKiezSMVrpbK40Yp8tAxhPnCYAW9avnn/PQF+DG22XVGmkFbtt9EWkFtiyfvWCEn/Rc3XH3PdkWDwU0WZZ4i0eQbv70WNWKg3xYJq3RVcyT2PrKYoh0CHoDQCvxT6fBPCRxD17AHyLiH4XwAoAH+c9iIjuBHAnAKxZs8a0rQvRbAKSPOhFQ2dzpXU3s+HiNdDXwMgrbyzwmGEAHj86uSAPjWxhRhuI109d42bWvvAOhS3rV3cUTogOGN0x9yUZlgmKaLPonaJ8KZEve16utjbytqQhzkUxBDo6dOJcSzpa3A7gK4yxSwB8EsDXiKjj2YyxBxlj/Yyx/tWrTUKBOJiZ8Zqg66gTdNU3rnD4+Gmh+2ME0QJcWa/Ni6Ei5KWeSKoFDh8/LeWodMY8y9ykUT/Z0DfnvZ6iA326OYvudpH26J06Tgo2VRAiLjpqVxImhDWNraeowhk6FPEkgEtjny9Bp0rlMwBuAADG2N8Q0TIAFwF43UYjufCcoOuqE4oMo9fhPLasX42vH3m1457m7Jw07SlQnHpC1S/dMU8zN2ndKG15imRZTyZ1PwF0RJEmozBlHHIEW4e+y7wtad2cgfxzHulQxGcBXEFEawFMomX0/HTinlcB/CqArxDRlQCWAThts6EdmJkBVqxw+oqs8D3niY5YePg4fxrPnJNv1CLVE7J+uXYlS6Nv9SHIxaTu5z1PPIdltS5pm5METZYHxgZc521Z2nO+v120ULoQPauI/a8k6IyxGSL6LIBhAN0AHmKMjRPR/QBGGGMHAdwN4L8S0U601DG/xVzn5Z2ZwY/OzuAWjuErQA86nEcaDkpU4DcvyAylrn2m0+hbffCokdX9TBLi6easlm46TtBEKX5tHfqqvC1ZpJbkc6MyCD763GvpLBhjTwN4OnHt3tjfLwDYbLdpckz9dBovzNL8qezj4PoOHbFQxPmIarj21muF55RO9mtlvQYicFVHtjlh0Xh1EWHtrqeMxjhPlVXaup9J6KTPcMGAiZ4PgOvtpAtRPpsIvqUL8FcJrcCPp36Gc+++YME13wa3DFBxL6ZeDG9NN623MQ3iofAqXa5NTliUyC3p3THyyhs4fPz0/IFT66b5UoZA/ior07qfvfUa3pmZM9Yru9ybyefbsE3orA2fYhO0IkV9BGvOYIa6O677NLhVgMhzQuTF4JuftorDAuy2OTlePC+L6eYsHj7y6nygytR0E2DAquW1QrydALFXxu0fuZR7fc9NGwr10NKBbnCPzMNIZ22I0n4UESVbWg59KeYw29VJ0H0jKD4hrUFQxFnZ0om6NFSqDngXnHB8vNbueop7T5LnzTNjJw8ylUj/ZRcK58cnAp6Ejm1CxcXLUmcD8rQfReS3KS1Bf8+y7g63Rd8DP4qE7UVmSyfqevHLipPkYbyVvT+JoqVL0cHtu7eWCDq2CZWHkcgeI6t9UKTXUmkJ+oouhk1rL0Kjtx68XDTgYpHZ2OiuF7/IBpCXemBw6zoMPnYMzTm1cTFNxs4AMbJ4cYm8dXRQpNdSaQk6ZmZw2c+vLMyjwuVGc/FsH1zjTN5vq11FBXjE33/fk+N486zcWOyb6F4FZPHiyqK6LdJrqbwEvcBcLi43mqtnuyxWkOUAymPxu1QZ6PR9SkLMCchVdC+a48/7/UmingwGcpH/RpT4Lko/7bLPpfVyKTL0P2tqTJkF3FXaTZ3cEmnShGbN+1xUzgsb0O276HBq9Nal+fhtSy9F5egu8v2qd7rIfxN/JqCXftoWAkFPAVkO6KyJ712pIHQWrs5hkjyM7ntyPNMBlHdCKZvQPXzTHlqig6CLKBVBKLpoh+j9ew6O5/7OeJ8H+tIVu5Ehemajt65MgGcT5VW5FEjQZZ4LKvWISox2qYJQqR5UhwlPHWT6rDTt8hW6h29aPb4sSCmNGk7XjS9vF9Kp6SaGRidTvUfV3qJtR3m/v7IcukvHfh7HFUF1+qom2KYKwnQMVGlCdYJ0VM8yhS9lzHgwSauahguMpBdRcJIpl6dqr2uViGxNpOFYddpro8xhFuT9/nISdMakBN10YZoSjWijiZCmkPPKeg2b9x7CzgNjWNrTlTlqMM3mVB0mulyFLR14HjrXLAdGHvr/gb4G5gT5VEy5PFV7s6hkdMZRNi5pOFad9hZto8n7/eVUuczNtf6t1bhfiyb6rgNj2HNwfEFggKy6jYyIytzRRER7aHQSZ96Z6bhe6yKcOTfTCgFHSwSt17qxf8em1OJuGg8JlWpApA7qrdewYmmPdTHdtY96Vo+ivFwibanhVAmsRCo0FbHVHcc0e0YGXR9yoFi31TzfX06C3mwvCAGHLluAU7HkUZNT0wtKsEWQEQ1Z2TUAqHUT9/QVJYlatbx1KCUXeVbClVZ3J9Nni1y89ty0wckCda1/tHFg5KH/37J+dcc6TcvlqRJY8RAntqZ1O5Njs3vbBmtugkXmvY9D9a48bUTlJOgzbS5XQNBNwq1F8Xs8oqGz+Fcs6eFOnkj/vHxJjxPC5cK4mje34dpHvWiDmQ6GRifx+NHJBeuUANxytR0iobKLxImtjbqdNteQad77eIZLW2vXt+CvShJ0VUIdHTCgIwhAxygoSh8rIx4uCJergsFpuY00HFPaPui+y4c85Crw1hzDwkpSWbhR2eGVzHUjq9tpUo1ItYZ0+yM6HETtjEs5tgivD9Wm4qgkQY9PtA6nLirWkJx0Hc5NZtUWEQ8XxLdo3WEcWbiYeOmvVctr2L1Nrt4xeVfWcc9DrE/jSmpCqGSl25JpNVzW7Yxg2h/e4bDzwBj32Saq1agtqvn1Tcorp5eLgqAD593EvrRjE2rd/MrfQGvh3XHtGmF+77jVXMW5yRaxzNrtKrjGRcBEGqTxnog2dtzm8XZzLvW77n70WIfnRZZxzyvqMY0rqYlLo4kXhiziVXccVd4wNoKfTCQsEeHNGgVclJRXOQ49eapuWb9aqChPipRrdz3FvTXuI87L0cA4z0pCxTGXNbhGB2m4mLSirIyLFHlepBn3vERtlRSRlUM0keSy1u3U4b5tcLyyfZqEiPDqzq8r1WZalJOgC7xceAuG58UCdIqUQ6OTysrkWdUYVSbaMqTRVafd2DKDOG9DplWb5CVqp3UlNeEQdddl1vWvQyRt9SfZzqR7MiAnvK6jgF2hnAQ94tATfugiAxIPvKolPGKenPTFSpSzIA0Xk3ZjqwziJtVqZMjToJrGlXRw6zonOv4s61+HSNrieHntlFVeSsJkfn2iCZXSoZtwR6qqJUCrHmQUEepr+HkZkEZXnTbCThYuD+hXq1GB1754itS81ohobAEUmlmRB51kY67sSdGzdW1KLiI880hjUW4OPUHQRadqUn+mW7VkljGMvPJGqkjSgIUw5WKyiLLRPTaq1ei0b3JqmpsiNX6fS/DGdvPeQ1650wH6ycZ84Hhtq1Ly8lcvFUGPRMjeiefxFICdjz+HoWeXzQ+2SFy75eqGNKBApnf9+pFXO64VvTGyoOgCB6btSts2nQ2ZVW0StY8XNu9yjdhypyuq6tbdjx7rUG/6uKdsHix5GdFLQ9DjJ9x75loDM9XEAnHyge0b8cD2jdbSlMrgUzShLmRcAlCcYcdVu1Qb0pa+Nk9fZF1OT3VYfX7oOSuBNjzCDUDaRpGfeBn3lC7yWiOlIejRCdd463X0nWrpOGfpvAkgOu3S+FtH998lWGg8+BRNqAtZgYF3ZuYKUysV1S5bYnWeBlIb7nRDo5PGOYx4EB0uy2pd3DbedWCsJWEvr1lL0FUW5LVGSmMUjU6yG4//Ne77338GAHhr2QXce9JgoK8hDC5Koiwl0pKQFRjIGsyRBpGRSKTuyqNdNoKv8kyRauJOJzIu7hueMMphJILocJEVxJ6cmsbP3p7pCPYr657SRV5rpDQcenTCHbzyoxi9eB3OLqlj/L3v77gnC3RUL6oAIp9hkrQMcCsC6yQ6E8E30TxPX2Qb7nRp8vXzkHYemnPMWcplX5HXGikNQY+I7Y/efRF+9O6LOr63lXgKAPYcHF8Qcg60cpZfsKy1ACMO0eZk5GGsFInhy2pdXK5qZZ2fb94GdLL8idolyzfPG8M8xjYvzwwben+ZN5iN5/TWawtUZTy8Nd3E2O7rtd9VBeSxRkpD0JMn3Mp6bUGhClubNBr0OBFYWa/hzLmZeeJiW5ebl7FSxCUAwOBjx9CcWyiInzk3k7rWowo6Wf4AtethBNEYltntVHYQZVkTotD4O65dk/k5UX78qI0iibDK+vIiQUxQ3so1+vv72cjISCHvNoVIz8vLSGfz+TxOp17rthZoEUff/d/icsO2+piE7pjqctei54lSu7rqly3wVFKyuTeVQnTut3GPaT8C1CCio4yxft53Whw6Ed0A4I8BdAP4c8bYXs49twLYg5Yn4THG2KdTtzhH5J0ik5c8TGYUTMKVv+6UwJDlSl+tqzrQFVNlwWEm9/sCmecPj6iaBq2oxtWkrJzsOb7lOqk6lASdiLoBfBnAJwCcBPAsER1kjL0Qu+cKAPcA2MwYe5OI3uuqwTZhy6c3y/t4gUsquCBGeRd7sL3RRe03Lb6QN0QMhcwjKakGU7kyprEh2AyE8SHyswgUEcSn47Z4DYATjLEXGWPnADwC4ObEPf8KwJcZY28CAGPsdbvNdAPdXB62XI50Kh7Fnx/VG03CBTHK0/Uugg2XwQii9t/+kUsLrfougyzntmyOk+tTJkGmzdvuW+GGsiGvfPlJ6BD0BoDXYp9Ptq/F8UEAHySiZ4joSFtF0wEiupOIRoho5PTp07xbcoUNn14b7+Phge0bsXvbhtyIkcukSHlA1P4vDGz0tl8yhkI2x8l1JCuykDYBWdrCDXkkoCoDbBTqSAMdHTovbV1Shu0BcAWAjwG4BMBfE9EvMsamFvyIsQcBPAi0jKLGrbWMvFNk6vqBN3rrHSJ1HmJb2UVjUft97ZeMoRjoa+C+J8e13DZl9oi0YfZp3CN9K5gsgg+lA11Bh6CfBHBp7PMlAE5x7jnCGGsCeImIJtAi8M9aaaUl8AySJknvs0IncImAVpWlNnwlRkXB1+RiaaBiKHZv26BtOAb4B7/IdVDFaaexcfhWMJmHvA6dvG1SEXRULs8CuIKI1hLREgC3ATiYuGcIwBYAIKKL0FLBvGizoVnB02k9fnQSt1zdsCaOq8RNnlpg8wcuXCACMQCPH51ctKKqDLw5HHzsGPru/1YhIn5W9YLKbmGiBhPZI3xMS1Ak8lKFFJUvX8mhM8ZmiOizAIbRclt8iDE2TkT3AxhhjB1sf3c9Eb0AYBbAIGPsJ05anBL3PTnOncjDx09b8UdWZa9Lcpb7d2yaT72aNUnSYgFvMzbnmLOALxlscHo6XLBMQtORVtJ6E6XpX1FcqQmyHjq6EmJ83KMxidOGwceOLbjPFhZFYNHQ6KQwkyIBeGnvjZmfv/PAmLB2qUgf+cD2jcLfZWmXa7WE7Pku3y0q4p2ETtBQ1na6DjZTwXXATpr+8doUpcywHdGdFlnmLe2Yb7rvW9yYkt56LVX6g8yBRWWHTJyywT2ostfJxDzbXI1rHaEqTYHLd+salVXclo0xKlq94FpfreqfTloC1ykz0iBLLpy0Y84j5rLrWVCa9LlZINtkNnSJqux1ss1hW8dpQ0co0w3Lnu9aP8kbKx5Uh6GNdqZ167MF1weKrH8yH+u4Ln/F0h40Z/mViUxhyx1S1y7Be1/Rh7gOFgVBFy3O3nrNGufIQ5S9TrY5bPt/29ARygIiZM/XfXfazZkcq956LVVebRsbs4hArDhcHyiy/ukeiLYIoO0gHVVAm+h9vSkD/UQBgqLrWbAoCLpocUZZ4bIi7mYYx6984EIM9DW0vBlsRUxm3eiqzSp7vs67s27O+CbSbxAAACAASURBVFiN7b4e+z51lfFhqNtOU4+lPAOWXB8osv7pEmpbh07eQTqi9zGGVGO+e9uGDsaj1k3Yvc0O/YljUejQXScIOnycH/X68k+mc3l/HFnzZYs2a6S7Vj1f9W7but80fvqqPthKTJWETYNxHmtK1D9du48oTa+IARIhb1WH6LlvTTexf8cm4zHPc/8vCoIOuA3QkRHBzXsPZa5gb4Ksi0dW/CCeFEr2fNl3PughVX1wYXB0YawuKujMJFPmyCtvLHDnjeIs+i+7MPOadGWvkL0v7ZjnNVeLhqC7hMz7Im/LflYuMAoXT3rtMGCeoMkWp2rh+uKrLGun7UNnaHQSdz96rCProy/xBkOjkwuqdK1aXsPubRuk7Vrac74QtOz+w8dPZ46zsFGlyQR5v88mFoUO3TVU3hd5JOUB7BiPBvoaVgoIi9p35p2ZjutZN4vthFA2DY7RnPial31odBKDjx1b4EL35tkmBv/iGHcco/7E73+7OSd8vioTpM685W2vKNo+kgWBQ7cAXlRYEpH6xaUOTUdVoMPBNxxw0clI2gg63KAMLlQZNjk0Vcpk3TF1FbC1b3iio/QgADRnGZeLNlVHiSSylfWa0bzZVFnoRtiWgYAnEQi6JUQLQBSJFuVxANypYXSCQUR1Nw8fP+0sadnQ6CSXmAPA8iU9mcbAhb7bphFLxoHrjqnLYDFZ+3jfmaqjRIcjEbTmzeQg0y2ZZ2MsfU0SF1QuliFKyiPSI9qESlUgIn5fP/Kq06RlqkjaLHBlZLXlSiqak24i7TEVzdvdj/LVIjbaJ/rOVB0lUl/olDw0USHq3msr8K6I4hU6CATdMngL2BUxS0Llm6z7vnjSMhu+8apI2iRMdOJFR2yqIJqTL956VeY0A7OMcQmJyfgNbl2HWldnyYNaN3GlhzT+77zDUWfeTIhvnsFORRWv0EEg6A6QXMANQ6JjK5IyyVmbEDmbh40qkjYOU+6nqDSlurBhYJPNW5KQmI7fQF8D+379KvTWz0ctLq914YKlPdh5YKxjDG0ZDHUOBhPim2ewkw+utyIEHXoOMDGyZdXxyYw5OgU2IpgscJU+URRgcse1azIb3ZIG6bh6y4dkUNG7s7xfNW9xQpLGphC1L+6+eLbtucIbQxsGQx07hYmLa5ZgJ1P7kC+utzwEDr0Nl7UQTbga0Ya8i8MtpW1HN/GqCp6HyQLX4Qh5/d+/YxO+MLCx43mqIC3e/EQSEU+95YsonAWqeYsTkrTcI88dMYKrMVTZKUzUO7r32pAwis7jI0Pg0JFPWSpdrka28WxFFwKdIfoRZ9vgcEoy6HKEuv2XRaqqvIR8FoWzQjRvSUKSlntUuVcWMYYm3kam92bZ1za9oGwjEHT4VQtRlfPbRrvycMtLSwBE6hmdaEMXueV92rQ687Zl/Wp8/cirHb9V5U9RzdfKuv3MgDowIb55+o776qdeGYKeZfPZJEo2Qu9Vem4b3JKtBWlKRFXjwyNaogMuOQ42A4J8rWCfHJ9IFRJdFyWKE12PoGIkzpybWZDLJ8BPVEKHntUv1Jbrm63Q+0jHJ0IXUeHeGxFM9Im645PWS8hmyLavrmlZ8tXLoEpfEUWOBviNSnDoWVUmtjg7W6qb6N7Bx45xw7Ij/+P4vUVBpgZIcuNnz82kGh+T+bEleaiMs0WpYVRrLK3aKektxEMVbBFVRyUIelaViS2dsk3VjSjHRgRfMvUBfCLKU1mIoBofF0Yoleoni3HWJVRrLAtzokpfUYRbnm92DBl8aGslCLoNY5gNzs6mUU7nEPCZY1J5TcShMz42jVA6+vEsxtn4e2xvcNUas3H4+ZI+1lc7Bg++tLUSOnRXfqGmvuk226FD5HwIZBBB97ApglDo+PpnTeHgKt+HzhpT+XerYNMWkQW+2jF48KWtleDQXYnkpieuzXaovF18CWQQQcRJ9tZrWLG0B6emprGyXgMRsPPAGPYNT+Qmopr4+sfbI1JFREbq+L2uXGHz8oH2wS2vTHEFvrSVmCDxvmv09/ezkZGRQt6tA9HmbfTW8cyu63JpQ1xkj4jf1Nmm97pEoPNABFqHUMTpqb531SaZ0S8O3jzz2hyhXuvGLVc35tMQi3YVAXhp740pWs+HD3pbV23IsgfzHpc86QURHWWM9fO+qwSH7gI+nLg+cElpoeIk73tyPNdgLhkx5oE3z1G7ROXkRDnf47CpJvNBb+uyDWl1+UWMiy92h0DQBfA5AU9ZIDqQhkYn8eZZdT5smzAx0gLy/N47D4xxv1MRc9sbPO8IZx7X67INKqZAxIUXEfntSzqARUXQTcQwX05cE/ggfutAZihydWCqKgeZzLMqqjIJav/G9nzkKUWKuF6dDJBZIGMKRFx4UdK1DxL1oiHopmKYLyeuLlyIma4OCNnGcnVgiohwI8Zp6vZT16Uxen5ch6o7pjr35SlFirjebiJuAWzXkqyMC89rXHxkoBYNQVe5FfEmxocTVxe2xUyXekiZB4xNX/P4nMrqpJrOM++w16nDqjumRdV+lUFWNclUwnHZnlNT09i/Y5Nz6doH+wUPWn7oRHQDEU0Q0Qki2iW571NExIiIa4EtErJQbl/rA5rAtpjp0q9W5Eu956YNmZ8N8H3AbddJTfp6f2Fgo9J3W3dMRfc97Lj2qwwi7jZ6Z94+66L2dLVzxrtuky9+50koOXQi6gbwZQCfAHASwLNEdJAx9kLivncB+LcAvueioVkh03v6kjpXBZmI17u8xjU0phUzXeohbamzTI1iUZ1UV1Bx+rpjKrqPF6Xquk8RZDalIiRZUZzGLGPYeWBsPrf//h2bnLTNBy84HnRULtcAOMEYexEAiOgRADcDeCFx378H8EcAfs9qCy3BpPwaUPzEJCET8QDgZ2/PdPxGVOg3/kwRUXWth8xKBHw0iqkgGlOGlh9zNP4mRte8+uSbTUnmQppHCUJfveB0VC4NAK/FPp9sX5sHEfUBuJQx9k2LbbOKKJxZVX4tQtETk4RMxBMl8lqxpEe4kFWh6T6X2QLURjEeXM2pbooIWZGJ+PiLCl/zYFr7NUuZxawpBWxjoK+BOUVgpCs1iK/7Q4eg89bS/CgSUReA/QDuVj6I6E4iGiGikdOn5Qn3XUBnAQB+TEwSMq5T9N1bnPqQEVQ6QF/yeYggG488N5tJzhZVkYm4qi859ndcuyZTn1zllikaOgeaCynG1/2ho3I5CeDS2OdLAJyKfX4XgF8E8FfU4n5/HsBBIrqJMbYgtp8x9iCAB4FW6H+GdqeGSFTqJsIcY1xR0gf3JJWIZyr+6aglfPbykY1HnuoBE+8ikwyavLHvv+zC1H3yqcyiTeioUl1JZj7uDx2C/iyAK4hoLYBJALcB+HT0JWPsLQAXRZ+J6K8A/F6SmPsCkXFHdLr64p40uHUdBv/iGJqz58/BuI7c1E3LVx2gLlSBX7Y3m+hQN9HX6+jGZeOfpU++2hWyIn54T05Nd8QD+Chtu4SSoDPGZojoswCGAXQDeIgxNk5E9wMYYYwddN3ItJBx1rqcjlecTVKmaX9OLupuog71SRJljISNI08uXHaomxyMRWbQLPsBLgNv/c8yNh805pu07RJagUWMsacBPJ24dq/g3o9lb1Z2qDhr3UnMm7ORueIlDZ/NOTZ/sET9MZEmlvZ0zd+7ankNu7dtKNXizkvklR3qpuXxouflnUGz7Ae4DMm9HgU78Yi5D9K2S1Q2UtQWZ50nZ5PVFU/U57sfPQbg/KLlZR58uzln3NaqcDqqvsjG3lRSyOsQ4vXpge0bKzNncejuda+kbUeoLEG3xVnnydlkzU8hC8+OcyJZF3aVOB2dvuiUffOp36I+3XK1P220iawBW2W3I8RRiRJ0PNjyRc7TPSmrK56sb3GdetaF7UPYc1af6gg6ffHV51gE3dQBvrotiuZWdF13r+cdn1AEKsuh2+Ss8+LAsrriqYxuEcHOqkYqmtOxKSHoum8C/kRJqmCSOsA3dYMsMVk8EVl8znX3epXtCBEqS9DLtgmB7K54snBo4DzBzrqwTQ4EF7p2m7pQ3b64PNRtj5GPqQN0IZrbb3zvNW6VqH3DE/O5bFRjWEaaYIrKEnTAP92mCjYWHM/bBeg8GLK8RyQJnD03s6BYsitdu4hYmRSdiJDmcLNJgF2MkUm+dt/UDTI7kOx+3b1eNppgikoT9LLBFqHQIdhZFnb0uz0HxzEVSy/w5tmmVeOrCKKiCrp5euIwPdxsE2AXY8TrU56507NAFsldRCENUxTt/RUIuiewTShccyIRwZ5K5IuJEyNXunYRtzbLGDbvPWS8mUzGyjYBdjVGtlMH5AWRxHTL1Q3vDyQfvL8CQfcEZfSRVREjVz78DcFzCefVLq42k20CnGecQxnUDTKJKY8DKQuH7cMeDgTdMXQXSNGeI2mgIkam+mndsdLVEbvYTLYJsE3Pi/j45RmFahuig8f1gZSVw9bZw65VMqUk6EXrqXRhskB8z7XBG3MdrxxATz9tMla854oMorYPRNuub7Y8L5LjF1eFuZJWyrIPI6jam5XDVu3hPFQyxDTyg7tAf38/GxkxT8jIC1uXZUssEpv3HuJOcLISPOB3v2RtA+y4gZmMlYvfm0CHkOVN7ET9j0O07tK00+f1yoNOe9fueorrCUQAXtp7Y+Z32FqjRHSUMcat21w6Dt0HPZUuTNQoPvvIysbcVuUaWRHveHk2EfIMGlGJ/kUYx0xyrUfI0s4y7UNAr70iDruLCGt3PdWxJ03z5eShVi0dQS+TrtlUjeKr0SqPMZepTXQIjU8HYhHELk2u9Szt9GEfmkgXOu2VFZ4GOuv48g7DB7ZvFHLbeahVS5fLpUz5GHzJAZI174numGd5D2+s4tDJFTPQ16p5uX/HJgDAzgNjmfK8pEURxE41frx1l6WdRe9D05J6Ou1N5m3ixTVE6zBNPqM86EHpCLqogK6sAG9RSC6QIuoO2qglqbMQs74nPlYi6BAaH2pnFkHskmutt17DquU16brL0s6imRVTgqrb3ogpeGnvjcL6w7I6vrI1mgc9KJ3KZaCvgZFX3sDDR16dN2AwAI8fnUT/ZRd6p7IoWo1iQ/zXUWfYeo/MeKRDaHzQ7RaVBMp0rWVpZ9EqLlOCmqa9KhVJmjXqmh6UjqADrerpPmSOK4PbVhax2raOUhdZCI0Put2iiZ0usrYzTpyitbLzwFgu/U2jj7Z94PmYubGUBN2HTetDmK8O0hpiTPtn0+CThdDYaIeNg7poySyCqi822lnEXshDCtKVTH06tEtJ0H0IwvFBtNdB2oVv2j8XATdpxjFrO8pyUOsgr74UtRfyqIkrW4e+HNpxlM4oChRvkAH8kBJ0kNYQk0ZHWbQB2EY7fKjGZAt59SXvvRAdVPFoWNOauFVFKTl0H3SUPkgJukjDSeSho3SFLO0oy0Gtg7zsJ3nvhbJIx0WglAQdkG/aPIyVVS9nlXf/fDEwl+mgViEv+0nea6VKh65tlFLlIkMWP2STwBhfVAyukGf/fPAdj+CDOs8W0vbFVFWT915Q+c/bKiBeRpQuOZcKqgQ4Ik6wbMmGdOEL5ytDnom1dFCGMdOFaV+GRidx14Ex7ne6SapcQ5Usror7OI5KJedSQSaOyUTJKurlyuKx4ZsInUYH7+shYNKXaL2I4IvaSWZD27z3UOX2sQkqR9BlekMZ0faNqOhCRkiKPKRk7Up+17u8hjfPNjueERehfSCWutLd5NQ07jowhvueHE/lSldUf3nrJYKp2sl1H0QHVVn3sS1UTocu0xvKJrvoZENpoNI9q1LSutItytrF++5nb8+g1r0wEVI0Z77o12XtEBHCqGi2SVuL7K+M6JmoLIrsQxn3sU1UjqDLDDSyyS6jMUxlvJIt4iSRtWlEkrWL911zjmHFkh7unPniF55GuovfY+M9riFaL43euhF3XWQfyriPbaJ0KhcdUU4kjsncq1z4trsWO1Xi5eDWdRh87Biac3zD93RzFnsOjuOdmTmrevY0Yu9b002M7b7eyrNcQCXdyXKRq/KU677HNWy5HxbZB9v7OOsezlt9ViqCntXIp5psm4ExWduqsxBUfsYDfQ3c9+Q4Vz8dIR5tFyGrnt1mljpf/MJl7RAVRohjaHQyt1w0aWGLGBY9Z7b2sY09nLdTgpbKhYhuIKIJIjpBRLs43/87InqBiL5PRN8hosvsN9WOKDfQdz7fsa3yaTxkaatKBxmpSCanppFMwZ/kqKYkxFyGLNyUTOw1FYltidBZ1UqydkRqPk49hHnE513WFt3+uvK1trE/qqL2yEpvilA9KTl0IuoG8GUAnwBwEsCzRHSQMfZC7LZRAP2MsbNE9K8B/BGAHbYb64v4rYMsbVUthPipz9DyD2Zo6TqTHJVMHVCvdWNZrUvqYZIGSU5vZb0GIsynVr3l6gYOHz+txQXa4BptcEo60t1Ogf82cH7eVW3R6a9Nzs+FSsAWp5+lgLWNPmWlN0XQKx2VyzUATjDGXgQAInoEwM0A5gk6Y+xw7P4jAH7DZiMjFC3KmSBLW2ULgUfsI2LOC8IRqQOi7HSAm7zOEXHiEZ/Hj04aeU1kFaFtuW+q2iE7PKN512mL6j0m/VG5j7pSCejMmYu2mbqQytqQld4UQa90VC4NAK/FPp9sXxPhMwD+F+8LIrqTiEaIaOT06dP6rWxDtxSaD2G/sraq2ijTJ9vIgvilHZsweu/185vOZdi2D14qJmPGmxvdNTW4dR1qXZ16l1o3za9RG1yb7jNUqrsi58ZV20xcSFVtyKo6KkL1pMOh8zSDXLcJIvoNAP0APsr7njH2IIAHgVbov2Yb5yES5QAs0ClHD7bBcaQV32RtVXEeMm+DfcMT1rMg2jQGJ+GDmkyXU+Jxd4OPHQMIaM52Vn7neVcBwJ6D4/PG5mSebhtcm+4zRETx7kePASh2blRSRtq26biQ6gbemaqOeLTige0bvfNyOQng0tjnSwCcSt5ERB8H8AcAPsoYe8dO8zqRJD7JTWizNJ0Nr5rkfTqhyaqFVKYsjz6oyXTd8UQ+8knI1pTqcLThGqj7DBFxm2UM9zzxnDJC1yVUBDvtulG5kMbfq3No6DI7IlrxwPaNueYj0lG5PAvgCiJaS0RLANwG4GD8BiLqA/BnAG5ijL1uv5liyMKVI6TlOGyLpEOjk8LFlmyjyNvAtYrENnzweNAdMxtqD1ttsfEMGfGbbs6CMRQ2N6qIzrTrhvc70XttRpX6oFoENDh0xtgMEX0WwDCAbgAPMcbGieh+ACOMsYMA9gG4AMBj1PLdepUxdpPDds9DZ2OlmSAT4qv7PFuJj1yqSGzDRcBW2nao3qni7pL3umyLDaj849+abmL/jk2FzI1Kyki7bnhqL97zddpgAh9Ui4BmYBFj7GkATyeu3Rv7++OW26UN1SZM67NsO+uczcRHZUNZDiDeBu8CwCtutmX96tzalYSuKjD6++5Hj2GWkyb74nZIfxFzEyfYk1PT6CZawNFG7UrTtriHlexA0GmDLnxQLQIlixTlgbcJZX7ZOnBBfG0lPgpwBx5XePbcDFfPfPi4uZeWLZi4Lfpsd+G1LW/XSVtt8KWCWekJuguR3gXxFZ3gpomPAsw8j0y9lJJEYO2up7j3FRnMlsZ1FShe7cWDD3UIbLTBlzEuPUEH7Iv0LoivLyd42TE0Orkg4di8ayE6uSkbgTO+iNLJd/tcwNvkELWle84SHWri169S4fDemWeCrsqlz7UBF54ZZfNO8RV7Do53uBI25xj2HBzvuFfleaATMOSDl04SNtrkKgDPNBe6DU+TrHWEuwRJeOJtSPuOvHPDV4JDtw1X4lNZjIM+g5cdUnQ9TTnCkVfe6Mgzk3dwiApZ16fLkH9T9YUNyTWtyiQaB57BONmGtO/IW6UUCLoAgfjmD9uiqUw1IdpoDx95tSPSOO/gEB1kWZ8uiUwR+v20ahuR80M3UYf0bDty1ZUNJhD0AuFLrUwfoMs1rhJEN65aXuu4JuP+RJkRbUYa+wqXRKYI/X5aO4eov3OMcYO0bEauurLBBB16QbCpW/MlIVkW6Eba7d62oaP+aK2b5jNHxiGzW5hsKNvcVNHzZTNCMokibA5p32kyDqII1LPnZqTzl/d4BA69INgSe4uoiuICulyjqYgu4v5k8QtJ2OSmfJgvlx5Xrt33ZFKt6TtF47Bl/Wps3nuI+6xkBGqUxTHe9zjydmcMBL0g2BJ7ffDjtQET0dSGfYO30basX43Hj046dS31Yb5cExlX9ifd4iAm7QTkayD5jn3DEx0GeNX85WmPCwTdEVT6cVu6NV9ySGRFEX76vI3Wf9mFWoQurf3Dl/kqo9HfxWEYH4eh0UlumgQbaX3zQiDoDqAjVtsiYD4GvqSBL5F2OoQui9qkKvOVB5KHps1kebx3iVwY4+/wff4CQXcAVWGBONHISsB8jEDNUhSkDFxjFk5RNV++ej7l3S7eoenSxqFKwx1P6+vbfosjEHQHUBUWALJlk4vDF842gg9GP9fIInbL5uvzQ89x/eDjvysCLudUdFCIaucmibotYiqbOxtpffNCIOgOIBMPbRnAfOXkfDD6uUZWsZt3kA+NTi4g5hF8GDtXcyo7KEQElqEV+DPLWOpsqjyI5pQXZOSzJBn80DnI6iesqpqSVeeXd34IE/huNLIBF77F+4Yn+IV6UfzYuZpT2UEhOxxnGZsfb1uEVTSnX7z1Km+JNw+BQ0/AhngZ3ScrLJAFLjmmrFy/70YjG310IXbLiGPRY2cyp6Lx5V2XHRT7d2ySVluyLbn4rkrRRSDoCegYNHUQ3efCgOKCY9I5yHSIoc9GI5u64KS7277hCew8MGb9ICSg8LHTnVNZwjOeb7esSHWcwLr0bonDZ1WKLoLKJQGVQdNErSELPc8CF6HbOqlmddQ8rvpsAy4K+dpSf/FEfgJwx7VrCh873TkVje83vvca97qqSPVAX6tQesNhqgIRik7PkBaBQ0/AtkHTxamvwzGZqhZUXL9p2TOb3KstpJFsVONoS/3lu8ivs45lzBAPukWq85b6yuypVXmCbkrYVJXSizZQAfzK5stq54Ut3oK868AY7ntyHLu3beD2X3SQrazXsHnvoVRib9qN4cqDx1S/r9N+Hyru5IW00c+RV0oSukWq8z7syuypVWmCLtuQAH+BuDZo2sQ7M+fr0ceTBImCJGSJhHgHWa2LcObcjLCoBCAfjzQbwyV3ZMLp6YSBA/JDQpdIl4EjzBL9fMvVjcw5cvLUb5fZU6vSOnQRQdlzcFyq9xzoa+CLt17lXemxOER9u+vAmJCbju7h6Yx5etILlvWgOStyplOPR5qN4ULPHUFXF6wbBg6I3d22rF+trVt32Wdb0GmjaHy/MLDRW7sKD1lsVEXr3ivNoYsIB4/jTHJevus0s3ALsuox8f6JKt4D0ArqkHkxmLbNFnekw+nphoFHz4t+E18nJtJJGThCk/TGvPEtkweJKLXy5NQ0Nu895LWkVWmCLiIoIuguTh8gM97q/DbLOxq9dWVJtqHRSfzs7ZmO67VuknL1afTcskM3jW5aNww8Am+diCoi8Z7tu+8+UI422kLSZTKebkBGpH3QvVdW5SIjKLxyZUC5FqcqGlUEE7VRlojIfcMTaM51qixWLOlRGqV1K8OoXAbTuhSK1gEvDNz0GbrVcHxS7wHlaKNNxF0mRekYkvBB0qosQZcRlN3bNpR+cUb6ym4i6X299Vpq3WUWn3LRIn5LYmCNv7O3vvDQjQy6cWKs0uum1U3bCAM3IYA+++5HKEMbXcCESLss7aeLyqpcZATFd/24LkTRqBHqtW7suYnvpmjyDp3fx1UbK+s1EAE8m6LO4h7o06sMo9psaTkmG+tD9xlJldD+HZu8XYd5qCBNVWRp3T11f2eiavIhSrqyBF01ET7rx02Q1Pe5yESnQtIYJHJzNFncOsRYNcdZ9L421ofqGT4Y0XyC6XhkiXMQpSg4fPz0AiJvQqR9YBSJCVyzXKO/v5+NjIw4e35y0oDWRJRFTCxDoEkEWeBRhG4iI5WF6Jlxg6xqjn1fAzp9XExQjUdyT5w9N8N1elCNn+g9vFzrD2zfCKCTSPOu5bWmiOgoY6yf910lOHSblcBdtsXkGWXi3HSMPnOMGbVdhzNSzXHy+0gVtPPAGPYNTxR+SPpgRPMJsvHg7QnT56i+Fxk/n9l1XYfnlI1Edi6gRdCJ6AYAfwygG8CfM8b2Jr5fCuCrAK4G8BMAOxhjL9ttKh9ZKoHbHnRbhNgH9ycT6LhQmhqGdA9k1RxH3/t4SC4mV0AdyMZDFRuQvD/Ne3iIE/+IXvB+G9+fRa41pZcLEXUD+DKAXwPwIQC3E9GHErd9BsCbjLFfALAfwB/abqgIaT0ZXBSJsBXxVzbOTeVCmTYFbOQ69tLeGzu4JFP4GI252FwBVZCNh+7a1xk/UWZLHqLDIU4vRNBJZOcaOm6L1wA4wRh7kTF2DsAjAG5O3HMzgP/R/vsvAPwqkcKfzhLSEj8Xg26LEPvg/mQCkash4E8KWB8PycXqCiiCbDxEaz+NWy7vPXdcu0Z6uOpICFEbi1xrOiqXBoDXYp9PAviI6B7G2AwRvQXgPQB+HL+JiO4EcCcArFmzJmWTFyKt2Opi0G2J0D64P5kirtrw0Zjrq3qjKt5WtiAaD9GeSOuWy3tP/2UXCteuii7E92eRa02HoPM47aT9QOceMMYeBPAg0PJy0Xi3EmmJn4tBt0WIfXB/SgtfCVQZD8mA88hjT8jWrkzvnnQRLnKt6RD0kwAujX2+BMApwT0niagHwEoAb1hpoQJpJ9rFoNtcdL4SxrKizIdkQAtF7gkRveCpeIpca0o/9DaB/nsAvwpgEsCzAD7NGBuP3fM7ADYyxn6biG4DsJ0xdqvsua790HXgq3ogICDAP/hCL2R+6FqBRUT0SQBfQsttm/3eEgAABJhJREFU8SHG2H8govsBjDDGDhLRMgBfA9CHFmd+G2PsRdkzfSDoAQEBAWVD5sAixtjTAJ5OXLs39vfbAH49SyMDAgICArKhstkWAwICAhYbAkEPCAgIqAgCQQ8ICAioCAJBDwgICKgICkufS0SnAbyS8ucXIRGFuggQ+rw4EPq8OJClz5cxxlbzviiMoGcBEY2I3HaqitDnxYHQ58UBV30OKpeAgICAiiAQ9ICAgICKoKwE/cGiG1AAQp8XB0KfFwec9LmUOvSAgICAgE6UlUMPCAgICEggEPSAgICAiqB0BJ2IbiCiCSI6QUS7im6PLRDRQ0T0OhE9H7t2IRF9m4j+of3vqvZ1IqI/aY/B94nol4preXoQ0aVEdJiIfkBE40T0ufb1yvabiJYR0d8S0bF2n+9rX19LRN9r9/kAES1pX1/a/nyi/f3lRbY/LYiom4hGieib7c+V7i8AENHLRPQcEY0R0Uj7mtO1XSqCrlmwuqz4CoAbEtd2AfgOY+wKAN9pfwZa/b+i/d+dAP40pzbaxgyAuxljVwK4FsDvtOezyv1+B8B1jLGrAGwCcAMRXYtWYfX97T6/iVbhdaDAAuyW8TkAP4h9rnp/I2xhjG2K+Zy7XduMsdL8B+CXAQzHPt8D4J6i22Wxf5cDeD72eQLA+9p/vw/ARPvvPwNwO+++Mv8H4C8BfGKx9BvAcgB/h1aN3h8D6Glfn1/nAIYB/HL77572fVR02w37eUmbeF0H4JtolaysbH9j/X4ZwEWJa07Xdqk4dPALVle5xNDPMcZ+CADtf9/bvl65cWiL1n0AvoeK97utfhgD8DqAbwP4RwBTjLGZ9i3xfi0owA4gKsBeJnwJwO8DmGt/fg+q3d8IDMC3iOgoEd3ZvuZ0bWsVuPAIWsWoFwEqNQ5EdAGAxwHcxRj7KRGve61bOddK12/G2CyATUTUC+B/AriSd1v731L3mYj+OYDXGWNHiehj0WXOrZXobwKbGWOniOi9AL5NRMcl91rpd9k4dJ2C1VXC/yWi9wFA+9/X29crMw5EVEOLmD/MGHuifbny/QYAxtgUgL9Cy37Q267fCyzs13yf8y7AbgmbAdxERC8DeAQttcuXUN3+zoMxdqr97+toHdzXwPHaLhtBfxbAFW0L+RIAtwE4WHCbXOIggN9s//2baOmYo+v/om0ZvxbAW5EYVyZQixX/bwB+wBj7j7GvKttvIlrd5sxBRHUAH0fLWHgYwKfatyX7HI3FpwAcYm0laxnAGLuHMXYJY+xytPbrIcbYHahofyMQ0Qoielf0N4DrATwP12u7aMNBCkPDJwH8PVp6xz8ouj0W+/UNAD8E0ETrtP4MWrrD7wD4h/a/F7bvJbS8ff4RwHMA+otuf8o+/1O0xMrvAxhr//fJKvcbwD8BMNru8/MA7m1ffz+AvwVwAsBjAJa2ry9rfz7R/v79RfchQ98/BuCbi6G/7f4da/83HtEq12s7hP4HBAQEVARlU7kEBAQEBAgQCHpAQEBARRAIekBAQEBFEAh6QEBAQEUQCHpAQEBARRAIekBAQEBFEAh6QEBAQEXw/wHdZASi+iH2gAAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plt.plot(randu,'o')\n", "plt.plot(records,'r')" ] }, { "cell_type": "code", "execution_count": 14, "metadata": {}, "outputs": [], "source": [ "length=500\n", "randu=np.random.exponential(1,length)\n", "records=np.zeros(length)\n", "records[0]=randu[0]\n", "k=1\n", "while (k]" ] }, "execution_count": 15, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAWoAAAD4CAYAAADFAawfAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nO2df5CcxXnnv8+uBmkksFaAcMwaWbjMiZhgdoPOJlEuZekchI3x6QyxzrFzqYrvVHV21QWFW5eoUEHy5QrlVLaI71K54CSVSpkQgcF7GJIIXyRXyjjClrwrg2Lk+Bfgle8kGy0E7aCd3e37Y6ZX77zTP9+335l+330+VarVzsy+091vv08//fxqEkKAYRiGiZeBfjeAYRiGMcOCmmEYJnJYUDMMw0QOC2qGYZjIYUHNMAwTOcuKuOjll18u1q9fX8SlGYZhKsmxY8d+IoRYq3qvEEG9fv16HD16tIhLMwzDVBIiekH3Hps+GIZhIocFNcMwTOSwoGYYhokcFtQMwzCRw4KaYRgmcgqJ+mB6z/jEFPYdPIlT0w1cOVTH2NYN2DY63O9mMQwTABbUFWB8Ygp3P/YsGs15AMDUdAN3P/YsALCwZpgKwIK6Auw7eBKN5jw+dPwpDL96evH1s19bBmy6uo8tY5glxsUXA5/8ZPDLsqCuAKemG7jk/Dn897/9LABgAXThzcN9ahTDLEXe+EYW1IyaK4fqOD91FgBwz80fx+dH3wcAGB6q4+ldW/rZNIZhAsBRHxVgbOsGvIEWAACzA621t14bxNjWDf1sFsMwgWCNugJsGx3GqhffCvwPoLmshmGO+mCYSsGCuiL8ytvWAAD2//o7gQ+xuYNhqgSbPqrC7Gzr5/Ll/W0HwzDBYUFdFc6fb/286KL+toNhmOCwoK4KrFEzTGVxslET0Q8B/DOAeQBzQoiNRTaKyQBr1AxTWXyciZuFED8prCVMPlijZpjKwqaPqsAaNcNUFldBLQA8RUTHiGiH6gNEtIOIjhLR0TNnzoRrIeMGa9QMU1lcBfUmIcTPA3gvgE8Q0S+nPyCEeEAIsVEIsXHtWuVBukyRsEbNMJXFyUYthDjV/nmaiL4I4J0A/r7IhjGeRKxRc61shsmHVaMmolVEdIn8P4CbATxXdMMYTyLVqGWt7KnpBgQu1Moen5jqd9MYpjS4mD7eCOCrRHQcwNcBPCmE+Ntim8V4IwV1ZBq1rJWdpNGcx76DJ/vUIoYpH1bThxDi+wBu6EFbmDxI00dkGvWp6YbX6wzDdMPheVUhUtPHlUN1r9cZhumGBXVVmJ0FajVgIK5bOrZ1A+q1wY7XuFY2w/jBZU6rwvnz0WnTwIXDdTnqg2Gyw4K6KszORudIlGwbHWbBzDA5qLygXjIxvJFq1AzD5KfSglrG8MrwMBnDC6B6wjpijZphmHxUWlDLGN7//PRDeNtPX1p8vf74IPD2N/axZQXwta+xoGaYilJpQX1quoHlc7P47a8+iJ/W34Dp+iUAAAKA1yuWGbdiBXDbbf1uBcMwBVBpQX3lUB1zL7VKaH/6l38dfznyXgDA8FAdT+/iA2AZhikHcQXdBmZs6wb8THMGAHB2RUub5hhehmHKRqU16m2jw7h842XA54Dp+hswXOWoD4ZhKkulBTUA/NJlray4h+6+FXjHO/rcGoZhGH9KK6id46N/+tPWz0sv7W0DGYZhAlFKQe0VHy0F9WWX9bKJDMMwwSiloJbx0cOvnMY7fvydxde/8f1/wLZ/e33nh595BqjXW/8YhmFKSCkFtaxlvPdvPot/9cJk55t/qfiDt7+9+EYxDMMURCkF9ZVDdUxNN3DJ7Ay+Mfx2/M7WjwMArrhkBT7/H97V/QdvfnOPW8gwDBOOUgrqsa0bcPdjz6I2P4czq4bwnbXrUa8N4uMfvB74OQ69YximWpRSUEuHYf1zC5gbWMbx0QzDVJpSCmqgLazfUMNbf34d3luhdPAlU5aVYRhnSiuoAbRKe1aoBvOSKsvKMIwz5a710WxWSlDLsMMkjeY89h082acWMQwTA+UW1PJA14ogww5dX2cYZmlQfkFdIY36yiF1Uo7udYZhlgYsqCNibOsG1GuDHa+FKss6PjGFTXsP4epdT2LT3kMYn6jYwQkMU2HK7UysmI1aOgxDR32wk5Jhyk15BfX8fOtfhWzUQEtwhhaeJiclC2qGiZ/yCupms/WzQhp1UfTbScmx4QyTD2cbNRENEtEEET1RZIOcYUHtTD+dlNLsMjXdgMAFswvbyBnGHR9n4m8B+HZRDfFmdrb1s2KmjyIo0klpg2PDGSY/ToKaiN4M4FYAf1JsczyQgpo1aivbRodx3wevx/BQHYTWKez3ffD6npgf+m12YZgq4Gqjvh/AJwFcUmBb/GDThxdFOCldkCVpVa8zDOOGVaMmovcDOC2EOGb53A4iOkpER8+cOROsgVpYoy4F/TS7MExVcDF9bALwASL6IYC/ArCFiD6f/pAQ4gEhxEYhxMa1a9cGbqYCtlGXgn6aXRimKlhNH0KIuwHcDQBE9G4A/0UI8dGC22WHNerS0C+zC8NUBY6jZpiKwvHr1cFLUAshvgLgK4W0xBc2fTCMFi4bUC3Kq1Gz6YNhtHDZgAtUYWdRXkGtMH1U4YYwTAg4fr1FVXYW5S1zmjJ9cKoyw1yAa5u3qEpmbLSC2lo/OWX6qMoNYZgQcPx6i6rsLKI0fThtV1KCuio3hGFCUFRt87JRlczYKAW11I7v/OqD2PH1xxZfp98noNbeBMzNtX4uXw6gOjeEYULB8eutnUVS6QPKubOIUlBLLXjk1Hfw2kUr8cXrNgMACMCOX37rhQ9ecQVw9dUAqnNDmGrAju04qMrOIkpBLbXjFXPn8YM1V+K+zb8JoJV+vGPXFuXfVOWGMOWnKpEGVaEKO4soBbXUjlfMzeLV5asAuGnHVbghTPnhGOb+UdWdTJSCWg7sqj+Zw+nacgxXaMCZ6sOO7f5Q5Z1MlIIaaA/sJYO45sb1uFlj7qgqVdUKlgrs2O4PVd7JRBtHDQBoNIAVK/rdip7CiTvlh2OY+0OVdzLxC+p6PFqINQknAJy4U364Bnd/qHI2ZrSmDwDA669HI6h7Zf+qslawlGDHdu+pcohuvBq1EFFp1L3SdKusFQC92ZUwS5Mq72Ti1aibTWBhIRobdRGarsppWGWtoMpeeSYOqrqTiVejfv311s9INOrQmq7OaQigsloB29+rD++YiiFejbrR1lQjEdShNV2T0Hp615ZKCOY0bH+vNrxjKo54NWopqCMxfYS2fy1FoVV1+/tSh3dMxcEatQch7V9LMSmin/Z3TiIqnqWofPSKeDXqyGzUoVmKSRH98spzElFv4B1TcUSnUUvN52ee+yYeBfC1qXP4xX43qgCWarW/fnjlq5xaHBNVjljqN1EJaqn5zJ6fxdtmZwAA/+uZKZz+xalKPlBVDSWKDd6S94alqnz0gqgE9b6DJ7HstVdx5I9+E6vPnwMAvDywnDUfJhdL0R8Qgix2fVY+iiEqG/Wp6QbWzLyK1efP4YkNv4R7bv44/vGKq1nzYXKxFP0BeWG7flxEpVFfOVTHwMsLAID/c827MN4+gms4ofmw957xhbfk/rBdPy6iEtRjWzfgjz/3fQDAAhGATs2HA+qZrPCW3I9+2PVZCdMTlaDeNjqMi99zDfBHAEBdJ7vwKs8wvaEou75OGLMSZiYqQQ0A77l2LQDgsx+5EdjeebILe+8ZpjcUEWpnEsashJmxOhOJaAURfZ2IjhPRCSLaU2iLhGi3rLtpHFDPML2hiOQkkzBmJcyMi0Z9HsAWIcRrRFQD8FUi+hshxJFCWrTQciaibaNOwgH1DNM7Qtv1TcK4nyGUZbCNWzVq0eK19q+19j9RWIukoFZo1FUuDO4Kl5FkyoppR9yvEMqyhCE62aiJaBDAMQBvA/CHQohnFJ/ZAWAHAKxbty57iwymD2Bpe+/Z4cKUGdOOuF8hlGWxjTsJaiHEPIARIhoC8EUi+jkhxHOpzzwA4AEA2LhxY3aN22D6WOqUZVIxjAqbMO6HElYW27hX1IcQYpqIvgLgFgDPWT6eDYtGHYIy2KRUlGVSxURZ77UPZepjbDvispQXsApqIloLoNkW0nUA7wHw+4W1SGOjDjUZy2w+KMukioHxiSnsfvwEphvNxdfKdK9dKfN8NtGrxacsAQouauubABwmom8B+AaALwshniisRQrTR0iDf5lPoeCaFW7I+ZIU0pKy3GtXyjyfdfTSwVeWAAWrRi2E+BaA0R60RX5h62dbox6fmMJdDx/HvOg0e2e1zZbZfMA1K9xQCa8kZbjXrpR5PuvotS8mNnOMiugyE5OmD7mypoW0JMtkLLv5oAyTqt/Y5kVZ7rULZZ/PKqq4+OQlqjKnADpMHzbNKMtkXGrmg6UYd22aF1nvdazjWMX5zBnI3cQnqBOmD9MKmnUylsUmFYKyBPOHRiW8AGDNylqmex3zOFZxPldx8clL1KYP3bZukCjXZFwq5oOlGncd2pYf+zhWbT6zL6abeAU1Eca2/gtl6EzZNYZesZRtfSGF11Iex35RtcUnL/EJ6oTpg1fWfFTR0dQPYhnHMiW2MGGJT1CnEl54Zc1OWYL5YyeGcaxqYkuvKetiF58zkWt9BKOKjqZ+EMM4VjGxpdfE7BS2EZ9G3YNaH0uJmHYkaW1m87Vrcfj5M6XQbvo9jmwnz0/sTmET8QlqQz1qpryotu6fP/Li4vu93Mr3e/ub5ftjsZOXmTIvdvFJQzZ9VBJb8hLQm618v7e/Wb+fY4vzU+ZEmvgENZs+Komr1jI13Sg0+6/ftt6s3x+DnbzslHmxY9MH0xN0W3cVSU0TCGsK6ff2N8/3h7CT99vs00/KHO4br6Bm00elUIW42SjC0dNvW2+/D3F1CfGrsjDvt1M4K/GprWz6qCSqrftHb1q3+LuO0Jpuv7e//fx+F7NLv234jJp4Neo+CuoqaxT9xKTNbNp7KIimabt3/d7+9vP7XcwuZQ5hqzLxCuqcpo+swpYzwPpDiOw/13tn2/4WvVD3a/vtYnbptw3flaWmTMVnXwhg+sizfet3VMBSJURUQ4h7V+Wtv4vZpQwhbFW+Rzri1ahzCOo827eyaBRVJK+mGeLeVXnr72J2iaGuiY2Y7lGvNPv4BLXUqHOYPvI8sP2OCmCyE+LeVX2hti2G/bbhuxDLPeqlmTQ+00cAjTrP9q3fUQFMdkLcuzJs/Ytm2+gwnt61BT/Yeyue3rUlKiENxHOPemkmLZ2gdjm7Ls8Dyxlg5SXEveOFOn5iuUe91OxLZfpQbTXuPDCJOw9MAmidiXfvbdd5bd90NiYWzOUk770rw9Y/FGWNnIjlHvXSTBqfoDZo1LbCPmdnmhj7wnEAbg8sh+IxKpbCQl32uR/DPeql47VUpg+XLUVzXjjbiDgUj1mq8NzPTy/NpPFp1AbTh2thH1cbUSzeY4ZRUaRpoui5X1azii+90uzjE9QGjdq1sI+rjSiEjanqE7Lq/YuV8YkpjD1yHM2FluIyNd3A2CMXzHp5KdK+WnazSoyUyvSR3GroqA2Ss40or/e46hlSVe9fP3CJWgKA3Y+fWBTSkuaCwO7HTwRpR5GRE2xWCY9VoyaiqwD8BYCfAbAA4AEhxB8U1iJLwktyqzE+MYXdj5/AdKMJoDvqw0Ze73FMGVJFUPX+9RofTVPO6TTJ1/PsdoqMnGCTYnhcTB9zAO4SQnyTiC4BcIyIviyE+MdCWuSR8BLCPpTnGlWfkFXvX68JufCFMC8UZV+tWnZvDOY/qzQUQvxYCPHN9v//GcC3ARTXygjKnLoSS4ZUUVStf65mh6LwWfjWrKwpPytfj9m8EEtCig6feRCL+c9LGhLRegCjAJ5RvLeDiI4S0dEzZ85kb1GAWh+9IvYJmZes/eu3QNS1qd8PnM/Cd+9t16E22PkM1AYJ9952HQBoo59cjzsrkqLC1kLMK995EMuC6CyoiehiAI8CuFMI8Wr6fSHEA0KIjUKIjWvXrs3eohJo1HLC7DwwieXLBrBmZa2S6eZZHrgYBKKKGB44n4Vv2+gw9t1xQ8fY77vjhsWxH9QoMgREsUjKeiH7t48AAHYemMzVnlDzyncexGL+cwrPI6IaWkL6QSHEY4W2KBJBrbNLpW2D040m6rVB7N8+Ep2ADmFb09kxddeO1QEZwwPn68Az2ZDnhVC+LoUYcKHEwp4vnfBysociZJheqHnlOw9isbe7RH0QgD8F8G0hxGcKb1EEpg/TBItVEKUpMpbVdO0YBKKKXj5wpgUylANv2ONU97Mzzb7EMYd8VkLNK995EEt9bhe1dROAXwewhYgm2//eV1iLItCoTRPMd8L0y17repBplraZrh2rA1JldgCAmdm5oPekV6YfXX909MOuGnLRDjWvfP0usVTTtGrUQoivAsaDosMSgaA2TTCfFbmfGVq2hyRP20zX3r99JAoNJI3sUzLuHgivbfZqx7VtdBhHX3gZDx55EWojSDe93tWE3MWoNNvaAGFmdg5X73rS2bSXJX48hgJQ8XnsIjB9mFZvnxVZ99DemdOx4oJNA8njXDNdOxYNRMW20WGsWt6tm4TUNntp+jn8/BlnIQ30flcTMioqPa/qtQE0FwTOzjS9dy6xH4ygIj5BHYFGbZpgPoLI9HAWHQ1he0jyCBTbtWN+EIoWpL00/fi0uR+7mtCLdjKS5PXmQtf7scSRF0GpijL1Ctv2yHUrZKv2V6QT0taHPNvSWAq3u5B27A2trOHsTHd6tqnfPtEzvXQ+6e7hcLuNMdyfIswG+w6e1O4kQiy4MWQipolPUEdg+gDCTDCXan8ht8TJCba6XgMRMD3TLOS0adv4xDDZVXb42gChNkhozl941E399rXlh1rEXMbPdA9jsKsWhemZybtzibXyX3yCOgKNOhTJh1anWYfaEqviuyWqyVakVhzLZFfZ4ZsLAkP1GlYtX+bU7yzOwbxC0nX8yrSzCYluJ0FA7p1LrOG38QrqEqSQuyAf2vTDB6g1uayaqO2YMtVkK0rrimWy6zSvVxpNTN57c65rFBlB4TN+Vdacdah2EgTgIzetyz0WseYBxCeoIzF9hMZF+ykiZM73My7YFpNYJnuI8LB+ZKbFMn6xUuROIpZMxDTxCeqFhUqYPVTYtJ88mqjLMWW9Or2jV5PdtmCEcOz1IzMtVmERE6ZnKY9/JJZMxDTxSUQhKiOofTP/XJJUdNezZar5nHxjwiX+2jcLMEuGpEsGYIjwsH7EhZepKmPIzNt+VMdLE2seQJwatafZo6gIgzzXzWLGMGlStuvJa9718HFlwZ5VFy0LMiYu23KfLMCs5h7X3UcIG24I56B0KA8SYV6IxRC6IiNHiiak0zjUtUL4R2K0+8enunqaPlQr6M4Dk7hn/Nlczci7MmfJ/DNpUrrr3fXw8cU2bRsdxoKmqtormqOdfHFN6HDNAsyaIdkrO25eLS85j4ALVe9iKf+aB597ZxvHUGVoQ82L2GqqxyeohfDSqFU3WAB48MiLuQY378TJMmFM2y7d380L0fHAF50Z57MtdxmDrEWudAkPtsQVn4cvRIElUzSObj7dM/4sdh6YjK6mdxrXe+cyjiGr47m+rpsPMdZUj9P04aFR626kAHKFg+WdOFkdQrptl8lZmNzaFe0M8dmWu4xBniJXabIkrhx94WUcfv6Msi8httG2+aISaqpCS1nCG0OaBFXXcr13LuMYyoHqOv/LVso4Po3aU1CbbmSebfCQ5sw614mj0zw3X7s205bK5iyUfe2FM8S1loeL9p23yJVkeKiO229sHVygGlvdw/fgkRe1mlMILc82X1RCLUR6dEitUHetzdeudbp3LuMYyoHqOv9DljLuBfFp1J6mj7GtG7DzwKRycmfd7o9PTOG11+e6XveJnFBpnusvq3doSz4OE5uzMNnXWJwhLtq3j4aue1BkRprJGWXaeSVJak4htDxTGQEfoeb7vSG1Qt21Dj9/Bvd98HrrvXMZx5AOVJf5H6qUca+IT1B7atTbRtV1efNs9/cdPInmQv7IieSEGZ+YUi4oPg+P/ExscZ55TzRxXVhMD5BNMLnEmUvkQxzCjJQuI2CL+giVHh1SKzRdy+Xebb52rdPz6XKtUOYc01yKMZa69IIaAH5v2/XY+JZLg9njTKnHaVwnTqgtbWyhW72s62F6gHYemFT+jUnoEro1auCC5mQaax+B4bPDcU2Ptn1/SK3Q5VqmM0YfPTbVMc4E4PYbh63jmH5v87Vr8eixqSBzzVbQCojnGQNiFNSepg9JyO2+6yT3EVIhK37JvsqJvPPAJPYdPNmXydRLx4vpAdIVvjIJ3fSDD3RrTqp5VeTi5CIkXL4/pFZou5avY06gdeiB6e8AdL0XyskKhCtl3CtIaOJu87Bx40Zx9OjRbH/8iU8ADz8MnDkTtlEaVCs6oDYvpJ0Sm/Ye0tYDfnrXlo7XdJ8lYPEEcx8tTVfkqddZVFfvelKplRKAH+y9dfH3osueZh0Pkyaoa6/rfS+qz/34/qzjcartgExDMNfTBqB8T0V6rpUVIjomhNioei8+jbqHtT50K/p9H7y+w0kiazsnNVf5eRUq7dm2pfXV0nqlyYbYYvfCPJJ1u5pFY3ax/xbZZ1f7s04rtN1T3fu6dmd1zIWyow8QeZ2b6EMMddWBGMPzMpo+smATdvLYn/NzCx1ns409chxjXziuva6cOMkQMVXY0P7tI/i9bddb26KiFyFEqrCssUeOY/RTTy32TxWiRe3Pyv6Hyjqz4Ro2aMPWXpekiiL7nCepyRa2lyWsz9QeU9id6e98zIHzQhSSmBJT4suS1qhdhJ2u+LyJdJowcEG7yaKVqOhFCJGu7/Ioq6npBh49NoXbbxzG4efPYGq60eGgk/3XxT7HWrbTdi90uyO5OI1t3eB1P31rgeSxP9uUE937ux8/YYyXz+qYM/Vj7JHjXc/aAAGr6zVMzzQx0B4rXV+S3DP+LB565iXMC4FBInz4XVctKklZx6qXLGlBrRN2A0QYn5gybnNdyVumVOXAlA91Omohi7PItLVz6buMp3161xalrbLRnF8UPra+2doF9MYTb7sX6ZA71eLkejZj2kSiW+ST+Jh50uNoM9fp7vl0o7n4TKTJ6piz/d2eL53oGsMFAay8aBkmfvdmXL3rSWNfJPeMP4vPH3lx8fd5IRZ/NwnrmBJf4hPUPTR96JIRZP0MwK3Osw2XG+uiJaUfaoELIWYmDUyHzY6qEza6/pnqkdRrg84aoKpdY48cBwiLZx0WYef2WQSl8NEtTsuXDTj12aUWiE7IucQcp8fRFpJomu8mhSNrlITp76Y1c0/OM1fl5qFnXlJe56FnXjIKat/yBkUqEfHZqB006lCVraTdeFCxMMiHRGVjkwekJqnXBrEmR9q5yoadjljQhTpJb7/vxDBt7XTZmSpk/+o19X0bqte80tp1JpfkgbTJtoYgXeVOLoKwtNcUc+/SZ99aID7o5kt6tneYGww7sl5rkjZbvGvauWo3Z3pd4nL98YkpjOx5CncWXEQrPo3aIqhDe9O3jQ4bkyV02zPda3liV21aSeitmO7vpqYb2lT1NLJ/4xNTmGkuKD9D5Kdx+fQnlPCwLYI6TFqXS59tO7Y8PgdT2rwMnVOZKVQmh7xtyYJtl+lqAtKZ3lQKWhLb9U1FwkLbsuMT1BbTRxEGfps90WZjU7WxiC2Qi009xPUAu7YhkVripr2HtJ/RbWGztEv1WRdsGXA+oZZJ8iaW+NYC8cEUp2xafO697booUqhdBLHLYvjhd13VYaNOvq7CNlc27T2EU9MNpTMzScgdSHyC2qJRh9YqQxRgSlJkRpOLTd3FbpmMD68NUpdJwZXhttYIhM28VPWzNkAdNmrAXXi4ZMBlbXvW+G3V37tGfdjI63B2zY7shWM3xPMk7dAuUR8+2ZI2ZSbkDsQqqInozwC8H8BpIcTPBftmHRZBHTosTVeAadkA5Z4goSez/FuVWcJlV5GehNONJmoDhAFqedN9SD/wPsWEbOPiY25yGU9bTHOW+tZJ8giT5FjkEc7J64VwOJv61Mv6Lunvzfo8/d62663heIB+rtz18HG8ob5MO1fShN6BuGjUfw7gfwL4i2DfakII/PPsPG5pby/SNyR0ZSudJthoLmQyJ0iKmsw2m7qpPSoBb4sJlxBaJqLpmaYy1VqlvemKCbmMS1pYS8euacuuI+suzCUdP4/wKGKO+Nras7Q/tPnRpQ33jD+buUSwy3cl57GKeSGcIqAAYM3KGu697bqgi5ZVUAsh/p6I1gf7Rgs/+uk5nH+tuThgqqQRIJwdOGs4ko0ig+V9dxVSILjandPoamaYtLehdtr9g0dexOHnz3QUT9JpLDsPTGods1PTrbMwj77wspNmlMQ2Xjo7ru0+qYTHzgOTuPPAZIf2qhMORcwR30Qb20KhantI86NLG9Jx0BLfsTKd8pMuzuXDIBEWhCjUBBSdjfrk/30VV6UCiNI3RArsENXjxrZuwJ0ZNFQbWTLTXAsD+e4qTLG6QEuonp9bUJYANW2Xddpb+nrJh88Ua5387IragPLaDx55ERvfcqnXfbaNV5Ydmu7IrLTGlxYCLmORZ975LOK2hUIn2FyTeVxwacODCiEt8XmedN8lbddZ6FUhtGCCmoh2ANgBAOvWrct8ndnzTQhF0Ef6hoTaNhYVjpS3VKrpAffdVZge/HptELs/cJ3T9Vyz3KYVdbvlw+cS0dFozmsXlixnYbqMV5btv+3R1gkB21isrqvj8V1wKUkq+2qrj64TbKpkHgCYmZ3zNheaQkQ37T2Emdk54zi7Pk+m7/IR0kP1GlYtX9bzIk3BBLUQ4gEADwCtMqdZr7OyNoAF6nYmpm+IbhLd6ahd26IfdFqVq03PVev1WeWTmoaPA0snEAaJOrQBk/No9+MnOgSwKctNx6npBvZvHzHW/3C9ji+m8criDHRtg04IyLFQ1bM45ynw0nNS1l5R7dBcxl4+a6Zknv3bR7rmxNmZpreyZFq4bQu6ylFt0tB9s4xVETO7PxDW9uxKdJmJ11yxqivqw/dsOVtmULoq1nSjCYiWEy/EgBgAABlCSURBVMCUReZSeUxmTO47eBK33zjclZkGoCOr0jeOOYuQ0mVYffpDNzg5y+5+7FmllqzLcjNlaG4b7czA1CUdDNVrXddOXqff5G2DHIuLV3TrSs15gTsPTHZk3eqycVVz8tFjUxjbuqGriqDNBAZ0PmumzMBto8NYtby77b7Zoqq56UrSUS3HxxQP7/Nd9dogPnLTukIPifbBJTzvIQDvBnA5Ef0IwL1CiD8tqkFXXrIcK4fq2sypxc9ZVkeTo0GXoiyLveiwhXmlt1yPHpvquLk+tReyFDLSkccBa3u4VVluQLfdF+jcGpuyu6TmojoLk9A6g8+ETzW6rFEbpkSVZFt1RfNlNqcpksDFDObjkDQt8rKQv0+EVQgbezqO3AUZTSSdyi47Bbm4AND6pOS1e2nScMUl6uPDvWjIIgsLGFq13BqG5fKgTE03lFtI3wlmC92R4WO2B8ZUeyG9xbr9xmHrMVE+uG7vXe3QElOWm8vW2LSIyPeSwloAePTYlNah6FONLo+fQ9Xuzdeu7TA56MZO9sWUaCOxmcF85rJvpqJtgQ+V0yDvtU4jttmFbcpEOu1c9yzbMjb7SXRRHxDCqcyp60qsevB8q2K5rNamB8Ym6HW1F1wO7HXRCF21Rh+NHzAvHPKBSJtMkjsQlzYdfv6M1zl5PtXo9nzphHGHZGufbfEzHVHlYoaQmMxgLnM5T6aiqY+ufpi8fh2bXdikwQ8SddxTW9RULCe6pIlPUC8sOJc5lQOY1tySqB5q3Y3afO3axTx+eZNcV2udIF5dr1kFvUmjsTlEXeJgXbVGV40fcAvqN3n0VWVM93zpRFdCjekaPt+Zft9kdlC1L0tEkUkg6JKWVJjMYC5RHqZMxc3Xrl0Mcc2T/m5KNXcdy6wmutX1mvb5T+6o5KlMuu8Bus2Xybb2U4jHd7jtLbcA09PAkSPWj7p6sQG1HTW9bVWZGmwC1uRRr9cGsaI2YLRD5onDdDnk1OcgVJPtzuYzkNdIjunM7Jyy7zrBk0Y6JVXXSB4KnMTkUJL90B1yYGtflq1xOrqICNrTSWS/kq/KkrrpyoTJeWMSIKb7rxPyuvnoexKN7ftDmRlGP/WUc9bgAAGf+VD3vAH8xwoIm4VYisNt5SS47+RpXNxs4IO7nrROAtftozwmCeg8wDY5UTbtPaTcBrs+tLpV2qQ55XVw2eyTrlXh5CJjaqct/Xh1vYZzs3Mdhf2lkEmHPbpu+c/OtGqRqNDFU7tWozNp3iEjbuSuSGc7T7cvGVonxzQtpNPCwbTzMu1IVAuzzqyU5SQa0/eHrCznU51xQajNoaY26XxQQLaQxCxEIajlJFg5/VMsn29ioR2YlXUSpHGxcZqC4V1PJ1E9ML6OC9VWUZWWDJht7Tbhu7pecyrXaIonTxd4StNcEEpHkI+H31SLRHXP0r4LneanGzvZ3pCFv2R7VA+6Kf14095DynFdedGyIHHKOlTj6mL7l59L3uvQRdRU+PZRtxiZ2mqSNXnT/l2IQlDLSXDsjz+Glc3z+Nq6dyy+ZxqELJNQkh5404Mr2wGYtzoqTThEureuEI3p2qYHqzZAODc7tygETGYI3TbYdTfzSqOJyXu7Qx7zJr0ArXFZv+tJLw1TYnJcqdqXtxqa7kFfEAI/2Hur19/4aKMu0VFpVELU9p06u37o6CUVWfqYdPLbzJ8uykXIHYKKKBJeZCf/65b/iHt+5T9hz3t2KN9Powpgdz1tMT0ZVddKCzQAeF1ziokuGQaA1zFUthue9mDrrm26zkXLBpxqUJsKE7lOTNVDn273UDsz1ITp3bMzTYx94bjX0UemsTO9lxVT8kjIv0mT7osNnRC1faeMrkjSaLYOPg49lmlU9+ujN63Tms2AC07+dKKQKklNKkWmZJmik7Ci0KilNvvQyC3a91Wo7MIuGrYqaUJ1LZUzTKfhm+Kofc4zdOlD8n2d9mi6zrlZu+Zh03pc2qk6X05ne1elqSexLSvNeYHdj5/wEgAmzTv9XvJkD5vHP8TOCghX0jfZF5MT1eQTGtu6QZnuDsB4+MSp6YZxnENEUuiusfEtlyrnVG2A8Orrza4a7HJhMcWUq67Xi9Nvooj6MEVv+EZF2Lz+ktoA4eIVy7rCwZJcvetJbWZZcrtqiphIfzaNy/ZLdU1VxEP6uqYoDhv357w+EbA/4V1X3WNVhT5bzHnedmdBF9HjUvoV6C79appzqu8OERJmi6V2ecZ00RUmu77JF6MTej7Pu8u9MTm909ie1/T1QobpRR/14eoAcsHVXtVcuFAIXOe0dE0mMDntTFsileNQbr8OP3/GmCCT1upVk2eNphylqqxpEpdazNtG9VUHgdYJOUl8bO/bRoe1i6SNux7ujJUNgU+atqmf040m6rXBxUXWRUt3sbfbGJ+Y6tCGk+OqSgjRoYuueKXRxO4PuJ+zaFLMGs15r52Ry71J7yh0uzbArXJhiHviSxSCGgjXeXkN11O0JT6JMcnJZ3Kq2bZEukkmC+3b0tYluqQCnSNHOszybuNUh6BKmvOiYzxdbe+2qAwbqvMj82ZvumSdupreTLVhigrz2v34CW30TDrM7ugLLysr7wH6eyJw4XlzUbJsjujpRtO5eqCvw9U2D6cbTYx+6ilrbHSvk1+icCYWwSWKqmQ20jfRxankErajc3LZMvdMD31y5TcJfJPDbPLem3H/9pEOp96K2gB2piq36ZDjoyPZPxdnS/LzOueuzekIdApDW8VDl8/o2q5ySLlgqg1z18PHuyrk+aCqsmfSINPf/+CRF7XjYHKoSYEvw1lNgsvFEe1agc/X4eoyD2VsdHqOyHEd2fMUxr5w3DinQlM5QW0qy2lDF6Hw9K4tXSUjTX+TxHQThzTlQAH9gasSWbMYMGsPtvbL9/dvH8H5uQWcnWl6Tb5toy0vuYrk2LiUmEx+XrVI7vvVG7DvjhsWv88ksuWY2CoeunxGVyaWyH6fVJjicueFyPzwmyKPXNHlHACd98REUkFRlWb1XbRN6ObVTOL5cPl8GtNiP91odtm403MqNJUT1LptVb02sFj7mAAMpmyoWT23LjdedRPHJ6bw2utz3t8nkaYFIEwYl4tA06ETZMnxTD/kqjrW6fFXLTLytR/uvRX7t49o61nbit8nX7d9RhVOaCsNoEP20+Xe+D78unvoWDpHS3J85PjbLpncFaYXHpdnxvWUG3lvhlKfV2nFyc/bFhvAvNibPl8ElRPUusFqNBcWt2cCrY6vrF3o/oqa31BIbWHngUksXzaweOiAa7v2HTzpfAK47ZougtL1Wq6vuxySoNPepZDNG1u7bXQYn/7QDcpY+qnp1lFOul1LUlC6LHSqnYcrWeJyJenxH5+Ywsiep7B+15NYv+tJjH7qqUVhpHU+C3SZjAao+zUdqvHJGlct/RA2YXlOoxGr2Dbqd4iBvJf3bx9xio3OkzMQimiciSEYn5gypkMnaS4IzCUEpU/Ovkv6dJpkynaejMokcmKoYsB9nRt5Sr+qDkmwEdp5rAo9M9UbSS5iPvHKPuVJAfdaz7p5m44wSscyy2Sfoy+8bDyoYPu/vKrLSZj8/tX1mjK2uDZIHZX4kp/XxU+b6rkkdylyDFRhf2lntA2Tv2fT3kPaiBrA7lTPkjMQmijiqEPgU0nPhMvR766x2pLaAAGEjkntct6grcpc1phhlcca6E6b1p1ErounDVERLY83PWvh+eT32sJDfcIGa4OEfXfYjzuT368a/+RJJlkq/kls90Z3T4laGrlqvspchLMzza4x86lx45qvYML2TNris3XPhC4JyyUPw5fo46htuDy8vpqODpeqYL62qHkhsJDKPNfVeZbIiaWb8EP1WmYhPfaF4x1V7sa+cBz77rih4/vSmmnSKaXb9ue10eU5ccX0/bp6I0m2jQ4vHv1lmgM+u6FVHsWTWoiu3x491tr+m+LqAftJ2qZ7Y6rNLS+rurru+LrxiSmcO9/tf9Flqupa7mNKMGVOAvbCSapMVN31QpY2dSVqQa3KXtI9vL5CwkWj9a2ypUNnik6e7LLakLWm0rTef8ObnL9fMj4xhZ0PTyL9TDfnBX774Un82rvWdbQtiYtjK6+NziexJI3J7KUz36QzQtPnM6a/f3xiCi+fO+/cn1ccI48uLFDddWRkyJxtrkrNV4fp3uSJVlDZ0E11mwFgZM9TVnOhrynBloClaqsJkw/p7EzTKUEoJNEKalv2Uvrh9RGeQ/UaZufmu+r8qlDd3CzVulQMElm3+ElNz/XcQBVyPHUP84IAPn/kReM1bBM9r40ua5KC7JuuxnO6XSrN3dR3mdiiE6Y6XBYIl1OEnEwthg/ZhF6enZAAOmzAur6svKglakxaryRLRjJgr0vto0i4VAvUnUpUBFEJ6uQktjkF0wPpKjwJauffAKk1X11sNeB3crIK1+LrPucG6sxEIUxDcixCm2J8HGoqTLWeVXZJ37FIZt654rpA5KnHIjHtDl2EXl7ndnL+2orv24Q0Ac5+jrRfwXRlXw3dZUxcylCEIprwvHRQue2hSD+8qlhXVfiR7qpvWFHzCnFLhprdv32kI45TV11xkKgVw60IbDWZFlw1TVOGXV77sRwLXSigTEv3weWeuzxgplrPqocmy1jY5uNQvWYNN8yyWNoC6Oq1QaOAkhq7KdtxbOsGY0lQF+T81S2qrpmbrlpvcu4A5vuzZmVtMYNWlYCjIsuYFJn0Eo1G7TOJXU9Y8dHQX2k0sX/7SOaaEOnvNVX0unrXk8o26ASIa+icycabR2tSaWUh6hzo7nlSQ3SJbx/SFJ8yxUeHCI9MYjspG/BfINasrOHWd7ypq15LOhrHpJW7OGdd7LsuTE03uhJPfEifBJ6MwFFFldjkhelMU1ttE1Ponomikl6iEdSuHfTxuPqYKK5sV4yzFWJxnfhHX3gZDz3z0uIku/3GC9f2PZ7INc7XpHnv3z7ibVeXiwvQGsPkSdUhDibVtTe5nMr4dt0DpcvwTMb/JtFFJOShXhvoGh/VPHJZIAjousbGt1xqXBh1wkSV4q4zmfmcO2hqe5bSDcCF5zrtj0nWEAG6T5IxtSU5R3VKTPK7VMfeTd57s1dIZlFJL9EIap1WJMniYPCJrdZtr5MaNBSeddXEH5+YwqPHpjomWdL5p4owsJlZALsWa1oAVNfYfO1aPHH8x+pz+WoDWF4bwJ0HJrWhenltca6abaM53+HsS7Zhz5fUleFUoXG2+eASCZSmNkCYWxCL/TBpajY/iik5xjTWuz9wXZeTrjZAWnuwaoHU3QvXMckydkCn4jU+MeUU4eLyPEuBaatvnv4u1Tx3nadFJr1EkfBiilkE/IuJS1wTU4bqNWWcraugTwfm+x47n05s8CVvUfi0OWf9ZXV87XsvGx+YEMkt94w/6/Rg6hiq14waXFo7Nc0HF0GTPCU8uS13Ibk7CVEwP40umsQ16URnrrv9xmFrNJCvkNbNHd9EMh3Jsc4bnTVUrylrbavIe2iFKeElCkHte4NczR8uWxZV9liWU0aSGr8p00q3OusWCxu2xcQ3OH98Ygo7D0za43bhnjWm+54QIY6u6A5RsKHKUsza9qSA8s3AzJKxafKVAObdVXLeuDyfrouWaUHKelhEkmHHhRlwX2Du3z4CwGxCTV4ra0JM9II6yw1ySc91mWBpAZlHy7NlE8rkFt21XVdkmzlG1y7Abj5xXTQHifDpD7mlR6vaHtqZ54Kv5qdbjPJofiobtA2VwHVNYU7X5yBq2f1tY5EUqK4Lk6nGh2TNypq2zXnHNX2vTHJleMjt2Dv52eQC6zIWPuUDJNEL6jw3yGS7dhnU5A121SZNSK93emsrHwyXmgzp8wNtNTlc25U+fkul3fgsmr71E1wfjCJRmYZ0JUuTEQZpQZcXH1OHrg6Hz/Wy7AJUOwDbwbi+URLJvwX0c1uX56Bqq8RkgnTtF6A+I9Wln77mQZOgdoqjJqJbiOgkEX2XiHY5f7MjriUfVZgKrCdjq3UkvbSmugOuTDeayhuY9mLrkJ7nj3zuH5Qx0Xu+dCKToJtuNJ3qTft4rU1xo6qY7gePvNhXIQ1cCGtLxjvfe9t1ytjwzdeu7SoYH0JIA+6nuZjqcKSvZ4rhzRLDrapDrSoNKsdq38GTmaI+ko47VW1pAFi+bBAfvWmdMrZZF+HjWidd5kOscSiJKzk/ly2rOStWQU1EgwD+EMB7AbwdwIeJ6O3BWgB1sopPrLlpktomWPKmFVn42wcB4OnvvawUrKEEhUSV4amq7+z69xKVYAi/d/NHajnpwwhUR5Ydfv6MV2z/R29a56VwuJzm4pNAYZq/Wea2Lis3PVbybM485qxk9JSutvTh589g36/e0CHI16ysaU0Muvuq23XoFuz0IuC66IUM1XMJz3sngO8KIb4PAET0VwD+DYB/DNYKqJNVfLZqtonoEuJWRCJELOicaaoMT6B7nHTbQ91k9BUMQ4miVKvrNZybneuqIW2y/7tgC4FMP8A7HdO7k+YqVcyzS5t18c0+42gSDL5zm6APWU2P1aa9h4LslGRfbUfL+dh9fT7vGgbrck90Wn5WXAT1MICXEr//CMC70h8ioh0AdgDAunXr0m9741tPw2X1st00XejcyosGcW62eyKm6xzPzM45a7zyb32FjsrW7PI3qhPDXTM8Ja5/D7jH5epsq6YoB58FXDpwQx6mkL5+0g7pOnYqfOKb09hieHVzW7fLEXCPk7cJrnptEARhLYImn2HfhLCQuAh22z0pogyqi6BW7Xy77q8Q4gEADwAtZ2LOdgHoHDSThh0q0Fy3ogJqIZVOG9Z559OHBiT/1uTAVAk1WVND58wYHCDML6i/S9W30NqGRJdNKeOQXaoG6l5Pt0O3QOaN9bYlqLjOu3SbfYpP5R1HXRtcYq1dMQmupKPQliux+DmPE3f6ga59eeLgbVijPojoFwDsFkJsbf9+NwAIIe7T/U1RJ7y4nsJRBK5xrLpIDdPfqkICXR7GLN/VS/Kc1uL7PUU9OKrwtrxlLX3bW+Q4hhg712v4jGWv5k5WimhfrvA8IloG4DsA/jWAKQDfAPBrQogTur/px1FcZSf2iRk7ZRu/mNoboi0x9aes5I6jJqL3AbgfwCCAPxNC/DfT51lQMwzD+JH7zEQhxF8D+OugrWIYhmGciObgAIZhGEYNC2qGYZjIYUHNMAwTOSyoGYZhIqeQ6nlEdAbACxn//HIAPwnYnDLAfV4acJ+XBln7/BYhxFrVG4UI6jwQ0VFdiEpV4T4vDbjPS4Mi+symD4ZhmMhhQc0wDBM5MQrqB/rdgD7AfV4acJ+XBsH7HJ2NmmEYhukkRo2aYRiGScCCmmEYJnKiEdRFH6DbL4joz4joNBE9l3jtUiL6MhH9U/vnmvbrRESfbY/Bt4jo5/vX8uwQ0VVEdJiIvk1EJ4jot9qvV7bfRLSCiL5ORMfbfd7Tfv1qInqm3ecDRHRR+/Xl7d+/235/fT/bnwciGiSiCSJ6ov17pftMRD8komeJaJKIjrZfK3RuRyGoe3GAbh/5cwC3pF7bBeDvhBDXAPi79u9Aq//XtP/tAPBHPWpjaOYA3CWE+FkANwH4RPt+Vrnf5wFsEULcAGAEwC1EdBOA3wewv93nswA+1v78xwCcFUK8DcD+9ufKym8B+Hbi96XQ581CiJFEvHSxc1sI0fd/AH4BwMHE73cDuLvf7QrYv/UAnkv8fhLAm9r/fxOAk+3//zGAD6s+V+Z/AP43gF9ZKv0GsBLAN9E6W/QnAJa1X1+c5wAOAviF9v+XtT9H/W57hr6+uS2YtgB4Aq0T5Kre5x8CuDz1WqFzOwqNGuoDdKt8PMQbhRA/BoD2zyvar1duHNrb21EAz6Di/W6bACYBnAbwZQDfAzAthJhrfyTZr8U+t99/BcBlvW1xEO4H8EkA8uTay1D9PgsATxHRsfah3kDBc9vp4IAe4HSA7hKgUuNARBcDeBTAnUKIV4lU3Wt9VPFa6fothJgHMEJEQwC+COBnVR9r/yx9n4no/QBOCyGOEdG75cuKj1amz202CSFOEdEVAL5MRM8bPhukz7Fo1D8CcFXi9zcDONWntvSC/0dEbwKA9s/T7dcrMw5EVENLSD8ohHis/XLl+w0AQohpAF9Byz4/1D53FOjs12Kf2++vBvByb1uam00APkBEPwTwV2iZP+5HtfsMIcSp9s/TaC3I70TBczsWQf0NANe0vcUXAfh3AB7vc5uK5HEAv9H+/2+gZcOVr//7tqf4JgCvyO1UmaCW6vynAL4thPhM4q3K9puI1rY1aRBRHcB70HKwHQZwR/tj6T7LsbgDwCHRNmKWBSHE3UKINwsh1qP1zB4SQnwEFe4zEa0iokvk/wHcDOA5FD23+22YTxjZ34fWaeffA/A7/W5PwH49BODHAJpora4fQ8su93cA/qn989L2Zwmt6JfvAXgWwMZ+tz9jn38Jre3dtwBMtv+9r8r9BvAOABPtPj8H4Hfbr78VwNcBfBfAIwCWt19f0f79u+3339rvPuTs/7sBPFH1Prf7drz974SUVUXPbU4hZxiGiZxYTB8MwzCMBhbUDMMwkcOCmmEYJnJYUDMMw0QOC2qGYZjIYUHNMAwTOSyoGYZhIuf/A5Umvwx9CALNAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plt.plot(randu,'o')\n", "plt.plot(records,'r')" ] }, { "cell_type": "code", "execution_count": 46, "metadata": {}, "outputs": [], "source": [ "length=500\n", "randu=np.random.pareto(2.5,length)\n", "records=np.zeros(length)\n", "records[0]=randu[0]\n", "k=1\n", "while (k]" ] }, "execution_count": 47, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXAAAAD4CAYAAAD1jb0+AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nO3df5RcZZ3n8fe3K0VS4UcqkSCkIQZ3YnAgmjg9EIzjSBgJCmJvxjHDiIO77GRn17MK62ZMZjyCc3RgTnYH3B2Pu1nF0ZGDUcxkGHA3MIZZz0QFEzohMBAV+ZVOlFbSUZM2/evZP6pup7r63qpbVfepqlv1eZ2T0+lb1VXP/fF87/N8n+fea845REQkfXpaXQAREamPAriISEopgIuIpJQCuIhISimAi4ik1KxmftnZZ5/tlixZ0syvFBFJvb179/7UObewfHlTA/iSJUvYs2dPM79SRCT1zOyFsOVKoYiIpJQCuIhISimAi4iklAK4iEhKKYCLiKRUU2ehiKTNjoFBtuw8yOHhERblc2xcu4z+lb2tLpYIoAAuEmnHwCCbtx9gZGwCgMHhETZvPwCgIC5tQQFcJMKWnQcZGZtg1YtPcPkLT0wtP/rtWbD6whaWTFLp/e+HpUsT/UgFcJEIh4dHALj1H7fy+qHnmcROvfhIiwol6fXmNyuAS/qkNY+8KJ/j8NHjXHj0MP/7N/v51Jp/B0BvPsfuTWtaXDoRzUIRz4I88uDwCI5TeeQdA4OtLlpVG9cuY8nIMHPGR3l+/iIActkMG9cua3HJRAqqtsDN7G7gWuBl59wlxWVbgHcBo8CzwL9xzg37LKikU5BHPucXP+OSnzw7tXz3XXvp/703trBk1fUDS+1pAJ5b0EtvinoP0h2s2jMxzeytwC+BL5UE8KuAXc65cTP7SwDn3EerfVlfX5/Tzay6y4WbHsQBd3/tNtb8KKX7PpOBQ4fg3HNbXRLpUma21znXV768agvcOfctM1tStuyhkl+/C7yn0QJKZ1qUzzE4PMIZoyfYf+5SPnbVfwTgnDNn8/kP/GaLSxfTggUK3tKWkhjE/LfAtqgXzWwDsAFg8eLFCXydpMnGtcvYvP0A2YkJjs05gwPnLSWXzXDTuuWgVIRIQxoK4Gb2Z8A4cE/Ue5xzW4GtUEihNPJ9kj5Bvjj3BcfRnozyyCIJqjuAm9mNFAY3r3TVEunS1fpX9sLZOS5a0ssaTb8TSUxdAdzMrgY+Cvy2c+5EskWSjjQ2BrN02YFIkqrOAzeze4HvAMvM7JCZ3QT8NXAm8LCZ7TOz/+m5nJJ24+OQzba6FCIdJc4slOtDFn/eQ1mkk6kFLpI4XYkpzTE+rgAukjAFcGkOpVBEEqcALs2hFIpI4hTApTnUAhdJnAK4NIda4CKJUwCX5tAgpkjiVKOkOZRCaZq0PkBDaqcALv45BxMTaoE3gR7E3F2UQhH/xscLPxXAvQseoFFqZGyCLTsPtqhE4pMCuPgXBHClULwLHsQcd7mkmwK4+Dc2VvipFrh3i/K5mpZLuimAi39qgTfNxrXLyGUz05bpQcydS00i8U8t8KYJBio1C6U7qEaJfxrEbKr+lb0K2F1CKRTxTykUES8UwMU/pVBEvFAAF//UAhfxQgFc/FMLXMQLBXDxT4OYIl4ogIt/SqGIeKEALv4phSLihQK4+KcWuIgXCuDin3LgIl5UDeBmdreZvWxmT5YsW2BmD5vZD4o/5/stpqSaUigiXsRpgf8NcHXZsk3AN51zS4FvFn8XCacUiogXVQO4c+5bwCtli98NfLH4/y8C/QmXSzrEjoFBNn1lLwAf+NvH2TEw2OISiXSOenPgr3bOHQEo/jwn6o1mtsHM9pjZnqGhoTq/TtIoeLzXsV8UHiZw+Pg4m7cfUBAXSYj3QUzn3FbnXJ9zrm/hwoW+v07aSPB4r1mThUd8TfRk9HgvkQTVG8B/YmbnARR/vpxckaRTBI/x6nGTAExaz7TlItKYegP4/cCNxf/fCPx9MsWRThI8xqvHOQAmigFcj/cSSUacaYT3At8BlpnZITO7CbgDeLuZ/QB4e/F3kWmCx3sFAdyZ6fFeIgmqOjHXOXd9xEtXJlwW6TDBU2GeemYXAOfMy/GRdcv1tBiRhOjKCvGqf2Uv/e+4CO6Dr3/wLbBYwVskKbqUXvybLAxiYtbacoh0GAVw8a+YA6dHh5tIklSjxL+gBa4ALpIo1SjxTykUES8UwMU/pVBEvFCNEv+UQhHxQjVK/FMKRcQLBXDxTykUES9Uo8Q/tcBFvFAAF/+UAxfxQjVK/FMKRcQL1SjxTykUES8UwMU/tcBFvFCNEv+UAxfxQjVK/FMKRcQLBXDxTykUES9Uo8Q/pVBEvFCNEv+UQhHxQgFc/AtSKArgIolSABf/JieVPhHxQLVK/JucVOtbxAMFcPHPObXARTxoqFaZ2S1m9pSZPWlm95rZnKQKJh1ELXARL+oO4GbWC3wI6HPOXQJkgN9PqmDSQZQDF/Gi0Vo1C8iZ2SxgLnC48SJJx1EKRcSLumuVc24Q+K/Ai8AR4Jhz7qHy95nZBjPbY2Z7hoaG6i+ppJdSKCJeNJJCmQ+8G7gQWAScbmY3lL/PObfVOdfnnOtbuHBh/SWV9FILXMSLRmrV7wDPOeeGnHNjwHbgzckUSzqKcuAiXjRSq14EVpnZXDMz4Erg6WSKJR1FKRQRLxrJgT8K3Ac8DhwoftbWhMolnUQpFBEvZjXyx865W4FbEyqLdCqlUES8UK0S/5RCEfFCAVz8UwpFxAvVKvFPKRQRL1SrxD+lUES8UAAX/5RCEfGioVkoIrGoBS4ptWNgkC07D3J4eIRF+Rwb1y6jf2Vvq4s1RQFc/FMOXFJox8Agm7cfYGRsAoDB4RE2bz8A0DZBXLVK/FMKRVJoy86DU8E7MDI2wZadB1tUoplUq8Q/pVAkhQ4Pj9S0vBUUwMU/pVAkhRblczUtbwXVKvFPKRRJoY1rl5HLZqYty2UzbFy7rEUlmkmDmOKfUiiSQsFApWahSHdTC1xSqn9lb1sF7HKqVeKfcuAiXqhWiX9KoYh4oQAu/imFIuKFapX4pxSKiBeqVeKfUigiXiiAi39KoYh4oVol/qkFLuKFArj4pxy4iBeqVeKfUigiXqhWiX9KoYh40VAAN7O8md1nZs+Y2dNmdnlSBZMOohSKiBeN3gvl08D/dc69x8xOA+YmUCbpNEqhiHhRdwA3s7OAtwIfAHDOjQKjyRRLOkrCKZR2f06hSLM00ix6LTAEfMHMBszsc2Z2evmbzGyDme0xsz1DQ0MNfJ2kVoIt8OA5hYPDIzhOPadwx8BgIp8vkiaN1KpZwJuAzzrnVgLHgU3lb3LObXXO9Tnn+hYuXNjA10lqJZgDT8NzCkWapZFadQg45Jx7tPj7fRQCush0CaZQ0vCcQpFmqTuAO+d+DLxkZsHzha4E/iWRUklnSTCFkobnFIo0S6O16j8B95jZE8AK4C8aL5J0nARTKGl4TqFIszQ0jdA5tw/oS6gs0qkSTKGk4TmFIs2iZ2KKfwnPA2/35xSKNIuurhD/dCm9iBcK4OKfLqUX8UK1SvzTpfQiXqhWiX9KoYh4oQAu/imFIuKFapX4pxSKiBeqVeKfUigiXiiAi39KoYh4oVol/imFIuKFapX4pxSKiBcK4OKfWuAiXqhWiX/KgYt4oVol/imFIuKFArj4pxSKiBeqVeKfWuAiXiiAi3/KgYt4oVol/imFIuKFapX4pxSKiBcK4OKfUigiXqhWiX9KoYh4oYcai38pSqHsGBjUE+8lNRTAxb+UpFB2DAyyefsBRsYmABgcHmHz9gMACuLSlhquVWaWMbMBM3sgiQJJB0pJCmXLzoNTwTswMjbBlp0HW1QikcqSaIF/GHgaOCuBz2oadZWbKCUplMPDIzUtF2m1hppFZnY+cA3wuWSK0xxBV3lweATHqa7yjoHBVhetM6WkBb4on6tpuUirNdoCvwv4E+DMqDeY2QZgA8DixYsb/LpkBF3ldz7zz/zRY383tfy0L/bAolR1JNLhxIlUtMA3rl02LQcOkMtm2Lh2WQtLJc2Sxl553QHczK4FXnbO7TWzt0W9zzm3FdgK0NfX5+r9viQFXeKrfvAdLhp6nscuuPjUi/l8i0rVwa6+Gt797laXoqqgsqatEkvj0jqA3UgLfDVwnZm9E5gDnGVmX3bO3ZBM0fxZlM8xODzCaeNjvDTv1dz43j8HoDefY/emNS0unbRS/8retq6w4kelAex2Ph7qDuDOuc3AZoBiC/y/pCF4w6mu8mkTY4zOygLp6SqnsZsn6dKNx1haB7C7ch54cDCe9ZVJRjOz6PV0kCZdEdLazZP06NZjLOiVhy1vZ4lMDXDO/ZNz7tokPqtZ+lf2cumi03nT0nPZvWmNl+Cd9EwXzVMW37r1GNu4dhm5bGbasjT0yruyBT7l5EmYN8/LR/vIqSXZzevGbrJUl9ZUQqPSOoCtAD57tpeP9lERkurmdWs3WapLayohCWkcwG7/qyt88hjAfVwUklQ3r1u7yVJdWlMJ3UoB3FMA91ER+lf2cvu65fTmcxiFaY+3r1tec6uhW7vJUl1Sx5g0h1IongK4r5xaEt28bu4mS3VpTCV0KwVwTwEc2rci6JJxkc6gAO4xgLertI64i8h0CuBdGMChfXsHIhJf9w5iOgejo10bwEUk/bq3BT46WvipAC4x6MInaUfdG8BPniz8VACXKnThk7QrBXAF8Kq6vfWZ1luNSudTAFcAr0itT134JO2rewcxFcBj0WX3elamtK+uC+A7BgZZfccu3n77QwB87/DxFpeovan1mextEYLj78JND7L6jl16kLY0pKtSKKXpgIvHC7NQvvj4EQYHBrsmHVArXXaf3IVPSkd1hnYaE+qqAB6kA84/9hP+/WPbAfily2gwqgJddl+QxIVPGgxNv3Y7CXdVAA+6/f/6yV1c9/S3OHTWOTz7qvO7Jh1QT8tBl90nR+mo9Gu3k3BXBfAgHbDw+DBH55zJW/7D3UDhlpmdrpGWgy67T4bSUenXbifhrhrEDAajzj5+lKHT5wPdkw6IO5tEg2z+6GEJ6dduM5K6qgUetCJ7//bnDJ2R9/Y0+nYUp+XQbvm9TqN0VG3aabAw0G5jQh0VwOPs8P6VvTBrBC69lN2b1rSopM0Xp/uedH6vHStgqykdFU+7Niba7STcMQH8YzsOcM93X8QVf6+4w3/8Yzj33OYWsMXitBySfup9O1ZASYd2Gyws1U4n4boDuJldAHwJOBeYBLY65z6dVMFqsWNgkHu++yKvG3qe+768kdNHf3XqxU+F/IFzcN55TStfO4jTckhykK2dK6C0v1oaE93c02ukBT4OfMQ597iZnQnsNbOHnXP/klDZYtuy8yAO+NDue3HWw/948/qp1wz40JVLp/9BNgvvf39Ty9gOqrUckszvtdtovaRL3MZEt/f06g7gzrkjwJHi/39hZk8DvUDTA3gQFN40+AwPL72MO3/rhqnXevM5PpTiXHczWxdJ5vfSNGWu2S24bm4xxhW3MdHtPb1EcuBmtgRYCTwa8toGYAPA4sWLk/i6GRblc/z4lV9yzvGjHDrrnFPfDameotWK1kVS+b12G62P0uxt3O0txrjiNibizq7q1BNmwwHczM4Avg7c7Jz7efnrzrmtwFaAvr4+V/56LaJ2xMa1y/j0F3aRcZMcOfPsQrmA961anOodlebWRbuN1kdp9jZO8z6NK6mAGacxUa2n1+knzIYCuJllKQTve5xz25MpUrhqO2LB/rMAOHLWwprnd7frGTrteeR2Gq2P0uxtnPZ9Wk2zA2a1nl6nnzAbmYViwOeBp51zf5VckcIFO2LeyC/47ef20uMKjfn93/9/9L/rYt66bw8AX/z4Orjkktif285n6DTlkdOqnm3cyAm/0/dpswNmtZ5ep58wG2mBrwbeDxwws33FZX/qnPtG48WaKdjgN31vBx/6zrbpL361+HPuXFiypKbPbeczdDPzyO3aC/Gt1m0cdr3BLdv2cfO2fbF6fmkZGygX9/hoRcCs1NPr9BNmI7NQ/plCqrkpgh1x5ugJfnFajms/UJhy/uoz5/DVP7688KYFC+CMM2r63HY+Qzcrj9zOvRDfatnGwfUG5QM5sS4eq+P7ktLoybmW4yPJgBlV7lrWJ60nzLhScyVmsCNOmxjjV7Nm88L8ReSyGW5Ztxx+rf6DPz83y9ETYzOWt8sZuhl55Di9kE5uocfdxsH1BpXE6b01c2wgiZNzLb3UpAJmVLn3vPAKX987GHt9mtkIakX9SE0ADzZG9v84RjPZRG5EtWNgkF/+anzG8mzGKh5w7RzM6ilbtV5IK1ro7biN4/bK2qH3FkgiRVhLLzWpgBlV7nsffYkJ52Ysr7Q+vk+YrezBpiaAQ3FjvG4BHJuXyI2otuw8yNjkzDbV+KTjlm372LLz4IyDr53TDbWWLQiSUa3KoBfS7HGCZmzjek4QUemBsPf5OAH5ODnHUWtapDyIB7csrmX9o8pXHryrvb8ZWjmOlr77gZ88OeNJ8vXewzpqpztXyGsGgaP085r5lPZa16uWsgVBMioglXZ7mz1O4Hsbf2zHAW7Zto/B4ZHI/Rwm7H7e5XLZDFdctHBq29by+ZWU7q9aPjOJ+1fXeh/zessap3wZCx92a2XKs5XjaKkP4I0cLHF2enngaNbOqme9ailbWJAM9OZz3L5u+VTrIYkgUMvJyOc2jhqIjHOC6F/Zy+3rltObz2EUttMNqxZP+/32dct55JmhxE9A1U5qUds3iYdIhK136fFRa1njiCr39Zdd0HYPxahWP3w+JCVVKRRgRgBvpPsSNuASpjRwNGtaUi3rFTcVUioqGBrMSE81OjAVlhK5eds+PvEPT3Hruy72OpOhXKXtFOcEEZVPDfbBLdv2RX5+nPRLlEontTgpp0bTObXkkZM4AVcqd99rFrTV+MgVFy3ky999MXS573Rg6gN4I/dCKD1IKlWu0sDhY1pSWPmi1muwWGFLZ4dUOglZ8W9W37Fr2oFeS5BsNAhEtfaPnhgLPZh9Tv2qFETqPUFU2wcBK7436ZuD+crB1pvHT+oEHHXSaLcrfB95ZihyeaXeWBLrkPoUSpzuS6VURHAvlUq5zSsuWjj1/1q7k9VElS8/Nxv5N6Xlr5QKMWbOUa63a92/spfdm9bw3B3XsHvTmkQGpCC8a530Ni4VtV0bufFZpX1QyhXfG6ila11pf8VpndeaYmwkNdltz/6stP19p1xT3wJP4l4I1Spg+Rk2iRZA0LoJa6mMjE0we1YPuWwmtFyl5a90IETlect7H767otVmb0RNRyu9aKN8VlDc1mHp++blsvz8VzPn/ENjNz6rpTKWTs3c+LX9U7OgBodH2Pi1/UDt85mjjqNF+Ryf+Ien6moB1tqqL98fv/sbvTzyzFDbpDlK61vGjAnnEnsmbrUeh8+Ua+oDeBL3QqhWAX0NUFY6aRwbGePO9Su4edu+0NeDMsWd2lb+dzDzRBS0CIPteMVFC2dUQqg96Fcba4iadgc0dDFH+XYeHgkP3vlclk/2L6+4DpXUsg+Cinvb/U/NmMI6Num47f6nap7PHNWIicrNQv3HfNjysNsLfH3vYMVeU5IBNerYKT1xHx8dZ2yiUMJgKmJS+ehqjUifV4KmPoBD5YGlnuLBUW5e7lRXuloF7DHjwk0Phgasai3BsAM1E1GmUovyuaqtK4g+eOZke0KvMI1al7DBltLKP9VCNKYqQtwKELx22/1PzQiipdPuygNyj9HQxRxxUxvHIgJ7XHEHw0srbtTJZHhkrObcc1QjptKsj2otwFqeiFNpVk9UvSzdXo0E1LDjtvw4jdrW1cpZ7XvDehxBPQ8+d+PaZdy+brm3nq65KoEkSX19fW7Pnj2NfUhvL7zjHfC5z00titN6K5fNGFve88ap7nicCgiFShi0LMpbHqWvQ3jAiqP0O8LKVvp6pfUv7aJX+57Vd+yqe5ZEbz4X+8Kqek9ocRjw3B3XTP1+4aYHq176DvHKX8uJOqps+blZhk+MVW0wlKfO4uzvsIBQaf3vWr+i6vETdtyVp0ZOjI6HNhSCdS7dH4Fqx1ot+6ORmT2BqHJW+u6obVPaMwzMn5sNnW1VUxnN9jrn+mYsT10AP/tsWL8ePvMZYGb3DSq3QEuVHii1BJagqxc1ZSyfy3JyfDLWCSHqs2tp5YfZMTDIxvv2T7VCKn3f7k1rYge7KAY1lS3uCTNM1H7J57Lsu/Wqqd/jnpRuWLW4YgolrLzZHuOMObOmAnKw3o1uxx6DqHNubzGtVR4kygN8IGr9S7dTpQYCTG/Vh313JVGBuNo2igqopXW0dIC+UbU2QD7y1f2hx1+lhkjUPoorKoCnOoVSqfsW5yCrlA+GwoEW9XeV5hM32uouV8+g6ZadB6sGb6g/l16udKYCVO4Chw2s1eL6yy5g22MvzehdHB8dnzZNL25q44H9RyoOuIWlYsYm3VQDofSWstZAZMlmrOI+K09rBWq9sdRt11089Xulwcry2Uar79gVe79FzeqplNYMhKV3yk80SQXveq5niCp7pXXydWl9qgN4nLvDVdJIHjCJgc2MGZPOxZ5FUUv+LG75KuXS6zEyNsHNEfeRgcL6VOsZAZEtrGDA8cEnjsz4nLEJx+btT0y1kDJmrHrtfL79o1eo1NEcHhmbOukGOdQ/3f4EJ8Ymq5YzEHx8rR3a4BiYl8tiRqxtEybquod5uSxzsj0zegphfxf1edXeWy7qcYbVAiBEB9S44xlQ6B1NAhMlJ/geK4x7HT0xVvegabUyVEsF+ri0Pl0BfHISxsamAnilDZLL9jBSoQIGg2fBrIugApUe6JVGl6vlO6vV4zhdqkau4orTojZOzXEPGwgbPjHK8dH6AnpYWYPuZ5TSE1pUqiBoPQ5HBLrSfT7hHLuffaXmso9NuopjB0madI47169o+OQZ9QzI4ZGxqYAaliaq5aKbOMdUkAJ65JmhGRctxQnCc7I9U+tRehKK26udPzfLNW84j22PvUTpN2XMGs5DV4o3Bqx67Xwef/FYxdlWSUtXAB8dLfycPbtiV8yAOdlMZAAPyyWWHiBBl9hRaPFFtWDCKl2l4B0c3MGZH5h2Ahkdn5hq9c0vXnAS1r39yFf3c8u2fZFT/YKLk6oFBQfcU+ySf7J/+YxUTVQKKa6grIFqra//9t43Tvv+SpdMN5ryaRf5udmaWpdRKl33EOznvtcsmDG2cvzkzNspR7WC4xxTwXkv7AQepwV69MRYTbNIyjlH6CylYIpm1PFUqadb7VYVUNjGj794jN/9jV4e2H8kdLaVjwuZ0jWIeewY5PMc+M8f571zLw89kILWRlhuPHj9uTuuqWnWRVQLpnynVxyRL0b2uLNk6pXtMdZfesHUlKa45mZ7pk4ec7M9zM5m6u7Ol4szqBwMrIVVpD0vvDJVKTNmXH/ZBfS9ZoGX7RdXpsemddEryWUzTExOMhqS387nshwbGWtsALnk2Kq0z/O5LKfPnhU6NzpQ3sgob7F+bMeB0AAZR6UB2lYIGlthja7SgdxajrNgu5XOQOvqWSg7BganNsarjg+z969v4ONv/2O+9KZrQ98fzCioFKB762y9zS+ZBhZ2cDc6A6Hb5SOCSpgbVi2e0UI/cmykaQEiCIZRMyJOPy3DidGJqWMlasaSQU0pgmqSnJ2RxJTYblQ+GN3oDBSIDuBtfS+U4HJj+9lPuWD4xyw5ehiAk5no+4Q8sP8IKz7xUMUAXW/X++iJsanZFhu/tn/afSGClI7Ub3hkLFbwhkJKoLyl/geXLY79XT0N7qpjI2NsXLuM3nwuNGAeH50gPzc7daKvdM+eJA+bJM9fI2MTU7NrFLzjKz+Gg4H9pG8lC23eAg9a0X/+0Gf5w4FT+dgPXvdRHnz9b/koYs168zmWvCrHt599Ra3vLmLArCrT/uBUSissLwqFXl1SqSppf/W2xlM5DzwY9Nh+yRr2LXodACczp/GPSy9rZbGmGRwe6YjBNKmNY2ZLK8zYpIu8HwnUP22wkqSubJXkJT0fvK0DeDAos2/RMvYt6sxbUYokTcG7vSU5H7yhHLiZXW1mB83sh2a2KalCBTr1/sEi0r2SnA9edwA3swzwGeAdwK8D15vZrydVMCjMH71hVfyBKRGRdpdkw7SRFvilwA+dcz9yzo0CXwHenUyxTvlk/3LuWr8i6Y8VT3o9XG0mp+Rz0TOwpP3lsj2J3g+lkQDeC7xU8vuh4rJpzGyDme0xsz1DQ+HPjqumf2Wv98Awf26Wdp8E2Gj5wv6+3s/M57KRj82qtq968zlW/6sFbb+9201vPse+W69qeAqkJKeWfZHtMW5f94Zkv7+Bvw0r+ozRE+fcVudcn3Oub+HChSF/Es/GtcsqFnZutqeuA9soXBQy8PGreO6Oa9q2BZnLZrhz/QruWr+CbB0rWvr3pc+afN+qxRWfBxr1Wbddd3Hkcys3rl0WWsZsxrhr/Qp2b1rDPX90OXeuXxGrRZntMWbPautLFsgUJ3MHP6utV7bHmJuNv06ll2LXMt+9FrWWSeeRwr6IU3/mz82y5ffemGjrGxqbhXIIuKDk9/OBw40VJ9qp+488MXWPkx4rbMDgEvfSqzYrCa5WC7tcuNL9HrI9helZYVf7BfdXiboPQvkN+csvtb3mDefNuHlTpXKGXaoL0x8jVX5zruDvyw+i0isaa73EutLzG6tdTlz63Muo/VZ635iw+5uflrGqU/qiLpkuv6d33PtdV9ovpQo3lnpixj15Sv8u7Fgo35fl+y843itNTyx36oEDh0LvEVS+f6L2Sem613p/cChs89L7nES+r/jAFWjsKtDgjpTP/2yk7vuIh5W59PYa5VcER92fyIe6L+Qxs1nA94ErgUHge8AfOOeeivqbRB7oEFP53cyiglmlvw17Xh9UfyZkvbeArffvktYu5YhbrvLlURUoznpVe75iO20PiD7OKx3ztezfuE8iitr2cZ+rGraslu9Jso614zHg5V4oZvZO4C4gA9ztnPtUpfc3M4CLiHQKL1diOue+AXyjkc8QEZH6tPfIkIiIRFIAF585708AAAQPSURBVBFJKQVwEZGUUgAXEUmppt4P3MyGgBfq/POzgZ8mWJw00Dp3B61zd2hknV/jnJtxJWRTA3gjzGxP2DSaTqZ17g5a5+7gY52VQhERSSkFcBGRlEpTAN/a6gK0gNa5O2idu0Pi65yaHLiIiEyXpha4iIiUUAAXEUmpVARw3w9PbhUzu9vMXjazJ0uWLTCzh83sB8Wf84vLzcz+e3EbPGFmb2pdyetjZheY2SNm9rSZPWVmHy4u79h1BjCzOWb2mJntL673J4rLLzSzR4vrvc3MTisun138/YfF15e0svz1MrOMmQ2Y2QPF3zt6fQHM7HkzO2Bm+8xsT3GZt+O77QN4Mx6e3EJ/A1xdtmwT8E3n3FLgm8XfobD+S4v/NgCfbVIZkzQOfMQ593pgFfDB4r7s5HUGOAmscc69EVgBXG1mq4C/BO4srvdR4Kbi+28Cjjrnfg24s/i+NPow8HTJ752+voErnHMrSuZ8+zu+nXNt/Q+4HNhZ8vtmYHOry5Xg+i0Bniz5/SBwXvH/5wEHi///X8D1Ye9L6z/g74G3d9k6zwUeBy6jcFXerOLyqeMc2AlcXvz/rOL7rNVlr3E9zy8GqzXAAxQeYtOx61uy3s8DZ5ct83Z8t30LnJgPT+4gr3bOHQEo/jynuLyjtkOxm7wSeJQuWOdiOmEf8DLwMPAsMOycGy++pXTdpta7+Pox4FXNLXHD7gL+BAie3/YqOnt9Aw54yMz2mtmG4jJvx3dDD3RoklgPT+4CHbMdzOwM4OvAzc65n5tFPh63Y9bZOTcBrDCzPPB3wOvD3lb8mer1NrNrgZedc3vN7G3B4pC3dsT6llntnDtsZucAD5vZMxXe2/B6p6EF3tSHJ7eBn5jZeQDFny8Xl3fEdjCzLIXgfY9zbntxcUevcynn3DDwTxTGAPLFZ8vC9HWbWu/i6/OAV5pb0oasBq4zs+eBr1BIo9xF567vFOfc4eLPlymcqC/F4/GdhgD+PWBpcQT7NOD3gftbXCaf7gduLP7/Rgp54mD5HxZHrlcBx4JuWVpYoan9eeBp59xflbzUsesMYGYLiy1vzCwH/A6Fwb1HgPcU31a+3sH2eA+wyxWTpGngnNvsnDvfObeEQn3d5Zx7Hx26vgEzO93Mzgz+D1wFPInP47vVSf+YAwPvBL5PIW/4Z60uT4LrdS9wBBijcDa+iULu75vAD4o/FxTfaxRm4zwLHAD6Wl3+Otb3LRS6iE8A+4r/3tnJ61xcjzcAA8X1fhL4eHH5a4HHgB8CXwNmF5fPKf7+w+Lrr231OjSw7m8DHuiG9S2u3/7iv6eCWOXz+Nal9CIiKZWGFIqIiIRQABcRSSkFcBGRlFIAFxFJKQVwEZGUUgAXEUkpBXARkZT6/1oVgaFJDMPiAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plt.plot(randu,'o')\n", "plt.plot(records,'r')" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "from statsmodels.distributions.empirical_distribution import ECDF" ] }, { "cell_type": "code", "execution_count": 17, "metadata": {}, "outputs": [], "source": [ "def gumb(x):\n", " return 1-np.exp(-np.exp(-x))" ] }, { "cell_type": "code", "execution_count": 20, "metadata": {}, "outputs": [], "source": [ "def fre(x):\n", " return 1-np.exp(-x**(-alpha))" ] }, { "cell_type": "code", "execution_count": 25, "metadata": {}, "outputs": [], "source": [ "def weib(x):\n", " return 1-np.exp(-(-x))" ] }, { "cell_type": "code", "execution_count": 16, "metadata": {}, "outputs": [], "source": [ "length=500\n", "maxlist=np.zeros(200)\n", "\n", "k=0\n", "while (k<200):\n", " randu=np.random.exponential(1,length)\n", " maxlist[k]=max(randu)\n", " k+=1\n", "\n", "maxlist=maxlist-np.log(length)" ] }, { "cell_type": "code", "execution_count": 18, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAD4CAYAAADvsV2wAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nO3de3RU5b3/8fc3QBASREXqOopcWls9qD+rREVFICQoIEKltgpTxQsiWtpKz+9XUarWtohSD/VSQk/qDTFeUBRECWKCSvVYFS9VudiFKJFIRVAwCZKY5Pn9sRMMYSaZkJnZe2Y+r7VmDbNnZ+bLLPjMzrOf/X3MOYeIiKS+DL8LEBGRxFDgi4ikCQW+iEiaUOCLiKQJBb6ISJro6HcBLTn00ENd3759/S5DRCRpvPnmm9uccz3DPRfowO/bty+rV6/2uwwRkaRhZpsiPachHRGRNKHAFxFJEwp8EZE0ocAXEUkTCnwRkTSRsMA3sywzm29mfzOzUKLed4+iIujbFzIyvPuiooSXICLip3YFvpndZ2Zbzez9ZttHmNkHZrbBzKY3bB4HPOGcuwIY0573bbOiIpg8GTZtAue8+5/9DA49tPXgv/pq6NgRzLz7/Px9vzj0ZSIiScDa0x7ZzAYDlcCDzrnjGrZ1AP4FDAc2A28A44GxQLFz7h0ze9g5N6G118/JyXH7Mw//rrvuYvv27ZgZGRkZZMyZQ8bOnXTAu/CgI9Cp4ZaZmUnnK67ggPx8OnfuTNeuXcnKyiI7O5vsW27hoAULyAIs0pt16gT19VBXF7kgM5gyBQoK2vx3ERFpCzN70zmXE/a59vbDN7O+wDNNAv804HfOubMbHl/XsOtm4Evn3DNm9qhz7sIIrzcZmAzQu3fvAZs2RbyGIKL+/fuzbt26Nv9cJB2Bg4AewGHAdxrujwR6N9z6AofTwhfDnhfr6H059O4NM2dCKPGjWyKSuhId+OcDI5xzkxoeXwScClwL/AXYDbzsnGt13GN/j/AbOedwzlHXrx/1ZWXUAbUNt28abjVANVD99tt8/fXXfP3111RVVVFZWUnFhAnsAHYAXwLbgc+ArcC/G7Y3lQV8H/gBcDxwQsPtSFr5IsjIgCuv1G8AItJuLQV+PForhMs255yrAi6Nw/tFLsTMG9a55RZvDH/XrvA79ukDP/zhvtsvuqjFoZoK4BNgE/AR3jjWv/DGsBY22a8HMBA4DTi94c9dmr5QfT3Mm+fdAHr0gDvv1NG/iMRUPGbpbMY7qG3UC/g0Du8TvVAICgu9IG2ua1dvaCWcyZNbfNluQH9gJHA1cAewDNgI7AReBubinaH+EPgtMAw4GMgDZgFvAvv8jrV9u3dSuUsXnQAWkZiJR+C/AXzfzPqZWSZwIfB0W17AzM41s8KdO3fGrqpQCLZtg4ce8o7ozbz7wsLIR9IFBXDVVdChg/e4QwfIy9v756+6CjIz9/nRA4Ez8L4I7gPW4Q0JPQtMbfjz9UAO3jfiVcByvCGmPXbv9oK/WzcFv4i0W3tn6TwCDAUOxRvevsk5d6+ZjcI74O0A3Oeci3AI3bL2juEnTFERzJgBZWXeydijjoLS0lZ/bCteyC8BngOq8I7+zwcmAGfifYB7ZGdDVZVO+IpIRHE9aRtPSRP4LTn2WFi7ttXddgMlwKPAYrzw74V30uNyoE+kH7zqKp3sFZE9Wgp8tVaItzVrvIu9Gm9XXRV2twOA0cBDeEf+jwLHAX8E+gEj8MbF9jmFPG+eN7x09dVx+guISKpQ4CdaQYEX/A89FP4kMtAVuAAoxpv9cwPwPt6Va0cDd+HNENrLvHneHH+N9YtIBIEM/LictA2axpPIjUf+Dz3kHak30we4GS/4H8O74OtXeMM9M4DPm+5cV+ed5DXz5vbrqF9EmtAYftDk57d6wvd14E/AIrz5/JOB/4d3pe8+Dj8cystjXKSIBJXG8JNJScm3Y/1hjvgBTgEeB9YAPwHuBr6HF/rbmu/86afe65hF1yxORFKWAj+oCgq8K3CbDvlkZe21y38CD+Bd3ftTYA7wXeB3eB3t9tF4QZdCXyQtKfCTRSgElZVhZ/p8F5gPvAechTfm/33gXsLM6gEv9EUk7QQy8NPipG17NM70ycvba3N/4AngVbypnJOAk4CV4V7j2GPjXKSIBE0gA985t9Q5N7l79+5+lxJsjeP9zY76BwKv4M3q+Qqvb88FeE2O9li7VrN4RNJMIANf9kPjUf/h3lwdwxvXX4s3pv80cAwwG68tNODN3Vfoi6QNBX6qKS/fq+FbF+AmvODPw1uUYADwWuP+8+apK6dImlDgp6KCAqit9Wb2NOiH16RtCfAFXm/+X9JwxW5jV06FvkhKC2Tg66RtjIRC+5zYHYN3tH813vJjx+E1bQNg4sREViciCRbIwNdJ2xgqKdkzrt/oQLywfxmvadtwYApQUVen2TsiKSyQgS8xVl4O/fvvs/l04B3gv4BCvHV4X1q7Fg4+OLH1iUhCKPDTxZo1YVszdwFuxzva7wTkAtfu2EH1sGGJrU9E4k6Bn04KCvY6kdvU6cDbwBV4UzdPfeEF1tx2WwKLE5F4U+Cnm1Ao4iIs2cD/4M3Z3wLkTJ/O/wwZQpA7qopI9BT46ahxcfYIzgXeBQYDU1at4ifdu/Pll18mqjoRiZNABr6mZSZA45W5ERyGt+LWbGBJRQU/PPxw/vGPfySqOhGJg0AGvqZlJlCzefpNZeD12P9foMPu3Zw5aBB//vOfNcQjkqQCGfiSQGHm6Td3MvAWMLqujl//+teMGzeOHTt2JKQ8EYkdBb548/RbONIHOAh4Em+RlWeWLGHAgAH885//TER1IhIjCnzxlJREnLLZyIBpwCrn2F1WxmmnncZDrfyMiASHAl++FQpBdnaru50GvFVbyyl1dVx00UVMnTqVmpqa+NcnIu2iwJe9/fWvERdPb+owoKSmhv/q0IG5c+eSl5fHZ599Fv/6RGS/KfBlb6EQLFgAPXq0umtH4Pa6Oh4B3nz9dQYMGMDq1avjXqKI7J9ABr7m4fssFIJt27wx/SiO9i8E/remho7l5QwaOJAFCxbEv0YRabNABr7m4QdE49F+w+pZLfkhsBo4va6Oiy++mGuvvZa6urq4lygi0Qtk4EuAhELe6lmtzNUHOBR4DrgKmD17Nj/60Y+oqKiId4UiEiUFvkQnirn64LVYLgDmAsXLlnHaaafx0Ucfxbs6EYmCAl+iV1ISVeiDt4TicxkZlJeXc+qpp/Lqq6/GtzYRaZUCX9qmpMTrtBnFuH5ebS3/uP56DjzwQHJzc3nkkUcSUKCIRKLAl7YrKPDG9Z2Djh1b3PXo3/yG16qrObVfPyZMmMDvf/97NV8T8YkCX9rngQda3aXH5s2sWL+eizt35qabbuLSSy/VlbkiPlDgS/uEQlGN63cGHqiu5mZg/vz5jBw5Uh03RRJMgS/tV1IS1W4G3Ag8CPx95UrOOO44Nm3aFM/KRKQJBb7ERgtLJjZ3EbAC+LS8nIEnnMA777wTt7JE5FuBDHy1VkhCBQVRT9kEGAq8DHTauZMzzzyTFStWxKsyEWkQyMBXa4Uk1dhTv08f73ErfXiOBV4FvltZyTlnn838K6+Me4ki6SyQgS9JLBSCjz/2pmzW17d61H8EsAoYAlxSWMitAwdq2qZInCjwJb6iuDq3O7AMGA9c99prXDNiBPX19YmoTiStKPAl/hqvzm1BJvAQcA1w14oVjB8/nurq6kRUJ5I2FPiSGFGc1M3AWyT9T8DChQs555xz1G1TJIYU+JI4UfThMeD/AvOvvJIXX3yR3NxcPv/884SVKJLKFPiSWI19eB56qMXdLv7b31hyzTWsXbuWQYMG8fHHHyemPpEUpsAXf4RCLT9fX885//3fPN+pE1s3b+aMM85gzZo1ialNJEUp8MU/USyUfsZXX7Gqpob6qioGDx7M66+/noDCRFKTAl/8c+edUe12fG0tr+zcyUFffcWwwYMpibJ3j4jsTYEv/gmFou7B813g5dpa+lVXc87ZZ/Pkk0/GtzaRFKTAF38VFHgncKMY3vkPvKtyB9TX85Pzz2f+/PlxL08klSjwxX+hEGzb5gV/Rsv/JA/G67Q5zIxLLrmEu+++OyEliqQCBb4ERygEDz7Y6m7ZwNL6en7Urx+//OUv+eMf/6j+OyJRCGTgqz1yGmsc12+l0+YBwOMffcTFgwZxww03MH36dIW+SCsCGfhqj5zmCgpgwYJv2yxH0BG4/5VXuDo/n9mzZ/Pzn/9cTddEWhDIwBfZ02a5lStyM5zjLyUl/OaHP2TevHlceuml1NbWJqZGkSSjwJdgi2LqpgG3vvMOfzjlFB588EEuvPBCampqElOfSBJR4EvwFRRAdnaLuxjw29dfZw6waNEifjxwILt3705IeSLJQoEvyeGvf4XMzFZ3mwbMA555+23OPflkqqqq4l6aSLJQ4EtyCIXgvvuiukBrCjAfWPn++4wcOZKvvvoq7uWJJAMFviSPxgu0omjHcDHwCPDq3//OWT17sqOwMO7liQSdAl+ST0FBVHP1fwo8AbxVU0PelVey/bLLElKeSFAp8CU5RTlXfyywBFgD5N5/P58VFCSiOpFAUuBL8mqcq+9ci2P7I4FngQ+B3F/+ki1btiSoQJFgUeBLarjzTujaNeLTeUAxUFZXx9ChQykvL09YaSJBocCX1BAKQWFhi0M8g4HngC1btjBkyBDKysoSVp5IECjwJXU0HeKJ4Azg+eefZ9u2bQwZMkSLo0taUeBL2jl1wwZKSkrYuXMnQ4YMYePGjX6XJJIQCnxJPxMnkvPBB5SWllJZWcmQIUPYsGGD31WJxJ0CX1JTS1fk1tXBRRdx4kknsbJzZ3Y3HOl/8MEHiatPxAcKfElNd97Z8vMN4/wnbNnCC998Q21lJbm5uaxfvz4BxYn4Q4EvqSnKlbMAjtu9mxeysqirq1PoS0pT4Evqarwat0OHVnftv2ULL7zwAs45hg4dyrp16xJQoEhiKfAltYVCMH9+ixdlNer/9tu8+OKLmBlDhw5l7dq1CShQJHEU+JL6orgoC4BLL+WYQYN44d//JmP7dnJPO02hLylFgS/pofGirJZC/5tvYPt2jgFerKsj46uvFPqSUhIW+Gb2XTO718yeSNR7iuxj5syohneOBl4EOlRUkJuby5o1a+JdmUjcRRX4ZnafmW01s/ebbR9hZh+Y2QYzm97SazjnNjrnLm9PsSLt1ji8E8WJ3KOBF5yjQ00Nw4YN05G+JL1oj/AfAEY03WBmHYC5eN1n+wPjzay/mR1vZs80u30nplWLtEcbTuQeDbwAdNi+ndxjj2Xt4YdDUVHcSxSJh6gC3zm3Cvii2eZTgA0NR+41wKPAWOfce8650c1uW6MtyMwmm9lqM1v9+eefR/0XEWmTpidyzVq8MvfoHTt4oa6ODGDYli2smzRJoS9JqT1j+EcAnzR5vLlhW1hm1sPM/gqcaGbXRdrPOVfonMtxzuX07NmzHeWJtKLxRG59vbdWbkuhj3ekb0Du7t2s/81vElSkSOy0J/DDXcIYsS+tc267c26Kc+57zrlZ7XhfkfhoZRGVY/BCHyD300/Ve0eSTnsCfzNwZJPHvYBP21eOiI+aD/P06QMZe/8XOQZYCdQDubm5/Otf//KjUpH90p7AfwP4vpn1M7NM4ELg6VgUZWbnmlnhzp07Y/FyItFrOszTeN9Mf7zQr62tJTc3V62VJWlEOy3zEeBV4Ggz22xmlzvnaoGpeKvGrQMWOudiMlnZObfUOTe5e/fusXg5kf0XYfrmsWaUlpZSU1NDbm4uH374YYILE2m7aGfpjHfO/YdzrpNzrpdz7t6G7cuccz9oGJefGd9SRXxQVxd+u3Mc/+67lJSUsGvXLnJzc/noo48SW5tIG6m1gkhLWmrFMGMGJ5xwAiUlJVQ29NPftGlT4moTaaNABr7G8CUwZrbwi2tZGQAnnngizz//PDt37iQ3N5dPPvkk8s+I+CiQga8xfAmMUAiyssI/d8ghe/44YMAAVqxYwfbt28nNzaW8vDxBBYpEL5CBLxIoBxwQfvuXX+51xe3JJ5/Mc889x9atWxk2bBhbtmxJUIEi0VHgi7Tmi+ZdRRrU18PkyXuF/sCBAykuLqa8vJxhw4bx2WefJahIkdYp8EVa07t35Od27YIZM/badMYZZ7Bs2TLKysrIy8tDPaEkKAIZ+DppK4HSWg/9MDNzBg8ezNKlS/nwww8ZPnw4X0T6LUEkgQIZ+DppK4HSWg99s7DdM4cNG8aSJUtYv349w4cPZ8eOHXEuVKRlgQx8kcBp7KFvYXoGOrfPsE6js846iyeffJL33nuPs88+G/3WKn5S4ItEKxTywj2csjLvKL9vX6/hWt++e476R40axeOPP85bb73FqFGjqKioSFjJIk0p8EXaItKVt4cc4s3Y2bTJ+1LYtGmvGTxjx47l0Ucf5bXXXmP06NFUVVUlsGgRjwJfpC3CncBtfLxr197bm83g+fGPf8yCBQt4+eWXGTt2LF9//XWcixXZWyADX7N0JLDC9cwvLIw8V3/Tpr1O6I4fP57777+flStXMm7cOKqrqxNUuAiYizQmGQA5OTlu9erVfpch0rq+fcNOzwS83wAKC70viwb33nsvkyZNYsyYMTz++ONkZmYmpk5JeWb2pnMuJ9xzgTzCF0k6Lc3VD3Nx1uWXX87cuXN5+umnmTBhArW1tQkoUtKdAl8kFhqHeiJp6KzZ1NVXX82cOXNYtGgRF198MXWReu+LxIgCXyRWQqHIs3gyMvaZrgkwbdo0Zs2axSOPPMKkSZOoD7OkokisdPS7AJGUMnOmNx2z+YydxqP3xumasGdMf/r06VRXV/O73/2Ozp07M2/ePCzcBV4i7aTAF4mlxhOzM2Z4wzgZGfsuk9g4pt/kJO6NN97I7t27ufXWWznggAP485//rNCXmAtk4JvZucC5Rx11lN+liLRdKPRtmGdEGDVtNqZvZtxyyy3s3r2bO+64g86dO3Prrbcq9CWmAjmGr+ZpkjIitVYOs90efpg5Tz7JVcDs2bO5+cc/jm9tknYCGfgiKSPSlbnN18otKoLJk7GyMv4CXArc/NRT3HrBBYmqVNKAAl8kniJdmdtk/B7wxvQbTvRmAH8DJgDXLVzIHXfckeiqJUUFcgxfJKU0HdOPpNmYfgdgPlCNN3WzS5cuXHnllfGqUNKEjvBFgiDMmH5H4OFDDmH06NFMmTKF+fPnJ74uSSkKfJEgmDkTOnXaZ3NmZSWPjxvH8OHDueyyy3jsscd8KE5ShQJfJAhCITjwwH2319RwwM03s3jxYgYNGkQoFGLx4sWJr09SQiADX+2RJS1FarFcVkbXrl155plnOPnkk/npT39KcXFxYmuTlBDIwNc8fElLrczZ79atG8XFxRx//PGMGzeOlStXJrA4SQWBDHyRtBRuzn6nTlBZuafx2kHPPsuKFSs46qijGDNmDK+88oo/tUpSUuCLBEXzOfs9enj327fvtU5uj+XLKSkp4YgjjmDUqFG88cYbflcuSUKBLxIkoRB8/DHU10N2NtTU7P18Q+O1ww47jNLSUnr06MHZZ5/Nu+++60u5klwU+CJBFWbRlKbbe/XqRWlpKVlZWeTn57N+/foEFifJSIEvElRRNF7r168fpaWlZGRkkJeXx4cffpig4iQZKfBFgirKxms/+MEPKCkpobq6mry8PMoi/WYgaU+BLxJU4RqvTZzoNVprtlzicccdx4oVK/jyyy/Jy8tjy5Yt/tYugaTAFwmypidxZ86E+fO92TpNZu00hv5JJ51EcXExW7ZsIT8/n23btvlbuwSOAl8kWTRpobxH43KJDU4//XSWLl3Kxo0bOeuss9ixY0eCi5QgC2Tgq7WCSBitzNpplJuby1NPPcX777/PyJEjqaioSEBxkgwCGfhqrSASRhuWSxwxYgSPPfYYb7zxBmPGjOHrr7+Oc3GSDAIZ+CISRrTLJTY477zzWLBgAS+99BLjxo2juro6AUVKkCnwRZJFtMslNjF+/Hjuueceli9fzoUXXsg333yTwIIlaBT4IskkFPKO6Hv39sbuZ8zYM0snkssuu4y7776bxYsXc8kll1BXV5egYiVotKatSDIpKvKmYjbO1mmcmgktHulPnTqVqqoqpk+fTpcuXSgsLCQjQ8d76UaBL5JMWpqa2cpC6ddeey1VVVX84Q9/oGvXrtx5552YWRyLlaBR4IskkyinZkZy8803U1VVxZw5c8jKymLWrFkxLE6CToEvkkx69/aGccJtj4KZcfvtt1NVVcWtt95KdnY2M5pcuCWpTYEvkkxmztx7DB9anJoZjplRUFDArl27+O1vf0tWVhbXXHNNHIqVoFHgiySTxnH6GTO8YZzevb2wb2X8vrmMjAzuu+8+du3axbRp08jKyuKKK66IQ8ESJOac87uGiHJyctzq1av9LkMkZdXU1HDeeedRXFzMggULCLXxi0OCx8zedM7lhHtO87JEklVRkdciuVmr5LbIzMzkiSeeYOjQoUycOJGnnnoq5mVKcCjwRZJR43z8CK2S26JLly48/fTTnHLKKVxwwQUsX748DgVLECjwRZJRFK2S2yI7O5tly5Zx7LHHct555/HSSy/FoEgJmkAGvtoji7SinfPxwznooINYsWIF/fr1Y/To0bz22mv7/VoSTIEMfLVHFmlFG1olt0XPnj0pKSnhO9/5DiNGjOCf//xnu15PgiWQgS8irWhjq+S2OPzwwyktLSU7O5vhw4ezfv36dr+mBIMCXyQZ7Uer5Lbo27cvpaWlZGRkkJeXx8aNG2PyuuIvzcMXkYjee+89hg4dSvfu3Vm1ahW9evXyuyRphebhi8h+Of7441m+fDnbtm0jPz+frVu3+l2StIMCX0RadPLJJ/Pss89SVlbG8OHD+eKLL/wuSfaTAl9EWnXmmWeyePFi1q9fz6hRo6ioqPC7JNkPCnyRZBaD9grROuuss1i4cCGrV6/m3HPPZVfzC78k8BT4Iskqhu0VojV27FgWLFjAqlWrOP/886mpqYnbe0nsKfBFklWM2ytEa/z48RQWFlJcXMyECROora2N6/tJ7CjwRZJVHNorRGvSpEnccccdLFq0iMsuu4z6+vq4v6e0nxZAEUlW7VzusL1+9atfUVFRwQ033EB2djZz587VougBp8AXSVYxWO6wvWbMmEFlZSW33XYb2dnZ3HbbbQr9AFPgiySrGC132B5mxqxZs6isrORPf/oT2dnZ3HjjjQl7f2kbBb5IMguFEhrw4ZgZd911F5WVldx0001069aNadOm+VqThKfAF5F2y8jI4J577qGyspJf//rXZGVlMXnyZL/LkmY0S0dEYqJjx448/PDDjBw5kilTplAUx+sBZP8o8EUkZjIzM1m0aBGDBw9m4sSJLF682O+SpAkFvojEVJcuXVi6dCk5OTlccMEFPP/8836XJA0U+CISc926daO4uJhjjjmGsWPH8vLLL/tdkqDAF5E4Ofjgg1mxYgVHHnkk55xzDlrMyH8KfBGJm8MOO4zS0lIOOeQQRowYwZo1a/wuKa0lLPDN7Edm9jczW2JmZyXqfUVSWgLbI++vXr16UVJSQmZmJvn5+WzYsMHvktJWVIFvZveZ2VYze7/Z9hFm9oGZbTCz6S29hnNusXPuCuAS4IL9rlhEPD60R95f3/ve9ygpKaG2tpa8vDw++eQTv0tKS9Ee4T8AjGi6wcw6AHOBkUB/YLyZ9Tez483smWa37zT50d82/JyItIdP7ZH3V//+/XnuuefYsWMHeXl5fPbZZ36XlHaiCnzn3Cqg+UKWpwAbnHMbnXM1wKPAWOfce8650c1uW81zG1DsnHsr0nuZ2WQzW21mqz///PP9/XuJpD4f2yPvr5NOOoni4mLKy8u1Pq4P2jOGfwTQ9PeyzQ3bIvkFkA+cb2ZTIu3knCt0zuU453J69uzZjvJEUlykNsgJao+8v04//XSWLFnCBx98wMiRI7U+bgK1J/DD9UB1kXZ2zt3lnBvgnJvinPtrO95XRMDrjNm1697bEtweeX/l5+fz+OOP8+abb2p93ARqT+BvBo5s8rgX8Gn7yhGRqIVCUFgIffqAmXdfWOh798xojRkzRuvjJlh7Av8N4Ptm1s/MMoELgadjUZSZnWtmhTt37ozFy4mkrlAIPv4Y6uu9+yQJ+0ZN18cNhUJaHzfOop2W+QjwKnC0mW02s8udc7XAVOA5YB2w0DkXk6sqnHNLnXOTu3fvHouXE5EAmzRpEnPmzOGJJ57giiuu0Pq4cRRVP3zn3PgI25cBy2JakYiknWnTplFRUcFNN91EdnY2d911l5ZKjAMtgCIigXDDDTdQUVHB7bffTrdu3bjlllv8LinlBDLwzexc4NyjjjrK71JEJEHMjNmzZ1NRUcGsWbPo1q0b1113nd9lpZRABr5zbimwNCcn5wq/axGRxDEzCgoKqKys5Prrr6dbt25MnTrV77JSRiADX0TSV0ZGBg888ABVVVX84he/oFu3bkycONHvslKC2iOLSOB07NiRRx99lOHDh7Nw4UKci3hNp7SBjvBFJJA6d+7M4sWL6dChg2bsxEggj/B14ZWIAHTt2pXOnTv7XUbKCGTg68IrEZHYC2Tgi4hI7CnwRUTShAJfRCRNKPBFRNJEIANfs3RERGIvkIGvWToiIrEXyMAXEZHYU+CLiKQJBb5Iqigqgr59ISPDuy8q8rsiCRj10hFJBUVFMHky7NrlPd60yXsMSbfOrcSPjvBFUsGMGd+GfaNdu7ztIg0CGfialinSRmVlbdsuaSmQga9pmSJt1Lt327ZLWgpk4ItIG82cCV277r2ta1dvu0gDBb5IKgiFoLAQ+vQBM+++sFAnbGUvmqUjkipCIQW8tEhH+CIiaUKBLyKSJhT4IiJpIpCBr3n4IiKxF8jA1zx8EZHYM+ec3zVEZGafA5v288cPBbbFsJxUoc8lMn024elziSyIn00f51zPcE8EOvDbw8xWO+dy/K4jaPS5RKbPJjx9LpEl22cTyCEdERGJPQW+iEiaSOXAL/S7gIDS5xKZPpvw9LlEllSfTcqO4YuIyL5BUzAAAAKMSURBVN5S+QhfRESaUOCLiKSJlA58M/uTma03s3fN7CkzO8jvmvxkZiPM7AMz22Bm0/2uJwjM7Egze8HM1pnZGjP7ld81BY2ZdTCzt83sGb9rCQozO8jMnmjIl3VmdprfNUUjpQMfeB44zjn3f4B/Adf5XI9vzKwDMBcYCfQHxptZf3+rCoRa4L+cc/8JDAR+rs9lH78C1vldRMDcCSx3zh0DnECSfD4pHfjOuRXOudqGh/8AevlZj89OATY45zY652qAR4GxPtfkO+fcFufcWw1/rsD7j3uEv1UFh5n1As4B7vG7lqAwswOBwcC9AM65GufcDn+rik5KB34zlwHFfhfhoyOAT5o83oyCbS9m1hc4EXjN30oC5Q7gN0C934UEyHeBz4H7G4a67jGzLL+LikbSB76ZlZjZ+2FuY5vsMwPvV/ci/yr1nYXZpjm5DcwsG1gEXOOc+8rveoLAzEYDW51zb/pdS8B0BE4C5jnnTgSqgKQ4J5b0Sxw65/Jbet7MJgKjgTyX3hcdbAaObPK4F/CpT7UEipl1wgv7Iufck37XEyBnAGPMbBRwAHCgmT3knPuZz3X5bTOw2TnX+JvgEyRJ4Cf9EX5LzGwEcC0wxjm3y+96fPYG8H0z62dmmcCFwNM+1+Q7MzO8sdh1zrk5ftcTJM6565xzvZxzffH+vaxU2INz7t/AJ2Z2dMOmPGCtjyVFLemP8FvxF6Az8Lz3/5p/OOem+FuSP5xztWY2FXgO6ADc55xb43NZQXAGcBHwnpm907DteufcMh9rkuD7BVDUcPC0EbjU53qiotYKIiJpIqWHdERE5FsKfBGRNKHAFxFJEwp8EZE0ocAXEUkTCnwRkTShwBcRSRP/HxrGusAR5aQFAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "ecdf_max = ECDF(maxlist)\n", "plt.plot(ecdf_max.x, np.ones(len(ecdf_max.y))-ecdf_max.y, 'ro')\n", "\n", "t=np.arange(-2, 5.0, 0.1)\n", "plt.plot(t,gumb(t),'k')\n", "\n", "plt.yscale('log')" ] }, { "cell_type": "code", "execution_count": 19, "metadata": {}, "outputs": [], "source": [ "length=5000\n", "maxlist=np.zeros(200)\n", "alpha=2.5\n", "\n", "k=0\n", "while (k<200):\n", " randu=np.random.pareto(alpha,length)\n", " maxlist[k]=max(randu)\n", " k+=1\n", "\n", "maxlist=(maxlist+1)/length**(1/alpha)" ] }, { "cell_type": "code", "execution_count": 23, "metadata": { "scrolled": true }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAD8CAYAAAB0IB+mAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nO3de3yU1b3v8c9KAkICRQVEC2VijVajVkUEEeOVIigQtFYrQ2uVbbgogqW0emirZ2+px4oVEBIIStUm1VYOWyWRi9IKiBhEfW0RORaKTblVBCkUMAUm6/yxMmEYJvfMzDMz3/frNa+QZ55nZtE633lYl98y1lpERCT5pcW7ASIiEhsKfBGRFKHAFxFJEQp8EZEUocAXEUkRCnwRkRSREe8G1KdLly42Ozs73s0QEUkY77///m5rbddIz3k68LOzs1m3bl28myEikjCMMZV1PacuHRGRFKHAFxFJEQp8EZEUocAXEUkRMQt8Y0yWMeZ5Y8w8Y4w/am9UWgrZ2ZCW5n6WlkbtrUREEkmLAt8YM98Ys8sY83HY8UHGmE+NMZuNMQ/WHL4FWGCtvQcY1pL3rVNpKRQUQGUlWOt+FhS44+PGQUYGGON+jhvX8Gt16eLONwbS091PfYmISIJq6R3+c8Cg0APGmHRgNjAYyAXuMMbkAj2ArTWnBVr4vpFNmQKHDh1/7NAhGD2ad4uKeD0QYCXwQSDAX4qK2HHnnezfv59AIKw5paVw112wZ8+xY9XV7mdlJYwceeyLoEuX478Awr8o9GUhIh5hWloP3xiTDZRZay+o+b0f8Ii19oaa3x+qOXUbsNdaW2aMecla+/06Xq8AKADo2bPnpZWVdU4pPVFamruzjyAfeK2eS9u1a0eHDh3cY/t2Ohw5QgegA5BV8zP4OBk4C8gBzgTatW0L8+e7F7rrLjhypHFtHT0aCgsb+ZcTEWmYMeZ9a23viM9FIfBvBQZZa/+j5vcfAH2BnwGzgCrgbWttg7e6vXv3tk1aeJWd7e7AI/gbsAs4EP548kkOHDjAwYMHOXDggHu88MKJ59U8wv79gAG+AeS0a0dORgZnHzhADu7L4JtAZkNtTk+H558Hf/SGNUQkddQX+NFYaWsiHLPW2oPAXVF4v2OmTnV99qHdOpmZUFVFdnU12eHnp6fDj3984uusWFHnF0c18CWwBdgc8thUVcVCYHfY+d2h9gsg+LgM8AVPCARcF9HIkce3q6BAd/8i0qqiMUtnG+6mN6gHsCMK73Mivx+Ki8Hnc33mPp/7ffToyOcXFEQ+PnUqtGkT8ak0oAvQBxgB/BJ4AVjj8/GFz8deYB3wEvAo8B3cgEU58BDwPSAbF/wFNed9Hv4mgQAUFbm/Q8eOmnEkIq0iGl06GcBfgOuB7cB7wAhr7YYmvOZQYGhOTs49mzZtalH7ao0b58I/EGjcHXRpKUyYcGzgNi3t2MBtuEb24f8L9z/MauBPwFvAvprnzgeuq3lcgxsnqNP118Obb9Z3hoikqKj14RtjXsTlUxfcjerD1tpnjTE3AtOBdGC+tXZqc16/yX34sRD+RdC5M8yYcawPPvz5egSAD3Dh/ydgFfAV7l8RvTj2BXAlbuA4IoW/iISI6qBtNHky8Jsq9F8WDfg3sBZYjvsCeBc4ArTBhf4PgFuBjpEuVvCLCPUHvkorRFthIRw96qaLlpS47qQ6nATkAY8AK4G9wFLgAdzAyN1AN2AksIywxQzLlx+b9z9gQDT+JiKS4DwZ+MaYocaY4n379jV8ciLx+134l5S4AWVwAV2HLGAg8DjwKfAOcCduAPgGoCdurusn4ReGhn96esOrikUkJahLxyua0PdfBZQBzwOLcXf6lwI/BO4AIm51c/LJsHdvqzVXRLxJXTqJwO+H3btd10/wMXZsxC6gdri+/EW4+a7TcesDJgBfx60qLqs5Vuuf/1R5B5EUp8D3stD+f2vdwGyY03BB/wHwEa6//z1gKHAhrtjR4dALgrWAzj8/yo0XEa9R4CeSN988vv8/zIXAr4FKoBQ3u+cuXL2fJzg25x+ATz5xawt0ty+SMjwZ+Ek7aNsa/H7429+OzfrJOnGGfhvcKuAPcbN8zgN+yrFB3tplz9a6u33N6hFJCZ4MfGvtImttQadOneLdFG/z++HAgWN3/ca4hWA1DG6Wz5u4cg+DgWm40g6jgI3BE4OzesJLPYtIUvFk4EsTBe/6q6vdwG9JyQmnXIqr27MJV8PnRdxmBfm4fwkAbobQyJGayimSpBT4ycjvr7O755u4GtWVwMPA27gyDncAfw2eVF3tire1b687fpEkosBPVsHunmBff5iuuBW9W4ApuM1hzgXuBf4RPKmqyt3xp6Xpjl8kCXgy8DVo28r8fhf8ubknPNUJV8Z5M3APUIzbzevnhMzqsdbd8bdpozt+kQTmycDXoG2UbNhQZz2fM4BC3EDuMGAqLvh/g1vZC7g1ASNHKvRFEpQnA1+iKFjPZ+zYiE/n4AZ038cN9E4CzsGVcahduaupnCIJSYGfqgoL653L3ws3h385cDrwI6A/bnon4KZytmmj3bhEEogCP9UFB3fHjo1YufM6XF3+54HPcFs7FlCzd2+w7EOwXIMGdkU8TYEvTmGhm44ZoW5/Gq4S56fARGA+rptnNnA09DWKihT6Ih7mycDXLJ04C/bzf/3rxx3uhBvE/QjX5XMf0Bu3NWMthb6IZ3ky8DVLxyO2b49YoTMXeAN4GfgSuAq3C9fO4AlFRerTF/EgTwa+eEiwQmcYg6vJ//9wc/YX4Iq0FVMzm0czeUQ8R4EvDQuWaggpzBaUCfwXsB7XzTMauAb3RcDy5dC9e+zaKSL1UuBL44TvyBW2avds3BTO+cDHwEW4L4LDO3a42T8ZGdpxSyTOFPjSPBs2nBD6BrfhykbgFuCXuLv+NQCBgDspOIVThdlEYk6BL80XIfQBuuFW65YB+3ELtu6r+XOtYGE2zegRiRkFvrTMhg11lmm4CdgAjMfV6bkQ1+1znKIi6NhRd/siMeDJwNc8/AQTLNMQNm8foCMwA3gHaA8MAMYBB0JPOnDA3e2rj18kqjwZ+JqHn6C2b3d3+xGqcV6O21nrx8Ac4NvAikivEezj112/SKvzZOBLAissPFZjJ7jXbo32wJO4oE/DTd+cCByK9DrBu3718Yu0GgW+RE9wr11rj1uxmwf8D65vfwZuCufqul6jqEibq4u0EgW+xEZwxW5NKeYsYCbwZ1wBtjzgfwGHI10b3Fxdd/siLaLAl9gJlmIOCf5rcMXY7gYeA/rh5vFHVFTkxgcU/CLNosCX2AsNfp+PjsAzwH8DlbjFWrMBG+na6mpN5RRpJgW+xE9oH//YsQzH1eS5BrdQ60ZCKnCGC53KqT5+kUZR4Is3FBbC2LGcAbyOu8N/C7dY678bunbPHrjzToW+SAMU+OIdhYVQUoLp3JlxuHn72bi6PGOoY/pmUCCggV2RBngy8LXSNoUFq3KWlHBu5868A/wUmAtchhvgrVdREZx/frRbKZKQPBn4WmkrweBvW1LC4xkZLAX24DZRr3NAN+iTT7T5ikgEngx8kVp+Pzz3HAN9Pj4CrsMN6A4Hdtd33fLlqs0jEkaBL95XM5vnNGspCwR46oorWIJboftWQ9dWVsJddyn0RVDgS4JJS0tj4urVvPvoo3QwhuuB/wQC9V105AhMmBCbBop4mAJfEtIlU6bw/v79jOjUiYeBgdQzZx/c1E2RFKfAl4TVoUMHXti7l/nnncca4GLgjfou0ECupDgFviQ0Ywx3ffIJ626/na7ADcAUXEG2Eyxfrr58SWkKfEkKuS+9xNpnn2UU8CvcbJ4dkU4cORIyMrRAS1KSAl+SRubddzOvpITfAe/junjejHRiIOAWaCn0JcUo8CW5+P2MvP563gO64gZz/zd1zOIpKtI8fUkpCnxJPm++SW67dqwFRgKPAIOBXZHO1Tx9SSEKfElOzzxDVtu2PI+rtb8KuIQ6tlI8cgR++EOFviQ9Bb4kJ78f5s/H+HyMAt4FMnG19p8iQi2e6mo3oKva+pLEPBn4qpYprSK4wUpJCRdlZLAOGAL8GLgN2B/pmuD+udpRS5KQJwNf1TKlVdUUYOvk87EQ+DVuU5XLgI/ruia4o5Zm8kgS8WTgi7S6mrt9U1LC5LZtWQ7sA/oC9d7Ha/qmJBEFvqSWmr79qzt35kPgUtxMnvuBw3Vdo9CXJKHAl9RTs7nKGSUlLD/1VB4AngaupY7VuQBz5qhPXxKeAl9Sl99Pmz17+M3YsbwE/A/QC1gZ6VxrYcqUmDZPpLUp8EUKC7m9pISK9u3phKvDE3HqZmVlzJsm0poU+CIAfj/nHzrEe6NGMQw3ddMPHAw/T335ksAU+CIhvvbMM/xfXMXNl4ArgL+GnlBUpL58SVgKfJEwxufjIWAxsBXoXfPnWtouURKUAl8k3NSpYAw34MosZwM3AY8C1aDtEiVhKfBFwvn9MGYMAGfiCq75gV8AtwL/AnXrSEJS4ItEUlhY+8dM4AVgOvAabnXuX1RdUxKQAl+kLp071/7RABNwm6R/AVxWXU3ZqFFu1k52NqSlaTMV8TwFvkhdZsw44dC1uH79HGDov//NfxUVUV1Z6RZmaTMV8TgFvkhd/H4YOxaMOe5wT+Bt4AfAL4HvUdOvD24zFc3iEY9S4IvUp7AQfvc7SE8/7nB74HncitxXgcuBzcEnNYtHPEqBL9IQvx+efx4yM487bICJwFLgH7j6+kti3zqRRlPgizSG3w/FxeDzuS6e4E/gemAd4MPN138CsCUl8WurSB0U+CKNFdwysbra/bTHyqsF5+t/F/gpMHLUKL7q2VOzd8RTYhb4xphvGmOeNcYsiNV7ikSVz3fcr1nAH3B1eF48fJgrt27l78HZOwUFCn2Ju0YFvjFmvjFmlzHm47Djg4wxnxpjNhtjHqzvNay1W6y1o1rSWBFPqSnBEMoAD+EWaG3C9eu/DXDokOrpS9w19g7/OWBQ6AFjTDowGxgM5AJ3GGNyjTEXGmPKwh6ntWqrRbwgWIIhLPQBhgAVUFtffx6onr7EXaMC31q7Evgy7HAfYHPNnfthXDXZfGvtemvtkLDHrlZut4g3BKdthg7mprmP1Xm40L8OKADuM4YjR47EsbGS6lrSh98dVz02aFvNsYiMMZ2NMXOAS4wxD9VzXoExZp0xZt0XX3zRguaJxEj4YG51de1TpwDlwE+A2dZyw8UXs3v37vi0U1JeSwL/xH/HRtgVrvYJa/dYa8dYa8+y1j5Wz3nF1tre1treXbt2bUHzROIkbJFWOm6q5gvAO598Qp/cXNavXx+PlkmKa0ngbwO+EfJ7D2BHy5ojkgQCgYiHf4DbIL1q9276XXQRrxijKZsSUy0J/PeAs40xZxpj2gLfx01OaDFjzFBjTPG+ffta4+VEYitsumaoPsA6a8m1lpuBRysrsffco9CXmGjstMwXgTXAt4wx24wxo6y1R4H7cCvLNwJ/tNZuaI1GWWsXWWsLOnXq1BovJxJbU6eeUIYh1Ndxd/o/wG2q8v2vvuLQD3+oRVoSdRmNOclae0cdx18HXm/VFokkOr/f/Zwwoc5Cau1wxdcuBH4GbK6u5hXgG8FFWqGvI9JKVFpBJBr8fti9G0pKjp+yGbapymRgEccWaa0BLdKSqPFk4KsPX5JG+JTNGTNO6O65CXgX6ABcg5vNw9//Htt2SkrwZOCrD1+SVnjVzZpVurm4RVpXAncCPz3pJAJ1zPYRaS5PBr5IUgu96z/11NrDnXH19McCT1RVkZ+fz/79+49dV1qq/XOlRRT4IvH05fEVS9oAhbgiVUuWLKFfv35s2bLFhXtBgavHowqc0kwKfJF46tkz4uFxPh/Lli1j586d9OnThxWTJrnB3FAa3JUm8mTga9BWUkakOfuZmTB1Ktdddx1r166la9euDPj8c1dxM5wGd6UJPBn4GrSVlBFp68Ti4to5+Dk5OaxZs4YBbdpQAEwAjoZeHzIGINIQTwa+SEoJn7oZtuDq5JNPZlHHjkwEZuJq7f8z+GRVlQZypdEU+CIJIGPvXp7CbaSyHOgHbAY4eFADudJoCnyRRFAzuPsfwJvAF0Bf4M/h52kgV+rhycDXoK1ImJDB3atxi7S6AQOB4vBzKyt1ly8ReTLwNWgrEiZscPcsn481p5zCAGA0MJGwwVx17UgEngx8EYkgbHC309NPs6h9eyYAM4BhQO2/idW1IxEo8EUSld9Pxrx5TAfmAG8AVwBbgs9rjr6EUeCLJDK/H3w+RuN2ItqJ21VrFbg5+pqyKSEU+CKJrmZA9zrcYG4X4HrguX/+U1M25TgKfJFEFzKge7YxrOnRg6vbtOGuQICfArVFltWvn/I8GfialinSRCEDuqds3crrR44wDngCuAU4EDyvslJdPCnMk4GvaZkiLdPG52M28DRQhttYpXYIV108KcuTgS8iLVTTr38fUA58hhvMrQg9R108KUeBL5KMQvr1BxnDGiATt2fuH0LP09TNlKLAF0lWIf36uT4fFUBv4PvAfwIWIm/Aoq0Uk5YCXyQVTJ1K18xM3sRtkv4w4E9Pp+rhh48/T1spJjUFvkgqqOniOcnn47fA/zn5ZF4MBLi2uJjPP//82HlTpmgrxSTmycDXtEyRKKjp4jHW8rO9e1m4cCEfffQRffr0Yf369e6cuvr01defFDwZ+JqWKRJ9N998M6tWreLo0aNcccUVlJeX17mpep3HJaF4MvBFJDZ69erF2rVrOeeccxg2bBgz8vKw7dufeOKBA+rHTwIKfJEU1717d1auXEl+fj4TS0oYd8UVHAnfHH3PHg3eJgEFvoiQlZXFggULePDBB5mzfDk3Hjp0bKP0IA3eJjwFvogAkJaWxmOPPcZvf/tbVlRV0Q/4a/hJGrxNaAp8ETnOj370I97o1o1duI3S3w59UoO3CU2BLyInuPrJJ3m3XTs642rrvwBgjFuIpdW3CUuBLyIn8vs5+5lnWNOjB/1xq3N/bi3VoNW3CUyBLyKR+f2cunUrS3r2ZBQwFVeH5yvQAG6C8mTga6WtiHe03bqVebjNVBbgKm7+AzSAm4A8GfhaaSviIT17YoCfAAuBj3GDuetPPz2uzZKm82Tgi4iH1GymAjAcWAUcNYb+e/eyePHiuDZNmkaBLyL1C9lMBWPo5fOxdsYMcs47jyFDhvD000/Hu4XSSAp8EWlYyGYqTJ1K9yefZOWHHzLkpJO4//77GT9+PEePHo13K6UBCnwRabyQDVI6AAu/+opJGRnMmjWLYcOGsX///ni3UOqhwBeRxgvbICUdmHb0KHNPPZVly5bRv39/Kisr49c+qZcCX0Qar46pmAV797JkyRK2bt1K3759qaioiHHDpDEU+CLSePVskDJgwADWrFlDZmYm11xzDS+//HJs2yYNUuCLSOOFTNGslZnpjgPnnXceFRUV9OrVi9tuu41f/epXWGvj0FCJRIEvIo0XNkUTn8/97vfXntK1a1eWL1/OiBEjmDJlCnfffTeHDx+OY6MlKCPeDRCRBBMM9ylTXJ9+sKZOSOi3a9eOkpISzjnnHB555BE+++wzFi5cyKnhO2lJTOkOX0SaJmRqJtbWWT3TGMPDDz9MaWkpa9as4fLLL2fTpk1xarSAAl9EmipsaiZQb/XMESNG8Kc//Ym9e/dy+eWXs3Llyhg0UiLxZOCrWqaIh9VVJbOe6pn9+/fn3Xff5bTTTmPAgAG88MILUWqc1MeTga9qmSIeVs/UzPqcddZZvPPOO+Tl5XHnnXfyi1/8gurq6ig0UOriycAXEQ9rYGpmfU455RSWLFnCqFGjePTRRxkxYgRVVVVRaqiEU+CLSNM0Ympmfdq0acO8efN4/PHH+cMf/sC1117Lrl27otxoATBeXhTRu3dvu27dung3Q0SiZOHChYwcOZJu3bpRXl5Obm5uvJuU8Iwx71tre0d6Tnf4IhI3t9xyCytWrKCqqop+/frxxhtvxLtJSU2BLyJxddlll1FRUYHP52Pw4MHMmzcv3k1KWgp8EYme0lLIzoa0NPczbHFWUM+ePXn77bf5zne+Q0FBAZMnT9YMnihQ4ItIdDRyRW7Q1772NRYtWsS4ceOYNm0at956KwcPHoxxo5ObAl9EoqOJK3IBMmp2z5o+fTqvvPIK11xzDTt37oxyQ1OHAl9EoqMZK3LB1eCZMGECr776Khs3bqRv376sX78+Cg1MPQp8EYmOZq7IDRo6dCirVq0iEAjQv39/Fi9e3IqNS00KfBGJjhasyA265JJLWLt2LTk5OQwZMoTCwsJWbmRqUeCLSHS0cEVuUPfu3Vm5ciU33XQT9957Lw888ACBQCBKjU5uWmkrIgkhEAgwadIkZsyYwdChQ/n9739Phw4d4t0sz9FKWxFJeOnp6UyfPp3Zs2dTXl7OVVddxfbt2+PdrISiwBeRhDJu3DjKysrYtGkTffr04cMPP4x3kxKGAl9EEs7gwYNZvXo1aWlp5OXlUVZWFu8mJQQFvogkpG9/+9usXbuWc889l/z8fGbOnBnvJnmeAl9EEtYZZ5zBihUrGDZsGBMmTGD8+PEcPXo03s3yLAW+iCS0rKwsFixYwKRJk5g1axb5+fn861//inezPClmgW+MGW6MmWeMedUYMzBW7ysiCaqRlTbBzeCZNm0aRUVFLF26lLy8PLZt2xazpiaKRgW+MWa+MWaXMebjsOODjDGfGmM2G2MerO81rLWvWGvvAX4E3N7sFotI8mtipc2gMWPGUF5ezpYtW+jTpw8ffPBBjBqcGBp7h/8cMCj0gDEmHZgNDAZygTuMMbnGmAuNMWVhj9NCLv15zXUiIpE1o9Jm0A033MDq1atp06YNeXl5vPbaa1FqZOJpVOBba1cCX4Yd7gNsttZusdYeBl4C8q216621Q8Ieu4zzOLDYWlvn164xpsAYs84Ys+6LL75o7t9LRBJZMyttBl144YVUVFRw/vnnM3z4cGbMmIGXqwrESkv68LsDW0N+31ZzrC7jgQHArcaYMXWdZK0tttb2ttb27tq1awuaJyIJq4WVNgFOP/103nrrLYYPH87EiRM1g4eWBb6JcKzOr1Br7Uxr7aXW2jHW2jkteF8RSXatUGnTXZLJggUL+MlPfsLs2bMZNmxYSs/gaUngbwO+EfJ7D2BHy5ojIkKrVdoESEtL44knnmDu3LksW7aMK6+8kq1btzZ8YRJqdLVMY0w2UGatvaDm9wzgL8D1wHbgPWCEtXZDixtlzFBgaE5Ozj2bNm1q6cuJiACwbNkyvve975GVlcWiRYu49NJL492kVtfiapnGmBeBNcC3jDHbjDGjrLVHgfuApcBG4I+tEfYA1tpF1tqCTp06tcbLiYgAMHDgQFavXk3btm256qqrUm4GT2Nn6dxhrT3DWtvGWtvDWvtszfHXrbXnWGvPstY2rXNNRCQOLrjgAioqKrjgggsYPnw406dPT5kZPCqtICIpp1u3bvz5z3/m5ptv5oEHHkiZGTyeDHxjzFBjTPG+ffvi3RQRSVKZmZm8/PLLTJ48OWVm8Hgy8NWHLyKxkJaWxq9//euUmcHjycAXEYmlgoICysvL+eyzz+jbt2/S1uBR4IuI4GrwvPPOO0ldg0eBLyJSIziDJ1iDJ9lm8Hgy8DVoKyLxElqDJ9lm8Hgy8DVoKyLxFF6DJ1l20fJk4IuIxFuwBs+cOXOSZhctBb6ISD1Gjx5du4tWos/gUeCLiDQguItWRkYGeXl5LFq0KN5NahZPBr4GbUXEa4K7aOXm5pKfn8+MGTPi3aQm82Tga9BWRLwoOIMnPz8/IXfR8mTgi4h4VVZWFgsWLGDSpEnMmjWL4cOHJ8wMHgW+iCSu0lLIzoa0NPeztDQmb5uens60adMoKipiyZIlCTODR4EvIomptBQKCqCyEqx1PwsKYhb6AGPGjKGsrKx2Bs+HH34Ys/duDgW+iCSmKVPg0KHjjx065I7H0KBBg1i9ejXp6enk5eVRXl4e0/dvCk8GvmbpiEiD/v73ph2PouAMnnPPPZdhw4bx9NNPx7wNjeHJwNcsHRFpUM+eTTseZWeccQYrVqxgyJAh3H///UyYMIFAIBCXttTFk4EvItKgqVMhM/P4Y5mZ7nicZGVlsXDhQiZOnMjMmTO5+eabOXDgQNzaE06BLyKJye+H4mLw+cAY97O42B2Po/T0dJ566ilmzZpFeXk5V111FTt27Ihrm4KMl2s99+7d265bty7ezRARaZbXX3+d22+/nZNPPpmysjIuuuiiqL+nMeZ9a23vSM/pDl9EJEpuvPFGVq1ahbWWK6+8ksWLF8e1PQp8EZEouvjii6moqCAnJ4chQ4ZQWFgYt7Yo8EVEoqx79+6sWrWKG2+8kXvvvZdJkybFZQaPJwNf8/BFJNl06NCBV155hfHjx/Ob3/yG7373uxw8eDCmbfBk4Gsevogko/T0dGbOnMnMmTNZtGgRV199NTt37ozZ+3sy8EVEktn48eN59dVX2bhxI3379mX9+vUxeV8FvohIHAwZMoRVq1YRCATo378/S5cujfp7KvBFROKkV69eVFRUcOaZZ3LTTTcxd+7cqL6fAl9EJI569OjB22+/zcCBAxkzZgyTJ0+muro6Ku+lwBcRibOOHTvy2muvMXbsWKZNm8Ztt93GofDSz60go9VfUUREmiwjI4PZs2eTk5PD3LlzOXToEJnhxeFaSLV0REQ8pqqqinbt2jXrWtXSERFJIM0N+4Z4MvC10lZEpPV5MvC10lZEpPV5MvBFRKT1KfBFRFKEAl9EJEUo8EVEUoQCX0QkRXh64ZUx5gugshVfshPQGnM9W/o6zbm+qdc05fwuwO4mtifZtdZ/K9EU6zZG8/302TxRcz+XPmtt14jPWGtT5gEUe+F1mnN9U69pyvnAunj/f+O1R2v9t5JMbYzm++mzGfG8Vv9cplqXziKPvE5zrm/qNa31d01VifC/X6zbGM3302czBjzdpSOxYYxZZ+uovSEi8RGNz2Wq3eFLZMXxboCInKDVP5e6wxcRSRG6wxcRSREKfBGRFKHAFxFJEQp8OY4xJssY87wxZp4xxh/v9oiIY4z5pjHmWWPMgkptfV8AAAFfSURBVOa+hgI/BRhj5htjdhljPg47PsgY86kxZrMx5sGaw7cAC6y19wDDYt5YkRTSlM+mtXaLtXZUS95PgZ8angMGhR4wxqQDs4HBQC5whzEmF+gBbK05LRDDNoqkoudo/GezxRT4KcBauxL4MuxwH2BzzV3DYeAlIB/Yhgt90H8fIlHVxM9mi+kDnbq6c+xOHlzQdwcWAt81xhSRGOUFRJJNxM+mMaazMWYOcIkx5qHmvHBGa7ROEpKJcMxaaw8Cd8W6MSJSq67P5h5gTEteWHf4qWsb8I2Q33sAO+LUFhE5JmqfTQV+6noPONsYc6Yxpi3wfeC1OLdJRKL42VTgpwBjzIvAGuBbxphtxphR1tqjwH3AUmAj8Edr7YZ4tlMk1cT6s6niaSIiKUJ3+CIiKUKBLyKSIhT4IiIpQoEvIpIiFPgiIilCgS8ikiIU+CIiKUKBLyKSIhT4IiIp4v8DMc8IIQLZEZsAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "ecdf_max = ECDF(maxlist)\n", "plt.plot(ecdf_max.x, np.ones(len(ecdf_max.y))-ecdf_max.y, 'ro')\n", "\n", "t=np.arange(0.5, 10.0, 0.1)\n", "plt.plot(t,fre(t),'k')\n", "\n", "plt.xscale('log')\n", "plt.yscale('log')" ] }, { "cell_type": "code", "execution_count": 24, "metadata": {}, "outputs": [], "source": [ "length=500\n", "maxlist=np.zeros(200)\n", "\n", "k=0\n", "while (k<200):\n", " randu=np.random.random_sample(length)\n", " maxlist[k]=max(randu)\n", " k+=1\n", "\n", "maxlist=(maxlist-1)*length" ] }, { "cell_type": "code", "execution_count": 28, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "[]" ] }, "execution_count": 28, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXQAAAD4CAYAAAD8Zh1EAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nO3de3xU9Z3/8dcnNy5ykxBEBRLaghatWhtZ13b703q3FuturdrghVK5WNSqvdDi1gvlp6X1p3QxICpegGqtv211hRZ7c7tu1RJqdYuKohJAVEIotwRIQj77x5ngZDKTmSRzz/v5eMwjOXO+mfkeLu/55nO+53vM3RERkdxXkOkOiIhIcijQRUTyhAJdRCRPKNBFRPKEAl1EJE8UZeqNhw0b5hUVFZl6exGRnLRmzZpt7l4WbV/GAr2iooKamppMvb2ISE4ys9pY+1RyERHJEwp0EZE8oUAXEckTCnQRkTyhQBcRyRNxA93MlpjZVjP7W4z9ZmY/MbP1ZvaKmZ2Y/G6KiEg8iYzQHwLO6WT/ucDY0GMqsLDn3eolli+HigooKAi+Ll+e6R6JSA6LOw/d3f9oZhWdNLkAeMSDdXhfMLMhZna4u7+XpD6289xzz/HMM89QUFBAYWHhwUfkdvhzRUVF7Z5v2w7/2tmjuLi4w/fFxcUHH4WFhZhZ1w5k+XKYOhUaG4Pt2tpgG6CqKvbPzJ4NGzfC6NEwd27stiLS6yTjwqIjgU1h25tDz3UIdDObSjCKZ/To0d16s+eff545c+Z062dTqS3cS0pKDn6N/L7t0adPH0qee44+e/dSAvRpezQ20ufqq+nz2mv07du3/WPNGvotWULfpib6Af1qa+n3ta/Rb8sW+l9yCf3796dfv3707duXgoI0nhrRh4xI1rBEbnARGqE/7e7HRtm3Arjd3Z8Lbf8O+La7r+nsNSsrK70nV4q6OwcOHGj3aG1tjfr9gQMHaGlpibrd0tLS4fuWlhaam5s5cOAAzc3N7Z5r2277vrNHU1PTwcf+/ftpbm5m//79wXPPP08TsA/YH/HYZ0ZPbjzSr18/DjnkEPr373/w64ABA6I+Bg4cyMCBAxkwYACDBg1q9xg8eDCDBg2ib9++0X8DifwtA6B/f1i8uPuhrg8IkU6Z2Rp3r4y2Lxkj9M3AqLDtkcCWJLxup8zsYBkkJ1VUBGWWSOXl+Dvv0NLSwr59+9i3bx/79+9n76hR7CX4ANgb+bj/fhobG9m7dy8NDQ00NjbS2NhIQ0NDu8d7771HQ0MDe/bsYc+ePezevZvW1ta4XS0pKWHIkCEMHjyYIUOGMGTIEA499FAOXbmSQxsbGQocCpQCpY2NlH7rWww74wyGDh1KcXFx4n8m3SlDRf68PgykF0tGGj4FzDSzx4B/AHamqn6eV+bOjT66nTsXMztYwhk4cGCwr7w85gcAU6Z0qwvuzt69e9m9eze7d+9m165dB7/u2rWLnTt3snPnTnbt2sWOHTvaPTZu3Mjf9+xhO9AS+cLvvQcjRgAwePBghg8fTllZ2cHHYYcdxogRI9p9Pfzwwxk4e3b7Pw8ItmfPjh/MOichEr/kYmaPAqcCw4APgJuBYgB3X2TB7+ILCGbCNAKT3T1uLaWnJZe80JVASUV5o6cqKvDaWhqA7aHHNqB+2DDqb72Vbdu2sW3bNrZu3UpdXR11dXUHv4/2m8EA4Iiwx5EEv/qNBEatXs2oUaMoKyuLfo6gk9942LCh4/PR/jyLi2HQINi+XQEvWauzkktCNfRUUKB3Q7aNKLv5IXPgwAHq6+v54IMPeP/99w8+tsyZw5bdu9kCvBt6NEX8bJ8+fRg9ejQVFRVUVFRQXl7OmDFj+GhVFR8lKPu0q/abQbSyUqwPgHDFxfDggwp1ySoKdEmdZH7IRHxAOFDXrx+bvvc9Nn/iE2zatImNGzdSW1vLhg0b2LBhA1u3bm33EoOAj4Ue44CjSks56le/Yty4cQwePPjDhgUFkMi//UMOgT17unc8IimgQJfc0cUPiMbGRt555x3eeuAB3lqwgLeam3kLeBN4Bwgfm48YMYLx48dzzDHHcMzy5RyzfTvHEJzQ7dSMGVBd3dMjE0kKBbr0DhEfBk233spbEybwxhtvsG7dOl5//XXWrl3Lq6++yp6wUfco4HjghLCvHyHsMmozWLo0+GDJtrKX9DoKdJEwra2tbNq0ibULFvC3++/n5R07eBl4HTgQajMYqAx/HHkk5XfcgU2b1v6cgRlMn64RvKSNAl0knuXL2TdpEmuBl4Ca0OMVoDnUZERhIaccOMApwCnAiQRX+GoEL+mkQBdJxNVXw8L2a8vtJwj1Pw8dygvbt/Mn4O3Qvj4EF16cCpx62GGcfPvt9Js5s+OsnyuugJUrFfKSFAp0kURdfTUsWtR+BkzbVMzZs6G2lveB54H/Bv4T+AvBydcS4GTgDOAsglJNYbT30HRI6YHOAl03uBAJV10dlE/Ky4NSSnn5h/Pq584FM0YAFwI/BlYTXFD1dFkZ1wINBFfenQwMBy4GHiBiLYzmZpg2LZ1HJb2ERugiXZHACL4O+C3wDLCKD5cdPQmYSLDe9LGALVumUbp0mUboIskSbwTfvz9lwKXAgwRXu74MzCX4z/avwHEE0yJvuPJKXnzxxR6trCkSTiN0kWSKnOVSX9/uStP3gKeBJwlG8M1AeXk5F110ERdffDGf+tSnun6zFOlVNEIXSZeqqmAxsNbW4OuiRe12Hw5cRRDqWwnu73jMkCHcfffdnHTSSRx11FHM/fKX2ThypG5NKF2mQBdJpaoqGDAg6q4hwBXAih07+OCDD7j//vs5orCQm37+cyrefZfT3VlaW0vDpEkwbJiCXeJSoIukWsQovYONGxk6dChTpkzh2b17eRu4BdgAXE6wlPC19fW8pmCXOBToIqlWVRUs8BVL+P11N25kDPB9YD3BPPfzgXuB8cBp9fU8PmkSTWYqx0gHCnSRdKiuDkI98oRn6C5VB4WFuwGfBZYT3IX9doJR+8XAGOCHtbXsmDQpmEopggJdJH06m/LYJjT1MdJwYBbBqP1p4OOh7dHAjQsXsunQQzVaF01bFMk6y5fDddcFUx478RfgTuBnBKP5S4Gbhgxh3IIFumApj2naokguqaqCbdtg2TIoLY3Z7ESCcsxbwEzg/wMf37GDyydN4k3V2HslBbpItgoP9vLymM3KgbsI7tB0PfAEcDRwRW0t6zUzpldRoItku7aLldyjn1gNGU6wYNjbwHXA4wS19mvq69l25ZUK9V5AgS6SS9pOrHZSihkB/D+CYJ8CVAMfbWlh3uTJ7Nu3Lz39lIxQoIvkmgRr7IcDi4D/Af4J+E5zM0cPH86jjz6qBcHylAJdJFe1BXsnZRgILkh6mmBJ30N37+YrX/kKp59+Oq+//nq6eippokAXyXXh89s7cTrBfVIXlpTwl2ef5biPf5ybBg9m74MPpqWbknoKdJF8EH7idNmymM0KgelNTaxz52Jg7q5dHDtlCr/6/OeDaY5a4TGnKdBF8k28tWOAw4ClwO+AYnfOW7mSybW17HSH2lqYPDmY7qiAzykKdJF8VF0dc9necJ8juKPSbOAR4BMEtXaam4MrVdsC/rLLgjq9wj2rKdBF8tWiRVBSErdZH+AHwJ+A/sCZwNeBPeGN2mbFtIW7FgTLSgp0kXxVVQVLlrRfDGzGjI6Lf4VmyPwD8BLB1aYLgROA1dFe1x0WLlSoZ6GEAt3MzjGzdWa23sxmRdk/2sz+YGYvmdkrZnZe8rsqIl0WeUu86upghcfwkJ8+/WDI9yO4KOlZgvudfhr4CRB11rpCPevEDXQzKwTuAc4lmNJ6qZmNj2h2E/C4u38SuITg4jQRyUbxQr60lM+WlPAScA7BMgL/AuyI9loLF2qtmCySyAh9ArDe3d929ybgMeCCiDYODAp9PxjYkrwuikjKhYf8tm2wZAlDy8t5kmCJ3v8gWN0xagmmvl519SyRSKAfSXDDlDabQ8+FuwWYZGabgZXANdFeyMymmlmNmdXU1dV1o7sikhahgDd3bpgxgz8CLQQlmAeitW+rq2uknlGJBHq0a4ojS2qXAg+5+0jgPGCpmXV4bXdf7O6V7l5ZVlbW9d6KSPpVV/OPM2bwV+BU4GvADQQB38G0aWnsmERKJNA3A6PCtkfSsaQyhWC1Ttz9eaAvMCwZHRSRLFBdzdAZM1gJXEuw/voXgJ2R7RoaNErPoEQCfTUw1szGmFkJwUnPpyLabCRYKgIz+zhBoKumIpJPqqspWraM+aWl3EtwAdLJBPc5bee669LeNQnEDXR3byG4w9Uq4DWC2Sxrzew2M5sYanYjcJWZvQw8ClzpWp9TJP+EVnicOmAAvwG2Esya+K/wNvX1WjIgQ3STaBHpuuXLYdIk3iI4abaRoOb6hch2xcUwaBBs3w6jR8PcubqBdQ/pJtEiklyhUP4o8BxwLHAh8HBku2hrwmh6Y8oo0EWke0J3SyoDfk8wA+ZKgvuaxqTpjSmlQBeR7pk//+DiXwOBFcCXgG8B3yHGcgFtLr9coZ4CCnQR6Z6Ixb/6lJbyWHEx04F5BFcXxgz11laVX1JAgS4i3RexZEDhgw9SPXo03yRYAOobdBLqWrUx6RToIpI8VVVYbS3zWlu5/vrr+QnBnOZOyy8K9aQpynQHRCT/mBl33nknB9au5a5nnqGQoAwTbR0RIAh1CFZ+lG7TCF1EUsLMuPvXv+brZ57Jj4HvopF6qinQRSRlzIx/W7WK6dOn80NgTrwf0JTGHlGgi0hKmRn33HMPV1xxBTcD98f7Aa0F020KdBFJuYKCAu677z7OPvtsppvxdGeN6+vT1a28o0AXkbQoLi7miSee4IQTT+TLRUW82Flj1dK7RYEuImkzYMAAVqxYweGjR/N54I1YDVVL7xYFuoik1WGHHcavf/1rzIyzgfdjNVQtvcsU6CKSdmPHjmXFLbewlWCVxv3RGqmW3mUKdBHJiAnf/z6PnHUWLwBXE2OOumrpXaJAF5GM+ZdVq7gJWEKw9ksHqqV3iQJdRDLq1qFD+QLBQl7PRmugWnrCFOgiklEFP/kJy4CxwEVAbWSD+nqN0hOkQBeRzKqqYtCMGTwJNBOcJG2MbKO10xOiQBeRzKuuZhzwU+CvwPTI/bp1XUIU6CKSHUpLOQ+4GVgKPBKtjW5d1ykFuohkh9A9Sm8C/g/BVMYOV5Lq1nWdUqCLSHYI3aO0sLSUZUAf4BKiXHSkW9fFpEAXkexRVQXbtjFywAAeBF4CZsVqq5p6Bwp0Eck+ixYxsaSEa4C7IfZyu6qpt6NAF5HsEyq/zBs6lOOBK4F3o7VrbYXJkxXqIQp0EclOVVX0ra/nZ5dcwl6CUI+63ktzM0ybltauZSsFuohktaMefZQ7+/Tht3Ry+7qGBo3SUaCLSA6Yet99nGbGjcDGWI1UT1egi0j2K7jsMh64805aganEKL2onp5YoJvZOWa2zszWm1nUWURm9mUze9XM1prZT5PbTRHp7cZcfz13XH45q4CHYjVqbu7VqzOae9TPug8bmBUSXLB1JrAZWA1c6u6vhrUZCzwOfM7d/25mw919a2evW1lZ6TU1NT3tv4j0Iq2trZxWXMzLra2sBY6M1TBOruUyM1vj7pXR9iUyQp8ArHf3t929CXgMuCCizVXAPe7+d4B4YS4i0h0FBQU8MG8eTcA0YpReoNeWXRIJ9COBTWHbm+n4wTgOGGdm/21mL5jZOdFeyMymmlmNmdXU1dV1r8ci0qt97MYbuX3SJFYQLOIVVS89QZpIoFuU5yI/GIsI1qc/FbgUuN/MhnT4IffF7l7p7pVlZWVd7auICADXPPwwny4q4gZge7QGvXQRr0QCfTMwKmx7JLAlSpsn3b3Z3d8B1hEEvIhI0hUUFLBwzhx2ALNjNeqFa6gnEuirgbFmNsbMSggWQHsqos0vgdMAzGwYQQnm7WR2VEQk3CdmzWJmcTH3Ams6a9iLriKNG+ju3gLMBFYBrwGPu/taM7vNzCaGmq0C6s3sVeAPwLfcvT5VnRYRAbj1nnsYDnwdaI3VqKGh15Re4k5bTBVNWxSRZHhk+nSuuPde7gemdNZw2bJg0a8c19NpiyIiWeuyhQv5zGc+wyxinCBt0wtKLwp0EclpZsaCBQvYbsa/dtawFyzgpUAXkZx3/PHH8/WZM1kI/KWzhnm+LIACXUTywm233UbZ8OFcU1AQ+wrS+vyeq6FAF5G8MGTIEH7wgx/wp9ZWftlZwzye8aJAF5G8MXnyZI4++mi+C7TEarRwYd6GugJdRPJGUVERd9xxB+uABzprmKdXkCrQRSSvTJw4kU+PG8ctQENnDfPwBKkCXUTyipkx78EHeR+4q7OGeXiCVIEuInnnlFNO4Ytf/CLziovpTQt1K9BFJC/dfvvtNLa28oNMdySNFOgikpeOPvpopkyZwkJ6z9KvCnQRyVs333wzRXSyZnqeUaCLSN464ogjuB74GfBqvMZ5QIEuInntBqA/MDfazjybi65AF5G8VkpwA4zHgDcid+bZXHQFuojkt9JSbgT6AP83cl+ezUVXoItIfps/n+HAdGAZ+T3jRYEuIvktdNu5bwFFwO2R+/Oojq5AF5Fe4XDgKuAhoDZ8Rx7V0RXoIpL/SksB+DZgwA/D9+VRHV2BLiL5b/58AEYBXyVYWvfd8P15UnZRoItI/gvV0QFmAa3AvPD9eVJ2UaCLSO8QKrtUAJcBi+HDlRjzpOyiQBeR3iFUdgH4JrAPuDdjnUkNBbqI9A5hZZfxwLnAAmB/25N5UEdXoItIr3Q98AHwaNsTeVBHV6CLSO8RqqMDnAEcS3CbOoe8qKMr0EWk9wiroxvBKP0V4PdtT+Z42UWBLiK9R1gdHeArwHDCbiad42WXhALdzM4xs3Vmtt7MZnXS7ktm5mZWmbwuiogkUVjZpS9wNbACeB1yvuwSN9DNrBC4h+Ck8HjgUjMbH6XdQOBa4MVkd1JEJGnCyi4AMwiW1r07I51JrkRG6BOA9e7+trs3EawTf0GUdnMILr7al8T+iYgkV0TZZTgwCXgE2JaJ/iRRIoF+JLApbHtz6LmDzOyTwCh3fzqJfRMRSYvrgb3k/oVGiQS6RXnOD+40KyA4p3Bj3Bcym2pmNWZWU1dXF6+5iEhaHAOcRXChUdPDD2e4N92XSKBvJlikrM1IYEvY9kCC6ZzPmtkG4GTgqWgnRt19sbtXuntlWVlZ93stIpJk1wHvA0/OnJnprnRbIoG+GhhrZmPMrAS4BHiqbae773T3Ye5e4e4VwAvARHevSUmPRUR6KmymS5uzgdHAvXv2pL07yRI30N29BZgJrAJeAx5397VmdpuZTUx1B0VEki5ipgtAIcEdjX4HrL/zznT3KCnM3eO3SoHKykqvqdEgXkQyxDqeHtxCMEq/sW9ffrh3b9q7lAgzW+PuUa/10ZWiItI7RSm7HAF8AXhw3z6amprS3qWeUqCLSO8UpewCMI3gxhe/uDHuxL2so5KLiPReUcourcBHgI8WFfG75ua0dykelVxERKKJUnYpIDg5+vuWFt588820d6knFOgi0nvFKLt8FSgCFn/jG2ntTk8p0EWk94pY16XN4cBE4KGVK9n/0EPp7FGPKNBFRKKYRrBY1y+mT890VxKmQBeR3i1KHR2CW9SNAe7dvz9n7mSkQBeR3m3+fCgp6fB028nRZ4E3r7oqzZ3qHgW6iPRuVVWwZAkUdIzDywiWm122d29OjNIV6CIiVVXwyCMdnh4JfA5YBvi116a7V12mQBcRgZgzXi4D3gb+tH171o/SFegiIm2inCD9Z6AfsBRg8uSsDnUFuohImygXGg0ELgQeB/Y3N8N116W7VwlToIuItKmqggEDOjx9GfB3YAVAfX2aO5U4BbqISLhFizpMYzwDOIxQ2SWLKdBFRMK1TWMMUwR8hWCEXg9ZW0dXoIuIRIoy4+UyoJmgls60aWnuUGIU6CIi0UTMeDkBOIZQ2aWhIStH6Qp0EZFoIma8GMEo/XlgPWTlKF2BLiISTZQZL1UEwb4csnKUrkAXEYll0aJ2myOB0wjKLg5ZNyddgS4iEkuUUfplwFvAnyHr5qQr0EVEOhMxSv8iUAw80faEGVRUZEX5RYEuItKZiCmMQ4AzgZ8TKrsA1NZmxTovCnQRkXgipjBeBNQCa8KfzIJ1XhToIiLxRExhnEhw9egTke0yXFNXoIuIxFNVBTNmBPVyYChwOkGge2TbDJZdFOgiIomoroalS6G8HIAvEcx2eTmyXQbLLgp0EZFEVVXBhg1QWsoXgUKyq+yiQBcR6ar58xkGnErEbJc2GSq7JBToZnaOma0zs/VmNivK/hvM7FUze8XMfmdm5cnvqohIlghNZfwS8AawNnJ/hsoucQPdzAqBe4BzgfHApWY2PqLZS0Clux9H8BvIvGR3VEQkq5SWciHB2i7ZUnZJZIQ+AVjv7m+7exPwGHBBeAN3/4O7N4Y2XyBY8kBEJH/Nn89hwGeJEuiQkbJLIoF+JLApbHtz6LlYpgC/irbDzKaaWY2Z1dTV1SXeSxGRbBMqu1xEUHJ5LXJ/BsouiQS6RXmuwzkAADObBFQCP4q2390Xu3ulu1eWlZUl3ksRkWyUZWWXRAJ9MzAqbHsksCWykZmdAcwGJrr7/uR0T0Qki82fzxHAp4lRdkmzRAJ9NTDWzMaYWQlwCfBUeAMz+yRwL0GYb01+N0VEslDYbJdXCGa8tHP11WntTtxAd/cWYCawiqBM9Li7rzWz28xsYqjZj4ABwM/N7K9m9lSMlxMRyTsXhr7+R+SOhQvTenLU3KOWw1OusrLSa2pqMvLeIiJJM2wY1NdzPMEaL3+I3F9aCtu2Je3tzGyNu1dG26crRUVEeiK0EuP5wH8BOyL3p/HkqAJdRKQnQisxng8cIKhNZ4oCXUSkp6qrmQAMI0odPY0U6CIiSVAInEdwVWVL5M40nRhVoIuIJMn5wHaC9U/aSdNVowp0EZFkKC3lLIJb0z0duS9NJ0YV6CIiyTB/PoMJFuvqEOiQlrKLAl1EJBlCV42eT7BY1zuR+9NQdlGgi4gkS2kpXwh9uyJyXxrKLgp0EZFkmT+fjwFHEaPskmIKdBGRZAkru/wB2JPmt1egi4gk2flAE/DbNL+vAl1EJMk+DQwm/VeNKtBFRJKsGDiH4MRoaxrfV4EuIpIC5wMfAGvCn0zxXHQFuohICpwLfAyoC38yxXPRi1L66iIivU1pKdTXUwq8GbkvxXPRNUIXEUmm0A0vYkph2UWBLiKSTKG56DGlsOyiQBcRSbbS0tj7Ulh2UaCLiCRbvLJLiijQRUSSLV7ZJUV1dAW6iEi6paiOrkAXEUmFDNTRFegiIqmQgTq6Al1EJBXi1dFTQIEuIpIJKTgxqkAXEcmEFJwYVaCLiKRKmk+MKtBFRFIlzSdGFegiIqmS5guMEgp0MzvHzNaZ2XozmxVlfx8z+1lo/4tmVpHUXoqI5KNp05L6cnED3cwKgXsI1msfD1xqZuMjmk0B/u7uHwPuAn6Y1F6KiOSqzuroDQ1JHaUnMkKfAKx397fdvQl4DLggos0FwMOh758ATjczS1ovRURyVbw6+uzZSXurRAL9SGBT2Pbm0HNR27h7C7AT6PCxZGZTzazGzGrq6uoid4uI5J+qKhgwIPb+jRuT9laJBHq0kbZ3ow3uvtjdK929sqysLJH+iYjkvkWLYu8bPTppb5NIoG8GRoVtjwS2xGpjZkXAYGB7MjooIpLzqqpgxgyIrET37w9z5ybtbRIJ9NXAWDMbY2YlwCXAUxFtngKuCH3/JeD37t5hhC4i0mtVV8PSpVBeHgR7eTksXpzUNV+K4jVw9xYzmwmsAgqBJe6+1sxuA2rc/SngAWCpma0nGJlfkrQeiojki6qqlC7aFTfQAdx9JbAy4rnvh32/D7gouV0TEZGu0JWiIiJ5QoEuIpInFOgiInlCgS4ikicsU7MLzawOqO3mjw8DtiWxO5mkY8k++XIcoGPJVj05lnJ3j3plZsYCvSfMrMbdKzPdj2TQsWSffDkO0LFkq1Qdi0ouIiJ5QoEuIpIncjXQF2e6A0mkY8k++XIcoGPJVik5lpysoYuISEe5OkIXEZEICnQRkTyR04FuZteEbl691szmZbo/3WVmt5jZu2b219DjvEz3qafM7Jtm5mY2LNN96Q4zm2Nmr4T+Pp4xsyMy3afuMrMfmdnroeP5hZkNyXSfusvMLgr9f281s5ybwmhm54Qya72ZzUr26+dsoJvZaQT3Mj3O3Y8BfpzhLvXUXe5+QuixMn7z7GVmo4AzgeTdWyv9fuTux7n7CcDTwPfj/UAW+w1wrLsfB7wBfDfD/emJvwH/DPwx0x3pKjMrBO4BzgXGA5ea2fhkvkfOBjowA7jD3fcDuPvWDPdHPnQX8G2i3IYwV7j7rrDNQ8jtY3kmdK9fgBcI7jqWk9z9NXdfl+l+dNMEYL27v+3uTcBjBIPSpMnlQB8H/JOZvWhm/2lmJ2W6Qz00M/Qr8RIzOzTTnekuM5sIvOvuL2e6Lz1lZnPNbBNQRW6P0MN9FfhVpjvRSx0JbArb3hx6LmkSusFFppjZb4ERUXbNJuj7ocDJwEnA42b2kWy99V2cY1kIzCEYBc4B7iT4j5eV4hzL94Cz0tuj7unsONz9SXefDcw2s+8CM4Gb09rBLoh3LKE2s4EWYHk6+9ZViRxLjrIozyU1r7I60N39jFj7zGwG8O+hAP+zmbUSLHhTl67+dUVnxxLOzO4jqNlmrVjHYmafAMYAL1twM9yRwF/MbIK7v5/GLiYk0b8T4KfACrI40OMdi5ldAZwPnJ6tg542Xfh7yTWbgVFh2yOBLcl8g1wuufwS+ByAmY0DSsjRldjM7PCwzQsJTvzkHHf/H3cf7u4V7l5B8A/4xGwM83jMbGzY5kTg9Uz1pafM7BzgO8BEd2/MdH96sdXAWDMbY2YlBPdefiqZb5DVI/Q4lgBLzOxvQBNwRbaPPDoxz8xOIPj1awMwLbPdEeAOMzsKaCVY5nl6hvvTEwuAPsBvQr85veDuOXk8ZnYh8G9AGbDCzP7q7mdnuFsJcfcWM5sJrAIKgSXuvhYbF+cAAAA4SURBVDaZ76FL/0VE8kQul1xERCSMAl1EJE8o0EVE8oQCXUQkTyjQRUTyhAJdRCRPKNBFRPLE/wKY5+FUc6WkpAAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "ecdf_max = ECDF(maxlist)\n", "plt.plot(ecdf_max.x, np.ones(len(ecdf_max.y))-ecdf_max.y, 'ro')\n", "\n", "t=np.arange(-6, 0, 0.1)\n", "plt.plot(t,weib(t),'k')\n", "\n", "#plt.yscale('log')" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "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": 2 }