Skip to content

Instantly share code, notes, and snippets.

@SusanEAllen
Last active June 9, 2020 00:58
Show Gist options
  • Save SusanEAllen/968d6d54e1a375f804a10cf360ca7ed4 to your computer and use it in GitHub Desktop.
Save SusanEAllen/968d6d54e1a375f804a10cf360ca7ed4 to your computer and use it in GitHub Desktop.
Test notebook for ChoseFractionSpilled
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [],
"source": [
"import importlib\n",
"import matplotlib.pyplot as plt\n",
"import numpy as np\n",
"\n",
"import moad_tools.midoss.random_oil_spills as ros\n",
"\n",
"%matplotlib inline"
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<module 'moad_tools.midoss.random_oil_spills' from '/Users/sallen/Documents/MIDOSS/moad_tools/moad_tools/midoss/random_oil_spills.py'>"
]
},
"execution_count": 7,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"importlib.reload(ros)"
]
},
{
"cell_type": "code",
"execution_count": 21,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"[0. 0.24266816 0.35966341 0.42851504 0.47787627 0.51841066\n",
" 0.55415534 0.58670183 0.61673461 0.64459766 0.67050337 0.69460976\n",
" 0.71704939 0.73794026 0.75739031 0.77549928 0.79235978 0.80805793\n",
" 0.82267388 0.83628224 0.84895248 0.86074926 0.87173281 0.88195918\n",
" 0.89148058 0.90034559 0.90859948 0.91628437 0.92343948 0.93010134\n",
" 0.93630395 0.94207896 0.94745587 0.9524621 0.95712322 0.96146302\n",
" 0.96550364 0.96926571 0.97276844 0.97602969 0.97906613 0.98189324\n",
" 0.98452546 0.98697622 0.98925803 0.99138254 0.9933606 0.99520229\n",
" 0.99691702 0.99851354 1. ]\n"
]
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXQAAAD4CAYAAAD8Zh1EAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+j8jraAAAevklEQVR4nO3deXyV1b3v8c+PzPMcwAwMMiuIEsCpVmsdwPZqLY61g7aX2mrvua/r6dHT9rSn1Xuq3vZe66lKOdZSq621xx6LiNrJAQeUQYgIAiEM2RDIPM/Z6/6RQEMayAZ28uzh+3698sp+9rPY+7dM+LpYez3rMeccIiIS/sZ4XYCIiASHAl1EJEIo0EVEIoQCXUQkQijQRUQiRKxXb5ybm+smTpzo1duLiISlDRs21Djn8oY651mgT5w4kfXr13v19iIiYcnM9h7rnKZcREQihAJdRCRCKNBFRCKEAl1EJEIo0EVEIsSwgW5mT5hZlZltOcZ5M7OHzazMzErN7JzglykiIsMJZIS+ArjyOOcXAVP7v5YCj516WSIicqKGXYfunHvDzCYep8nVwJOubx/etWaWaWbjnXOVQapRRCTk+f2Otu5eWjt7aOnsoaWjh9auHlo7e2nr/97a2ffcOcVZXDRtyGuDTkkwLiwqACoGHPv6n/u7QDezpfSN4ikuLg7CW4uInLpev6OpvZumjm6a2nto6uimecDjpo4emju6aenoD+vOHpoPP9fZH9ZdPQR6e4mvXXx6yAa6DfHckN1yzi0HlgOUlJTozhoiElQd3b3Ut3VR19pFQ1s39W193xvb+74a2rr6v/cdN3f00NjeF8rDSU2IJS0xltSEWFITY0lPiqMgM4mUhBhSE+JITYwlNSGGlIT+NgmxpCTEkhIfS3JCDKkJsSTHx5AcH0vMmKFi89QFI9B9QNGA40LgQBBeV0SimHOOpo4eals6qW3torali9rWTupauqht7ToS3HWtXdS3dlHX1kVHt/+Yr5cYN4aMpDgyk+LJSI6jMCuZjKQ4MpLiSE+K7fueGEd6UhzpibGkJcaR1h/cqQkjF8LBFIxAXwncaWbPAAuBRs2fi8ixtHf1UtXcQVVzJ1VNnVQ3d1Dd0kl1cyc1LV1UN/c9rm3tpLt36H/IpyXEkpUST3ZKPGPTE5kxLp3slDgyk+PJSo4nK7nvcXZKPJnJfaGdGBczyj0dfcMGupn9BrgYyDUzH/A9IA7AObcMWA0sBsqANuDWkSpWREJXr99R09JJZWMHBxs7ONTUQWX/94ONHRxq7qC6qZPmIaY3YsYYuanx5KYmkJeWwPRxaeSkxpOXmkBOajw5KQlkp/Sdz0qJIyE28sP5ZASyyuWmYc474I6gVSQiIam5o5v9De346trZ39DOgcZ2Khs6qGxs50BDX3D3+I8eUcfFGPlpiYzLSGTGuDQumppHfnoC+WmJ5KclkJ+eQF5qAlnJ8YwJgymNUOfZ9rkiElrau3rx1bexr66Niro2Kurb2VfXxv76dnz1bTR1HD2yjo8Zw7iMRE7LTGThpGzGZyYyPiOJ8RmJjE3vC/FsBfWoUqCLRJGmjm721rSxp7aVvbWt7KltO/K9urnzqLZJcTEUZSdRmJXMvAlZFGYlUZCVREFm3/fclASFdYhRoItEmO5eP/vq2iivbqW8uoXdNa19j2taqGnpOqrt2PQEJuSkcPG0PIqzkynOSaYoO5mirGRyU+MxU2CHEwW6SJjq7OmlvLqVnVUtlFW1UFbVzM5DLeypbT1qdUhOSjyT81K4dMZYJuWlMDEnhYm5yRRnJ5McrwiIJPppioQ4v9+xr66N7Yea2X6w+cj33TWt9PZ/CDnGoDg7mSn5aVw6cyxT8lOZnJfC6bmpZCTHedwDGS0KdJEQ0tHdy/aDzWytbGLrgSa2VjaxrbKJtq7eI22Ks5OZNjaNK88Yx9SxqUwbm8ak3JSoWGctx6dAF/FIR3cvWyub+MDXSKmvkS37G9lZ1czhlX+pCbHMGp/O9SVFzByfxvRx6UzNTyUlQX9tZWj6zRAZBb1+x45DzWyqaGDTvgY2+xrYWdVyZMokNzWe2QUZXH7GWM44LZ1Z4zMozErSKhI5IQp0kRFQ29LJ+r31vL+vgU0V9ZT6Go9Mm2QlxzGnMJPLZo1ldkEGswszGJeeqBUlcsoU6CKnyDlHeU0rG/bUs25PHRv21lNe0wr0XSk5a3w6180r5OziLOYWZTIhJ1nhLSNCgS5ygpxzlFW1sLa8lrXldby7u/bI+u7M5DhKJmRxXUkR8ydmcWZBhj6slFGjQBcZhnOO3TWtvFVWw9ryOtaW11Lb2hfgp2UkctHUPBZMyqZkYhaTc1M17y2eUaCLDKGmpZO3ymp4q6yGN3fWcKCxA+gL8I9Pz+PcyTmcNzmHwqwkTZ9IyFCgiwA9vX427mvgte1VvLa9mq2VTQCkJ8ZywZRcvn5JLhdOydX8t4Q0BbpErermTl7fUc2r26tYs6Oapo4eYscY8yZk8c0rpnPhlFzOLMgIizvViIACXaLI4Q8z/7j1EH/aeohNFQ0A5KclcOWZ47hkej4XTM0lPVGXykt4UqBLROv1Ozbuq+dP/SG+u3854ZzCDO66bBqfmJnPrPHpmkaRiKBAl4jT63es21PH6g8qeWnLQaqbO4mLMc6dnMNtF07ikzPzGZ+R5HWZIkGnQJeIcDjEXyyt5OUP+0I8IXYMl0zPZ9HscVwyI19TKRLxFOgStpxzbK1s4vn397Ny8wEONXWSGNcX4otnj+cTM/K1kZVEFf22S9ipqGtj5eYDPP/+fnZWtRA7xrh4eh7fuapAIS5RTb/5Ehbau3p5aUslz66vYG15HQDzJ2Zx3zVnctXs8WSlxHtcoYj3FOgSspxzbKpo4Nn1Pl7YfICWzh4m5CTzj5dP4+q5BRRlJ3tdokhIUaBLyGls6+Y/N/p45r197KxqISkuhqvmjOe6eYUsmJStJYYix6BAl5BR6mvgV+/s5YXSA3R0+zm7OJP7r53NVXPGk6YVKiLDUqCLp9q7enlh8wGeencvpb5GkuNjuPacQm5ZOIFZp6V7XZ5IWFGgiycONXXw5Dt7ePrdfTS0dTM1P5UfXH0G15xdoPXiIidJgS6jasv+Rn7+5m5WlR6gx++4fNZYbr1gEgs1Ny5yyhToMuL8fsdfPqriP9aU897uOlLiY7jl3Ancev4kinO0UkUkWBToMmJ6ev2sKq3k0dfK2HGohYLMJL5z1Uyun1+kaRWREaBAl6Dr6O7ldxt8LH9jFxV17Uwbm8pDN8zlU3PGExszxuvyRCJWQIFuZlcCPwFigMedc/cPOp8BPAUU97/mj5xzvwhyrRLi2rt6+dXaPfzHmt1UN3cytyiT737qDC6dka/7bIqMgmED3cxigEeAywAfsM7MVjrntg5odgew1Tn3aTPLA7ab2dPOua4RqVpCSkd3L0+/u4/HXttFTUsnF0zJ4Sc3zuW8yTn6oFNkFAUyQl8AlDnnygHM7BngamBgoDsgzfr+9qYCdUBPkGuVENPZ08tv11XwyKtlHGrq5LzJOTx2yznMn5jtdWkiUSmQQC8AKgYc+4CFg9r8FFgJHADSgBucc/7BL2RmS4GlAMXFxSdTr4SAnl4/v9vg49//spMDjR3Mn5jF/7thLuefnut1aSJRLZBAH+rfzG7Q8RXAJuATwOnAn8xsjXOu6ag/5NxyYDlASUnJ4NeQEOec48/bqrj/pW3sqm5lblEmDyyZw4VTcjW1IhICAgl0H1A04LiQvpH4QLcC9zvnHFBmZruBGcB7QalSPLepooF/W72N93bXMTk3hZ99fh6XzxqrIBcJIYEE+jpgqplNAvYDNwI3D2qzD7gUWGNmY4HpQHkwCxVv7K1t5cFXtvNiaSW5qfHce82Z3Di/iDgtPxQJOcMGunOux8zuBF6hb9niE865D83s9v7zy4B7gRVm9gF9UzR3O+dqRrBuGWEtnT389K9l/PzNcmLHjOF/XDqVpRdNJlV3AxIJWQH97XTOrQZWD3pu2YDHB4DLg1uaeME5xx82HeDfVm+jqrmTJfMK+acrppOfnuh1aSIyDA235Igt+xv515Ufsn5vPXMKM1j2+XmcU5zldVkiEiAFutDY1s2Dr3zEr9/bR3ZyPA98djbXzSvS1Z0iYUaBHsWcc6wqreT7L2ylvq2LL50/kf/5yWlkJGnjLJFwpECPUr76Nv7l+S28ur2a2QUZrLh1PmcWZHhdloicAgV6lOnp9bPi7T38+I87MIPvfmoWXzx/IjGaXhEJewr0KPLRwSa++btSPtjfyCdm5HPvNWdSkJnkdVkiEiQK9CjQ0+vnZ2+U89Cfd5CRFMcjN5/D4tnjdJWnSIRRoEe4XdUt3PXsZjZVNLB49jjuu2Y22SnxXpclIiNAgR6h/H7Hirf38MDLH5EUH8PDN53Np+eM16hcJIIp0COQr76Nu57dzLu767h0Rj4/vHa2rvQUiQIK9Ajz0geV3P1cKX4HDy6Zw3XzCjUqF4kSCvQI0d7Vyw9WbeU37+3jrMIMHr7pbCbkpHhdloiMIgV6BNh+sJlv/GYjOw618NWPT+auy6YTH6vtbUWijQI9jDnneOrdfdy3aitpiXE8edsCLpqW53VZIuIRBXqYau3s4e7nSllVWslF0/L48XVnkZeW4HVZIuIhBXoY2lXdwu2/2sCu6ha+ecV0vvbx07Uzoogo0MPNy1sO8o+/20x87BievG0hF07N9bokEQkRCvQw0dPr58d/2sFjr+3irMIMHr1lnvZhEZGjKNDDQF1rF9/4zUbeKqvl5oXFfO/Ts0iIjfG6LBEJMQr0ELfzUDO3/XIdh5o6eXDJHK4vKfK6JBEJUQr0EPb6jmrufHojCXExPPvV85hblOl1SSISwhToIeqXb+/h+y98yPRx6Tz+xRLNl4vIsBToIaan188PVm3lyXf28smZY/nJjXNJSdCPSUSGp6QIIU0d3dzx9EbW7Kxh6UWTufvKGbo1nIgETIEeIg41dfDFJ96jrKqFBz47mxvmF3tdkoiEGQV6CNhV3cIXfv4eDW1drLh1gS4WEpGTokD32Pv76rltxTpixhi//ep5nFmQ4XVJIhKmFOgeenV7FV9/aiN5aQn86ssLtH+5iJwSBbpHntvg4+7nSpk+Lo0Vty7QTokicsoU6B54fE059724jQum5LDslnmkJcZ5XZKIRICAbmtjZlea2XYzKzOze47R5mIz22RmH5rZ68EtM3I88moZ9724jcWzx/HEl+YrzEUkaIYdoZtZDPAIcBngA9aZ2Urn3NYBbTKBR4ErnXP7zCx/pAoOV845fvKXnTz0551cM/c0fnTdWcTG6DZxIhI8gSTKAqDMOVfunOsCngGuHtTmZuD3zrl9AM65quCWGd6cc/zoj9t56M87WTKvkB9fP1dhLiJBF0iqFAAVA459/c8NNA3IMrPXzGyDmX1hqBcys6Vmtt7M1ldXV59cxWHGOccPX/qIR17dxU0Linnws3N09aeIjIhAPhQdKn3cEK8zD7gUSALeMbO1zrkdR/0h55YDywFKSkoGv0bEcc7x/Re2suLtPXzxvAn86387AzOFuYiMjEAC3QcM3IS7EDgwRJsa51wr0GpmbwBnATuIUs45vrfyQ558Zy9fvnAS37lqpsJcREZUIFMu64CpZjbJzOKBG4GVg9r8AfiYmcWaWTKwENgW3FLDh3OO+1/+iCff2cvSiyYrzEVkVAw7QnfO9ZjZncArQAzwhHPuQzO7vf/8MufcNjN7GSgF/MDjzrktI1l4KHvk1TJ+9no5t5xbzD8vmqEwF5FRYc55M5VdUlLi1q9f78l7j6RfvLWb77+wlWvPLuBH153FGH0AKiJBZGYbnHMlQ53T2rkgenZ9Bd9/YStXnDGWB5fMUZiLyKhSoAfJi6WV3PNcKR+bmsvDN52tdeYiMuqUOkHw6kdV/MMz7zNvQhY/+/w8EmJjvC5JRKKQAv0Ulfoa+PrTG5kxPo2ff2k+yfHa70xEvKFAPwW++ja+/Mv1ZKfE84svLSBdG22JiIc0nDxJTR3d3LZiHR3dvfz6Kwu1n7mIeE4j9JPQ3evnjqc3Ul7dyrJb5jF1bJrXJYmIaIR+opxz/MvzW1izs4YHl8zhgim6obOIhAaN0E/QY6/v4pl1Fdx5yRSuLyka/g+IiIwSBfoJWFV6gAdf3s6nzzqN/3XZNK/LERE5igI9QNsqm7jr2c2UTMji/+gqUBEJQQr0ALR09nDH0xtJT4rj0VvOITFOFw6JSOhRoA/DOcc9z5Wyp7aVf7/pbPLTEr0uSURkSAr0YTy1di+rSiu56/LpnDs5x+tyRESOSYF+HKW+Bu5dtY1LpufxtY+f7nU5IiLHpUA/hsa2br7+9Eby0hL4v9fP1YegIhLydGHREJxz3PW7zRxq6uDZr55HVkq81yWJiAxLI/QhPL5mN3/edoh/XjSTs4uzvC5HRCQgCvRBth9s5oGXP2LRmeO49YKJXpcjIhIwBfoAfr/jW//1AWmJsfzvz8zWzZ1FJKwo0Af47foKNuyt51uLZ5KteXMRCTMK9H7VzZ38cPU2Fk7KZsm8Qq/LERE5YQr0fve9uJWObr+mWkQkbCnQgTd2VPOHTQf42sWnMyU/1etyREROStQHekd3L//yhy1Myk3haxfralARCV9Rf2HRT/9axt7aNn79lYXaRVFEwlpUj9B3HGrmZ2/s4tpzCjhft5ITkTAXtYHunOPb//UBKQmxfHvxTK/LERE5ZVEb6Gt21rBuTz3/dMUMclITvC5HROSURW2gP/7mbvLTErTmXEQiRlQG+vaDzbyxo5ovnj+R+Nio/E8gIhEooDQzsyvNbLuZlZnZPcdpN9/Mes1sSfBKDL6fv1lOUlwMn1tY7HUpIiJBM2ygm1kM8AiwCJgF3GRms47R7gHglWAXGUxVzR08//4BlswrJDNZ+7WISOQIZIS+AChzzpU757qAZ4Crh2j3DeA5oCqI9QXdr97ZS7ffz20XTvK6FBGRoAok0AuAigHHvv7njjCzAuAzwLLjvZCZLTWz9Wa2vrq6+kRrPWXtXb08tXYvn5w5lkm5KaP+/iIiIymQQB9qpyo36Pgh4G7nXO/xXsg5t9w5V+KcK8nLywu0xqB5bqOP+rZu/vvHJo/6e4uIjLRALv33AUUDjguBA4PalADP9O9SmAssNrMe59zzQakyCPx+xxNv7mZOYQbzJ+q2ciISeQIJ9HXAVDObBOwHbgRuHtjAOXdkQtrMVgCrQinMAf76URXlNa08fNPZ2h5XRCLSsIHunOsxszvpW70SAzzhnPvQzG7vP3/cefNQ8fib5ZyWkciiM8d5XYqIyIgIaLdF59xqYPWg54YMcufcl069rODasr+RteV1fHvxTOJidCGRiESmqEi3x9eUk5oQyw0LioZvLCISpiI+0KuaO1hVWskN84tIT4zzuhwRkRET8YH+5s4aevyOa88pGL6xiEgYi/xAL6shOyWemePSvS5FRGRERXSgO+d4u6yW807PYcwYLVUUkcgW0YG+q7qVg00dXHC6bi8nIpEvogP97V01AFyo+4WKSBSI6EB/q6yGwqwkinOSvS5FRGTERWyg9/od7+yq1XSLiESNiA30Lfsbaero4YKpCnQRiQ4RG+hv9c+fn396jseViIiMjsgN9LIaZoxLIzc1wetSRERGRUQGekd3L+v31HOBVreISBSJyEDfuLeezh4/F0zRdIuIRI+IDPQ3y2qIHWMsmKRAF5HoEZGB/tauWuYWZZKaENB27yIiESHiAr2xvZsPfA2cr/lzEYkyERfoa8tr8Ttd7i8i0SfiAv3tshqS4mKYW5TpdSkiIqMq4gL9zbIaFk7OJj424romInJcEZV6Bxs72FXdqv1bRCQqRVSgH94u93ytPxeRKBRRga7bzYlINIuYQNft5kQk2kVMoJfX6HZzIhLdIibQdxxsBmBOYYbHlYiIeCNiAt1X3w5AUbZuNyci0SliAr2ivo30xFgykuK8LkVExBMRE+i++nYKszQ6F5HoFTGBXlHXRlF2ktdliIh4JqBAN7MrzWy7mZWZ2T1DnP+cmZX2f71tZmcFv9Rjc85phC4iUW/YQDezGOARYBEwC7jJzGYNarYb+Lhzbg5wL7A82IUeT21rF+3dvRRlaYQuItErkBH6AqDMOVfunOsCngGuHtjAOfe2c66+/3AtUBjcMo/v8AoXjdBFJJoFEugFQMWAY1//c8fyZeCloU6Y2VIzW29m66urqwOvchgVdW0AFGoOXUSiWCCBPtR19G7IhmaX0Bfodw913jm33DlX4pwrycvLC7zKYWiELiICgdx00wcUDTguBA4MbmRmc4DHgUXOudrglBeYivo2spLjdA9REYlqgYzQ1wFTzWySmcUDNwIrBzYws2Lg98DnnXM7gl/m8fnq23WFqIhEvWGHtM65HjO7E3gFiAGecM59aGa3959fBnwXyAEeNTOAHudcyciVfTRfXRszxqeN1tuJiISkgOYonHOrgdWDnls24PFXgK8Et7TA+P0OX0M7l80a68Xbi4iEjLC/UrS6pZOuHj+FWoMuIlEu7APdV394yaLm0EUkuoV9oFfU9W+bqxG6iES5sA/0IyN0rUEXkSgX9oFeUddObmoCiXExXpciIuKpsA90X4O2zRURgQgI9Io6bZsrIgJhHui9fseBhnZ9ICoiQpgH+sGmDnr8TiN0ERHCPNB9/dvmag5dRCTMA71C2+aKiBwR1oHuq2/DDE7LTPS6FBERz4V1oFfUtTM2LZGEWK1BFxEJ60D31WsNuojIYWEe6FqDLiJyWNgGenevn8pGrUEXETksbAO9sqEDv9MKFxGRw8I20P+2D7pG6CIiEMaBXtEf6EUaoYuIAGEc6L76dmLGGOMztAZdRATCONAr6toYl55IbEzYdkFEJKjCNg199e1agy4iMkDYBnpFfZtWuIiIDBCWgd7Z08uhpk59ICoiMkBYBvr+I7ssaspFROSwsAx0X3+gF2VrhC4iclhYBvrhNegaoYuI/E1YBrqvvp24GGNsutagi4gcFpaBXlHXxmmZScSMMa9LEREJGWEZ6L76dq1wEREZJEwDvU3z5yIigwQU6GZ2pZltN7MyM7tniPNmZg/3ny81s3OCX2qf9q5ealq6tMJFRGSQYQPdzGKAR4BFwCzgJjObNajZImBq/9dS4LEg13mETytcRESGFMgIfQFQ5pwrd851Ac8AVw9qczXwpOuzFsg0s/FBrhX42xp0XfYvInK0QAK9AKgYcOzrf+5E22BmS81svZmtr66uPtFaAUhLjOWKM8YyIUeBLiIyUGwAbYZaG+hOog3OueXAcoCSkpK/Ox+IkonZlEzMPpk/KiIS0QIZofuAogHHhcCBk2gjIiIjKJBAXwdMNbNJZhYP3AisHNRmJfCF/tUu5wKNzrnKINcqIiLHMeyUi3Oux8zuBF4BYoAnnHMfmtnt/eeXAauBxUAZ0AbcOnIli4jIUAKZQ8c5t5q+0B743LIBjx1wR3BLExGRExGWV4qKiMjfU6CLiEQIBbqISIRQoIuIRAjr+zzTgzc2qwb2nuQfzwVqglhOOFCfo4P6HB1Opc8TnHN5Q53wLNBPhZmtd86VeF3HaFKfo4P6HB1Gqs+achERiRAKdBGRCBGugb7c6wI8oD5HB/U5OoxIn8NyDl1ERP5euI7QRURkEAW6iEiECOlAD6WbU4+WAPr8uf6+lprZ22Z2lhd1BtNwfR7Qbr6Z9ZrZktGsbyQE0mczu9jMNpnZh2b2+mjXGGwB/G5nmNkLZra5v89hvWurmT1hZlVmtuUY54OfX865kPyib6veXcBkIB7YDMwa1GYx8BJ9d0w6F3jX67pHoc/nA1n9jxdFQ58HtPsrfbt+LvG67lH4OWcCW4Hi/uN8r+sehT5/C3ig/3EeUAfEe137KfT5IuAcYMsxzgc9v0J5hB5SN6ceJcP22Tn3tnOuvv9wLX13hwpngfycAb4BPAdUjWZxIySQPt8M/N45tw/AORfu/Q6kzw5IMzMDUukL9J7RLTN4nHNv0NeHYwl6foVyoAft5tRh5ET782X6/g8fzobts5kVAJ8BlhEZAvk5TwOyzOw1M9tgZl8YtepGRiB9/ikwk77bV34A/INzzj865Xki6PkV0A0uPBK0m1OHkYD7Y2aX0BfoF45oRSMvkD4/BNztnOvtG7yFvUD6HAvMAy4FkoB3zGytc27HSBc3QgLp8xXAJuATwOnAn8xsjXOuaaSL80jQ8yuUAz0ab04dUH/MbA7wOLDIOVc7SrWNlED6XAI80x/mucBiM+txzj0/OiUGXaC/2zXOuVag1czeAM4CwjXQA+nzrcD9rm+CuczMdgMzgPdGp8RRF/T8CuUpl2i8OfWwfTazYuD3wOfDeLQ20LB9ds5Ncs5NdM5NBP4T+HoYhzkE9rv9B+BjZhZrZsnAQmDbKNcZTIH0eR99/yLBzMYC04HyUa1ydAU9v0J2hO6i8ObUAfb5u0AO8Gj/iLXHhfFOdQH2OaIE0mfn3DYzexkoBfzA4865IZe/hYMAf873AivM7AP6piPuds6F7ba6ZvYb4GIg18x8wPeAOBi5/NKl/yIiESKUp1xEROQEKNBFRCKEAl1EJEIo0EVEIoQCXUQkQijQRUQihAJdRCRC/H/v3mO1RLDoOAAAAABJRU5ErkJggg==\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"fraction = np.arange(0, 1+0.02, 0.02)\n",
"cumulative = ros.CumulativeSpillFraction(fraction)\n",
"print(cumulative)\n",
"plt.plot(fraction, cumulative);"
]
},
{
"cell_type": "code",
"execution_count": 18,
"metadata": {},
"outputs": [],
"source": [
"random_seed = 101\n",
"random_generator = np.random.default_rng(random_seed)"
]
},
{
"cell_type": "code",
"execution_count": 19,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"0.63\n",
"0.03\n",
"0.31\n",
"0.15000000000000002\n",
"0.03\n",
"0.55\n",
"0.43\n",
"0.05\n",
"0.77\n",
"0.01\n"
]
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXQAAAD4CAYAAAD8Zh1EAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+j8jraAAAgAElEQVR4nO2de5Bc1X3nv79u9aAeYTQjPLbFgKJHMCocAbJlI6Kt3YhsAOPYTBywjc0mm3JMUYm3NnFKG7FRMNhyQaK1nUrshMWJ81hYbEOoCQt4wWW0zha2MPIKIQOSeclCw2sASSBphHpmzv7RfVqnb5/H79x7bnffnvOpopC67z33vPqo+3t+5/sjIQQikUgkUnxK3a5AJBKJRMIQF/RIJBLpE+KCHolEIn1CXNAjkUikT4gLeiQSifQJ87r14Le//e1i6dKl3Xp8JBKJFJKf/OQnrwohRnTvdW1BX7p0KbZv396tx0cikUghIaKfm96LkkskEon0CXFBj0QikT4hLuiRSCTSJzgXdCL6JhG9QkQ/NbxPRPSXRPQ0ET1GRO8NX81IJBKJuOB8Q/8HAJdY3v8ggDMb/10N4G+yVysSiUQivjijXIQQ/0pESy2XXAbgn0Td5WsbEQ0R0WIhxIuB6hic8R0T2HL/HrxwcAqnDVWx4eKzMLZ6NHU5EwenUCbCjBAYzVBeCEK1rRMUqa4csrYn7f3jOyZw/d2P4+BUDQAwPFjB5z/8Huu9yWetXzmCrbsncx+L0GOepu151LFX5nKIsMVRAM8rf9/feK1tQSeiq1H/Fo8lS5YEeLQ/4zsmcO1duzBVmwEATBycwrV37QIA7w+fWs5Mw7UybXkhCNW2TlCkunLI2p6094/vmMCGO3aiNnvCNfXA0Ro23LnTeK/uWbdu29d8P6+xCD3madqeRx17aS6H2BQlzWtaT14hxC1CiDVCiDUjI9q4+NzZcv+eZsdLpmoz2HL/nszlZCkvBKHa1gmKVFcOWduT9v4t9+9pWdAktRlhvNc2d32e7UvoMU/T9jzq2EtzOcSCvh/AGcrfTwfwQoByc+GFg1Ner/uWk7a8EIRqWycoUl05ZG1P2vtt7+dVp7SEHvM0bU9bZp79HJIQksvdAD5LRN8CcD6AQ53Sz226lem904aqmNB09GlDVa9nm8pJW14IfNvmq/tl0QmT9w4NVnDgaI1d125harNPezg6b9p5aZuHpntdc1dSIsKyjfcG04RNzxUA1t30YNszXPMtTdtdZfvW0VoPqj9HN1/y0tg5YYu3A/gRgLOIaD8RfZqIriGiaxqX3AfgWQBPA/gGgN8LXksNUreaODgFgRO61fiOCet7Gy4+C9VKuaWsaqWMDRef5fV8XTlZyguBT9tsfaTD93rXvYePTaNSblXrutVvJkxt3jS+i92e9StHsOGOnc3FHDih86p9l3Zebrj4LFRK7apnpUzGe21zV2VGCO+xtrF+pVlmTT6DM9/StN1Vtq1vTP1gqocQwIY7d2rnS4j+1OFc0IUQVwohFgshKkKI04UQfyeEuFkIcXPjfSGE+H0hxAohxCohREcMWmy6le29sdWjuPGjqzA6VAUBGB2q4saPrvL+11ItBwDKVB/QtOWFwKdtvrpfFp1Qd29tVmDBwLzM45Anpjbf/vDz7PZs3T3J0nnTzsux1aPYcsW5GKpWmq8ND1aw5fJzjffqnnXV2iXNv8u5nGx3Vk146+5J6/vqMzjzLU3bXWUnP9e2Oqr1OHm+XuyozQjtfMlLY++aOVdWsmhdY6tHgywcocoJCbdOvv2XRSc0XXNoqoZHP3+R8/5uYar3jCEPr649f/jtR9nlp51Pae6z3bNs473a17Nqwj5zhTvf0rTdVbYsc9nGe7XRHbr7D2rkNolpvuShsRd2QXdpjiF0ci5cLX9htQKi+uBn1aCzanC+mm2WvQfTvQsb36xcMdHrV47gnp0vajXoPLVJU73lmQPd9dwyTNf3AmnGmjMOHO1ePiPUXpeuniXm+PnUwdY2n/mSlcJ6udg0x1A6OQcfLf/gVA0HjtaCaNBZNTjfPsrSpyaN8cjxaa2+eOu2fW1/12nQeWuTpjZfef4Z7L5Iq/N2E9+x5s5Pl3avPiOPz7Csp25x1ZXtUwfbOPvMl6yUr7/++uCFcrjllluuv/rqq1Pfv3LxKTh9uIpdE4dw+Ng0RoequO7DZ2Ns9aj1vdB8+h+34/Wjx1tem54V2DVxCD985rW293TXffrfLEv9DNe9Jnz7KEufrlx8Cv7+oedwrDbb8vqsAB6feAPHZ2YNd5ox3Zu1X1RMbf699b/I7ouVi0/BkkWD2Pbsazg2Xa/r8GAFX/qN3tovUPEda+78TJY7VK2gOlDGW7XZtmfk8RnW1ROof4P+88vPaSvbpw62cfaZLxxuuOGGF6+//vpbdO+RMOg7ebNmzRrRDwkuTDqb/Lfa1bsE4LmbPpT6Ga57ewVTG/KgSP3SDxRlfhalni6I6CdCiDW69wqrofuQp+dKGi1fd12aZ5SImnGuPnD0ztDatClO26QvckirTYZoW694d2Txf+HE13NeB0H7zSWtRpyXz0xaXb4XPZtMFFZy4SJ1M/lTS867N49N4wc/m8Tpw1WsXHxK6vJPXTCAH/xsEtNKaFq1UsZ1Hz4b637x7W3vqcjrXM/XPUO2xbcNyf7Q9QPnGh/Gd0zg3sdeRLIbKmXCJz5wBp56+bCxj2z88opFePXwcW3fm+oZom2h+yctaethuu+lN6bwVw8+7fX6V7/3M8xohq5SJlz/kfd494eubo/tP4Q3j017tVGH7bPKnS+h14802CSXwm6Kcsnbc8UWP5x8b6hawfBgxTvuWpYTIj6YE9/bKc+NBQPzsHlslTUmenSoigUD+o20va9Necduh2hbr3h3ZPF/4cbXu17XjStQH9s0317z9JlJE+vfi55NNvpecumE54otFjZkzLspptmnDZz43k55bhxqRK64+sgWF+3bvyHa1iveHaH9X0zSl+/rwImx9SVvn5lQ8yVrPfKisAv6pvFduP3h5zEjBMpUDw3aPLaq7TqTdispEWHT+C6tRhcqhjwUIbxaOGVwn8PVb33qrSuT4/uSR11M2MpwzRkAbM1bneMSVbsN7f9i2pPwfR2oSxMrrr3PqDf7eqno2tAJQng2dXK/pZAa+qbxXbh1276mniUAPLb/EF49/BYuXPnO5nUm7VZF3pvU6JK64bHp2WbYXbf0Mx8N0KSTXnrOu9o062QZz716GI/tP9T2/MvOO63Zvz76LbfeujK//+TLmKrNtO25qRptHnWxYSrj0nPeZZ0z33/yZXzvyZebMfW2eibnuES9R7dHw2mLqf4fe//p2rlhe/3JF94wfr5MerNtvFz7Ttw2hsK0f8WtRx77LX2nod/+8POs103arQuTbpi8ptP6WQivlq27J51lmDw31Nd99FtuvU2eL7ohVDXaPOpiw1TG1t2T1jlTmxWoJXYQTfU0zXH1niz+L7r7dPsZrteTXiquOgN8LxXTnkonPX+yejZ1er+lkJILV9PLom9xQum6oZ+F8GpxlZGHzs6pt09/qhptHnVxoSvD5ttiQ1dP1/zL6ktkui/L666zBlyfll7zSMpSn07vtxRyQbdpd6pnMVeP832GROpneWlk3HLTauUmTJr1wmoF6256EC8cnDL6YYDQ4qEN8DVjn/FSvboHB8o4crz9m3EoH3guaeebrp6u+aeeQUi2Z3CghKdeOdK8dt2KRbjtMxc465G1X1ztl/PH1CrpOe6KO+d4/6SJWw85L2RZprbmtQdQSA198vAxrcYLuPW4SolQLpFVVzfphipSw9390pu5xCRztbcsWrnpubp9hxLVj3NL/df1+8VXMwb0eqVpvFR9NiljAOY46DxjyG16qwnTmNjmOHDiDIIuRvz1I63/GD9/YAqPPPcafvN9Z2hKqhOiX05dMIDvP/my9rOVnD8mXHHnnDj1NHHrIedFsqwkWfcA+k5D3zy2CletXaKNywbsetyWK87FlsvPdWp0Ujc0PUNquHlpZNxys2jlpufq9h0EoF04bfhoxoBe202Ol2k8kpjioPPUNJP1N1W1TOQcE9ccl/V27fVIHnrmdev7IfplbLXZFzzN/NHVgxOnbrvfRMh5Yatj3nsAhZRcgPqE3zy2yulZbNP+XNhiv6WGm5dGxi03i1bu89yQlj+2vnGNlykmPYkpDjpvTTOpKeuYFYLlHSLnuCxLNwRpbROShOoXky941mq69Hfu/b7vp3me6R4C8NDGC73L86GwC7rElgPwvBseaMYBc+LIddqczTvZx1s5ZLs4+wRqfL2PB0WWfQcuCxlREYDfeCTJw9fdl5DPCj0u3JyovnlFQ8W468q1lc+93/Z+3mPVidj5QkouKjaPZdV/3OVFrvN0vnXbPqN38vqVI17eyiHbxcmBOCNE01dc/j15r+9zk1RK1JZHk8OR49NOz3Kf8UiSl6+7LyGf5TMuOs58x4Lmn7k5XgH/vKK+HvK615O4fNJdcHOz5jlWncqXW8hNURXVs1huhHBI+jWbvJKTSO/k2x7e5+Wt7IurXbL+W644F6cPV/G9J15m29PaPMPlc03lqXG413/kPbjo7He1+Fsfn5l11mNWwOlZ7jMen1q7BK8dPp67r7svIZ/lGhcX1co861yfFcAp8ys49eSTcPjYNMpEbc/heM37esjrXr/svNOM46krP3m97X7femcZq7zm2JzxQ/f13FZ9kLn3yns66a3MeVaWtqd9ps993Sqn30jrK8+ZK77XRLrDnPBD3zS+y3uiq7G8XG1O1eZ1/xbqNMcsntXX3/24sV3qs1yeNbp2qFp8ktA+IXmXk/RQOT49g6OJDElqLtJukTXW2XecJfMrJecZgjR+PiHx6ZtO+qP0mqeTjcJr6MAJ3wtfZoRwatE6Dk7VjHHsSc0xbd7L8R0T2HDHTmvcrvoskwZqw1aXtDrg+pUjzuem1TRN5ehytyYXc+BELtJQOUd9yZobdnzHBA57yIoqU7XZ5nM5+z6d1oF9+iaPHLvceqXNC9wp+mJBt/leSP/xoWoFmhyuLA+JtNi8pDlxsT4+NLVZgQUD87zrG9rzxOQDIykTpfYbMXl6+MQm12ZE1zyss8Y6p/UmMmGLhw/heeODT9900h/FNbd6zRO9LyQXW+TDjusuav7Z5qsNmGOgs+TDNNUtD5/lQ1M1PPr5i7zrG9LzxFXvWSHYZWb1rQl1fSiyxjqHrrcrHr6Tnio+fdNJfxROmb3kid4X39BNp+mSr9tik21wdENuHUI+03SP770hddE82pX1mZ2oQ5bncuvDvY57mrZb/aDDp2+y9qMPnDJ7qR/7YkG/8ny9R8Xa5cMtf0+jC47vmMCRt+y6ZaVEOGle+4fIFn/L0ZArOo0IdV+MpF6ejNU13avjaCMufNP4Lqy49j4s3XgvVlx7HzaN72KXodbbpH1XSpRJgx3fMYF1Nz2IZRvvxXk3PIDVX3gAyzbei9ePvMUuo1LOVgefOq676UGM75hovjZxcArJUfHRpbn7CpzYbqCuAS/deC9Wf+EBlg6sa1sofD6bumsJvP2bEPVK4tuPeVL4OHQAuHDlO/HIc6/h+QOtP31ePXy8xVzHNz5UbogctizoQ9UKarOzeGu6VeQYHqzgi2O/ZIy/5cTFLlk0iG3PvoZj0yc2+IYHK7jxo+e0xH/rYnV195o4VpvFA4+/hJ37DzmThriQfax7drlEuOjsd6VONK0aHqnJI3RmWANlaktePDxYwZd+Iz8d2JSgQzUoA9Bc1H3jkzkx2K7Ybt2ZhmO1WXx/98tYsmgw1+Tavm0z9c3KxafgpTemsCthXvbUy4eDJ51J1muoWkF1oNyceyqcfgzBnIhDl9+AkowOVVP7J5jKVMsGEPy5IXG1wUWZCM/ceGmw56btF992dKP/ferYrflhq6OtTnl8vrLQ7fqk7ccQzIk49Dw2SrJsXPbKRknWeqQ1f+pUoulQ14cgRKLpvEk7Z3slMbbruZ2qT69+9gu5oGdJ6KC7F9AnYeAmiA1xAIOb9FoeNpI/4V2HZbKaOuk22DgHLWzmYksb0UauuqvP4ZpySbqxUeXT1/Mr7u0r1+GZNIdrbAeT1AQmyYNxeZnQpSVk4vA02Ma6m5ukrE1RIrqEiPYQ0dNEtFHz/kIi+l9EtJOIHiei3wlf1TqmQwXrV444N1V09264Yyc23LlTe0jBtiEiyw5xAEMejJIfGGmsldyU1B02ch2W4W6QmiZCcsOZe9CCs5lkq3vyOT6LeaeMkJL4HE6bqs1aN51dh2fSHK6xHUwqUd00LVmePBiXlwldGkztUDe88z58ZNuAzWNzlotzQSeiMoCvA/gggLMBXElEZycu+30ATwghzgXwKwC+TEQDgesKIFtCB1MSYlMSBk6C2BAHMLImvbYdlhlbPepM5Fsmwlc+fl5LQoUyEa5au6TtVwL3oEWy70yY6m56jjwMM1SttBwak3/udBJhFdlmbtig7UCc6/BMmsM1pvlDVP92rvscmBJocA+I5YGpHWkTh6fBdoDOdbguTziSywcAPC2EeBYAiOhbAC4D8IRyjQDwNiIiACcDeB1AujPKDvJIfmy7lnO4IusBjBBJrzlJI0wHjuRhn7HVo1qZh/uc5DWu59rKND2HmxyiW9iSoiSx/epwacRpNGTje8KcmMJUR58DYqExtSNL4vBQdQj5jDRwFvRRAOpXif0Azk9c8zUAdwN4AcDbAHxcCOGOl0uBS5sdHqzgQ+cs1iaN9dnec2nvvols5fs6rdxm/s9pu2SpchJWp0+79hk4miNHJ07qmC79WzVJk3A0Ug4hdFTfMrhaOlnKdo2V7f20ZfrstXDHwbQ3pNsL+tA5i3HPzhed+0Octptmmy65Spo5UmQNXff7MdlfFwN4FMBpAM4D8DUiagvEJKKriWg7EW2fnEz3s8SlUx44WmsmdlATI5g6X5ekgaO9J5/B0Tc/9Y0fabXy5SOD2rol9esNF59lHLDkgOj0aZvez9UcXf2v9p0s06V/qyZp8j6XRsohhI6apgwfLd1n/8aV6EFNvKIr06Ttrl854p04gqMTm/aGPvWNH2n3gm7dto+1P8Rpu4lkcpW0c8S0N5X3wTUXnAV9PwB1ZTkd9W/iKr8D4C5R52kAzwFYmSxICHGLEGKNEGLNyEi6jQOuNsvBlDSao70n4eibpkS9z04eZenXY6tHsXCQl74NaNenbXo/V3NMlmHTr31Ms5L959JIOYTQUdOUwZ2jusTJNrM4tW9N72/dPWmsr0nb3bp70vtzxdGJTXsEDz3zOttkTLfH4tN2V3lp54hub2p4sIItl5/bVTtdjuTyCIAziWgZgAkAnwDwycQ1+wD8KoD/S0TvBHAWgGdDVlSFq83aSCZsDaG9p01kOyNES0JgGyat01UniUnv99Ec8zbN4mikPuVlqVfaMrLMUe7+je59k37P0Xx96szpwzwTWPu03VZeljnSSfMyLs4FXQgxTUSfBXA/gDKAbwohHieiaxrv3wzgiwD+gYh2ob5W/rEQ4tU8KpwlNlnFppEnY6u5SQWkPpcm9lvq3wsGytbj6b5l2zRVTrlZNGvf8ZFaOjdBtosQOnzWfkkzF7JosGl08uTzOHXW7Xsk4SaAdsHtD+7nlJPIg5OoxrWP1g1YcehCiPuEEO8WQqwQQnyp8drNjcUcQogXhBAXCSFWCSF+SQhxax6VzRKbrOLSyJOx1dzkEW++NW3VP9etWOQs48jxGfzRHfa4cnYS5zJZNVVXub5xxlnHR2rpujMFEp8EISF0+Kz94jIrc+3f+GKrL7ctnDmW3PfQYTLNI0Cbm0AHd7y4yT84iTwAXqIa2z5atyiU26JJk1XnxvBgpS0RgikxgqtciZo8gmC2J52ZFVb987bPXNCilZuQ5egwJX5IxplLPc+mqbrK9Y0z9o0dNyUcUc8U6MiSIMRXh8/aL7azDJz9G19s9eW2RbdPYksOY2Lz2CosGGhfLAXqCaltZyOAenw8V5PmJP/Qxc4n26r7bNri8ZPXdTvZRaHMufJKXMvRDLkJmX3qstSQcMO3HBu9lsza9/rQyapjkuN05DkOIcbK9zOcpZysz8hK35hzmfQukweFxOV/wtEMVe8Nm1bno3/aNEaOBw2njTaNMKmB+vrEJOtl0sxLRNg0vqvtbIDpenU8fT1EXHHI3PFJ66PCOX8gGc1Zd3Vpvj4acJq4d9t9ajx41n0KzlkHoHXO+8bqc/cDpPau7sPJPx84WmuWk9fYF8oP/blXD+OxhAcyUJdEZARE0qdZ+p8cVX4uJX2LT10wgB/8bFLrqy2ZnhV49fBbeGNqGvc+9qI2SXSlTLj+I+9heyFPHj6mbU+5RLhBKcflRW1r46+fuxhPvXy4rW0C8O6nJMl6mXpPeqtLL+43j03jMcV7XaVSItRmZ5v/qOiuqVbKuO7DZ7fVKVkf7n2udun6W/f+S29M4a8efLrt9R89M4l7d73U1pbQnuKuNujGQP27rS66z0i1Usal57xL22ZZzqkLBvD9J19u+7zMCGH9/PmO1fEZ9zlGOedN43T6cBXrfvHt2rp87P2naz9HumcArX796p/l+1nG3uaHXigN3RT7OpPoZE5MsxqPmtTRTNz+8PNWPwzfGNTNY6tw1dolLc9cMFDGl69oLYfj62Fqo9SjTdqgTz8l8Ykz51Amwsnz57XFZsv3XDqzrT4++nRaHxVTQnDT+YNkuSFJMza2uqSJe5f3nTy/XQiwff6yjhVg3nh1JW431WXz2Cpr0nKuf4/umSEplOSSJn6Y67mgxpSatO0ZIax+GGl+PnHiz9P6esj3bP4ivv3EeT0ts0IY4+w5Hi6m+iTPHLhI29+hveOzkIfvfNq4d9OYmj5/PtjqaxoOV+J2U11sdTQloHcReuwLtaD7xPS6PC/Ua5LY9DKbTqxqZ756WRod0hW3Dbhj46VebVuKkv0k9eDQ2+m2WGlC+z+0yX7l6LXceqTxUUkbd62LeTbB9R1J64OfJuG29jmEpl7t6yEUStsH9HPJNE55xP9z7gtJoSQXrrc3J2GyLb7VFD8L2J0R1Rh29VpXjKrLT8IWKyvjtk39Ir0rTN4bb741bZ2IyX5K+nOEQo2V1rVFp5Dq+kl3b9K/w0VaH5Urzz/DKzm3JBnznHaeqKTx5E4TA2/qcyHQ9GHx9RDyie/29afJkrjdhq8XjlrPkBRqQTfpcYBZY03juWCKn82CTS9zabZS1zPp4Ft3T2LLFedCJ+NJrZK7/6Ci6yebh3cSdUy4ZwNsY6zDV6/lkNZHZfPYKq+6u9qTxMd3xNeTO62/uW28ZL/7egglsfUJ15/GpYeHjv9Xz1zY5PXQ3umFklwAsx5n01jT6HNHj4fb7JP4atFJjdGmV3J0ci62WFqfb+Zpvcuz+NVw9FoOaXxUbM/3Ics8cV1rIou/ua3NLl3a1yfJ9/U0enhaTGXa9PU5raEDYTxHOHlFub4QPkgtNxnvXSL9Bg7XY0PGApv0fQE/fdfmmeJTDje/Y/I9377n+GCbYuE5Wq1PzHlaLTXZHtPrtlwA6p5Cmj5MS1pvcG78eLIcztkH9SyDKe+tDU7eXC6d9E4vVBw6YI6H5cSsAvr43O8/+TK+9+TLzQX2zWPTmJkRKJdIG2+elhkh8Mqbx/DV7/2sJd6bG2vtigW2xeL6NsMUJ2uKnbfV3xbXvfulN9ve8+l7XSy0DlMsvCsO2zfm/NJz3sWKV7a1xzSXXecl1Lrt+PlB9tz1+fzoMMWa285l+MSP2+YS5yyDGgfOif9OPsP3/iSnLhjAA4+3n0MAgMvOOw0XrnwnqxxJ38ShA9m9Nbh5RZP+LaNDVauunvQr0VGbEbj94eeNnhOuWGufWOAQ6LRLGTtvg+OLLss2jUey7898xwLjc0K1X9de35hzlw9Nsv4unyEVOf6uNtjmmO8zOaTZp7L5/tjq5/ILGh2qGs8ySDi5Vzl5c7mMrTbnMZjzGjqQTf/y0awOTdXw6Ocvav7dpoUl9WKTJ4Tt5yVHc87iA50GXX9tHluF27btM/pvJGO+08SxJ/veRsj2J+vkG3Mu9zPGVpt9xbP4fYytHsUfONprqptvPL4Pvp9JW85Y27kMTq5ZTkx41pygvtp3qL0dF4Vc0CXcmNy0Ht1JLZmrhdm0QR//Fl39dXqej27rGys9v1LCupsexMTBqZa4eh+vcde1Ln8el4YZQrduosRP28o29aNabxM+cec6XGNoer9EhKUb7wXhhFSh8+tJfq6WnlrFtmcPtOUFVdHNUdtZjLR7YVnyiXKfw5lPrr0B7t7KnI5DV+HG5HI8unW+1JJknkdXTLstj6aMgTV1ui4m1eXV7opD1z3fJ1Z6qjbbnIhqXP0hwzeOZBtcvuS6+N1KiXDk+LS1za5cqbYxtaHGT+vaI1m7fNhZb9PCwo07N2E7J1EpkzbOWj4XiXol83bqPlcPPfN6W17QTeO7mmWY5qj6zGRbbfHjNpaeql8ABwdKznyiElfMuU/e3CSmdUnn8Z819l1HYRd0bkwuR3NTfal1qD4PLq3Q9jwZA+ujp3H0vK27J9nx+SFipQH9QR+gvQ0uX3LdvoCvBqorQ45pGo8NNW7dpHHufW3Ku96ctnCwnZNYMDCvLc7a1Qdqe7keMOp5BO49altt8eM2tj17QPv6U68cMdbBlvdWhy2m3HW/aV1S91ZC7V3oKKzkwtVlOZqbxKZ7cvM82p4n7/PR07Lqebp2hoiV5taFkx802adpNFDTuKTV112+LapWLumkn4fpnITsV7VuPv3JrYv6CzSkx5KrLN8TygSw92JUQnvL6OZLHhRuQXfpZDIm1+UH7hvra/MD4cbFSgYHyjii+UDq6sTR8xZWK1hw0jxjjPLqLzzQ1El9Yn/ToPpzXH/3405f8iyx3NLbxeTbPr5jAi2CsQcuDxyfvQIXurJkvyT3Llz7Oao+L+PsfTRlbt+r3/pDeiy5chv47gGF8r/n3Gf1DCJk2jfhUijJRdWnXEht0Fe7snme6LRObh5N6SeyaXyXdjEvl/TeMhz/miPHp615OGVfyLyI1g21FF4kEqmLS5nkmOwAACAASURBVG/1g1PmxW39yhEvvTGJTQcG0KxD2n+35Hhx83Byc1omMZWlznOdDs3JhSk9UVyoc4/rSaLq+Fw/EnWvibN3ottjMO0frFuxKLVG7eORY7vv4FTN/EVToCO5Rwu1oPt6PKt+4FztyuZ5otM6uXWSOqXJC2V2Vn/0muNt4vI9l9dw8iK+7aTW+G8ZEwyc+FZmeobUxTn5HbfunmTrjaqGaWufOj6uOshyhqoVrQc+x4dEhdNm9dmuskzjZPLtTrNXIHnbSSfyrOrau27Fopaxv2rtkpYoF24stZrPlbt3ktxjkOcgkvW57TMXBPVU5+atTXP+Ic/co4WSXNJojb7alW98qq9+aJOKTHA0b5efC8DTHznx3yZNVuq3WXV/25hxfTFs5ROAZ2681Fkmd9/E9bwkrvMGrrJ09Uqr3wOt+xnJcjlw2+56jmscJKYcAnno3mnuy/LMrBRqQU8TbxzM3xkn9D1V1/T1Nnnp0DHt9abvV1wd2FcDNSFj76X+qmury0+aU4cSEeaVCW9Nt8fLpI0RljHzrvMGcp+lTIS1y4eteUtdvj8y5tpH2dH5kqj97GJhtdJ2X9YdkbP/9Ls42jjebtuTSPqW37PzRfazF1YrOO+GB7Q5a117O2f+13vRqB5KBHzy/CXOxDBcfP3aOedSOM/Mg0JJLr6ewzbPcxOmOFcATX0PQEtcLgep55k0wFKJtDH0HB1Y1QpD+Csn9ddkW01tln3H0f1nhNAu5qa9BImtbBkzb9vLSNZBjbFWqVbKWL9ypE1b3XDHTmy4c6c25pqD7sxCsp9dvHGs1qyDz30mBNBczAHznoTOt9y2T6JSKRHeOFZruZ67tyMAKNXDrEBbLHwWfP3apQae9rOWR/y5pFALutTdOHqhy0vChCnOFYB3jLFOKzXFEOs0epsuSwStVhjaG8IH2Xe6eH0upr0EiSw7g2TsRJ4Z0HnE6Hx/TGXoUHXktBrsrEg3F4H6ngFn31u3J+Fb16S/im4qc/d2dPh489vw9WuXGrjPZ80nDj4LhZJcALsveBaPDEnIcD6TVmqKIebGcAMAhL6teWlzHNS+0+mjnJ7lXOPaK8iKPDOQ9hnUKEOHqiN3aqy4ewZJuHsSJrj+Kmk/cyE/q75+7b79seM6/1j4NBRqQXfFoIfQpdLmhdRh8uzgxja7vGPqPwcfw1Tj92iJgIF5Ja2U0SnU2PAPnbO46TfOjX0vN/KkJtuV1EyD+rckEADOu+GBujaeYirYPGrUPYq0MfK+qO3xkYfU+Zsmxl71QrKNV9rPnPoryBVHnsxBYNonSGKrN3d++NQzK4WRXFwx6KF0KZPGneYXvsmzgxvbbPOOWb9yBJ/79qPNRQ+o/xTv5mKucuBorSU3JPcDu3b5sLZdSc00bQ5HLgenalqJoFIiq2Sh5kY11U/qz2m/N6TxqfHV+oHW+Xv42HQqbxyX35DNe8aF/Ky64sh15yJ0+wQ6TOMoALbfPLeeISjMgm7T8ELqUro41wUDZdYXKakZmnJ/JnOEumJmbd4xW3dPGv1UbPUDTsSXJzVu+XkzxZ+HQNVVdfHNe1+bMrZL1Ux1faiOWx6UibDlinONJwLVvJxq/UKgarBJ76E82yzRedTr5pAONe+tbi7rcnyqnuhD1QoqykpVIrTEwrviyE17UZxcsz77dia49QwBiZyOgLtYs2aN2L59O/v6PLylsz7bVJdO1JVbp5DP9n1mmnq4nrHX0YYQdTSRdnyz1ok7dkszxKJnqQenfXl+Tl3jYasft15ZxnCvso8QYl0gop8IIdbo3iuMhs7xflB9sxdWKzg+PaONreXobWlyXLrisG0+2C7vjuQ1vqj6/KbxXbj94eebMeVrlw9j72tTLf2x/eevt1xz5fln5JovU33fprWqaHOzZqqdHTnXTM8oNfR/dY6F0vmXbbwX8yv1/ZFZAa0veaj9H5M2rGriALzap/aNii3Of+LglNG73SevqK1HuL7mWfyPZLs74YlemG/oUn9Sf7JUSvXYPW4IV6VM+Pj7z8A//2SipZxqpdz8qZz2Oa4y0l4vrwPgLNOG/Im6aXwXbt22z3ptCXp73DPfsQBPvXIk1fOB1jabGN8xgc99+1Ht89Wf2VIX5R63D0G5RJhxPK9aKeM33zfaNsfyQu0TzthycO3V+n7uJMnxT/NZM32GdbjGq1wifPkKfWgz5zPMxfYZ5nwmkti+obM0dCK6hIj2ENHTRLTRcM2vENGjRPQ4Ef2AXTsmaXyzk5hiXpN6GyfHpS3vIcdnw/XM5HVpY5YlMmaWE7tr0rDTLOYu35IkY6tH8ZWPn4eqIpomNVPAzzvFVT/ufoFrMQfM+UbzQh3P5P4PAVgwUHbmuk3iaqUrFp8A7cZxUi/m5vdted8jbt01XqqHTZKsnzcVk/9OHjHpTsmFiMoAvg7g1wDsB/AIEd0thHhCuWYIwF8DuEQIsY+I3hGshgoh/KdtuSDV/yfxyXEJ8Hw2uB7UIeKVZRl52eaa4ORJTRLaO8WGLodlVs27k32cfJbJ50SSt84uMXVB1th2IFz/Jr1lVEKfEfDxBcoCR0P/AICnhRDPAgARfQvAZQCeUK75JIC7hBD7AEAI8UroiuoIGYvs8r92+TQn4WioUlPP059cIvXPEvHDrYJAdSngnp0vescAS3T7C6H0YqntAgiilwJhzzJwniXRea3IcwByznaibi5//qUb78203xF67JN7VCHmQJK8vFuScCSXUQDq7/T9jddU3g1gmIj+DxH9hIh+S1cQEV1NRNuJaPvkZPYj6iFjkW3+1xyfZhWub7uM83VNHIEwm30TB6c6cpBFRTRiyNPEAANmv5NQH7YZIdr8WXzzzqrIvK15xsirrF0+DMDstZKcs/J6E5x2uq5x+fMD5i8Vrjj/SpmcbeAyI0Tzc8zNa6CDE9AYwmOJA2dB19U32dp5AN4H4EMALgbwp0T07rabhLhFCLFGCLFmZCR7A025/9Jg87/m+DSrhNTfQjMLoFoptcR/r1uxyCteOkTUMycGGMjel5zwYZNua8o7m5xrany/zNuqxqD7+qj4sPe1+j90nH6aqs1g72tTRp1d107dfpErBy/Hn1+HLNsU508EbLn83GabfTD1fZA9KkYTO+WxxJFc9gNQj0+eDuAFzTWvCiGOADhCRP8K4FwAPwtSSws6TSqtBmrSubg+za7Xe4VjtVmtru2jr4Y4tZ7VN52FSF9XU95ZDmn83CVyfeDU2TcX6AsHp5w6O8BrpysHr48fDgF4aOOFACx5YEU6Hx9CPa+o7+eYC+fLfKfWBM6C/giAM4loGYAJAJ9AXTNX+RcAXyOieQAGAJwP4KshK5rEFLe9fuVIav3LFCtr1OrpRJ5AVa90Pd+kAXZKezXpeT5+3CaN1AeppwJoi7lX65rlOTZfFe69XP8Nk1+4unfAOdPgMwN892FUDXtWtPe7+rlK/kM4WCnhpEq5ec7D5gPv2pPS1cvl+yKv8fWVcfVRljkC8D43sj+67uUihJgG8FkA9wN4EsB3hBCPE9E1RHRN45onAfxvAI8B+DGAvxVC/DRYLRPYci7eum1fJvc2nTZu8lRR8wRyfUtMGmsntVeTnmfysUly5Pi01Tc+DaZ9iSz7JNJXxVXXEqFNE+b4Yatw/MIPHK3h0FQtlSeKCe4+TBKpYavtSX6ukiUerc3iwNEaywde3ZNyeeNLXL4v8hrf3K22PpLjnEXjXj4y6LzmyPHppu97171chBD3CSHeLYRYIYT4UuO1m4UQNyvXbBFCnC2E+CUhxF8Eq6GGEBq1qiOr6LRxTl5PDtLrQ+ddodNeTWRdDkx6XjKO2URtRlh949Ni6nudHq3Tdk1nA1x1PWV+pU035vhhq3Dn5KxAy5mGtPs+pjGy+QmZCHXWQUXdk/Lxxpe+L7bPG/f8genfEd3ZiCwa97OTR53XcM7AhKAwR/9VQuhRx2pmaytd+Zy8ni6kzzZg1ljl6zbfh6zY+k/qq659iLykIV3dssbuuup6aKqW2Q/bZ07qzjT47vuY/Nal7u97RiMPjdeU+9TVzqx12evwYNftjWR5Jvez4DoDE4LCLeiu3INcbHeXiPCpb/wIP3zm9eZ1oX4kq37hyVhsTux6Vr0PMOuowIlYbNdOYlqvcE7dpJY6tnq0zXfmyvPPwJpfWNSsp+pv0lZHexOalIiwaXxXW8y2zX9DAFh+7b2YFfVv2T79Ift/lPEcHWUizK+UcESTKCVtblmB8Hs4ur0abr7ZU6rzvO1+JSuuvc96VkF3piSN37svrly8ISiMlwsQ1l+hF6iUqZkmj9u2q9YuwZpfWBTcx8TXn0P+nM3rkFK1UsZ7lyzEQ8+83vYex1MlxPNv/OgqbP/560H8UVzPAcAeU5OnjupN0u3PismjhFsv23yslCjT3E/On6zlcTB5/HTFy6VX6OX4bqBVm1P1XJOeqcZic9u2dfekV15NAlgaLTdXpmRW5HtGaao2o13MAZ6nSojn++aNzPIcndac1ICl74xJs1XzsSbPU3TAMr2JzaNEd3ZEp3Wr/klA695J1v2s5PyxLebqZ1qXC9h0j/p/3R5Z17xceolej+82+ZbY9Mw0ccSAXzzujusuysUnvEs/7jpGp+abr8+H6RdDcjg4fkKhUePJTXDPeZj8kzrVFoCfF1Xi8uvvBS+XniHPPJJcbDqjALD6Cw+0aeO2eucRr5vk7D/9bqdP/HuRlx6flU547AAnNF3dWOrixE2UG3sBt23b1xzvgTJhwUnzOjb+qh6syw17wfJFbd77tn2K1V94oJmbVp456eRUUfV2zvkSFY63TujFvVCSS955JF2US/VNOVtcrc6nxBbjavOQMcGJ15UIoJnkoxeplEi74Vwp1y0JukknFnPVJ0hHMk58wx07jWUtHxms5ypVXjs+I3Lf7JNUytTcXJe+9sncsA8983pbHLZtHqu5aYHsY6LzoTEtggS0eDi5nq2e4+B664SMQQcKtqBzfMZNmDyauSwYKOPLV9TzH7r066RPiU2H1XnIcODE6/Y6ZSKcPH+edmN1wcA83PaZC3DV2iW5PDvrfIByP3cayuuTerBr70KNEzfpvdVKiRUPLdF50WTJmyl9VuQ3zi3372HlvO3kPJY5YZNnDhaa9pg8ggSqlVKbXz/HWydkDDpQMMkFyKYLpv3HPamLcfRrH9/npIbK1bt7fU/BxawQxvh+6VW9eWxVi4QQkixf9tQ5wZ2Hz97ot7+i4hrrY7VZdh9Jb5O0ddEiWv1ffOZmp+Zx8hyIxNRun/mRPNfiuycWisIt6ADwqW/8yBgBYSLLgqDLA+rSrxdWKzjvhgdajn7b6pYq8QABg5WyNh65CMyvlHB8WhhjhSV5xAifNlTF0ePTqcotNzx/fHOG6nyCfNpmm8Np/F/UOT2+YyKT41oyttunXVnGwof5Fb0gEWJvLnl+gltmaJ/0QkkuQLrFPCs6zcvmD0IA3jhWYy3mmeolUNjFHACmarNGXfLNt6aburGvd4ekRHpZpVImrF85gkMpF5C1y4dZfvdJkpqpT9tC/kKR3iZyTm8a34UNd+xM/YulRGjLF8BtV7VSzjQWPkzVZrFpfFfb66G8ytU1grsnFtonvXALumsxp4RO6cOZ71hgvU/VvKz+IIGzAmXRNrM+M+mdosZJDw9W2Bq3bwtmZoVTN5Y6sIxnHlS+gQ0PVvCVj52Hr3zsvLY6b7n8XGzdPcnSeFVkHPje16ZSnYfQ5dR0HWjJe+xlDtQsB2sE2rXm2qww7i0kvVTSjEVadDl1Q541MOUPNY1j6HMOhZRcrIhWfZMrZRCA733uV5p/t/k8A/Yd75DBEaovRSfDtWaF0MbU6ny0bScpCUjlLQLY9UWTDqxDFxqWxlP7mRsvBZBNa+bsrcg+y/osLpltNAy3m15PntfwHYss6NraifyheXmxJ+m7BT2pSXH9KeZXSixdVGq7nfAuV/9V73QMvi3fYpv2auG0oWp6/x2qf/PWaasyRvng0VqzPgBw/d2Ps3KXpvE6WXfTg5n89mU5Lq8R4MQXkVKPxuir+H4WVL0Z6Fy8v0T14Dk+PcP+osRtp9wTU+efab6F1tAL5+XyB5Z/zXW+CJvGd7G9ODgeIdJ/JW+PD6Du2yK/EcsY5Lw9J1RUn5Gk/4bJm0KlXCJc+YEzrNe44Ca1rpTqH7bktapfjko3+rMfyTLGvv5B3aRcIqxdNuy9fyfnH6D/DKU5+t9XXi4mTL4I0uObA8cjRMaNy3JVZSyk2qku5gBYXh9JuPqr6SqbT7bUXk0fYhm3v3X3ZCb/He56W5ttX8wBc+5S2Z+9QBoPc1M5neZtJ81j+/gn8fUPUul0S2dnRapcprZcxaF9XICCfUO3eYTrPFQ496bB9LxQz+C0x/Y8Vbd21cd1rU9+y2SZtjp2Elt/pgoXDYzPeLnKCTX/fGLan0vE5GetQ6h2hCZtvbifZ3Z5lm/ohVrQXX4XSf1UekBIf+/QTR01aLehkD7bB47Wmvqd1P1cx/m5ep/Mw2qSj4aqFbx5bFpbls2DRfaNb6x2Xsh2Jr00/ug7Ozuq39rq14lYbA5cmUuFAAwOhDkTYZtzrLr02L6D+g+BbV+HXV6/LOg6L2Wp5X77x9lCr4B0E9mk3RYBm+c44NY4CcC8Mhnf5+jsLiplwtJTB7X+31lxtb/TFHkuhaRSJnxgqb9eXRRM+zpc+kZDN+lQW3dPBtncOmV+xVuHNGm3vYgaty37zhZP7/IZEUCLZ3US6dORHLOr1i6B4dBeCzJm/OjxfKKUp2oz2PvaVEse1TK5TcGk5p02H6iJ2qzAKfMr2lh/9TVO3xWZ2oxoG5dOk9Xnx4ZpXycEhfqGbiKkdg30pn6XFZOOZ9OQOZohR4O36Yec+/LU4dP2Sxbd2DbPsu6f+NQhhG6fF93ch+nEs7Po6rZv6IWNQ1fjokPFsXY6HrajkP8GIKcnXD40Aq0fisFKCSdVys34cZPuqsYq5xmDn/Tv5mj+yXt8583CagWHjtW0HSzj1NcuH27zDVd/omf1t8l0PiAAnD0eec6g03VM7a3kgepVFJJC/nhLeg2HGuy+XczR3U0i9dFHa7M4cLTW9PywbaKpftl5+eBLTx51Ttkol1o9v6+9a5fXvKmUCG8cq1nHY0YIrW+4PMQ1vmMis/fJ9MyMd9190PmOqywfGbTeL4DmPAlZR1udOonMgxCaQi7oaXKLJr0+Biulpp7cKZ1Oaq/VfhdBA6LT4UOOl9xD4MwpGVuven77zEOb/7sL1QeG6zVu4+U3jwfNz0uo94/cJ5G+46ax8vFuDwUBTS/0bpOXjl5IycXX/8Dl+xHKL8OlSybzE4aOWe+FmOo8eOHgVG75MeW3P86cevwLl7TVyweb/zsH3/yzncKmB5t8Wrrxa1jAP+dAnuQxjoVb0NPofjq/BFWDD3WSQT7HliNRLrq+3zE5vh/9SlKHD+1vwu0/6QFSpnoqQpOObRqrrFWW/jUD80p4a7q30grqfHW23L+n64umSje9kXSE9nEBCia5uDRL3SKp5jlMliM1yhCLg/ocjheyzyOrlbJRc+ylD0yeqO3sVpiofO6MELh12z7jpqRLH87CgaO1nlvMVb174uAUNtyxExvu3Nn1BTOJmvMztA+5L9VKuW1dCkGhFnSXZrmw2h7DqwvgT6PBu1gwMK8lvdWNH12VSetNekZ3Q3OMpGOuj1UWj5a8SHojZfUhHx2qZkpinoePC1AwycWlOR2aqrFiOznaldQFuT/HDyWO/XPyjtpIekbbXCY7Sa/6bKh0u479HC1VRAjtPv5Z9GsC8NDGCwGk3wvLYzEHCnawyOTlIikT4csfa/1GrktZx4mBdXmccAixsKSxI5jrDOeQg3Qu0Qmv/26g+jtllYNUn6U0JH8x+NA3R/9d2vSMEC3xuqb8o67JWi4Rlp5azex3HuIjERdzf+Jinp5qpdyiNfcTB47WcOu2fUG0/YNTtUzz7NZt+7T5TbPCWtCJ6BIi2kNETxPRRst17yeiGSK6PFwVT8DRptV43bTmPrOzwp4ztEdRY+uHqhWrH4X0VOmWV0YWhgfDeqjkRRF6Vufvs3lsVa5eJjoGKyXv/irg1G1Bl980K04NnYjKAL4O4NcA7AfwCBHdLYR4QnPdnwG4P3gtFTjadNb4ztCn07hkkWgIwBNf/GDLa6Z4bVUDvC3nrEuhIQA7rruoI7k2Q7C3x/1STOczOvXLcG/CE4cLAb2/meMgjzWGsyn6AQBPCyGeBQAi+haAywA8kbjuPwH4ZwDvD1rDBJvGdznHcX6lhHU3PZjpOd3QrrM+Tt3ALZG5vE54VeRGCk+ablEiwrKN9/asR5Cs32mN/aJ7dr4Y3NPfhZq31sefRqAYv4Bs5PHrmLOgjwJQfxvsB3C+egERjQL4DQAXwrKgE9HVAK4GgCVLeGnhVLj5Qadqs5l1sqJp18nqdqL+8md5J/uqB9dFI3IR78XFHDhRr4mDU7nnxzVx7V0ndOTDx6a97u3NXuWTx14FR0PX/TOS7Mu/APDHQghrcLcQ4hYhxBohxJqREf/A/jw0p0h6Tplf8XaNUz3R5xqqj3rRv12GQs1bm3fCbukt3wteLtVKKXWUiw3ON/T9ANR/Sk4H8ELimjUAvkX1nxBvB3ApEU0LIcaD1LJBr37TAXovj2YnSMbec1Bj67v1rbBbJL18InVse16hchTsTZxP6fZn9JgjhWRaOAv6IwDOJKJlACYAfALAJ9ULhBDL5J+J6B8A3BN6MQd6Oz5W6tI27brfEPDzVUlqhr08njqyniuQ3h3jOyZyO/1UtD4F7N0gAJyU0buGUO9zNe9vaD8gX1S//5CHjJySixBiGsBnUY9eeRLAd4QQjxPRNUR0TbCaMChCfGzRtPes+LQ3OX5FGE+VrEO7fuUIxndMYMMdO3NbTNYuH67ngu0jsnrXCACf+86jLRu+vfA5Tfrch6BQJ0WB+sbo7Q8/X7hvIXOdaqWEJxNhlQB/oztPOpUlXmq3eZpWjQ5VcfT4tPehl27bJcxlRoeqzTBiDn1zUhSoezI8c+OlbZqYDkLxQ5s6AaGuMe696UO59ZdJM9w8tqqrY9TJeOYXDk7l7mX+wsGpVJ7rvjlR85wrc42Qc6Jw39BVH/O58I2iU98eJXnG31Oj/B4z4ot+OSkoolbfq4T8hl4ot0XpYx7a+raXIdQX9U4tOHk+p34CN7/y0xIXc3/iYh6G0L7ohZJc8vAx73VmRdQ2u0mf7S9GeoCkd07IKJdCfUPvtVyKnSJ+GeoOhNj3RSRU7Lqt/CyeSzuuM+c3zkqhvqHnkYMvErFRdEe/OQnVfWHywvd0dAuEoGGKSQq1oHc7D2BkbiEQ9fUiIkR+nviVEuHIcT/PGRUhgA137sxtUS/Ugp41D2ARiF8IwxK/YUdCMTpUxcnz52XOl1qbEc2cDaGJGnqPwZkqeWuE/UInYsxVDx+g+x4hRaMoc3lvYA+evNayQi3oPn7J/cxpHThx2Bd0IIa/0N7yPUCvL+SS0GNcImrxgg9WbtDScmR8xwQO9chi3u1f8etXjmDDxWf19OCd+Y4F3a5CjFCJ9CzJ/Meh6OU1oYUt9+9BPoaTfAj1bN1f/fh5XdVmt+6exNjqUXzl4+ehWunNITx6fLawOUs7QeyXiJr/OBSFkVy6qZ8ndVIAzrymeSL7Ymz1KMZWj+am22aJt33h4BQ2j61qmvhHbbmV2YL8fNg7xzz+O03oda0wXi7rbnqwq5qx9PuQHhZzwctiOO5ZRCK5MlStGBN1m+gLt8UNF5/VVZ9nGY/c63kiQ+Kb4zESifjx5lvTQXX0wizoY6tHseWKc2NccQepzYrY35FIjszMho1JL4yGDtQX9W5q13OROfBDJJIR3R5TDOXkE1JHL9SCDtRjsGP8dSTSO8RY/GyE9KgqjOQi6fX460gkEuFSKdPc9UMH6rJL10/2RCIF5aq1S7pdhYjClsvPDXpatHALOhAd8CKRNBDQPBcQ6Q3m7NF/lXjKLlIkemW+xnwCvUUe86KQC/qV55/R7SpEImzyOLOQZimYODiFFdfeF7wukXSsXT4cvMxCLuibx1bhqrVLopTeoKj90CvfXItI2n8i5sKBuKLw//YdCm7OVbiwRcnmsVXYuntyzocwjg5V8dDGCwF03x5B1gdwW/uq9V5x7X1xodEwF+wleg11Xkry+lxJc645myQ6yVxIeOFC7YNe6A9uHdTr4qKlJ/ZL59HN3zw/V6HLLuSCvml8F25/+Pno/oa6Uf6m8V3YunuyJ/pDgCcBLaxWsGl8F27dti/vKkUibDp9SCr0RnXhFvS4CLQyI0TP9QfnH5aDU7Weq3ck0mlCJ74v3Kbo7Q8/3+0qRJjELc9IxE7oxPeFW9Cjrlgs4qIeiZgJraEXbkHvoiV6xJP4T28kYmdhtRK0PNaCTkSXENEeInqaiDZq3v8UET3W+O+HRHRu0Fo2GN8xEY/9F4w4XJGImSPHO5zggojKAL4O4IMAzgZwJRGdnbjsOQD/TghxDoAvArglWA0VQidUjUQikW5Sm+l8gosPAHhaCPEsABDRtwBcBuAJeYEQ4ofK9dsAnB6shgq9EGcdifQqyUQTMbFzMeh0gotRAGpoyX4A51uu/zSA72aplImY3CISMSNjqGVC80gx6HSCC902pHa6ENF61Bf0Pza8fzURbSei7ZOT/uE6S0+NbnGRiIu4mBeLkLHonAV9PwDV3vB0AC8kLyKicwD8LYDLhBCv6QoSQtwihFgjhFgzMuLfiG3PHvC6PgbERCKRXidkLDpnQX8EwJlEtIyIBgB8AsDd6gVEtATAXQD+gxDiZ8Fql8A3Bv25mz4UF/VIJNLTdFRDF0JME9FnAdwPoAzgm0KIx4nomsb7NwO4DsCpAP6a6pao00KINcFq2cDXfS4m/MaJKwAAB+9JREFUro30AoQYvhkxE1JDZ3m5CCHuA3Bf4rWblT//LoDfDVYrA1eef0b0/4gUjriYR0yEThJdKHMumQ8xLuqRSKToDA9W8PkPv2duJ4nePLYq6uKRSKTw7LjuopgkenzHRAxfiUQihSd0+jmgYAv6+I4JbLhjJ6LhYiQSKTob7twZfFEv1IK+5f49qMVTE5FIpA8I7eMCFGxBj14ukUikn5jTfuih8+9FIpFINwm9phVqQY9eLpFIpJ+Y0zlFfb1cIpFIpJeZ0zlFO5VPdG/0gCkkccwiRWNOa+hl6sxHdmlMDFBI4phFisac1tCvPP8M90WRSCRSEELvCxZqQd88tgpXrV3S7WpEIpFIEELvCxZqQQdOGHRF5g5RG4/0K6H3BQu3oOfhfxDpbaI2HulXQu8LFmpBH98xgc99+9FuVyMSiUSCsHxkMGh5hVrQt9y/B7PdrkQkEokE4tnJo0HLK9SCHr1cIpFIPzGnNfTo5RKJRPqJ0Bv+hVrQBwcKVd1IJBKxIhA20KNQK+RTrxzpdhUikUgkKCE90Qu1oEcikUi/EXJvMC7okUgk0kVC7g0WakGfX45nBiORSH+x4eKzgpVVqAX92Ew8MxiJRPoHAjC2ejRYeYVa0CORSKSfCP0VNS7okUgk0iXmrJdLNOWKRCL9xttPrgQtrzALeshYzUgkEukFXn7zeNDyCrOgRx+XSCQSsVOYBT36uEQikYidwizoIWM1I5FIpFf41Dd+FKws1oJORJcQ0R4iepqINmreJyL6y8b7jxHRe4PVsEHIWM1IJBLpFR565vVgZTkXdCIqA/g6gA8COBvAlUR0duKyDwI4s/Hf1QD+JlgNI5FIJMKC8w39AwCeFkI8K4Q4DuBbAC5LXHMZgH8SdbYBGCKixYHrGolEIhELnAV9FMDzyt/3N17zvQZEdDURbSei7ZOTk751xboVi7zviUQikV4m5LrGWdB1R5mSJ1Y510AIcYsQYo0QYs3IyAinfi3c9pkL4qIeiUT6hnUrFuG2z1wQrLx5jGv2AzhD+fvpAF5IcU0QQjY+EolE+gnON/RHAJxJRMuIaADAJwDcnbjmbgC/1Yh2WQvgkBDixcB1jUQikYgF5zd0IcQ0EX0WwP0AygC+KYR4nIiuabx/M4D7AFwK4GkARwH8Tn5VjkQikYgOjuQCIcR9qC/a6ms3K38WAH4/bNUikUgk4kNhTopGIpFIxE5c0CORSKRPoLpa0oUHE00C+HnK298O4NWA1SkCsc1zg9jmuUGWNv+CEEIb9921BT0LRLRdCLGm2/XoJLHNc4PY5rlBXm2OkkskEon0CXFBj0QikT6hqAv6Ld2uQBeIbZ4bxDbPDXJpcyE19EgkEom0U9Rv6JFIJBJJEBf0SCQS6RMKt6C70uEVBSI6g4i2EtGTRPQ4Ef3nxuuLiOh7RPRU4//Dyj3XNtq9h4guVl5/HxHtarz3l0SkszPuGYioTEQ7iOiext/7us1ENEREdxLR7sZ4XzAH2vyHjXn9UyK6nYjm91ubieibRPQKEf1UeS1YG4noJCL6duP1h4loqbNSQojC/Ie6OdgzAJYDGACwE8DZ3a5XyrYsBvDexp/fBuBnqKf4+3MAGxuvbwTwZ40/n91o70kAljX6odx478cALkDdl/67AD7Y7fY52v45AP8TwD2Nv/d1mwH8I4Dfbfx5AMBQP7cZ9eQ2zwGoNv7+HQD/sd/aDODfAngvgJ8qrwVrI4DfA3Bz48+fAPBtZ5263SmeHXgBgPuVv18L4Npu1ytQ2/4FwK8B2ANgceO1xQD26NqKuvvlBY1rdiuvXwngv3e7PZZ2ng7g+wAuVBb0vm0zgFMaixslXu/nNssMZotQNwC8B8BF/dhmAEsTC3qwNsprGn+eh/rJUrLVp2iSCyvVXdFo/JRaDeBhAO8UDS/5xv/f0bjM1PbRxp+Tr/cqfwHgvwCYVV7r5zYvBzAJ4O8bMtPfEtEC9HGbhRATAP4bgH0AXkQ9P8ID6OM2K4RsY/MeIcQ0gEMATrU9vGgLOivVXZEgopMB/DOAPxBCvGG7VPOasLzecxDRrwN4RQjxE+4tmtcK1WbUv1m9F8DfCCFWAziC+k9xE4Vvc0M3vgx1aeE0AAuI6CrbLZrXCtVmBmna6N3+oi3oHUt11wmIqIL6Yn6bEOKuxssvE9HixvuLAbzSeN3U9v2NPydf70XWAfgIEe0F8C0AFxLRrejvNu8HsF8I8XDj73eivsD3c5v/PYDnhBCTQogagLsA/DL6u82SkG1s3kNE8wAsBPC67eFFW9A56fAKQWMn++8APCmE+Iry1t0Afrvx599GXVuXr3+isfO9DMCZAH7c+Fn3JhGtbZT5W8o9PYUQ4lohxOlCiKWoj92DQoir0N9tfgnA80R0VuOlXwXwBPq4zahLLWuJaLBR118F8CT6u82SkG1Uy7oc9c+L/RdKtzcVUmxCXIp6RMgzAP6k2/XJ0I5/g/rPp8cAPNr471LUNbLvA3iq8f9Fyj1/0mj3Hii7/QDWAPhp472vwbFx0gv/AfgVnNgU7es2AzgPwPbGWI8DGJ4Dbb4BwO5Gff8H6tEdfdVmALejvkdQQ/3b9KdDthHAfAB3oJ7a88cAlrvqFI/+RyKRSJ9QNMklEolEIgbigh6JRCJ9QlzQI5FIpE+IC3okEon0CXFBj0QikT4hLuiRSCTSJ8QFPRKJRPqE/w8d8jZ2ypXJUgAAAABJRU5ErkJggg==\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"npoints = 10000\n",
"spills = np.zeros(npoints)\n",
"for ii in range(npoints):\n",
" spills[ii] = ros.ChoseFractionSpilled(random_generator)\n",
" if ii < 10:\n",
" print(spills[ii])\n",
"plt.plot(spills, 'o');"
]
},
{
"cell_type": "code",
"execution_count": 20,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAX0AAAD4CAYAAAAAczaOAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+j8jraAAAQW0lEQVR4nO3df6zddX3H8edrVAmbMutaSNfWlZm6CWSi3HXN3BYc2aj4RzHRpGwRYkjqGC6a+IfFP6bJ0qQmUze2gKlKgERpmomji+DGmI4ZQbw1lVK6aicMrm3oVZfJXMLS8t4f59vlrJz2nvvr3J5+no/k5HzP+/v9nvP+5Dav8z2f8z3fpqqQJLXhZ5a6AUnS6Bj6ktQQQ1+SGmLoS1JDDH1JasiypW5gJitWrKh169YtdRuSNFb27t37w6paeWr9rA/9devWMTk5udRtSNJYSfLvg+pO70hSQwx9SWqIoS9JDZkx9JOsTfLVJAeTHEjyga7+sSQ/SLKvu13bt8+tSQ4nOZTkmr76lUn2d+tuS5LFGZYkaZBhvsg9Dnyoqr6d5NXA3iQPdes+VVV/3r9xkkuBLcBlwC8C/5jkDVV1ArgD2Ao8BjwAbAIeXJihSJJmMuORflUdrapvd8svAAeB1WfYZTOwq6perKqngcPAhiSrgAur6tHqXeXtHuC6eY9AkjS0Wc3pJ1kHvBn4Zld6f5InktyZZHlXWw0817fbVFdb3S2fWh/0OluTTCaZnJ6enk2LkqQzGDr0k7wK+CLwwar6Cb2pmtcDVwBHgU+c3HTA7nWG+suLVTuraqKqJlaufNlvCyRJczRU6Cd5Bb3A/3xV3QdQVc9X1Ymqegn4DLCh23wKWNu3+xrgSFdfM6AuSRqRGb/I7c6w+RxwsKo+2VdfVVVHu4fvBJ7slvcAX0jySXpf5K4HHq+qE0leSLKR3vTQDcBfLdxQXm7dti8PrD+z4x2L+bKSdNYa5uydtwLvAfYn2dfVPgJcn+QKelM0zwDvA6iqA0l2A0/RO/Pnlu7MHYCbgbuAC+idteOZO5I0QjOGflV9ncHz8Q+cYZ/twPYB9Ung8tk0KElaOP4iV5IaYuhLUkMMfUlqiKEvSQ0x9CWpIYa+JDXE0Jekhhj6ktQQQ1+SGmLoS1JDDH1JaoihL0kNMfQlqSGGviQ1xNCXpIYY+pLUEENfkhpi6EtSQwx9SWqIoS9JDTH0Jakhhr4kNcTQl6SGGPqS1BBDX5IaYuhLUkMMfUlqiKEvSQ0x9CWpIYa+JDXE0Jekhhj6ktQQQ1+SGjJj6CdZm+SrSQ4mOZDkA139tUkeSvK97n553z63Jjmc5FCSa/rqVybZ3627LUkWZ1iSpEGGOdI/Dnyoqt4IbARuSXIpsA14uKrWAw93j+nWbQEuAzYBtyc5r3uuO4CtwPrutmkBxyJJmsGMoV9VR6vq293yC8BBYDWwGbi72+xu4LpueTOwq6perKqngcPAhiSrgAur6tGqKuCevn0kSSMwqzn9JOuANwPfBC6uqqPQe2MALuo2Ww0817fbVFdb3S2fWh/0OluTTCaZnJ6enk2LkqQzGDr0k7wK+CLwwar6yZk2HVCrM9RfXqzaWVUTVTWxcuXKYVuUJM1gqNBP8gp6gf/5qrqvKz/fTdnQ3R/r6lPA2r7d1wBHuvqaAXVJ0ogMc/ZOgM8BB6vqk32r9gA3dss3Avf31bckOT/JJfS+sH28mwJ6IcnG7jlv6NtHkjQCy4bY5q3Ae4D9SfZ1tY8AO4DdSW4CngXeDVBVB5LsBp6id+bPLVV1otvvZuAu4ALgwe4mSRqRGUO/qr7O4Pl4gKtPs892YPuA+iRw+WwalCQtHH+RK0kNMfQlqSGGviQ1xNCXpIYY+pLUEENfkhpi6EtSQwx9SWqIoS9JDTH0Jakhhr4kNcTQl6SGGPqS1BBDX5IaYuhLUkMMfUlqiKEvSQ0x9CWpIYa+JDXE0Jekhhj6ktQQQ1+SGmLoS1JDDH1JaoihL0kNMfQlqSGGviQ1xNCXpIYY+pLUEENfkhpi6EtSQwx9SWqIoS9JDZkx9JPcmeRYkif7ah9L8oMk+7rbtX3rbk1yOMmhJNf01a9Msr9bd1uSLPxwJElnMsyR/l3ApgH1T1XVFd3tAYAklwJbgMu6fW5Pcl63/R3AVmB9dxv0nJKkRTRj6FfVI8CPh3y+zcCuqnqxqp4GDgMbkqwCLqyqR6uqgHuA6+batCRpbuYzp//+JE900z/Lu9pq4Lm+baa62upu+dT6QEm2JplMMjk9PT2PFiVJ/eYa+ncArweuAI4Cn+jqg+bp6wz1gapqZ1VNVNXEypUr59iiJOlUcwr9qnq+qk5U1UvAZ4AN3aopYG3fpmuAI119zYC6JGmE5hT63Rz9Se8ETp7ZswfYkuT8JJfQ+8L28ao6CryQZGN31s4NwP3z6FuSNAfLZtogyb3AVcCKJFPAR4GrklxBb4rmGeB9AFV1IMlu4CngOHBLVZ3onupmemcCXQA82N0kSSM0Y+hX1fUDyp87w/bbge0D6pPA5bPqTpK0oPxFriQ1xNCXpIYY+pLUEENfkhpi6EtSQwx9SWqIoS9JDTH0Jakhhr4kNcTQl6SGGPqS1BBDX5IaYuhLUkMMfUlqiKEvSQ0x9CWpIYa+JDXE0Jekhhj6ktQQQ1+SGmLoS1JDDH1JaoihL0kNMfQlqSGGviQ1ZNlSN7AU1m378sD6MzveMeJOJGm0PNKXpIYY+pLUEENfkhpi6EtSQwx9SWqIoS9JDTH0JakhM4Z+kjuTHEvyZF/ttUkeSvK97n5537pbkxxOcijJNX31K5Ps79bdliQLPxxJ0pkMc6R/F7DplNo24OGqWg883D0myaXAFuCybp/bk5zX7XMHsBVY391OfU5J0iKbMfSr6hHgx6eUNwN3d8t3A9f11XdV1YtV9TRwGNiQZBVwYVU9WlUF3NO3jyRpROY6p39xVR0F6O4v6uqrgef6tpvqaqu75VPrAyXZmmQyyeT09PQcW5QknWqhv8gdNE9fZ6gPVFU7q2qiqiZWrly5YM1JUuvmGvrPd1M2dPfHuvoUsLZvuzXAka6+ZkBdkjRCcw39PcCN3fKNwP199S1Jzk9yCb0vbB/vpoBeSLKxO2vnhr59JEkjMuOllZPcC1wFrEgyBXwU2AHsTnIT8CzwboCqOpBkN/AUcBy4papOdE91M70zgS4AHuxukqQRmjH0q+r606y6+jTbbwe2D6hPApfPqjtJ0oLyF7mS1BBDX5IaYuhLUkMMfUlqiKEvSQ0x9CWpIYa+JDXE0Jekhhj6ktQQQ1+SGmLoS1JDDH1JaoihL0kNMfQlqSGGviQ1xNCXpIYY+pLUEENfkhpi6EtSQwx9SWqIoS9JDVm21A2cTdZt+/LA+jM73jHiTiRpcXikL0kNMfQlqSGGviQ1xNCXpIYY+pLUEENfkhpi6EtSQwx9SWqIoS9JDTH0Jakhhr4kNWReoZ/kmST7k+xLMtnVXpvkoSTf6+6X921/a5LDSQ4luWa+zUuSZmchjvTfVlVXVNVE93gb8HBVrQce7h6T5FJgC3AZsAm4Pcl5C/D6kqQhLcZVNjcDV3XLdwNfAz7c1XdV1YvA00kOAxuARxehhwXl1TclnSvme6RfwD8k2Ztka1e7uKqOAnT3F3X11cBzfftOdbWXSbI1yWSSyenp6Xm2KEk6ab5H+m+tqiNJLgIeSvKvZ9g2A2o1aMOq2gnsBJiYmBi4jSRp9uZ1pF9VR7r7Y8CX6E3XPJ9kFUB3f6zbfApY27f7GuDIfF5fkjQ7cw79JD+X5NUnl4HfB54E9gA3dpvdCNzfLe8BtiQ5P8klwHrg8bm+viRp9uYzvXMx8KUkJ5/nC1X1lSTfAnYnuQl4Fng3QFUdSLIbeAo4DtxSVSfm1b0kaVbmHPpV9X3gTQPqPwKuPs0+24Htc31NSdL8+ItcSWqIoS9JDTH0Jakhi/GL3Gb4S11J48YjfUlqiKEvSQ1xemcROO0j6Wzlkb4kNcTQl6SGGPqS1BDn9M8Sfg8gaRQM/RE6XbBL0qg4vSNJDfFI/yzntI+kheSRviQ1xCP9c4yfDCSdiaE/pmb7pbBvBpLA6R1JaopH+poVPzFI483Q10D+pkA6Nxn6jTPcpbY4py9JDfFIXwvCuX5pPBj6WhK+SUhLw9DXovI7A+nsYujrrDKXNwk/HUjDM/Q19pwqkoZn6OucNdtPDb5JqAWGvjQDP0noXGLoS52FuojdbPnmoVEy9KUlNpdPEn760FwZ+tJZaiFPd/VNQicZ+tI5ZLGnqHyTGH8jD/0km4C/BM4DPltVO0bdg6S5WezvMRbyE8lsn6uVT0OpqtG9WHIe8F3g94Ap4FvA9VX11On2mZiYqMnJyTm9nr8GlTRqs31Tme3zDCvJ3qqaOLU+6iP9DcDhqvp+19QuYDNw2tCXpHFyth9sjjr0VwPP9T2eAn7j1I2SbAW2dg//K8mhWbzGCuCHc+5wfDnutjjuc1w+/v8ezmXcvzSoOOrQz4Day+aXqmonsHNOL5BMDvpIc65z3G1x3G1ZyHGP+j9RmQLW9j1eAxwZcQ+S1KxRh/63gPVJLknySmALsGfEPUhSs0Y6vVNVx5O8H/h7eqds3llVBxb4ZeY0LXQOcNxtcdxtWbBxj/SUTUnS0vI/Rpekhhj6ktSQsQz9JJuSHEpyOMm2AeuT5LZu/RNJ3rIUfS6GIcb+h92Yn0jyjSRvWoo+F9pM4+7b7teTnEjyrlH2t1iGGXeSq5LsS3IgyT+PusfFMMS/859P8ndJvtON+71L0edCS3JnkmNJnjzN+vlnW1WN1Y3eF8D/Bvwy8ErgO8Clp2xzLfAgvd8FbAS+udR9j3Dsvwks75bffi6MfZhx9233T8ADwLuWuu8R/b1fQ+8X7a/rHl+01H2PaNwfAT7eLa8Efgy8cql7X4Cx/w7wFuDJ06yfd7aN45H+/13Koar+Bzh5KYd+m4F7qucx4DVJVo260UUw49ir6htV9R/dw8fo/RZi3A3zNwf4E+CLwLFRNreIhhn3HwD3VdWzAFV1Lox9mHEX8OokAV5FL/SPj7bNhVdVj9Aby+nMO9vGMfQHXcph9Ry2GUezHddN9I4Kxt2M406yGngn8OkR9rXYhvl7vwFYnuRrSfYmuWFk3S2eYcb918Ab6f24cz/wgap6aTTtLal5Z9s4Xk9/mEs5DHW5hzE09LiSvI1e6P/WonY0GsOM+y+AD1fVid7B3zlhmHEvA64ErgYuAB5N8lhVfXexm1tEw4z7GmAf8LvA64GHkvxLVf1ksZtbYvPOtnEM/WEu5XCuXu5hqHEl+TXgs8Dbq+pHI+ptMQ0z7glgVxf4K4Brkxyvqr8dTYuLYth/6z+sqp8CP03yCPAmepcwH1fDjPu9wI7qTXQfTvI08KvA46NpccnMO9vGcXpnmEs57AFu6L7p3gj8Z1UdHXWji2DGsSd5HXAf8J4xP9rrN+O4q+qSqlpXVeuAvwH+eMwDH4b7t34/8NtJliX5WXpXrT044j4X2jDjfpbepxuSXAz8CvD9kXa5NOadbWN3pF+nuZRDkj/q1n+a3tkb1wKHgf+md1Qw9oYc+58CvwDc3h31Hq8xvyrhkOM+5wwz7qo6mOQrwBPAS/T+N7qBp/uNiyH/3n8G3JVkP70pjw9X1dhfcjnJvcBVwIokU8BHgVfAwmWbl2GQpIaM4/SOJGmODH1JaoihL0kNMfQlqSGGviQ1xNCXpIYY+pLUkP8Frde5HwzrVOsAAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"plt.hist(spills, bins=50);"
]
},
{
"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.8.3"
}
},
"nbformat": 4,
"nbformat_minor": 4
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment