diff --git a/Exercise sheet 1/exercise_sheet_01.ipynb b/Exercise sheet 1/exercise_sheet_01.ipynb index 0247783..8290ef4 100644 --- a/Exercise sheet 1/exercise_sheet_01.ipynb +++ b/Exercise sheet 1/exercise_sheet_01.ipynb @@ -199,14 +199,14 @@ { "data": { "text/plain": [ - "array([[1.60000000e+01, 4.26468053e-01],\n", - " [3.20000000e+01, 2.84594009e-01],\n", - " [6.40000000e+01, 1.98253852e-01],\n", - " [1.28000000e+02, 1.49882034e-01],\n", - " [2.56000000e+02, 1.01855227e-01],\n", - " [5.12000000e+02, 6.41574064e-02],\n", - " [1.02400000e+03, 4.84787461e-02],\n", - " [2.04800000e+03, 3.37458855e-02]])" + "array([[1.60000000e+01, 3.96949304e-01],\n", + " [3.20000000e+01, 2.84220425e-01],\n", + " [6.40000000e+01, 2.15515371e-01],\n", + " [1.28000000e+02, 1.59275704e-01],\n", + " [2.56000000e+02, 1.00798933e-01],\n", + " [5.12000000e+02, 6.87799562e-02],\n", + " [1.02400000e+03, 4.87684636e-02],\n", + " [2.04800000e+03, 3.85210905e-02]])" ] }, "execution_count": 3, @@ -368,7 +368,7 @@ "outputs": [ { "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAmsAAAEcCAYAAACYg/MAAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAA89UlEQVR4nO3dd3hUVf7H8fdJoTfpVZCydELviaCAgIACFhAVFEWavaC77uq6Cr9dUVABESxIExUUYUWliSIdBKRJFSEUaRJaQtr5/TGBHcJMMjMkuZPk83qeeZ7MLWe+95xz73xz5hZjrUVEREREglOI0wGIiIiIiHdK1kRERESCmJI1ERERkSCmZE1EREQkiClZExEREQliStZEREREgpiSNREREZEgpmRNREREJIgpWQsixpiaxpiNxpizxpjHnI7nEmPMFGPMq06Ul1GfbYzZZoxp59T6AX6mX/0hdV05EXMwMcbsN8Z0CHDdDO3zEhyupU+IuPPn+JoR/c7nZM0Yc48xZr0x5pwx5ogx5htjTNtr+XC5ynPAMmttYWvt276upANQ+qy1da21y3xZ1lN9+rN+BgqoP1ySWTGrv0l2oH6acxhj8hpjPjDG/J7yz+tGY0yXVMsUN8Z8aYw5n7LcPb7M8/J56fadrP5O8ClZM8Y8BYwFRgJlgOuBCcBtmRaZH4wxYU7HkEEqA9ucDiKj5KB2cUqm9Ae1iwQ7p/uo058vVwkDDgI3AkWBvwOfGWOquC0zHojHlaP0A941xtT1YZ7fHOkf1to0X7gq5hxwZxrL1AaWAadxfbn0cJu3H3gG+AWIAT4F8qXMex6Ynaqst4C3U/4uD8wBjgO/AY+lKndESrkXcTVmY2AjcBb4POWzXnVbJ73yPMaZMr8S8EXKuieBcemV6U9dAUuBJCAupb7/4mHdEcChlO3bCdwMTAOSgdiU9Z5zq9u9KctuB3r6sa2NgJ9T1v0UmJWqHtMrO3W7pFleqm1M77M91jfp96X9QIf0tiGN+ry8vrc29KVuM7g/pFdX7jF7apc0+y4e+ry3+vGlr6bXd9zifDYlzvPAB7gOsN+krLMYuM5t2RdSyvkT+Igr+7H79qe3rT71USAvcCalbc6lvC61UwcPy6d3TLJAdbf3Uy7N9yFmb3XsbXpA5fnRZwM+nnvro6nmDwa+xvXFewI4DHT0EmNa+7G376OrPj+NbX0AmO/2eXuAz9zeHwQa+lmvGblveD1G+VAPafZZD3E/BHwHvItrP9wF1AEeBw6ktFUvb+sH+kqJvXfK3wVxJWN/cZs/Dfi/tOYF0Hfc+0c0PnyneDgW+dQfrorLhwrpDCSSasdxmx+e0lH/CuQBbkoJoqZbkGtx7ajFgR3A4JR5lYELQJGU96HAEaAlrlG/DcA/UsqtCuwDbnErdxOuL5T8Kcv8ntJBwoFeKQ106cDnS3ne4gwFNgNjUho+H9A2vTIDqKtlwENe1q2J6wBQPuV9FaBa6o7gtvydKdsSAtyNa+cu58O2XqrHJ1PivQNI4MovmPTK9tQuXstzKzfNZdOqb9LoS57qyIdtSF2f+4EOPrSh17rN4P7gSztd3g4P7ZLe/uCxz3urHz/6qtd6dyt7Na4voQrAMVxJVCNcidJS4CW3ZbembFNxYIWn7fdhW33uoynLPwAsdHu/C4hKo408HpNSlvGYrPkQs8c6TmN6QOX502e5huO5pz7q4bMnAKdw7e8hwEvA4jT64X4878fejn1XfH4621oVVyIUApRLaedDKeVUxZW0hPharxm5b6QVd3r1gA991kPcb+P6Z+7mlDafk6qMx4GfPaz335Q69PT6r7fPS1m3DK5/kGqlvG8ExKZa5hlgflrzAug77v3jimXSaj/+dyzyuT9cFVO6C7iGDI+mMT8SOAqEuE37BHjZLch73eb9B5jo9v4n4P6UvzsCe1P+bgEcSPVZLwAfuZX7oNu8KFzZqklV9qt+lOcxTqAVrv8GU/+nl2aZAdTVMrx/OVfHtXN2AMLT61ge1t8E3ObDtkbh+o/VvR5XkvbOmrrs1O3iU3npLetDG3rsS77UkYdt8JasXVN/z8D+kG69cnWy5t4u6dWlxz7vY1167atp1btb2f3c3s8B3nV7/ygw123ZwW7zunpqcx+21a8+D7wJvJHyd0FcI2vFvbSR12NSyntvyVp6MXus4zSmB1ReAH02oOO5pz7q4bN/wm0kF9eXYyDJmrdj3xWf78O2HsQ1CtUHmIQr+amFK5mf5+++kFH7Rnpxp1UP+NBnPcT5A/CM2/t/4ZZs4UrifvFn29Opl3Bco4jvpW6rVMs9jOsY6nVeAH3nwbSW8dZ+/O9YFHB/8OWctZNAyTR+oy0PHLTWJrtN+x1X5n/JUbe/LwCF3N7PBPqm/H1Pyntw/ZdW3hhz+tIL138KZdzWPZgqjkM2pWY8zPelPG9xVgJ+t9YmciVfynTnS115ZK3dAzwBvAwcM8bMMsaU97a8MeZ+Y8wmt7jqASXdFvG2rZ7q8Xc/y06vXa4oz49l06tvb33pKj5sgzfX2t/9LSutOHyt10v82R+89fmrGGP6pVx4dM4Y801afdXHev/D7e9YD+/d69N9m37HVS+ppbet/tZlfWCL299HrLWnPCyX3jEpLWnG7K2O06j7gMrzsk1p9dlrOZ6nVz/1cY2UXFIP189N/kpr/0x97EprW38A2uFKcH7AlRjcmPL6Afw7bmfgvuHrccVTPQTSZxvgGiW7pI6H97+mU4ZPjDEhuH6mjAeGu806BxRJtXgRXCOKac3zl9e68KX9/P0ed+dLsrYK13Dj7V7mHwYqpVTiJdfjys598TnQzhhTEejJ/3bug8Bv1tpibq/C1tqubuu6d6gjQAVjjHGbVsntb1/K8+YgcL2HhNXfMq+prqy1M621bXEd+Czw70uz3JczxlQGJuPqzCWstcVw/VzkXjfeeKrH6/0sO712uR7P0ls2vfr21peu4MM2WE/rpbjW/p5RZflTr5ekPgCnVZfe+nzqcrDWzrDWFkp5dUmZdlVfvcZ+6Y37Pn49rjpNLb1t9bcu3ZO1CLe/U0vvmASuL8kCbu/L+hiz1+OBl+kBl5dKen32Wo7n4GXfM8bcgOscsp1ukxvhGr3wJq392Jd10tvWS8laZMrfP5AqWQPf6jWD942MPq6k7rOXpcSdB9epAJc05Mp2aYCHdjKuO0qc8/L6xsPyhv+dp9fbWpvgNnsXEGaMqeE2LQLX+XppzfPGW9/x1j99bj8f97OrpJusWWtjcJ1nMN4Yc7sxpoAxJtwY08UY8x9gDa7fZp9Lmd4O6I7rBN10WWuP4/qP5CNcO/OOlFlrgTPGmBHGmPzGmFBjTD1jTDMvRa3C9XPEcGNMmDHmNqC523x/y3O3Flcn/j9jTEFjTD5jTJsAygy4rozrnls3GWPy4kqeY1O2F1z/XVV1W7wgrk5wPGXdB3Bl+b5YhescxcdS6rEXV9ajv2WnV54/y6ZZ32n0pdTS24bU9enumvp7BpblT716kl7f9dbnIe36SauvXku/9GaYMaaiMaY4rpGaTwPYVp/r0hhTEijF/0Z0anHlF5W79I5J4PoSuyclps64vujTjdlbHadR9wGV52Gb0uyzGXg8T60BsCXVaFEjXOdVepNmP/VBevvnD0B7XOfXRQPLcZ3jXQLXCfr+1GtG7hvXelxJr8+6i8CtXYwxRXAlIb+kWuaqdrLWdnH7Jy/1q0vq5XFdwFAb6G6tjU1V1nlcF0O9knK8aoPrbhXT0pqXxnb523d8aj8/+sNVfLp1h7X2TeAp4MWUYA7iyiDnWmvjgR5AF1xXfUzAdc6CP8OeM3H9hnt5JMRam4SrgzXEdeXQCeB9XFeneooxHtfJkANxnaB4L66h2IuBlJeq7EvrVsd1dUs0cHeAMQZaV3lxXdlyAtfwdWlcX04Ao4AXjWv49Rlr7XbgDVw73h+4RgJW+PAZ7vU4ANdJsnfj6uiX5vtVdnrl+fnZvtT3VX3Jw+ektw1X1KeHGK+1v19zWf7Uq5f106xLb30+ZXWv9ZPCY1+9ln6ZhpnAQlwnq+/Ddb6Xv9vqT102wHUe1qUvi0PAXcaYFh4+N81jUorHU2I7jev84Lm+xIz344G3ug+0PE/blF6fvebjuQdXjM4YY0rgGoXcmsY66fXTNKW3rdbaXbh+Ylue8v4Mrj64ImV7wfd6zbB9I4OOK2n1WXcRXD2KtsdaewEu/2xZj7RHQNNlXCNXj+DqP0fN/0bg+rktNhTXif/HcJ2jN8Rau82HeZ741Xf8aD+f+oMnxl7x03TOYoxZg+vk0Y+cjkVEMpYxZj+uCzAWOx2Lr3RMkuxGfTY45KjHTRljbjTGlE0Zvu2PK8v/1um4RCR30jFJshv12eCU0+7SXBP4DNdVLXuBO6y1R5wNSURyMR2TJLtRnw1COfpnUBEREZHsLkf9DCoiIiKS0yhZExEREQliOe2cNckiJUuWtFWqVHE6DBGRbGXDhg0nrLWlnI5DshclaxKQKlWqsH79eqfDEBHJVowx6T0WTuQq+hlUREREJIgpWRO/GGO6G2MmxcTEOB2KiIhIrqBkTfxirZ1vrR1UtKivT4kRERGRa6Fz1kQk10hISCA6Opq4uDinQ5EcLl++fFSsWJHw8HCnQ5EcQMmaiOQa0dHRFC5cmCpVqmCMcTocyaGstZw8eZLo6GhuuOEGp8ORHEA/g4pIrhEXF0eJEiWUqEmmMsZQokQJjeBKhlGyJlkr/gL88jnoMWfiECVqkhXUzyQjKVmTrLVxOnzxEEy7Hf7c73Q0Io57+eWXGT16tNf5c+fOZfv27VkYkYgEGyVrkrWaPQS3vgHR62FCK1g9EZKTnI5KJGgpWRMRJWuStUJCXAnb0NVQuQ18OwI+7AzHdzodmUiWee2116hZsyYdOnRg505X3588eTLNmjUjIiKC3r17c+HCBVauXMm8efN49tlnadiwIXv37vW4nIjkbMbq3CEJQNOmTe01P27KWvjlM1fCFn8ebnwO2jwBobrUXTLHjh07qF27NgD/nL+N7YfPZGj5dcoX4aXuddNcZsOGDQwYMIA1a9aQmJhI48aNGTx4MA888AAlSpQA4MUXX6RMmTI8+uijDBgwgG7dunHHHXcAcPLkSY/LSfBx72+XGGM2WGubOhSSZFMaWRO/ZOgTDIyBiLth2Dqo1Q2WvgqT2sHhjddetkiQWr58OT179qRAgQIUKVKEHj16ALB161YiIyOpX78+M2bMYNu2bR7X93U5Eck5dJ818Yu1dj4wv2nTpg9nWKGFSsGdH0H9O+C/T8Hkm6H1cGj3AoTnz7CPEXGX3ghYZvJ0peCAAQOYO3cuERERTJkyhWXLlnlc19flRCTn0MiaBI9at8KwNdCoH6x4C95tA/tXOB2VSIaKioriyy+/JDY2lrNnzzJ//nwAzp49S7ly5UhISGDGjBmXly9cuDBnz569/N7bciKScylZk+CSvxj0eAfu/wqSE2FKV9doW1zGnlsk4pTGjRtz991307BhQ3r37k1kZCQA//rXv2jRogUdO3akVq1al5fv06cPr7/+Oo0aNWLv3r1elxORnEsXGEhAMuQCg/TEn4elr8HqCVCkAnQbA3/plLmfKTmapxO+RTKLLjCQjKKRNQleeQpC55EwcBHkLQQz74QvBsH5k05HJiIikmWUrEnwq9QMHvkRbhwBW+fA+Oaw9Qs9skpERHIFJWuSPYTlhfZ/dSVtxSrB7AdgVj84cwT9lC8iIjmZkjXJXsrUhYGLodOrsHcJjG/Bj7Pe4NnPNnH+YqLT0YmIiGQ4JWuS/YSGQetHYchKKFufG3f+i9u3DmXQW7PZeigDbtYrIiISRJSsSfZVohr0nw/dxtIi7+98cOFR5k38Gx8u36OfRkVEJMdQsibZW0gINH2AsOFrCal6I38NnUajRXfz4qTPOXHuotPRiYiIXDMla5IzFK1Anvs+x/Z6n9p5T/LS4cHMefNRVvx62OnIREREromSNfFLhj7IPaMZg2lwJ/meWE9sje48kvwpJWZ2YsrnX5CQlOx0dCKXvf3229SuXZt+/frRunVrAE6fPs2ECRMuL5P6va9efvllRo8enWGxZoRjx47RuHFjXnjhBXr16kVysvZHEX8oWRO/WGvnW2sHFS1a1OlQvCtYkqL3fszFO2dSLm8c9219kP+OHsiBo8edjkwEgAkTJrBgwQJmzJjBypUrgYxL1oLRunXr6Nu3L6NGjaJ06dKcPKkbW4v4Q8ma5Fh5695K0ac3cKjqXfSM/QLebcNPi750OizJ5QYPHsy+ffvo0aMHY8aMoVChQgA8//zz7N27l4YNG/Lss89e9R5g+vTpNG/enIYNG/LII4+QlJQEwGuvvUbNmjXp0KEDO3fu9Pi5mzdvJioqijp16hASEoIxhpdeeimgbZg9ezYtW7YkIiKCtm3bcvy46x+hnj178uKLLxIZGUnZsmVZvHgx4ErWIiIiAIiJiaFUqVIBfa5IrmWt1Usvv19NmjSx2cmxzQvt4VdqWvtSEbtyTD979vRJp0MSB2zfvt3pEKy11lauXNkeP37cWmttwYIFrbXW/vbbb7Zu3bqXl0n9fvv27bZbt242Pj7eWmvtkCFD7Mcff2zXr19v69WrZ8+fP29jYmJstWrV7Ouvv37F58XGxtqaNWvaNWvWWGutffHFF+0zzzxjk5OTr1iubdu2NiIi4qrXokWLrljuxIkTl/9++eWX7bhx46y11lavXv3yZ8+ZM8cOGDDAWmttnz597PDhw+3gwYPt4sWLA6y17MdTfwPW2yA4huuVvV5hTieLIlmhVIOOJP5lPeunjqD5oRmcGruSEx1fp0rr3k6HJuKTJUuWsGHDBpo1awZAbGwspUuX5tSpU/Ts2ZMCBQoA0KNHj6vWXbx4MY0bN6Z58+YANGjQgG+//RZjzBXLLV++3KdYpkyZwqeffsrFixc5evQoI0eO5MKFC8TExPDkk08CkJiYSLFixQAICwvjnXfeCWi7RQQla5J7hOUrRNNB49mytjf5v3mC6gsfZM/mz6h67zhCCutnGQlu1lr69+/PqFGjrpg+duzYq5Ku1LZu3Ur9+vUvv//5559p3LjxVctFRkZy9uzZq6aPHj2aDh06ADB16lTWrl3L0qVLKVSoEFFRUdStW5dt27bRpEkTQkNDAfjll1+oV68eANOmTfNvY0XkCjpnTXKd+s1vouRTK5l3XX+uP7qIc2OacGbtJ3owvDiqcOHCVyRKqd/ffPPNzJ49m2PHjgFw6tQpfv/9d6Kiovjyyy+JjY3l7NmzzJ8//6qyS5QowS+//ALArl27+OKLL+jTp89Vyy1fvpxNmzZd9bqUqAFs2bKF1q1bU6hQIebMmcPKlSupX78+W7dupWHDhpeX++WXX2jQoME114uIKFmTXKpY4UJ0f+wtvmv7OfuSSlFkwWBOTu4JMYecDk1yqRIlStCmTRvq1avHs88+e9X7OnXq8Oqrr9KpUycaNGhAx44dOXLkCI0bN+buu++mYcOG9O7dm8jIyKvK7tu3L+fOnaNevXoMGjSITz75hBIlSgQUZ//+/Xn77beJjIxk165dVK1alYIFC7Jly5YrkrWtW7deHlkTkWtjrEYTJABNmza169evdzqMDLHz8Gm+n/oK/WOnERIaRugtrxLW7AHX0xEkR9mxYwe1a9d2OgzJJTz1N2PMBmttU4dCkmxK30aS69UsX4z+T73O+NrTWJ9wA2HfPEXsB7fCyb1OhyYiIqJkTQQgf55QnunTmbN3zeFlBpMYvZGk8a1gxduQlOh0eCIikovpalARN53rl6NBpZd5ZkZb7jg6ho6L/k7S1i8IvX08lKnrdHiSkb55Ho5uydgyy9aHLv+XsWWKSK6nkTWRVMoXy8+EId3YHjWR4QmPcebIPuzEKPh+JCRedDo8ERHJZTSyJuJBaIjh8Y5/YV2NJ+k7swmDYydz+w//xm7/CtNjHFRq5nSIcq00AiYi2YRG1kTS0KxKcT594la+q/kKA+Kf5dSpE9gPOsK3f4X4806HJyIiuYCSNfGLMaa7MWZSTEyM06FkmaIFwpnQrzGdbrufTvH/YbbpBKvHw4RWsG+Z0+FJNhQbG8uNN954+UHsGeXBBx+kdOnSXu9vtnPnTho2bHj5VaRIEcaOHQvAmDFjqFu3LvXq1aNv377ExcURFxdH8+bNiYiIoG7dugE/+D2rfPvtt9SsWZPq1avzf//neeS0SpUq1K9fn4YNG9K06f/uoHHw4EHat29P7dq1qVu3Lm+99Vaa0+Pj44mKiiIxURcgSRZw+uGkemXPV3Z7kHtG2XX0jL1lzA/2zudftydG1rX2pSLWfjXc2gt/Oh2a+CBYHuQ+btw4O3bs2Awv94cffrAbNmy44gHw3iQmJtoyZcrY/fv32+joaFulShV74cIFa621d955p/3oo49scnKyPXv2rLXW2vj4eNu8eXO7atWqDI87IyQmJtqqVavavXv32osXL9oGDRrYbdu2XbVc5cqV7fHjx6+afvjwYbthwwZrrbVnzpyxNWrUsNu2bfM63VrXQ+ynT5/uNSY9yF2vjHppZE3EDzXKFGbusDbUatmZ1jH/Yna+3tiN02F8C/j1a6fDk2xixowZ3HbbbQDMnj2bli1bEhERQdu2bTl+/HjA5UZFRVG8eHGfll2yZAnVqlWjcuXKgOvB67GxsSQmJnLhwgXKly+PMYZChQoBkJCQQEJCQrrPIe3Tpw933303LVq0oHLlynz9ddbsF2vXrqV69epUrVqVPHny0KdPH7766iuf1y9Xrtzl56UWLlyY2rVrc+jQIa/TAW6//XZmzJiR8RsjkoqSNRE/5QsP5ZXb6vHOfa14Lb4PvRP+xQlbGGbdA58PgHPHnA5Rglh8fDz79u2jSpUqALRv357Vq1ezefNmOnbsyGeffXbF8pGRkVf8dHnptXjx4muKY9asWfTt2xeAChUq8Mwzz3D99ddTrlw5ihYtSqdOnQBISkqiYcOGlC5dmo4dO9KiRYs0y928eTNVq1ZlzZo1zJgxg3/+859Zsj2HDh2iUqVKl99XrFjxclLlzhhDp06daNKkCZMmTfJY1v79+9m4ceNV25p6er169Vi3bt01xS3iC10NKhKgTnXL0rjydbw8ryQtf6nE34st5L4dnxGybxl0/jc0uAvSGYWQ3OfEiRMUK1bs8vspU6bw6aefcvHiRY4ePcrIkSOvWH758uUZHkN8fDzz5s1j1KhRAPz555989dVX/PbbbxQrVow777yT6dOnc++99xIaGsqmTZs4ffo0PXv2TPOZn7GxsZw4ceLyuW116tThzz//vKbt6dChA0ePHr1q+muvvXZ5dBJcp/Sk5mkUcMWKFZQvX55jx47RsWNHatWqRVRU1OX5586do3fv3owdO5YiRYqkOT00NJQ8efJw9uxZChcu7Nd2ifhDyZrINShZKC/j7mnMwojyvDi3IDPjIvigwFQqfjkItnwO3cZAsUrpFyS5Rv78+YmLiwNg6tSprF27lqVLl1KoUCGioqKoW/fKmy9HRkZy9uzZq8oZPXo0HTp0CCiGb775hsaNG1OmTBkAFi9ezA033ECpUqUA6NWrFytXruTee++9vE6xYsVo164d3377rddkbevWrdSoUYN8+fIB8PPPPxMREXFN2+PriFvFihU5ePDg5ffR0dGUL1/+quUuTStdujQ9e/Zk7dq1l5O1hIQEevfuTb9+/ejVq9fldbxNB7h48eLl7RXJLErWRDJAp7plaVG1BKMWlCZqXTmeLLKMoftnEDqhJXT8JzR5UA+GFwCuu+46kpKSiIuLY8uWLbRu3ZpChQoxZ84cVq5cSf369a9YPjNG1j755JPLP4ECXH/99axevZoLFy6QP39+lixZQtOmTTl+/Djh4eEUK1aM2NhYFi9ezIgRIy6vd/PNNzN16lQqVKgAuH4CPXDgAHFxcSQlJfHSSy/xn//8J9O3B6BZs2bs3r2b3377jQoVKjBr1ixmzpx5xTLnz58nOTmZwoULc/78eRYuXMg//vEPwDUyN3DgQGrXrs1TTz11eR1v0wFOnjxJqVKlCA8Pz5RtErlE3x4iGaRo/nD+r3cDpj3Uis/DunHj+VHsyVMLvn4aptwKJ/Y4HaIEiU6dOvHTTz/Rv39/3n77bSIjI9m1axdVq1alYMGCAZfbt29fWrVqxc6dO6lYsSIffPABAF27duXw4cMAXLhwgUWLFl0xQtSiRQvuuOMOGjduTP369UlOTmbQoEEcOXKE9u3b06BBA5o1a0bHjh3p1q0bAMnJyezZs+eKCxo2b95Mv379aNeuHc2aNWPIkCG0adMm4O3xR1hYGOPGjeOWW26hdu3a3HXXXZdHKS9t/x9//EHbtm2JiIigefPm3HrrrXTu3Blw/Tw6bdo0li5devk8ugULFnidDvD999/TtWvXLNk+yd2Mp9/5RdLTtGlTu379eqfDCFoX4hN5c+EuPlyxjwcLrmSEmUa4jYd2L0Cr4RCqQW0n7Nixg9q1azsdBhs3buTNN99k2rRpTocSsK1bt/Lhhx/y5ptvXp4WFRXF5MmTqVmzpoORZZ1evXoxatQor9vrqb8ZYzZYa5t6XEHEC42siWSCAnnCeLFbHeYMacMPBW6h9dlRbM7XDBa/BO/flPEPEJdspVGjRrRv3z7Db4qblerVq3dFogawd+9eatSo4VBEWSs+Pp7bb7891ySm4iyNrElANLLmu4uJSUz4fi8Tlu2mZ94N/Ct8CnkSYjBtnoCoZyFcJydnlWAZWZPcQSNrklE0siaSyfKGhfJkx78w/9FIdha/iRYxI1lVoD0sHw3vRcKBNU6HKCIiQUzJmkgWqVW2CF8MbcOwrs15MGYgg+1fOX/uLPbDW+CbEXDxnNMhiohIEFKyJpKFQkMMD0dV5dvHozhdIYrmp19lYcHusGYivNsK9i51OsQcT6d+SFZQP5OMpGRNxAFVShZk5kMtebFXc545dy/3JL3M6fgQmNYT5g6D2D/TL0T8li9fPk6ePKkvUslU1lpOnjypm+VKhtEFBhIQXWCQcY7GxPHi3C0s3xHNq9ct4I64OZiCJaHraKjTw+nwcpSEhASio6MvP0FAJLPky5ePihUrXnXDXF1gIIFQsiYBUbKWsay1zP/lCC/P20aluF1MLjaF0ud3Qe0erqStcBmnQxSRDKBkTQKhn0FFgoAxhh4R5Vn81I3cUL81rU++yId57yd513cwvjlsmgn6x0pEJFdSsiYSRIoXzMPYPo2YNKAlk7mdTrGvcjDsepg7BKb3htMHnA5RRESymJI1kSB0U60yLHwyihbNWhF14jnGhj9M0u+rYHxLWDMJkpOdDlFERLKIkjWRIFU4Xziv9azPJ4NaMzfPrUSdG8WuvHXhm2fhoy5wYrfTIYqISBZQsiYS5FpWLcG3T0TR7cYWdD75OP8MHU7CHzvg3Taw/A1ISnA6RBERyURK1kSygXzhobzQpTZzh7VlVeFbaH1mFBvytYAlr8Dk9nBks9MhiohIJlGyJn4xxnQ3xkyKiYlxOpRcqUHFYswb3pb7OzbnvrPDGJzwBGeOR2MntYfFL0OC7h8mIpLT6D5rEhDdZ815J89dZNLyfcxduY1nmMqdoT8QX6wqeXpOgMqtnA5PRDzQfdYkEBpZE8mmShTKywtdavP1iO7sbvVvHkz6K3/8eRY+6kzMnMfh4lmnQxQRkQygkTUJiEbWgs/xsxf56PstlF3/Ovea7zidpzTxncdQtsmtTocmIik0siaB0MiaSA5RqnBenuvRlC7PTmVq7Yn8GR9K2fn3sG7MXRyIjnY6PBERCZCSNZEcplThvAzo04fCT6xiebkBNDy9mAKTWzH1/bf4/eR5p8MTERE/KVkTyaFKX1eMyEfe4ux9C4kvWI77o//Bjrdu55WZSzhw8oLT4YmIiI90zpoEROesZTNJiZxbNpa8P/2bC8nhjEy6FyL6MfzmGlQqXsDp6ERyDZ2zJoHQyJpIbhAaRqGbnyF82CryV6zPv8Peo8eWYdw3+jOen/MLB09ppE1EJFgpWRPJTUpWJ8/Ab+DWN2iVdx8L842g4Kb3uXn0Ul744hei/1TSJiISbJSsieQ2ISHQ7CFChq0hT9VI/h76MYuvG8WmDWtpP3oZL3yxhUOnY52OUkREUihZE8mtilWCfp9Dz0lcn3yIBfle4N1K3/PVhv20e/17/valkjYRkWCgCwwkILrAIIc5dxy+eQ62fUFCyTpMLv4UY7YVBODuZpUY2q465YvldzhIkexPFxhIIDSyJiJQqBTc+RH0mUl43CmG7n6En1v+xD2NS/PpuoO0e30Zf5+7lSMxGmkTEclqGlmTgGhkLQeLPQ2L/g4/T4Xi1Th20+uM3V2az9cfxGDo07wSj91cg5KF8jodqUi2o5E1CYRG1kTkSvmLQY934P6vIDmR0rN7MTL8I5Y92oTeTSoyc80BbhnzI4u2/+F0pCIiuYKSNRHxrGo7GLoKWg6D9R9SYeZNjKp3hAWPR1KmSD4enrqeEbN/4dzFRKcjFRHJ0ZSsiYh3eQpC55EwcBHkLQQz7+QvK55m7gO1GNquGp9vOEiXt35k3f5TTkcqIpJjKVkTkfRVagaP/Ag3Pg9b55BnYkueq7iNzwa1xGC4671V/PvbX4lPTHY6UhGRHEfJmoj4JiwvtH/BlbQVqwSzH6Tp6uF882B1+jSrxLvL9nLb+BXsPHrW6UhFRHIUJWsi4p8ydWHgYuj0KuxdSsHJbRhVeSPv39eE42fj6P7OT7y/fB/JybrSXEQkIyhZExH/hYZB60dhyEooWx/mP0aH9Q+zaEAl2tUsxatf7+Ce91frWaMiIhlAyZqIBK5ENeg/H7qNhcObuG5KO96rvorXe9dl66EzdBm7nDkbotH9HEVEAqdkTUSuTUgINH0Ahq6GqjdiFr7InZseZPG9pahdrghPf76ZIdN/5tT5eKcjFRHJlpSsiUjGKFoB+s6C3h/An/sp+0knZtVcxt9uqcbSX4/RacyPfP/rMaejFBHJdpSsiUjGMQbq3wHD1kLd2wn54f94eMcDfHdXAUoWysMDU9bx1y+3cF430hUR8ZmSNRHJeAVLQu/3XSNtsae5Ye7t/LfWtwxrU55P1h7g1reX8/OBP52OUkQkW1CyJiKZp2YXGLYaGvcnbPV4nt03gAXdLQlJljveXckbC3eSkKQb6YqIpEXJmohkrnxFoftY6P9fwFB7YT++r/Ul90QU452le+g5YQV7julGuiIi3ihZE5GscUOk675srR8lz+bpvHpoIHNuiuHw6ThuffsnPlrxm26kKyLigZI1Eck6eQq4nnzw0GLIX5wmK4ewsvp0Ot8Qxj/nb+e+D9dwJCbW6ShFRIKKkjURyXoVmsCgZdDur+Tb/TVjTwxiZosDbDzwJ7eM+ZGvNh3SjXRFRFIoWRMRZ4TlgXYjYPByTPGqtN78POurvk+LEnE8PmsTwz/ZyOkLupGuiIiSNRFxVunaMHAh3DKKAodWMunsMKbU38rCrYe5ZeyP/LDruNMRiog4SsmaAGCMqWqM+cAYM9vpWCQXCgmFVkNhyEpMhca02z2STZXfoXb4cfp/uJZ/fLWV2Pgkp6MUEXGEkrUcwBjzoTHmmDFma6rpnY0xO40xe4wxz6dVhrV2n7V2YOZGKpKO4jfA/V9Bj3coeGoHH118gsnVVjBj1T5ufXs5mw+edjpCEZEsp2QtZ5gCdHafYIwJBcYDXYA6QF9jTB1jTH1jzH9TvUpnfcgiXhgDje+HYWsw1W6m46HxbCr/H8pf3Evvd1cya+0BpyMUEclSYU4HINfOWvujMaZKqsnNgT3W2n0AxphZwG3W2lFAtywOUcR/RcpBnxmw7UsKL3iWaUnPMbfk3Yz4Ip5df5zjr11rERaq/zdFJOfTkS7nqgAcdHsfnTLNI2NMCWPMRKCRMeYFL8sMMsasN8asP35cJ31LFjAG6vWC4esw9e6g55kZ/FTsZTauXMjAj9dzJi7B6QhFRDKdkrWcy3iY5vXGVdbak9bawdbaaimjb56WmWStbWqtbVqqVKkMC1QkXQWKQ6/3oN9sSudJ4Iu8L9PutzfpO24J+0+cdzo6EZFMpWQt54oGKrm9rwgcdigWkYxRoyMMXYVpNpAHQr/hvXOP8tq4iazcc8LpyEREMo2StZxrHVDDGHODMSYP0AeY53BMItcuXxG49Q0YsIAyxQoxmVeI/vghPl2+xenIREQyhZK1IGGMeeka1v0EWAXUNMZEG2MGWmsTgeHAd8AO4DNr7baMiVYkCFRpQ/iwlVxs+Ri9Q3+g3eJuzPx4AolJyU5HJiKSoYyevxccjDHJwOtAceBnYJa19k9no7qaMaY70L169eoP79692+lwRABIit7IiRkPUyZ2N2vyR1H7wYkUKeX1ehoRxxhjNlhrmzodh2QvGlkLHhaIwzUSVglYaYyJcDakq1lr51trBxUtWtTpUEQuC63YiDLPrGJLzcdodGEldnwL/vhpCuifURHJATSyFiSMMdustXXd3v8FmGitvcnBsLxq2rSpXb9+vdNhiFxly6a1JM8dRgS7OFXuRorfPR6KVUp/RZEsoJE1CYRG1oLHCWNMk0tvrLW7AN0fQ8RP9Rs2p/jwpUzIP4i8h1eT8E5z7NrJkKxz2UQke1KyFjweA6YbY6YbY0YYY2YAvzkdlEh2VKlkYe5/YhSvVPqQVfHVMAueIXlKVzixx+nQRET8pmQtSFhrNwMNgU9SJn0P9HUsIJFsrlDeMEY9eCurWk/mmYRHuHBwC/bd1vDTGEhKdDo8ERGf6Zw1CYjOWZPs5MuN0Yye8yMj837MjUmroVwE9BgH5Ro4HZrkMjpnTQKhkTXxizGmuzFmUkxMjNOhiPisZ6OKvDOoC0+bZ3nCPsXFPw/BpHaw5BVIiHM6PBGRNClZE7/o1h2SXTW+/jrmDW/DruI30TJmJLvKdoXlb8B7kXBgjdPhiYh4pWRNRHKN8sXyM3tIK1rUqU6n3/ryfuXR2PgL8OEtsOA5uHjO6RBFRK6iZE1EcpUCecKY0K8xj91UnVd3luf+fG8R2+hBWDsJJrSCPUucDlFE5ApK1kQk1wkJMTzVqSbv9G3E2sMJdNjRjd9vmw1heWF6L5g7FC6ccjpMERFAyZqI5GLdI8rz2SOtSEhKpuuXiSxp9wVEPg2bZ8GElrB9ntMhiogoWROR3C2iUjHmDW9LtdKFeGjmFt4N7Ycd9D0UKgOf3Qef3gdn/3A6TBHJxZSsiV906w7JicoWzceng1pxa/1y/PvbX3n6R0vcgEVw80uw6zsY3xw2zdSD4UXEEboprgREN8WVnMhayztL9/Dmol00ur4Y793XhNIXD8K8R+HAKqh2E3QbC9dVdjpUyaZ0U1wJhEbWRERSGGN47OYavNuvMb8eOcvt41awLb40DFgAXUfDwbWuK0bXvKcHw4tIllGyJiKSSpf65fh8cCsscMe7q/h2+x/Q/GEYugoqt4JvnoOPusDxXU6HKiK5gJI1EREP6lUoylfD21CzbGEGT/+ZF+du4agpDf1mQ8/34MROmNgGfhwNSQlOhysiOZjOWZOA6Jw1yS3iEpIYuWAHM9ccICTEcE/z6xnSrhplQs7Agmdh+1woWx9uG+96QLxIGnTOmgRCyZoERMma5DYHT11g/Pd7mL0h+nLSNrRdNUofWgRfPw3nT0Cbx+DG5yE8n9PhSpBSsiaBULImAVGyJrnVgZMXGPf9bub8fIiwEEO/FpUZ0rI4pVb+CzZOhxLVocc417ltIqkoWZNAKFmTgChZk9zu95PneWfpHr7ceIjwUMO9LSozvEo0xRY/DacPQLOHocNLkLew06FKEFGyJoFQsiZ+McZ0B7pXr1794d27dzsdjojj9p+4lLRFkycshAebleJR8xn5N0yCohVd92Wr0cHpMCVIKFmTQChZk4BoZE3kSr+dOM87S3Yzd9Mh8oaF8kK9M9zzx+uEndoFEX3hlpFQoLjTYYrDlKxJIJSsSUCUrIl4tvf4Od5Zspt5mw9TODyZCZWW0vrIVEz+61w31q1zGxjjdJjiECVrEgjdZ01EJANVK1WIsX0asfDJG2lXpyL37utAr8TXOEIJ+Lw/fHovnD3qdJgiko0oWRMRyQTVSxfirT6NWPRkFBVrtyDy1N8YndyPhJ2LSB7X3HXlqH7ZEBEfKFkTEclE1UsX5p2+jVjwRHt+q/UQnS6OZENcefhqGAlTboM/9zsdoogEOZ2zJgHROWsigdl59CzvLN5J0R0zeCHsE/KGWhLb/Z38bYdASKjT4Ukm0zlrEgglaxIQJWsi1+bXo2eY+s0KOuz7P24K3cThwvUpfNdECleq53RokomUrEkg9DOoiIgDapUtwsgHulB28Dw+KP0C+c78Rt4PbmTlhyM4c/6C0+GJSBBRsiYi4qA6FYoycOjzHLv/RzYVjKT1gYkceb0ln8ydx9m4BKfDE5EgoJ9BxS96goFI5vp95ecUWTKCIomnmGp6cLHtczzQrjZ5w3Q+W06gn0ElEErWJCA6Z00kE8XFcOrLERTf+Qn7ksvyTqHHuPfue2hS+TqnI5NrpGRNAqGfQUVEgk2+ohTvOxHun0e5IuGMufBXtr//EKPmruVCfKLT0YlIFlOyJiISrKreSP7H1hLffCj9QpfSf+PdvPTGGFbsOeF0ZCKShZSsiYgEszwFyNN1FCEPLaLYdSV5/eKrHPv4fv4560diYnUBgkhuoGRNRCQ7qNiUAsNXkBj5HD3C1jB8Rz/+PXokC7cecToyEclkStZERLKLsDyE3fw3Qgf/SP7SNzAy6U3Mp/3428ffceLcRaejE5FMomRNRCS7KVOXAkO+J6njq7QL38rz+/oz4Y2/M/fnaHSFv0jOo2RNRCQ7CgkltM2jhA9fTWiFRvzDvkfpL+/k+ffncfh0rNPRiUgGUrImIpKdFa9KgYcXkNztLZrmOcDL0Q8xfcxzzFy9j+RkjbKJ5ARK1kREsjtjCGk6gDyPr8NWbcdzZip1FtzBcxNmsf/EeaejE5FrpGRNRCSnKFKeAvd/hu39AbXynWLUiWHMf/tR3l/2K4lJyU5HJyIBUrImfjHGdDfGTIqJiXE6FBHxxBhM/TvI9/gGEmvdzqMhc4hc2pvn3/6IX4+ecTo6EQmAng0qAdGzQUWyB7vzW+LmPk6e2GNMSerMhTYv8EiHeuQJ0//qTtCzQSUQ2ltFRHIwU7Mz+R9fR0LE/QwMXUCPlb3525jxbDp42unQRMRHStZERHK6fEXI1/MtGPA1pYrk5/XzL7Jj0gBGf7WW2Pgkp6MTkXQoWRMRyS2qtKXAY6uJbzGcu0N/4N6f7+KVN0azau9JpyMTkTQoWRMRyU3yFCBPl9cIeXgJha8rzaiLIzkx5R5e/exHzsTpwfAiwUjJmohIblShMQWH/0RC1F/pEraBodv68sbr/2LJ9qNORyYiqShZExHJrcLyEH7TCMKG/ETe0jX4Z9Jb8EkfXpr2HSf1YHiRoKFkTUQktytdi4JDlpDYaRRR4Tt4dk9/Jr7xN+ZvinY6MhFByZqIiACEhBLWeijhw1djKjbhb3YyDZbcCyf3Oh2ZSK6nZE1ERP6n+A0UfOi/JHd/h+vj98K7rWHFW5CU6HRkIrmWkjUREbmSMYQ0uR8zbC1U7wCL/gHv3wxHtzgdmUiupGRNREQ8K1IO7p4Od06BM4dgUjtY+iok6uIDkaykZE1ERLwzBur2hGFrof6d8OPrMDESDq51OjKRXEPJmoiIpK9Aceg5EfrNhvjz8EEn+OZ5198ikqmUrImIiO9qdIRhq6HZQ7DmXZjQEvZ+73RUIjmakjXxizGmuzFmUkxMjNOhiIhT8haGW0fDA99ASDhMux2+Ggaxp52OTCRHUrImfrHWzrfWDipatKjToYiI0yq3hiEroO2TsOkTGN8CdvzX6ahEchwlayIiErjw/NDhZXh4KRQqBZ/2g8/6w7ljTkcmkmMoWRMRkWtXviE8/D3c9HfYuQDGNXONtlnrdGQi2Z6SNRERyRih4RD1DAxeAaVqwtzBMOMOOH3Q6chEsjUlayIikrFK/QUe+Ba6/Ad+X+W6YnTtZEhOdjoykWxJyZqIiGS8kBBo8QgMXQUVm8GCZ2BKVzix2+nIRLIdJWsiIpJ5rqsM930Jt02AY9vh3Taw/E1ISnA6MpFsQ8maiIhkLmOgUT8Ytg7+0gmW/BMm3wRHNjsdmUi2oGRNRESyRuEyrgfD3zUVzh6FSe1hySuQEOd0ZCJBTcmaiIhkrTq3wbA1ENEHlr8BE9vCgdVORyUStJSsiYhI1itQHG6fAPd+AYkX4cPOsOA5uHjO6chEgo6SNRERcU71m11XjDYfBGsnwYRWsGeJ01GJBBUlayIi4qy8haDrf+DBbyEsL0zvBXOHwoVTTkcmEhSUrImISHC4viUM/gkin4bNs1wPht/+ldNRiThOyZqIiASP8Hxw8z9g0DIoXBY+ux8+vdd19ahILqVkTUREgk+5BvDwUujwMuxaCOObw8YZejC85EpK1kREJDiFhkPbJ2HICihdB74aCtN6wp+/Ox2ZSJZSsiYiIsGtZA0YsAC6joboda4rRte8B8lJTkcmkiWUrImISPALCYHmD7tu81G5FXzzHHzUBY7vdDoykUynZE1ERLKPYtdDv9nQ8z04scv19IMfX9eD4SVHU7ImIiLZizGuR1UNWws1u8LSV13PGT28yenIRDKFkjXxizGmuzFmUkxMjNOhiEhuV6g03PWx6+Hw54/B5Jtg0UuQEOt0ZCIZSsma+MVaO99aO6ho0aJOhyIi4lK7u+vB8A3vgRVjXT+N/r7S6ahEMoySNRERyf7yXwe3jYP75kJSvOvig6+fhrgzTkcmcs2UrImISM5RrT0MXQ0th8K6D1y3+di9yOmoRK6JkjUREclZ8hSEzqNg4ELX3zPugC8e0YPhJdtSsiYiIjlTpeYweDlEPQdbZ8O4ZrD1Cz2ySrIdJWsiIpJzheWFm/4Gg36AohVh9gOuB8OfOeJ0ZCI+U7ImIiI5X9l68NAS6PgK7FkM41vAz1M1yibZgpI1ERHJHULDoM3jMGSlK3mb9yhMvQ1O/eZ0ZCJpUrImIiK5S4lq0P+/cOubcOhneLc1rJqgB8NL0FKyJiIiuU9ICDQbCMNWQ5VI+O4F+KATHNvhdGQiV1GyJiIiuVfRinDPp9DrfTi1DyZGwg//gcR4pyMTuUzJmoiI5G7GQIM7Yfg6qNMDvn8NJrWDQxucjkwEULImIiLiUrAk3PEh9PkEYk/B+x1g4YsQf8HpyCSXU7ImIiLirlZX14PhG98PK9+BiW3gt+VORyW5mJI1ERGR1PIVhe5vQf/5rnuxfdwN5j8BcTFORya5kJI1ERERb26Ict2XrdVw+PljGN9SD4aXLKdkTUREJC15CsAtr8HAxZC/GJw55HREksuEOR2AiIhItlCxiesZo6HhTkciuYySNREREV+F5XE6AsmF9DOoiIiISBBTsiYiIiISxJSsiYiIiAQxJWsiIiIiQUzJmoiIiEgQU7ImIiIiEsSUrImIiIgEMWOtdToGyYaMMceB3z3MKgqkfnhe6mklgROZFFp6PMWXFeX4unx6y6U139s8X9oEnGsXp9rEn3Uyul18bSvtK4EvF6z7SmVrbakA15Xcylqrl14Z9gImpTcNWB9M8WVFOb4un95yac33Ns+XNnGyXZxqEyfbxde20r6SdW3iT1s52S565c6XfgaVjDbfx2lOyahY/C3H1+XTWy6t+d7mqU2ufZ2Mbhd/2sop2ld8+xyRTKefQSXLGWPWW2ubOh2HXEntEnzUJsFJ7SJZTSNr4oRJTgcgHqldgo/aJDipXSRLaWRNREREJIhpZE1EREQkiClZExEREQliStZEREREgpiSNXGcMaaqMeYDY8xsp2OR/zHG3G6MmWyM+coY08npeASMMbWNMRONMbONMUOcjkdcjDEFjTEbjDHdnI5FciYla5IpjDEfGmOOGWO2ppre2Riz0xizxxjzPIC1dp+1dqAzkeYufrbLXGvtw8AA4G4Hws0V/GyTHdbawcBdgG4dkUn8aZMUI4DPsjZKyU2UrElmmQJ0dp9gjAkFxgNdgDpAX2NMnawPLVebgv/t8mLKfMkcU/CjTYwxPYCfgCVZG2auMgUf28QY0wHYDvyR1UFK7qFkTTKFtfZH4FSqyc2BPSkjafHALOC2LA8uF/OnXYzLv4FvrLU/Z3WsuYW/+4q1dp61tjXQL2sjzT38bJP2QEvgHuBhY4y+VyXDhTkdgOQqFYCDbu+jgRbGmBLAa0AjY8wL1tpRjkSXe3lsF+BRoANQ1BhT3Vo70Yngcilv+0o7oBeQF1iQ9WHlah7bxFo7HMAYMwA4Ya1NdiA2yeGUrElWMh6mWWvtSWBwVgcjl3lrl7eBt7M6GAG8t8kyYFnWhiIpPLbJ5T+snZJ1oUhuo+FayUrRQCW39xWBww7FIv+jdgk+apPgozYRxyhZk6y0DqhhjLnBGJMH6APMczgmUbsEI7VJ8FGbiGOUrEmmMMZ8AqwCahpjoo0xA621icBw4DtgB/CZtXabk3HmNmqX4KM2CT5qEwk2epC7iIiISBDTyJqIiIhIEFOyJiIiIhLElKyJiIiIBDElayIiIiJBTMmaiIiISBBTsiYiIiISxJSsiYiIiAQxJWsiIiIiQUzJmojkGMaYL40xrxpjlhtjjhpjOjgdk4jItVKyJiI5ST3gtLU2EhgK9HM4HhGRa6ZkTURyBGNMAaAoMCZlUhhw2rGAREQyiJI1Eckp6gIbrLVJKe8bAFsdjEdEJEMoWRORnKIesMntfQPgF2dCERHJOErWRCSnqM+VyVo9NLImIjmAsdY6HYOIiIiIeKGRNREREZEgpmRNREREJIgpWRMREREJYkrWRERERIKYkjURERGRIKZkTURERCSIKVkTERERCWJK1kRERESC2P8DUvcZ5B6G/WIAAAAASUVORK5CYII=\n", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAmsAAAEcCAYAAACYg/MAAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAA+9ElEQVR4nO3dd3gU1f7H8fdJAqEHqQoBAcFQAqH3BJAiIB0UFBUQUZqiXqzXe/VexYoNBRUEsSDiD0U6KgJK70WKdBAQpElogYTk/P7YwA1hN9ksSWaTfF7Ps8+T2Zk5+51zzsx+c3aKsdYiIiIiIv4pwOkARERERMQzJWsiIiIifkzJmoiIiIgfU7ImIiIi4seUrImIiIj4MSVrIiIiIn5MyZqIiIiIH1OyJiIiIuLHlKz5EWNMmDFmvTHmjDHmUafjucwYM9EY87IT5aXXZxtjthhjmju1vo+fmab+kLyunIjZnxhj9hljWvm4brr2efEP19MnRJJKy/E1Pfqd18maMeYeY8waY8xZY8xhY8xcY0zT6/lwucZTwCJrbUFr7ShvV9IBKHXW2mrW2kXeLOuuPtOyfjryqT9cllExq79JVqB+mn0YY4KNMeONMfsT/3ldb4xpl2yZIsaYacaYc4nL3ePNPA+fl2rfyezvBK+SNWPME8C7wCtASaAsMAbonGGRpYExJsjpGNLJzcAWp4NIL9moXZySIf1B7SL+zuk+6vTnyzWCgANAMyAE+BfwjTGmXJJlRgOxuHKU3sCHxphqXsxLM0f6h7U2xReuijkL3JnCMlWARcApXF8unZLM2wcMBzYB0cAUIE/ivGeAqcnKeg8Ylfh3KeBb4BiwF3g0WblPJ5Z7EVdj1gbWA2eA/0v8rJeTrJNaeW7jTJxfBvgucd0TwAeplZmWugIWAPHAhcT6vtXNuk8DhxK3bzvQEvgCSABiEtd7Kknd7k5cdivQNQ3bWgtYl7juFODrZPWYWtnJ2yXF8pJtY2qf7ba+Sb0v7QNapbYNKdTnlfU9taE3dZvO/SG1ukoas7t2SbHv4qbPe6ofb/pqan0nSZxPJsZ5DhiP6wA7N3Gd+cANSZZ9NrGcv4FPubofJ93+1LbVqz4KBAOnE9vmbOLrcju1crN8asckC1RMMj3x8nwvYvZUx57e96m8NPRZn4/nnvposvkDgdm4vniPA38CrT3EmNJ+7On76JrPT2Fb+wEzk3zeLuCbJNMHgJpprNf03Dc8HqO8qIcU+6ybuB8EfgA+xLUf7gCqAsOAPxLbqpun9X19JcbePfHv/LiSsVuTzP8CeC2leT70naT94yBefKe4ORZ51R+uicuLCmkLXCLZjpNkfq7EjvockBu4LTGIsCRBrsK1oxYBtgEDE+fdDJwHCiVOBwKHgYa4Rv3WAv9OLLcCsAe4PUm5G3B9oeRNXGZ/YgfJBXRLbKDLBz5vyvMUZyCwEXgnseHzAE1TK9OHuloEPOhh3TBcB4BSidPlgFuSd4Qky9+ZuC0BQE9cO/dNXmzr5Xp8PDHeHkAcV3/BpFa2u3bxWF6SclNcNqX6JoW+5K6OvNiG5PW5D2jlRRt6rNt07g/etNOV7XDTLqntD277vKf6SUNf9VjvScpegetLqDRwFFcSVQtXorQAeCHJspsTt6kIsNTd9nuxrV730cTl+wE/JpneAUSl0EZuj0mJy7hN1ryI2W0dp/C+T+Wlpc9yHcdzd33UzWePAU7i2t8DgBeA+Sn0w3243489Hfuu+vxUtrUCrkQoALgpsZ0PJZZTAVfSEuBtvabnvpFS3KnVA170WTdxj8L1z1zLxDb/NlkZw4B1btablViH7l6zPH1e4rolcf2DVDlxuhYQk2yZ4cDMlOb50HeS9o+rlkmp/fjfscjr/nBNTKku4BoyPJLC/EjgCBCQ5L3JwItJgrw3ybw3gI+STC8B7k/8uzWwO/HvBsAfyT7rWeDTJOU+kGReFK5s1SQr++U0lOc2TqARrv8Gk/+nl2KZPtTVIjx/OVfEtXO2AnKl1rHcrL8B6OzFtkbh+o81aT0uI+WdNXnZydvFq/JSW9aLNnTbl7ypIzfb4ClZu67+no79IdV65dpkLWm7pFaXbvu8l3Xpsa+mVO9Jyu6dZPpb4MMk048A3ydZdmCSee3dtbkX25qmPg+8DbyV+Hd+XCNrRTy0kcdjUuK0p2QttZjd1nEK7/tUng991qfjubs+6uazl5BkJBfXl6MvyZqnY99Vn+/Fth7ANQrVCxiLK/mpjCuZn5HWfSG99o3U4k6pHvCiz7qJ8xdgeJLpl0iSbOFK4jalZdtTqZdcuEYRP07eVsmWG4DrGOpxng9954GUlvHUfvzvWORzf/DmnLUTQLEUfqMtBRyw1iYkeW8/rsz/siNJ/j4PFEgy/RVwd+Lf9yROg+u/tFLGmFOXX7j+UyiZZN0DyeI4ZBNrxs18b8rzFGcZYL+19hJX86bMpLypK7estbuAx4AXgaPGmK+NMaU8LW+Mud8YsyFJXOFAsSSLeNpWd/W4P41lp9YuV5WXhmVTq29PfekaXmyDJ9fb39NaVkpxeFuvl6Vlf/DU569hjOmdeOHRWWPM3JT6qpf1/leSv2PcTCetz6TbtB9XvSSX2ramtS6rA78l+fuwtfakm+VSOyalJMWYPdVxCnXvU3ketimlPns9x/PU6qc6rpGSy8Jx/dyUVintn8mPXSlt6y9Ac1wJzi+4EoNmia9fIG3H7XTcN7w9rrirB1/6bA1co2SXVXUz/XsqZXjFGBOA62fKWGBokllngULJFi+Ea0QxpXlp5bEuvGm/tH6PJ+VNsrYc13BjFw/z/wTKJFbiZWVxZefe+D+guTEmFOjK/3buA8Bea23hJK+C1tr2SdZN2qEOA6WNMSbJe2WS/O1NeZ4cAMq6SVjTWuZ11ZW19itrbVNcBz4LvH55VtLljDE3A+Nwdeai1trCuH4uSlo3nrirx7JpLDu1dimLe6ktm1p9e+pLV/FiG6y79RJdb39Pr7LSUq+XJT8Ap1SXnvp88nKw1k6y1hZIfLVLfO+avnqd/dKTpPt4WVx1mlxq25rWukyarEUk+Tu51I5J4PqSzJdk+kYvY/Z4PPDwvs/lJZNan72e4zl42PeMMeVxnUO2PcnbtXCNXniS0n7szTqpbevlZC0y8e9fSJasgXf1ms77RnofV5L32SsS486N61SAy2pydbvUwE07GdcdJc56eM11s7zhf+fpdbfWxiWZvQMIMsZUSvJeBK7z9VKa54mnvuOpf3rdfl7uZ9dINVmz1kbjOs9gtDGmizEmnzEmlzGmnTHmDWAlrt9mn0p8vznQEdcJuqmy1h7D9R/Jp7h25m2Js1YBp40xTxtj8hpjAo0x4caYeh6KWo7r54ihxpggY0xnoH6S+WktL6lVuDrxa8aY/MaYPMaYJj6U6XNdGdc9t24zxgTjSp5jErcXXP9dVUiyeH5cneBY4rr9cGX53liO6xzFRxPrsRtX12Nay06tvLQsm2J9p9CXkkttG5LXZ1LX1d/Tsay01Ks7qfVdT30eUq6flPrq9fRLT4YYY0KNMUVwjdRM8WFbva5LY0wxoDj/G9GpzNVfVEmldkwC15fYPYkxtcX1RZ9qzJ7qOIW696k8N9uUYp9Nx+N5cjWA35KNFtXCdV6lJyn2Uy+ktn/+ArTAdX7dQWAxrnO8i+I6QT8t9Zqe+8b1HldS67NJRZCkXYwxhXAlIZuSLXNNO1lr2yX5Jy/5q13y5XFdwFAF6GitjUlW1jlcF0P9N/F41QTX3Sq+SGleCtuV1r7jVfuloT9cw6tbd1hr3waeAJ5PDOYArgzye2ttLNAJaIfrqo8xuM5ZSMuw51e4fsO9MhJirY3H1cFq4rpy6DjwCa6rU93FGIvrZMj+uE5QvBfXUOxFX8pLVvbldSviurrlINDTxxh9ratgXFe2HMc1fF0C15cTwKvA88Y1/DrcWrsVeAvXjvcXrpGApV58RtJ67IvrJNmeuDr65flpKju18tL42d7U9zV9yc3npLYNV9Wnmxivt79fd1lpqVcP66dYl576fOLqHusnkdu+ej39MgVfAT/iOll9D67zvdK6rWmpyxq4zsO6/GVxCLjLGNPAzeemeExKNCwxtlO4zg/+3puY8Xw88FT3vpbnbptS67PXfTx346rRGWNMUVyjkJtTWCe1fpqi1LbVWrsD109sixOnT+Pqg0sTtxe8r9d02zfS6biSUp9NKoJrR9F2WWvPw5WfLcNJeQQ0VcY1cvUwrv5zxPxvBK53ksUG4zrx/yiuc/QGWWu3eDHPnTT1nTS0n1f9wR1jr/ppOnsxxqzEdfLop07HIiLpyxizD9cFGPOdjsVbOiZJVqM+6x+y1eOmjDHNjDE3Jg7f9sGV5c9zOi4RyZl0TJKsRn3WP2W3uzSHAd/guqplN9DDWnvY2ZBEJAfTMUmyGvVZP5StfwYVERERyeqy1c+gIiIiItmNkjURERERP5bdzlmTTFKsWDFbrlw5p8MQEclS1q5de9xaW9zpOCRrUbImaWKM6Qh0rFixImvWrHE6HBGRLMUYk9pj4USuoZ9BJU2stTOttQ+FhHh7L0sRERG5HkrWRERERPyYkjURERERP6Zz1kQkR4iLi+PgwYNcuHDB6VAkB8iTJw+hoaHkypXL6VAkG1CyJiI5wsGDBylYsCDlypXDGON0OJKNWWs5ceIEBw8epHz58k6HI9mAfgYVkRzhwoULFC1aVImaZDhjDEWLFtUorqQbJWuSuWLPw6b/Az3mTBygRE0yi/qapCcla5K51n8J3z0IX3SBv/c5HY2Io1588UVGjhzpcf7333/P1q1bMzEiEfFHStYkc9V7EO54Cw6ugTGNYMVHkBDvdFQifknJmoiAkjXJbAEBroRt8Aq4uQnMexomtIVj252OTCRTjBgxgrCwMFq1asX27a5+P27cOOrVq0dERATdu3fn/PnzLFu2jBkzZvDkk09Ss2ZNdu/e7XY5Ecn+jNW5Q5IGSR43NWDnzp3XV5i1sOkbV8IWew6aPQVNHoNAXeou6W/btm1UqVIFgP/M3MLWP0+na/lVSxXihY7VUlxm7dq19O3bl5UrV3Lp0iVq167NwIED6devH0WLFgXg+eefp2TJkjzyyCP07duXDh060KNHDwBOnDjhdjnxT0n73GXGmLXW2roOhSRZlEbWJE3S9XFTxkBETxiyGip3gAUvw9jm8Of66y9bxA8tXryYrl27ki9fPgoVKkSnTp0A2Lx5M5GRkVSvXp1JkyaxZcsWt+t7u5yIZC+6z5o4r0BxuPNTqN4DZj0B41pC46HQ/FnIldfp6CQbSm0ELCO5u0qwb9++fP/990RERDBx4kQWLVrkdl1vlxOR7EUja+I/Kt8BQ1ZCrd6w9D34sAnsW+p0VCLpJioqimnTphETE8OZM2eYOXMmAGfOnOGmm24iLi6OSZMmXVm+YMGCnDlz5sq0p+VEJHtTsib+JW9h6PQ+3D8dEi7BxPau0bYL6Xt+kYgTateuTc+ePalZsybdu3cnMjISgJdeeokGDRrQunVrKleufGX5Xr168eabb1KrVi12797tcTkRyd50gYH4pG7dunbNmjUZ+yGx52DBCFgxBgqVhg7vwK1tMvYzJdtyd7K3SEbSBQaSXjSyJv4rd35o+wr0/wmCC8BXd8J3D8G5E05HJiIikmmUrIn/K1MPHv4Vmj0Nm7+F0fVh83d6ZJWIiOQIStYkawgKhhbPuZK2wmVgaj/4ujecPux0ZCIiIhlKyZpkLSWrQf/50OZl2P0zjG4Aaz/TKJuIiGRbStYk6wkMgsaPwKBlcGN1mPkofN4JTu51OjIREZF0p2RNsq6it0CfmVxo+zb20HrXg+GXj9aD4UVEJFtRsiZpYozpaIwZGx0d7XQoLgEBvHCwHvcGj+LUjY3gh+dgfBs4us3pyERERNKFkjVJk3R9Nmg6aV21JHtjQ6i5qz9Tyr5Awsm98FEkLHodLsU6HZ6IiMh1UbImWV6rqiX58YlmPNCkAs/uDOP2uDc5VPp2WPQKjG0Gh9Y6HaLIFaNGjaJKlSr07t2bxo0bA3Dq1CnGjBlzZZnk09568cUXGTlyZLrFmh6OHj1K7dq1efbZZ+nWrRsJCQlOhySS5ShZk2yhQHAQ/+5YlelDmhIcUoImO+/h3eL/5dK5v+GTVvDDPyH2vNNhijBmzBjmzJnDpEmTWLZsGZB+yZo/Wr16NXfffTevvvoqJUqU4MQJ3dRaJK2UrEm2Uj00hO8HN+H5O6ow9q8wGp1+ha03dYXlH8CHjWHvYqdDlBxs4MCB7Nmzh06dOvHOO+9QoEABAJ555hl2795NzZo1efLJJ6+ZBvjyyy+pX78+NWvW5OGHHyY+3nUhzYgRIwgLC6NVq1Zs377d7edu3LiRqKgoqlatSkBAAMYYXnjhBZ+2YerUqTRs2JCIiAiaNm3KsWPHAOjatSvPP/88kZGR3HjjjcyfPx9wJWsREREAREdHU7x4cZ8+VyRHs9bqpVeaX3Xq1LH+7sDJc/aBT1fZm5+eZZ9+8wN7YWR1a18oZO2MR62NOeV0eJLJtm7d6nQI1lprb775Znvs2DFrrbX58+e31lq7d+9eW61atSvLJJ/eunWr7dChg42NjbXWWjto0CD72Wef2TVr1tjw8HB77tw5Gx0dbW+55Rb75ptvXvV5MTExNiwszK5cudJaa+3zzz9vhw8fbhMSEq5armnTpjYiIuKa108//XTVcsePH7/y94svvmg/+OADa621FStWvPLZ3377re3bt6+11tpevXrZoUOH2oEDB9r58+f7WGtZk7s+B6yxfnAM1ytrvYKcThZFMkroDfn4pE9d5m0+wgszgplx9kXGlfmJxus+x+z40fVg+LC2Tocpkqqff/6ZtWvXUq9ePQBiYmIoUaIEJ0+epGvXruTLlw+ATp06XbPu/PnzqV27NvXr1wegRo0azJs3D2PMVcstXuzdqPPEiROZMmUKFy9e5MiRI7zyyiucP3+e6OhoHn/8cQAuXbpE4cKFAQgKCuL999/3abtFxEXJmmRrxhjaVb+JJpWK8ea87dy7Mpjm+SMYFTCegpN7QngPaPc65C/mdKgiHllr6dOnD6+++upV77/77rvXJF3Jbd68merVq1+ZXrduHbVr175mucjISM6cOXPN+yNHjqRVq1YAfP7556xatYoFCxZQoEABoqKiqFatGlu2bKFOnToEBgYCsGnTJsLDwwH44osv0raxInINnbMmOUKhPLl4qUs43w5qzOECVan91/NMv6EPdut014Phf5uqR1aJYwoWLHhVopR8umXLlkydOpWjR48CcPLkSfbv309UVBTTpk0jJiaGM2fOMHPmzGvKLlq0KJs2bQJgx44dfPfdd/Tq1eua5RYvXsyGDRuueV1O1AB+++03GjduTIECBfj2229ZtmwZ1atXZ/PmzdSsWfPKcps2baJGjRrXXS8i4qJkTXKU2mVvYOYjTXmibThPH29H50uvcjRXKfi2P0zuBdGHnA5RcqCiRYvSpEkTwsPDefLJJ6+Zrlq1Ki+//DJt2rShRo0atG7dmsOHD1O7dm169uxJzZo16d69O5GRkdeUfffdd3P27FnCw8N56KGHmDx5MkWLFvUpzj59+jBq1CgiIyPZsWMHFSpUIH/+/Pz2229XJWubN2++MrImItfPWI0miA/q1q1r16xZ43QY1+WPE+f55/e/sXTnUf5ZdBH9LnxJQGAuaPNfqN0XAvS/THaybds2qlSp4nQYkoO463PGmLXW2roOhSRZlL6NJMcqWzQfnz9Qn3d61ebDi225LeY19uUJg1mPux4Mf2K30yGKiIgoWZOczRhD55qlmf9EMxrWqUPzo4/zWtBgLh1a77ov29JREH/J6TBFRCQH09WgIkDhfLl5rXsNutUO5blpBZl2tCrjin5FjZ/+BVumQecPoGQ1p8OU9DL3GTjyW/qWeWN1aPda+pYpIoJG1kSuUr98EWY/2pTerRvSI/oRhvMYF47vxX4cBQtfgUsXnQ5RRERyGI2siSQTHBTIoy0r0aHGTfxzWhEa7qnCe4Wn0OyX12HrdOj0AZSp53SYcj00AiYiWYhG1kQ8qFC8AF8NaMDzdzblsdhB9I97itPRJ7HjW8O85yD2nNMhiohIDqBkTdLEGNPRGDM2Ojra6VAyhTGGHnVC+fkfzSkc0YHGp19hetDtsGI0jGkEexY5HaJkMTExMTRr1uzKg9jTywMPPECJEiVSvL/ZqVOn6NGjB5UrV6ZKlSosX74cgAsXLlC/fn0iIiKoVq3aNQ95j4+Pp1atWnTo0CFdY05P8+bNIywsjIoVK/LaaymPnLrbnvfee4/w8HCqVavGu+++e9Xy7ubFxsYSFRXFpUu6AEkynpI1SRNr7Uxr7UMhISFOh5KpiuTPzVt3RTD2wRa8l2cQd138F0fPx8PnnWHGIxBzyukQJYuYMGEC3bp1u/JopvTSt29f5s2bl+Iyw4YNo23btvz+++9s3Ljxyj3AgoODWbBgARs3bmTDhg3MmzePFStWXFnvvffe8+t71MXHxzNkyBDmzp3L1q1bmTx5Mlu3bvW4fPLt2bx5M+PGjWPVqlVs3LiRWbNmsXPnzhTn5c6dm5YtWzJlypQM3z4RJWsiadC4YjHmDoukYYuO3HZuBJ/SmYR1X2JHN4DfZzsdnmQBkyZNonPnzgBMnTqVhg0bEhERQdOmTTl27JjP5UZFRVGkSBGP80+fPs2vv/5K//79AcidO/eVh60bYyhQoAAAcXFxxMXFXXnm6MGDB5k9ezYPPvigV3H06tWLnj170qBBA26++WZmz874/WLVqlVUrFiRChUqkDt3bnr16sX06dPdLutue7Zt20bDhg3Jly8fQUFBNGvWjGnTpqU6r0uXLkyaNCnDt09EyZpIGuXJFcgTbcKY9mhLfio9mM4X/8uuc8Hw9T0kfNMXzh51OkTxU7GxsezZs4dy5coB0KJFC1asWMHGjRtp3bo133zzzVXLR0ZGUrNmzWte8+fPT/Nn79mzh+LFi9OvXz9q1arFgw8+yLlz/zvvMj4+npo1a1KiRAlat25NgwYNAHjsscd44403CPDyiR4bN26kQoUKrFy5kkmTJvGf//wnw7bpskOHDlGmTJkr06GhoRw65P7Rce62Jzw8nF9//ZUTJ05w/vx55syZw4EDB7yat3r1ap/jFvGWrgYV8VGlkgX5akBDlu6qyDPzatDo8Bc8unUa7FxA0B1vEBDRExJHJ0QAjh8/fmU0C2DixIlMmTKFixcvcuTIEV555ZWrll+8eHG6ffalS5dYt24d77//Pg0aNGDYsGG89tprvPTSSwAEBgayYcMGTp06RdeuXdm8eTP79u2jRIkS1KlTh0WLFqX6GTExMRw/fvzKOW9Vq1bl77//vq5tatWqFUeOHLnm/REjRlwZoXT32ETjZt+bNWuW2+2pUqUKTz/9NK1bt6ZAgQJEREQQFBSU6rzAwEBy587NmTNnKFiwYJq2SyQtlKyJXKcmFYvReEgUP2+rzNC5zXk4+l3qfP8wx1ZMoliv0ZjCZZ0OUfxE3rx5uXDhAgCff/45q1atYsGCBRQoUICoqCiqVbv6xsuRkZGcOXPmmnJGjhxJq1at0vTZoaGhhIaGXhkx69Gjh9sT8QsXLkzz5s2ZN28eJ06cYMaMGcyZM4cLFy5w+vRp7r33Xr788ku3n7F582YqVapEnjx5AFi3bh0RERHXtU3ejLiFhoZeGe0C10+dpUqVuma5pUuXetye/v37X/mJ+LnnniM0NPTKeinNu3jx4pXtFckw1lq99Erzq06dOlauFR+fYGes/8OOGvGEPffv4vb8iyXtjlnv2IT4S06HluNt3brV6RCstdaGhobamJgYO3z4cPvuu+9aa62dOnWqDQwMtGfPnr2usvfu3WurVavmcX7Tpk3t77//bq219oUXXrDDhw+31lp79OhR+/fff1trrT1//rxt2rSpnTlz5lXrLly40N5xxx1Xpm+77TZ78ODBq5YZN26cLV26tI2JibFnz561jRs3tkuWLLmubfJGXFycLV++vN2zZ4+9ePGirVGjht28eXOK6yTfnr/++staa+3+/fttWFiYPXnyZKrzjh8/bitXruzxM9z1OWCN9YNjuF5Z66Vz1kTSUUCAoWPNMgx6+k0WtpzOZnMrlVa/wLbXoti0cY3T4YkfaNOmDUuWLKFPnz6MGjWKyMhIduzYQYUKFcifP7/P5d599900atSI7du3Exoayvjx4wFo3749f/75JwDvv/8+vXv3pkaNGmzYsIHnnnsOgMOHD9OiRQtq1KhBvXr1aN26dYq36UhISGDXrl3XXNCwceNGevfuTfPmzalXrx6DBg2iSZMmPm+Tt4KCgvjggw+4/fbbqVKlCnfddddVo5RJ68CT7t27U7VqVTp27Mjo0aO54YYbUp23cOFC2rdvnzEbJZKEsfba3/pFUlO3bl27Zo2Sj9RcjLvEmu8/oPqWN8htY5l+Qx+qdX+O8DJFnQ4tx9m2bZtf3H5i/fr1vP3223zxxRdOh+KzzZs3M2HCBN5+++2r3o+KimLcuHGEhYU5FFnm6tatG6+++qrH7XXX54wxa621dTMjPsk+NLImkoGCcwXR5M7HyPXoav4s3pSepz7BjruNl8d/w86/rj1vR7K/WrVq0aJFi3S/KW5mCg8PvyZRA9i9ezeVKlVyIKLMFxsbS5cuXXJMYirO0sia+EQjaz6wlvMbv8POHk7u2Gg+ju/IH9UGM7RNdcoWzed0dNmev4ysSc6hkTVJL7oaVCSzGEO+mt3h1hZcnP00Q7d8w+7fV/GPzQ9TqW5LHr2tEjeG6KoyERG5mn4GFcls+YoQfOc46P0t5QoFMCXXi4Ste5m2b87l5VlbOXH2otMRioiIH1GyJuKUSq0IHLqCgPoD6BM4jwV5n2HHshlEvbGQt37cTnRMnNMRZjs67UMyi/qapCclayJOCi4I7d+EfvMoUqggn+d+lU8Kf8pnCzYQ9cZCRi/cxfnYS05HmS3kyZOHEydO6EtUMpy1lhMnTuhmuZJudIGB+EQXGGSAuAvwy+uw9D3i8hRhbIHBvHkgjGIFghnS4hbuaVCW4KBAp6PMsuLi4jh48OCVJwiIZKQ8efIQGhpKrly5rnpfFxiIL5SsiU+UrGWgPzfAjKFw5Df+vrkdz164n3n7LaVC8vBoy0r0qBNKUKAGxUWyIiVr4gsd8UX8TamaMGAhtPw3NxxcwEfRg/ixxUFKFAzmme9+o9XbvzB9wyESEvSPlohITqBkTcQfBeaCyH/AwCVQvDK3Ln+KaSFv82X3m8iTK5BhX2+g/ajFLN113OlIRUQkgylZE/FnxW+FfnOh3ZuYP1bQ9McOzGm4jVG9IjgXe4nen6zkwc9Ws+fYWacjFRGRDKJz1sQnOmfNAaf+gJmPwe6foUxDLt7xHp9uz8UHC3ZxIS6e+xuVY1jLSoTky5VqUSLiDJ2zJr7QyJpIVlG4LNz7LXT5EI79TvC4KAYGfM/Cx5twZ90yTFy2l2YjF/LZsn3ExSc4Ha2IiKQTjayJTzSy5rAzf8HcJ2HrdLixOnQezVZbnpdnb2XZ7hNULFGAf95RhRZhJZyOVESS0Mia+EIja5ImxpiOxpix0dHRToeSsxUsCXd9Dnd94Urcxrag6ta3mdQ3gnH31yU+wdLv09X0mbCKnX+dcTpaERG5DhpZE59oZM2PxPwNPzwPG76EohWh0wfElm7A58v38d7POzkfG0/vBmV5rNWtFMmf2+loRXI0jayJL5SsiU+UrPmh3Qtg5jDXhQj1BkCrFzh5KZh35+9g0so/yJc7kGEtK3F/o3LkDtKguogTlKyJL5SsiU+UrPmpi2dhwUuw8mMICYWO70LFVuz86wwvz97GLzuOUb5Yfp5tV5nWVUtijHE6YpEcRcma+EL/XotkJ8EFoN3r8MAPkCsvfNkdpg2kUsE4PnugPp/2q0dggOGhL9bS+5OVbDt82umIRUQkFRpZE59oZC0LiLsAv74JS9+FvDdA+5FQrQtx8Ql8tfIP3pm/g9MxcfSsV4YnWodRvGCw0xGLZHsaWRNfKFkTnyhZy0IOb3I9GP7wRqjS0ZW0FbyR6PNxjFqwk8+W7SNPrkCGtKhIvyblyJMr0OmIRbItJWviCyVr4hMla1lM/CVY/gEsehWCguH2V6BmbzCGPcfO8sqcbczfdpQyRfLybLsqtAu/UeeziWQAJWviCyVr4hMla1nU8V0w4xH4YxlUaA4d34MbygGwZOdxXp69ld+PnKF+uSL8q0NVqoeGOBquSHajZE18oWRNfKJkLQtLSIC1E+CnF8AmQMsXoP4ACAgkPsEyZfUB3vpxOyfPx9KtVihPtQ2jZKE8Tkctki0oWRNfKFkTnyhZywZOHYBZj8OunyC0PnT+AIqHAXD6QhyjF+7i0yX7CAwwDGp+CwMiK5A3t85nE7keStbEF0rWxCdK1rIJa2HTNzDvaYg9B82egiaPQWAuAP44cZ5X525j7uYjlArJw9PtKtMpopTOZxPxkZI18YWSNfGJkrVs5uwxmPsUbPkOSoa7RtlK1boye8WeE7w0aytb/jxNzTKFebNHDSqVLOhgwCJZk5I18YVuiisiUKA43Pkp9PoKzh2HcS3hp39DXAwADSsUZcbQprzRowYH/z5PtzHLWLLzuMNBi4jkDErWROR/Kt8BQ1ZCrd6w9D34sAnsWwpAYIDhrrplmDG0KaVvyEvfT1fx9ao/HA5YRCT7U7ImIlfLWxg6vQ/3T4eESzCxPcx6Ai64Hk1VqnBe/m9gI5pULMYz3/3G6/N+JyFBp1OIiGQUJWsi4l6F5jB4OTQcAmsmwJhGsONHAArmycX4PnW5p0FZPly0m0cmr+dCXLyz8YqIZFNK1kTEs9z5oe0r0P8n10Piv7oTvnsIzp0gKDCAEV3C+Wf7KszZfJi7x63g+NmLTkcsIpLtKFkTkdSVqQcP/wrNnoHN38Lo+rD5WwwwIKoCH/auzbbDp+k6Zim7jp5xOloRkWxFyZqIeCcoGFo860raCpeBqQ/A1/fA6cO0Db+Jrx9qRExsPN3GLGPZbl0pKiKSXpSsiUjalKwG/edDm5dh9wIY3QDWfkbN0BCmDW5CyUJ5uH/8KqauPeh0pCIi2YKSNRFJu8AgaPwIDFoGN1aHmY/C550owxGmDmpMwwpFGf5/G3n7x+3oxtsiItdHyZqI+K7oLdBnJnR4F/7cAGMaE7L+Yz7tU5uedcswasEuHpuyQVeKiohcByVrInJ9AgKgbj8YvAIqNIMf/0muibfzWmQgT7UNY/qGP7lv/EpOnot1OlIRkSxJyZqIpI+Q0nD319B9PPy9D/NxMwYzldG9qrHxYDTdxixl7/FzTkcpIpLlKFkTkfRjDFTvAUNWQbUusOhV7lh2N9O7BHP6wiW6jlnKqr0nnY5SRCRLUbImIukvfzHo/olrpC3mFFVmd2dhjZ+4KZ/l3k9W8v36Q05HKCKSZShZE5GME9YOhqyA2n0IWf8xs4Ke4t4b9/PYlA28N3+nrhQVEfGCkjURyVh5QqDju9BnFoHG8O8TTzP5xsl8Mn8D//i/jcReSnA6QhERv6ZkTUQyR/lI133ZGj9Cw+jZLC/0LNEbZnD/hJWcOq8rRUVEPFGyJiKZJ3c+aPMy5sH5FChcgvG53+Leg/+l3+i57D+hK0VFRNxRsiYima90HXhoETR/jvZBq5lwbggff/A6a/edcDoyERG/o2RNRJwRlBuaP03AwMXkLVmRV+x7nJ7QnZ9WrHM6MhERv6JkTUScVaIKeR7+mXMtXqJR4FYazm3PokmvYxP0iCoREVCyJiL+ICCQ/M0exQxazp/5q9B85yvsGXkbcUd3Oh2ZiIjjlKyJiN8ILnELtw7/mR9v+SfFz+3AjmlMzKK3If6S06GJiDhGyZoAYIypYIwZb4yZ6nQskrOZgADa3PcUv7aZzS8JNci76D/EfnwbHNnsdGgiIo5QspYNGGMmGGOOGmM2J3u/rTFmuzFmlzHmmZTKsNbusdb2z9hIRbzXoUlt8t//NcN5jLNH92I/bgYLRsCli06HJiKSqZSsZQ8TgbZJ3zDGBAKjgXZAVeBuY0xVY0x1Y8ysZK8SmR+ySOoaVyzOwMFPcl/e95ke3wh+fQM+joIDq50OTUQk0yhZywastb8CJ5O9XR/YlThiFgt8DXS21v5mre2Q7HU004MW8VLFEgX4bEg7PrvxWfrGPkX0qb+x41vDvGchVjfSFZHsT8la9lUaOJBk+mDie24ZY4oaYz4CahljnvWwzEPGmDXGmDXHjh1L32hFUlCsQDCTBzSkXMMuNDkzglm528GKMTCmEexZ5HR4IiIZSsla9mXcvGc9LWytPWGtHWitvcVa+6qHZcZaa+taa+sWL1483QIV8UaeXIG82Kkao/pG8WL8A/S+9AKnY4HPO8P0oRBzyukQRUQyhJI1P2GMeSGdizwIlEkyHQr8mc6fIZLpbqtckrmPRRJQvgn1Tv6HuSG9sBu+gtENYNssp8MTEUl3Stb8xwvGmNeNMeOMMYOMMTdcZ3mrgUrGmPLGmNxAL2DG9Ycp4rwSBfPwWb/6PHlHBMOOd6FP4KuczXUDTOkN3/SBszoNU0SyDyVr/sMCF4AfcI2ILTPGRHizojFmMrAcCDPGHDTG9LfWXgKGJpa3DfjGWrslY0IXyXwBAYYHIyvw3eDGHMobRq0jz/JL6CDs9jkwuj5s/Bqsx1/+RUSyDGN1MPMLxpgt1tpqSaZvBT6y1t7mYFge1a1b165Zs8bpMEQAOB97iZdmbWPyqj9oVzKad/J+Qp4ja6Fia+jwDhQuk3ohIpnAGLPWWlvX6Tgka9HImv84boypc3nCWrsD8Luz+I0xHY0xY6Ojo50OReSKfLmDeLVbdT66tw7LzxSj7p/DWVftWez+ZTCmIawaBwkJTocpIuITjaz5icSfPL8G1gK/ATWAgtbaTo4G5oFG1sRfHY6O4fEpG1ix5yT3VTH8235Mrn2LoGxj6PQ+FKvodIiSg2lkTXyhkTU/Ya3dCNQEJie+tRC427GARLKom0LyMunBhjzVNozJ26H54UfZ0+RNOLoFPmwMS97Rg+FFJEvRyJr4RCNrkhVsOHCKYV+v58DJ8zzd9AYGnBlDwO8z4aYI6PQB3FTD6RAlh9HImvhCI2sikm3VLFOY2Y9G0rVWKK8u/pseJwdxvP04OH0YxjaHn/8LcRecDlNEJEVK1kQkWysQHMRbd0Uw6u5a7PzrLC1mhzCn2fdQoycsfgs+joQ/VjodpoiIR0rWRCRH6BRRijnDIrn1xoIM/m4fT8Q+zPm7voG4GJhwO8x5Ci6edTpMEZFrKFmTNNGtOyQrK1MkH1MeasiwlpX4fsMh2s7KzaZO86D+AFg11vVg+F0/Ox2miMhVlKxJmlhrZ1prHwoJCXE6FBGfBAUG8HjrW5nycCPiEyzdxm9idN6Hie87B4KC4ctu8P1gOH/S6VBFRAAlayKSQ9UrV4Q5wyJpG34jb/6wnXt+MPx5908Q+Q/Xo6rGNIStepyuiDhPyZqI5FgheXPx/t21GHlnBL8diqbd6NXMKzkAHloEBUrCN/fBlPvgzF9OhyoiOZiSNRHJ0Ywx9KgTypxHI7m5aD4GfrmOZ5fD+T4/QssXYMcPrgfDr5+kB8OLiCOUrImIAOWK5WfqwMYMan4LX68+QIcxK9lcoT8MWgolqsD0wa7z2f7e73SoIpLDKFkTEUmUOyiAp9tWZtKDDTh/MZ6uY5byybZAEvrMhvYj4cAq1xWjKz/Wg+FFJNMoWZM00a07JCdofEsx5g6LpEVYCV6evY0+E9dwtPJ9MHg5lG0Ic5+CT9vBsR1OhyoiOYCeDSo+0bNBJSew1jJ51QH+O2sLuQMDeP6OqtxZpzRm0xSY9wzEnYdmT0OTYRCYy+lwJQvQs0HFFxpZExHxwBjDPQ3KMndYFJVvKsRT327i3gmr+KNMZxi6GsLaw4KXYFwL+HOD0+GKSDalZE1EJBXli+Xn6wENGdE1nI0Homnz7i98sv4s8T0mQs8v4exRGHcbzH/R9fgqEZF0pGRNRMQLAQGG3g1u5qcnomhasRgvz95Gtw+X8fsNzWDISqh5Nyx5Bz5qCvuXOR2uiGQjStZERNLgppC8jLu/Lu/fXYuDJ8/TYdQS3l58lIt3jIL7vof4WNfFB7OHw8UzTocrItmAkjURkTQyxtAxohTzn2hGp4hSjFqwiztGLWFtUAQMWg4NBsHqT1y3+dg53+lwRSSLU7ImIuKjG/Ln5u2eNZnYrx4xsfH0+Gg5L/6wn3O3vQz9f4Rc+WBSd5g2UA+GFxGfKVkTEblOzcNK8MPjUdzf8GY+W76PNu/8yi8x5WHgYoh6En77P9cjq7ZM0yOrRCTNlKxJmuimuCLuFQgO4j+dw5k6sBF5cgXQZ8IqnvhuG383eMr1YPhCpeH/+sKUe+HMEafDFZEsRDfFFZ/oprginl2Ii2f0wl18uGg3hfPl4oWO1egQXhyzYgwsfAUCg+H2EVDrXjDG6XAlE+mmuOILjayJiKSzPLkC+UebMGY+0pRShfPyyOT1DPhyI0fCH4aBS+HGcJgxFL7oAn/vczpcEfFzStZERDJIlZsK8d2gxvyzfRWW7DpG67d/4avduUm4fybc8TYcXOu6YnTFh5AQ73S4IuKnlKyJiGSgoMAABkRV4IfHoggvHcJz037j7k9Wsbd8LxiyAm5u4nrO6ITb4ejvTocrIn5IyZqISCa4uWh+vhrQgNe7V2fr4dO0ffdXPtpwkUu9pkC3cXBiN3wcCb+8AZdinQ5XRPyIkjURkUxijKFnvbLMf6IZzW4tzmtzf6fLh8vYUux2GLIKKneAhSNcD4Y/tM7pcEXETyhZExHJZCUL5eHj++owpndtjkRfpNMHS3ljyQkudPkEen0F50/AJy3hx3/pwfAiomRNRMQJxhjaV7+J+U9E0bVWacYs2k379xazKrgRDF4Bte6DZaPgw8awb4nT4YqIg5SsiYg4qHC+3Iy8M4LPH6hPbHwCd328nH/9cJAzbd6C+2eATYCJd8Csx+HCaafDFREHKFkTEfEDUbcW54fHonigSXm+XLmfNu/8yoLYyjBoGTQaCmsnwpiGsOMHp0MVkUymZE3SRI+bEsk4+YOD+HfHqnw7qDEF8wTxwMQ1DJ26nR01n4H+P0FwQfjqLvh2AJw74XS4IpJJ9Lgp8YkeNyWSsWIvJTBmkeuRVRcvJRB1a3EebFSKyCOfYxa/BXlCoN0bEN5dj6zKQvS4KfGFkjXxiZI1kcxx8lwsX63cz2fL93PszEVuLVmAx6vHcfvulwk4vB7C2sMdb0GhUk6HKl5Qsia+ULImPlGyJpK5Ll6KZ9bGw4xfspeth09TIl8gb5ZdRtTBsZjAXNDmJajdR6Nsfk7JmvhCyZr4RMmaiDOstSzfc4IJS/Yyf9tRKgYe5aOQz6h4fj2Ui4ROo6BIBafDFA+UrIkvgpwOQEREvGeMofEtxWh8SzH2HDvLp0v30XFtCbok/My/9k8meHQjzG3PE9BoMAQEOh2uiKQDjayJTzSyJuI/Tp2P5atVfzB36VqGXfiQVoHrOR5SnYJ3fUhw6epOhydJaGRNfKFbd4iIZHGF8+VmcPOKfPv0nZzt+iVvFHgSc2ofAeOaseyTf3D0lG6mK5KVaWRNfKKRNRH/Za1l3badXJrzDA3O/swOG8qcCs/TpvUdVC1VyOnwcjSNrIkvlKyJT5SsiWQNf62eRp4fn6RA7HEmxLdjaZmHuS+qCi3CShAQoCtHM5uSNfGFkjXxiZI1kSzkQjQX5/6L4I2fcZCSPBn7IH8VqU+/puXpUTuUvLl1IUJmUbImvlCyJj5RsiaSBe1djJ3xCObvvfwQfDvDo+8kMF8I99QvS5/G5ShZKI/TEWZ7StbEF0rWxCdK1kSyqNjzsOgV7PLRxOUtzrhCjzByfwWCAgwdapSif9PyhJcOcTrKbEvJmvhCyZqkiTGmI9CxYsWKA3bu3Ol0OCLiq0NrYfojcHQL5yp1YUzeB5m44RznYuNpUL4ID0ZWoGVlndeW3pSsiS+UrIlPNLImkg1cioWl78Ivb0BwQc63GsGksw2YuHw/h07FcEeNmxh9T22no8xWlKyJL5SsiU+UrIlkI0e3wfShcGgNVLqdS+1GMu9gEIXz5qZppWJOR5etKFkTX+imuCIiOV2JKtD/R7j9Vdi3mKCPGtPh4lya3lLE6chEBCVrIiICrueINhoMg5ZB6dow+wn4rCOc2O10ZCI5npI1ERH5nyLl4f7p0Ol9OPIbfNgYlr4H8Zecjkwkx1KyJiIiVzMGat8PQ1bCLS3hp3/D+FZwZLPTkYnkSErWRETEvUI3Qa9JcOdEiD4IY5vBghFw6aLTkYnkKErWRETEM2OgWlcYsgrCe8Cvb8BHkXBgldORieQYStZERCR1+YpAt4+h91SIPQfj28DcZ1x/i0iGUrImIiLeq9QaBi+Hev1h5YcwpiHsXuh0VCLZmpI1ERFJmzyF4I63oO8cCMgFX3SB6UMg5pTTkYlkS0rWRETEN+WawKCl0PRx2DAZRjeAbbOcjkok21GyJiIivsuVF1q9CAN+hvzFYUpv+KYPnD3qdGQi2YaSNRERuX6lasFDC+G2f8H2OTC6Pmz8GvT8aZHrpmRNRETSR2AuiBoOA5dAsVth2sMwqQecOuB0ZCJZmpI1SRNjTEdjzNjo6GinQxERf1U8DPrNg3ZvwP7lritGV42DhASnIxPJkpSsSZpYa2daax8KCQlxOhQR8WcBAdDgYddtPkLrwZzhMLE9HN/pdGQiWY6SNRERyTg33Az3TYPOY+DoVviwCSx+Ww+GF0kDJWsiIpKxjIFavWHIari1Dfz8H/jkNji8yenIRLIEJWsiIpI5CpaEnl/CXZ/D6cMwtjn8/F+Iu+B0ZCJ+TcmaiIhkrqqdYchKiOgFi9+CjyPhj5VORyXit5SsiYhI5stXBLqMgXu/c42sTbgd5jwFF886HZmI31GyJiIizqnY0nXFaP2HYNVYGNMIdv3sdFQifkXJmoiIOCu4ALR/Ax6YB0HB8GU3+H4wnD/pdGQifkHJmoiI+IeyDV1PP4j8h+tRVaMbwNbpTkcl4jglayIi4j9y5YGW/4aHFkHBG+Gb+2HKfXDmL6cjE3GMkjUREfE/N9WAAQug1Yuw4wfXg+HXT9KD4SVHUrImIiL+KTAXNH0cBi2FElVg+mDX+Wx/73c6MpFMpWRNRET8W7FK0HcOtB8JB1a5rhhd+bEeDC85hpI1ERHxfwEBUH+A6zYfNzeCuU/Bp23h2HanIxPJcErWREQk6yhcFnpPha4fw/Ed8FFT+HUkxMc5HZlIhlGyJiIiWYsxrkdVDVkFle+ABS/BuBbw5wanIxPJEErWREQkaypQAu6cCD0nwdmjMO42mP8ixMU4HZlIulKyJiIiWVuVDq4Hw9e8B5a84/ppdP8yp6MSSTdK1kREJOvLewN0/gDu+x7iY+HTdjD7H3DxjNORiVw3JWsiIpJ93NICBq+AhoNh9XgY3RB2/uR0VCLXRcmaiIhkL7nzQ9tXof+ProfET+oB3z2sB8NLlqVkTUREsqcy9eHhXyHqKdg8FT6oB5u/0yOrJMtRsiZpYozpaIwZGx0d7XQoIiKpCwqG2/4JD/0CIaEwtR9MuRdOH3Y6MhGvKVmTNLHWzrTWPhQSEuJ0KCIi3rsxHB78GVr/F3bNh9ENYN3nGmWTLEHJmoiI5AyBQdBkGAxaBjdWhxmPwOed4eRepyMTSZGSNRERyVmK3gJ9ZkKHd+DQOviwMSwfAwnxTkcm4paSNRERyXkCAqDuA66b6ZaLhB+ehfFt4Og2pyMTuYaSNRERyblCSsM9U6DbJ3ByD3wUCb+8AZdinY5M5AolayIikrMZAzXuhKGroWonWDgCxjaHQ2udjkwEULImIiLikr8Y9JgAvSZDzEn4pBX8+C+IPe90ZJLDKVkTERFJqnJ717lste+HZaPgoyawd7HTUUkOpmRNREQkuTwh0PE911Wj1sJnHWDmY3BBNwSXzKdkTURExJPyUa77sjUaCus+04PhxRFK1kRERFKSOx/cPgL6z4e8heH0IacjkhwmyOkAREREsoTQOq5njAbmcjoSyWGUrImIiHgrKLfTEUgOpJ9BRURERPyYkjURERERP6ZkTURERMSPKVkTERER8WNK1kRERET8mJI1ERERET+mZE1ERETEjxlrrdMxSBZkjDkG7HczKwRI/vC85O8VA45nUGipcRdfZpTj7fKpLZfSfE/zvGkTcK5dnGqTtKyT3u3ibVtpX/F9OX/dV2621hb3cV3Jqay1eumVbi9gbGrvAWv8Kb7MKMfb5VNbLqX5nuZ50yZOtotTbeJku3jbVtpXMq9N0tJWTraLXjnzpZ9BJb3N9PI9p6RXLGktx9vlU1supfme5qlNrn+d9G6XtLSVU7SvePc5IhlOP4NKpjPGrLHW1nU6Drma2sX/qE38k9pFMptG1sQJY50OQNxSu/gftYl/UrtIptLImoiIiIgf08iaiIiIiB9TsiYiIiLix5SsiYiIiPgxJWviOGNMBWPMeGPMVKdjkf8xxnQxxowzxkw3xrRxOh4BY0wVY8xHxpipxphBTscjLsaY/MaYtcaYDk7HItmTkjXJEMaYCcaYo8aYzcneb2uM2W6M2WWMeQbAWrvHWtvfmUhzljS2y/fW2gFAX6CnA+HmCGlsk23W2oHAXYBuHZFB0tImiZ4GvsncKCUnUbImGWUi0DbpG8aYQGA00A6oCtxtjKma+aHlaBNJe7s8nzhfMsZE0tAmxphOwBLg58wNM0eZiJdtYoxpBWwF/srsICXnULImGcJa+ytwMtnb9YFdiSNpscDXQOdMDy4HS0u7GJfXgbnW2nWZHWtOkdZ9xVo7w1rbGOiduZHmHGlskxZAQ+AeYIAxRt+rku6CnA5AcpTSwIEk0weBBsaYosAIoJYx5llr7auORJdzuW0X4BGgFRBijKlorf3IieByKE/7SnOgGxAMzMn8sHI0t21irR0KYIzpCxy31iY4EJtkc0rWJDMZN+9Za+0JYGBmByNXeGqXUcCozA5GAM9tsghYlLmhSCK3bXLlD2snZl4oktNouFYy00GgTJLpUOBPh2KR/1G7+B+1if9Rm4hjlKxJZloNVDLGlDfG5AZ6ATMcjknULv5IbeJ/1CbiGCVrkiGMMZOB5UCYMeagMaa/tfYSMBT4AdgGfGOt3eJknDmN2sX/qE38j9pE/I0e5C4iIiLixzSyJiIiIuLHlKyJiIiI+DElayIiIiJ+TMmaiIiIiB9TsiYiIiLix5SsiYiIiPgxJWsiIiIifkzJmoiIiIgfU7ImItmGMWaaMeZlY8xiY8wRY0wrp2MSEbleStZEJDsJB05ZayOBwUBvh+MREbluStZEJFswxuQDQoB3Et8KAk45FpCISDpRsiYi2UU1YK21Nj5xugaw2cF4RETShZI1EckuwoENSaZrAJucCUVEJP0oWROR7KI6Vydr4WhkTUSyAWOtdToGEREREfFAI2siIiIifkzJmoiIiIgfU7ImIiIi4seUrImIiIj4MSVrIiIiIn5MyZqIiIiIH1OyJiIiIuLHlKyJiIiI+LH/B405I2YNW02wAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] @@ -499,7 +499,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 10, "metadata": { "deletable": false, "nbgrader": { @@ -514,10 +514,37 @@ "task": false } }, - "outputs": [], + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "N = 10000 \td = 1 \tV_d = 2.000 \tdirect-sampled V_d = 1.000\n", + "N = 10000 \td = 2 \tV_d = 3.142 \tdirect-sampled V_d = 0.788\n", + "N = 10000 \td = 3 \tV_d = 4.189 \tdirect-sampled V_d = 0.530\n", + "N = 10000 \td = 4 \tV_d = 4.935 \tdirect-sampled V_d = 0.304\n", + "N = 10000 \td = 5 \tV_d = 5.264 \tdirect-sampled V_d = 0.165\n", + "N = 10000 \td = 6 \tV_d = 5.168 \tdirect-sampled V_d = 0.083\n", + "N = 10000 \td = 7 \tV_d = 4.725 \tdirect-sampled V_d = 0.039\n" + ] + } + ], "source": [ - "# YOUR CODE HERE\n", - "raise NotImplementedError()" + "import scipy.special\n", + "\n", + "def volume_d_ball(d):\n", + " \"\"\"Returns the volume of a d-dimensional unit ball.\"\"\"\n", + " return np.pi**(d/2)/scipy.special.gamma(d/2 + 1)\n", + "\n", + "N = 10000\n", + "for d in range(1, 8):\n", + " print(\"N = {}\".format(N),\n", + " \"\\td = {}\".format(d),\n", + " \"\\tV_d = {:.3f}\".format(volume_d_ball(d)),\n", + " \"\\tdirect-sampled V_d = {:5.3f}\".format(number_of_hits_in_d_dimensions(N, d)/N)\n", + " )\n", + "\n", + "# TODO: Results not yet as expected." ] }, {