From 64ba3ea088eb335a2cc56d5b36ac4a971b588427 Mon Sep 17 00:00:00 2001 From: Kees van Kempen Date: Thu, 8 Sep 2022 12:00:48 +0200 Subject: [PATCH] 01: Extend direct-sampling to d-sphere --- Exercise sheet 1/exercise_sheet_01.ipynb | 36 ++++++++++++++---------- 1 file changed, 21 insertions(+), 15 deletions(-) diff --git a/Exercise sheet 1/exercise_sheet_01.ipynb b/Exercise sheet 1/exercise_sheet_01.ipynb index a7ef35f..0247783 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.10028581e-01],\n", - " [3.20000000e+01, 3.00020182e-01],\n", - " [6.40000000e+01, 1.96900792e-01],\n", - " [1.28000000e+02, 1.47571487e-01],\n", - " [2.56000000e+02, 1.03273922e-01],\n", - " [5.12000000e+02, 6.88982665e-02],\n", - " [1.02400000e+03, 4.86827171e-02],\n", - " [2.04800000e+03, 3.75673313e-02]])" + "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]])" ] }, "execution_count": 3, @@ -368,7 +368,7 @@ "outputs": [ { "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAmsAAAEcCAYAAACYg/MAAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAA8DElEQVR4nO3dd3gU1f7H8fdJI3SkWQDp0iH0mgAKCkgHBUQRBZGmohex/LwXG+JVREVEBRFUUEFUFCleadKkKiCCIKBcQJDeA2nn98cG7hJ2k90lyWySz+t59nkyOzNnv3POmdlvzszOGGstIiIiIhKcQpwOQERERES8U7ImIiIiEsSUrImIiIgEMSVrIiIiIkFMyZqIiIhIEFOyJiIiIhLElKyJiIiIBDElayIiIiJBTMlaEDHGVDLG/GyMOW2MedjpeC4yxkw1xrzoRHnp9dnGmF+NMS2cWj/Az/SrP6SsKydiDibGmD+NMa0CXDdd+7wEh6vpEyLu/Dm+pke/8zlZM8bcZYxZb4w5Y4w5YIyZb4xpdjUfLlcYASy11ua31o7zdSUdgNJmra1mrV3qy7Ke6tOf9dNRQP3hooyKWf1NsgL10+zDGJPLGDPZGLMn+Z/Xn40xbVMsU9gY85Ux5mzycnf5Ms/L56XZdzL7O8GnZM0Y8xjwBvAScC1wIzAB6JRhkfnBGBPmdAzppDTwq9NBpJds1C5OyZD+oHaRYOd0H3X68+UKYcBeoDlQEPgnMNMYU8ZtmbeBOFw5Sm/gHWNMNR/m+c2R/mGtTfWFq2LOAHekskwVYClwAteXS0e3eX8Cw4HNwElgBhCZPO9JYFaKst4ExiX/fQPwBXAY+AN4OEW5TySXewFXY9YBfgZOA58nf9aLbuukVZ7HOJPnlwK+TF73KDA+rTL9qStgMZAInE+u75s8rPsEsD95+7YDtwAfA0lAbPJ6I9zqdlfysluBLn5sa23gp+R1ZwCfpajHtMpO2S6plpdiG9P6bI/1Tdp96U+gVVrbkEp9XlrfWxv6Urfp3B/Sqiv3mD21S6p9Fw993lv9+NJX0+o7bnE+nhznWWAyrgPs/OR1FgLXuC37VHI5x4EpXN6P3bc/rW31qY8CuYBTyW1zJvl1sZ1aeVg+rWOSBSq4TU+9ON+HmL3Vsbf3AyrPjz4b8PHcWx9NMX8gMBfXF+8R4C+gtZcYU9uPvX0fXfH5qWzrfcAct8/bCcx0m94LRPlZr+m5b3g9RvlQD6n2WQ9x9we+A97BtR/uAKoCjwD/TW6rrt7WD/SVHHu35L/z4krGbnKb/zHwcmrzAug77v1jHz58p3g4FvnUH66Iy4cKaQMkkGLHcZsfntxRnwYigJuTg6jkFuRaXDtqYWAbMDB5XmngHFAgeToUOAA0wjXqtwH4V3K55YDdwG1u5W7E9YWSO3mZPckdJBzomtxAFw98vpTnLc5QYBPwenLDRwLN0iozgLpaCvT3sm4lXAeAG5KnywDlU3YEt+XvSN6WEKAHrp37eh+29WI9Ppocb3cgnsu/YNIq21O7eC3PrdxUl02tvkmlL3mqIx+2IWV9/gm08qENvdZtOvcHX9rp0nZ4aJe09gePfd5b/fjRV73Wu1vZq3F9CZUADuFKomrjSpQWAyPdlt2SvE2FgZWett+HbfW5jyYvfx/wH7fpHUBMKm3k8ZiUvIzHZM2HmD3WcSrvB1SeP32Wqziee+qjHj57AnAM1/4eAowEFqbSD//E837s7dh32eensa3lcCVCIcD1ye28P7mccriSlhBf6zU9943U4k6rHvChz3qIexyuf+ZuSW7zL1KU8Qjwk4f1vk2uQ0+vb719XvK61+L6B6ly8nRtIDbFMsOBOanNC6DvuPePy5ZJrf3437HI5/5wRUxpLuAaMjyYyvxo4CAQ4vbep8CzbkHe7TbvFeBdt+kVQJ/kv1sDu5L/bgj8N8VnPQVMcSv3frd5MbiyVZOi7Bf9KM9jnEBjXP8NpvxPL9UyA6irpXj/cq6Aa+dsBYSn1bE8rL8R6OTDtsbg+o/VvR5XkfrOmrLslO3iU3lpLetDG3rsS77UkYdt8JasXVV/T8f+kGa9cmWy5t4uadWlxz7vY1167aup1btb2b3dpr8A3nGbfgiY7bbsQLd57Ty1uQ/b6lefB8YCryX/nRfXyFphL23k9ZiUPO0tWUsrZo91nMr7AZUXQJ8N6HjuqY96+OwVuI3k4vpyDCRZ83bsu+zzfdjWvbhGoXoCE3ElP5VxJfPf+LsvpNe+kVbcqdUDPvRZD3H+AAx3m34Bt2QLVxK32Z9tT6NewnGNIr6Xsq1SLPcArmOo13kB9J37U1vGW/vxv2NRwP3Bl2vWjgJFUzlHewOw11qb5PbeHlyZ/0UH3f4+B+Rzm/4E6JX8913J0+D6L+0GY8yJiy9c/ylc67bu3hRx7LfJNeNhvi/leYuzFLDHWpvA5Xwp050vdeWRtXYnMAx4FjhkjPnMGHODt+WNMX2MMRvd4qoOFHVbxNu2eqrHPX6WnVa7XFaeH8umVd/e+tIVfNgGb662v/tbVmpx+FqvF/mzP3jr81cwxvRO/uHRGWPM/NT6qo/1/rfb37Eept3r032b9uCql5TS2lZ/67IG8Ivb3westcc8LJfWMSk1qcbsrY5TqfuAyvOyTan12as5nqdVPzVwjZRcVB3X6SZ/pbZ/pjx2pbatPwAtcCU4P+BKDJonv34A/47b6bhv+Hpc8VQPgfTZmrhGyS6q6mH6tzTK8IkxJgTXaco4YKjbrDNAgRSLF8A1opjaPH95rQtf2s/f73F3viRrP+IabuzsZf5fQKnkSrzoRlzZuS8+B1oYY0oCXfjfzr0X+MNaW8jtld9a285tXfcOdQAoYYwxbu+Vcvvbl/K82Qvc6CFh9bfMq6ora+0n1tpmuA58Fvj3xVnuyxljSgOTcHXmItbaQrhOF7nXjTee6vFGP8tOq11uxLO0lk2rvr31pcv4sA3W03rJrra/p1dZ/tTrRSkPwKnVpbc+n7IcrLXTrbX5kl9tk9+7oq9eZb/0xn0fvxFXnaaU1rb6W5fuyVott79TSuuYBK4vyTxu09f5GLPX44GX9wMuL4W0+uzVHM/By75njCmL6xqy7W5v18Y1euFNavuxL+ukta0Xk7Xo5L9/IEWyBr7VazrvG+l9XEnZZy9JjjsC16UAF0VxebvUxEM7GdcdJc54ec33sLzhf9fpdbPWxrvN3gGEGWMqur1XC9f1eqnN88Zb3/HWP31uPx/3syukmaxZa0/ius7gbWNMZ2NMHmNMuDGmrTHmFWANrnOzI5LfbwF0wHWBbpqstYdx/UcyBdfOvC151lrglDHmCWNMbmNMqDGmujGmvpeifsR1OmKoMSbMGNMJaOA239/y3K3F1YlfNsbkNcZEGmOaBlBmwHVlXPfcutkYkwtX8hybvL3g+u+qnNvieXF1gsPJ696HK8v3xY+4rlF8OLkeu3J5Pfpbdlrl+bNsqvWdSl9KKa1tSFmf7q6qv6djWf7Uqydp9V1vfR5Sr5/U+urV9EtvhhhjShpjCuMaqZkRwLb6XJfGmKJAMf43olOZy7+o3KV1TALXl9hdyTG1wfVFn2bM3uo4lboPqDwP25Rqn03H43lKNYFfUowW1cZ1XaU3qfZTH6S1f/4AtMR1fd0+YDmua7yL4LpA3596Tc9942qPK2n1WXe1cGsXY0wBXEnI5hTLXNFO1tq2bv/kpXy1Tbk8rh8wVAE6WGtjU5R1FtePoZ5PPl41xXW3io9Tm5fKdvnbd3xqPz/6wxV8unWHtXYs8BjwTHIwe3FlkLOttXFAR6Atrl99TMB1zYI/w56f4DqHe2kkxFqbiKuDReH65dAR4H1cv071FGMcrosh++G6QPFuXEOxFwIpL0XZF9etgOvXLfuAHgHGGGhd5cL1y5YjuIavi+P6cgIYDTxjXMOvw621W4HXcO14f+MaCVjpw2e412NfXBfJ9sDV0S/O96vstMrz87N9qe8r+pKHz0lrGy6rTw8xXm1/v+qy/KlXL+unWpfe+nzy6l7rJ5nHvno1/TIVnwD/wXWx+m5c13v5u63+1GVNXNdhXfyy2A/caYxp6OFzUz0mJXskObYTuK4Pnu1LzHg/Hnir+0DL87RNafXZqz6ee3DZ6IwxpgiuUcgtqayTVj9NVVrbaq3dgesU2/Lk6VO4+uDK5O0F3+s13faNdDqupNZn3dXiylG0ndbac3DptGV1Uh8BTZNxjVw9iKv/HDT/G4Hr7bbYYFwX/h/CdY3eIGvtrz7M88SvvuNH+/nUHzwx9rJT09mLMWYNrotHpzgdi4ikL2PMn7h+gLHQ6Vh8pWOSZDXqs8EhWz1uyhjT3BhzXfLw7b24svwFTsclIjmTjkmS1ajPBqfsdpfmSsBMXL9q2QV0t9YecDYkEcnBdEySrEZ9Nghl69OgIiIiIlldtjoNKiIiIpLdKFkTERERCWLZ7Zo1ySRFixa1ZcqUcToMEZEsZcOGDUestcWcjkOyFiVr4hdjTAegQ4UKFVi/fr3T4YiIZCnGmLQeCydyBZ0GFb9Ya+dYawcULOjrvSxFRETkaihZExEREQliStZEREREgpiuWRORHCE+Pp59+/Zx/vx5p0ORHCAyMpKSJUsSHh7udCiSDShZE5EcYd++feTPn58yZcpgjHE6HMnGrLUcPXqUffv2UbZsWafDkWxAp0FFJEc4f/48RYoUUaImGc4YQ5EiRTSKK+lGyZpkrrhzsPlz0GPOxAFK1CSzqK9JelKyJpnr52nwZX/4uDMc/9PpaEQc9eyzzzJmzBiv82fPns3WrVszMSIRCUZK1iRz1e8Pt78G+9bDhMaw+l1ISnQ6KpGgpGRNREDJmmS2kBBXwjZ4NZRuCguegA/awOHtTkcmkilGjRpFpUqVaNWqFdu3u/r9pEmTqF+/PrVq1aJbt26cO3eOVatW8c033/D4448TFRXFrl27PC4nItmfsbp2SPzg9ripB37//ferK8xa2DzTlbDFnYXmI6DpMAjVT90l/W3bto0qVaoA8NycX9n616l0Lb/qDQUY2aFaqsts2LCBvn37smbNGhISEqhTpw4DBw7kvvvuo0iRIgA888wzXHvttTz00EP07duX9u3b0717dwCOHj3qcTkJTu597iJjzAZrbT2HQpIsSiNr4pd0fdyUMVCrBwxZB5Xbw+IXYWIL+Ovnqy9bJAgtX76cLl26kCdPHgoUKEDHjh0B2LJlC9HR0dSoUYPp06fz66+/elzf1+VEJHvRfdbEefmKwR1ToEZ3+PYxmHQLNBkKLZ6C8NxORyfZUFojYBnJ068E+/bty+zZs6lVqxZTp05l6dKlHtf1dTkRyV40sibBo/LtMGQN1O4NK9+Ed5rCnyudjkok3cTExPDVV18RGxvL6dOnmTNnDgCnT5/m+uuvJz4+nunTp19aPn/+/Jw+ffrStLflRCR7U7ImwSV3Iej4FvT5GpISYGo712jb+fS9vkjECXXq1KFHjx5ERUXRrVs3oqOjAXjhhRdo2LAhrVu3pnLlypeW79mzJ6+++iq1a9dm165dXpcTkexNPzCQgNSrV8+uX78+Yz8k7iwsHgWrJ0CBEtD+dbjp1oz9TMm2PF3sLZKR9AMDSS8aWZPgFZEX2rwE/b6HXPngkzvgywFw9qjTkYmIiGQaJWsS/ErVhweXQfMnYMsX8HYD2PKlHlklIiI5gpI1yRrCckHLp11JW6FSMOs++Kw3nDrgdGQiIiIZSsmaZC3XVoN+C+HWF2HXIni7IWz4UKNsIiKSbSlZk6wnNAyaPASDVsF1NWDOw/wx9hZOHbjKJyqIiIgEISVrknUVKQ/3zmFemScpemor4e81YcNnL5KYkOB0ZCIiIulGyZr4xRjTwRgz8eTJk06H4hISQru+T7H/rqVszRVF3d9e5ffRTfjlp9VORyYiIpIulKyJX9L12aDpqHKlytR54js21BvDtYkHqPR1O+a+NYy/jgZJUikiIhIgJWuSbZiQEOq2f4DIR9azq1grbj86hTPjmvLZV7M5H5/odHgiAIwbN44qVarQu3dvmjRpAsCJEyeYMGHCpWVSTvvq2WefZcyYMekWa3o4dOgQderU4amnnqJr164kJSU5HZJIlqNkTbKd3NdcS5WhMznS4UOKhcdyx8a+fPnyfSz4eTd6Yoc4bcKECcybN4/p06ezatUqIP2StWC0bt06evXqxejRoylevDhHj+qm1iL+UrIm2VbRup25ZvhPHL6pJ3clfk3lr9rwwvj32HZAzxkVZwwcOJDdu3fTsWNHXn/9dfLlywfAk08+ya5du4iKiuLxxx+/Yhpg2rRpNGjQgKioKB588EESE12jxaNGjaJSpUq0atWK7du3e/zcTZs2ERMTQ9WqVQkJCcEYw8iRIwPahlmzZtGoUSNq1apFs2bNOHz4MABdunThmWeeITo6muuuu46FCxcCrmStVq1aAJw8eZJixYoF9LkiOZq1Vi+9/H7VrVvXZiUJO5faU/+uZu3IAnb6M53t87NW2aNnLjgdlmSirVu3Oh2Ctdba0qVL28OHD1trrc2bN6+11to//vjDVqtW7dIyKae3bt1q27dvb+Pi4qy11g4aNMh++OGHdv369bZ69er27Nmz9uTJk7Z8+fL21VdfvezzYmNjbaVKleyaNWustdY+88wzdvjw4TYpKemy5Zo1a2Zr1ap1xev777+/bLkjR45c+vvZZ5+148ePt9ZaW6FChUuf/cUXX9i+fftaa63t2bOnHTp0qB04cKBduHBhgLWWNXnqc8B6GwTHcL2y1ivM6WRRJDOElm9O/mFrOf/9i/Rc9w5/b97IyF8GULd1T+5uVJqwUA0yS/BatGgRGzZsoH79+gDExsZSvHhxjh07RpcuXciTJw8AHTt2vGLdhQsXUqdOHRo0aABAzZo1WbBgAcaYy5Zbvny5T7FMnTqVGTNmcOHCBQ4ePMhLL73EuXPnOHnyJI8++igACQkJFCpUCICwsDDeeuutgLZbRFyUrEnOEZGHyNtfgqhuFP5iMG8de5mvFyyj1+pBPNKxCc0qFnU6QhGPrLXce++9jB49+rL333jjjSuSrpS2bNlCjRo1Lk3/9NNP1KlT54rloqOjOX369BXvjxkzhlatWgHw0UcfsXbtWhYvXky+fPmIiYmhWrVq/Prrr9StW5fQ0FAANm/eTPXq1QH4+OOP/dtYEbmChhMk5ylRl1yDl2NbPEWHsHVMOj2EmVNeZ8CH6/jv0XNORyc5UP78+S9LlFJO33LLLcyaNYtDhw4BcOzYMfbs2UNMTAxfffUVsbGxnD59mjlz5lxRdpEiRdi8eTMAO3bs4Msvv6Rnz55XLLd8+XI2btx4xetiogbwyy+/0KRJE/Lly8cXX3zBqlWrqFGjBlu2bCEqKurScps3b6ZmzZpXXS8i4qJkTXKmsAhMiycJGbScAjdUZFzEeHrtHkHv17/k1e9+4+wFPQVBMk+RIkVo2rQp1atX5/HHH79iumrVqrz44ovceuut1KxZk9atW3PgwAHq1KlDjx49iIqKolu3bkRHR19Rdq9evThz5gzVq1dnwIABfPrppxQpUiSgOO+9917GjRtHdHQ0O3bsoFy5cuTNm5dffvnlsmRty5Ytl0bWROTqGWt1KwPxX7169ez69eudDiN9JCXCmnexi17gQpLh+Qu9WJynDU+0q0rnqBJpnmaSrGHbtm1UqVLF6TAkB/HU54wxG6y19RwKSbIojayJhIRC4yGYwT8SWboeL4VPZiLP88bM7+j2zio27zvhdIQiIpKDKVkTuahwWejzDXQYR42QPSzK/RQxRz6j69vLePzzTRw6fd7pCEVEJAfSr0FF3BkDde/FVGxN2Nx/MGz7R9xZeD0PbLyXm7cc5OFbKtC3SVkiwvR/TpY2/0k4+Ev6lnldDWj7cvqWKSKCRtZEPCtwA/T8BLpP4QZ7iG9z/R+jCs1hzLxfuO2NZSz+7W+nIxQRkRxCI2si3hgD1btCuRaYBU/SafPHtCr+IyPiB3D/1LO0qFSMf7avSvli+ZyOVPylETARyUI0siaSljyFoetEuOtz8hLL+NgnmFNxLlv/PMBtry/jvR92oV9Vi4hIRlGyJn4xxnQwxkw8efKk06FkvptuhcGrMfXup8be6fxY6BmGltnP6Pm/MXj6T5zRvdnEB7GxsTRv3vzSg9jTy/3330/x4sVTvb/ZiRMn6N69O5UrV6ZKlSr8+OOPqc47f/48DRo0oFatWlSrVi3gh79nhgULFlCpUiUqVKjAyy97HzktU6YMNWrUICoqinr16vm8fmJiIrVr16Z9+/YAxMXFERMTQ0KC9nvJeErWxC/W2jnW2gEFCxZ0OhRnRBaA9mOh7zxCQ8MZ9tdwFpT9nB+37qbT+BXsPHTl43pE3H3wwQd07dr10qOZ0kvfvn1ZsGBBqss88sgjtGnTht9++41NmzZddg8wT/Ny5crF4sWL2bRpExs3bmTBggWsXr06XeNOD4mJiQwZMoT58+ezdetWPv30U7Zu3ep1+SVLlrBx40Yu3ivSl/XffPPNy+orIiKCW265hRkzZmTMRom4UbImEogyTWHQSmj6CJUPfs3agk8TdXYlncavZN4vB5yOToLY9OnT6dSpEwCzZs2iUaNG1KpVi2bNmnH48OGAy42JiaFw4cJe5586dYply5bRr18/wJVsXHzYurd5xhjy5XNdkxkfH098fHyaN4nu2bMnPXr0oGHDhpQuXZq5c+cGvE2+Wrt2LRUqVKBcuXJERETQs2dPvv7663Rbf9++fcydO5f+/ftftl7nzp2ZPn16um2HiDdK1kQCFZ4bWj8P/RcRUaA4ryW9wnu5x/Ov6UsYPW8bCYlJTkcoQSYuLo7du3dTpkwZAFq2bMnq1avZtGkTrVu3ZubMmZctHx0dTVRU1BWvhQsX+v3Zu3fvplixYtx3333Url2b/v37c/bs2TTnJSYmEhUVRfHixWndujUNGzZM9XM2bdpEuXLlWLNmDdOnT+e5557LsG26aP/+/ZQqVerSdMmSJdm/f7/HZY0x3HrrrdStW5eJEyf6tP6wYcN45ZVXCAm5/CuzevXqrFu3LuC4RXylX4OKXK0SdWDAUljxBk2XvcKyvJt4esXd3LO3I2/1rkPRfLmcjlCCxJEjRy6NZgFMnTqVGTNmcOHCBQ4ePMhLL7102fLLly9Pt89OSEjgp59+4q233qJhw4Y88sgjvPzyy7zwwgupzgsNDWXjxo2cOHGCLl26pPrcz9jYWI4cOXLp2raqVaty/Pjxq9qmVq1acfDgwSveHzVq1KURSk8/8PE2Arhy5UpuuOEGDh06ROvWralcuXKq63/77bcUL16cunXrsnTp0suWCQ0NJSIigtOnT5M/f36/tkvEH0rWRNJDaDg0fxxTpQN5vnmIN/ZN4If9q+g/bhAj776N2jde43SEEgRy587N+fOuJ2F89NFHrF27lsWLF5MvXz5iYmKoVq3aZctHR0dz+vSV10GOGTOGVq1a+fXZJUuWpGTJkpdGxrp3737pQvrU5l1UqFAhWrRowYIFC7wma1u2bKFixYpERkYC8NNPP1GrVq2r2iZfRtxKlizJ3r17L03v27ePG264weOyF98vXrw4Xbp0Ye3atTRt2tTr+itXruSbb75h3rx5nD9/nlOnTnH33Xczbdo0AC5cuHBpe0UyjLVWL738ftWtW9eKF4kJ1v74jk184Vp7dmRx+69nhtmPV+22SUlJTkeWo23dutXpEKy11pYsWdLGxsba4cOH2zfeeMNaa+2sWbNsaGioPXPmzFWV/ccff9hq1ap5nd+sWTP722+/WWutHTlypB0+fHiq8w4dOmSPHz9urbX23LlztlmzZnbOnDnWWmtvvvlmu2/fvsvKnzRpki1RooSNjY21Z86csU2aNLErVqy4qm3yRXx8vC1btqzdvXu3vXDhgq1Zs6bdsmXLFcudOXPGnjp16tLfjRs3tvPnz/d5/SVLltjbb7/90vSRI0ds5cqVvcblqc8B620QHMP1ylovXbMmkt5CQqHRQEKGrCa8dEOeC/2AivN78vK0uZyPT9/bNUjWc+utt7JixQruvfdexo0bR3R0NDt27KBcuXLkzZs34HJ79epF48aN2b59OyVLlmTy5MkAtGvXjr/++guAt956i969e1OzZk02btzI008/fWl9T/MOHDhAy5YtqVmzJvXr16d169a0b9+epKQkdu7cecUPGjZt2kTv3r1p0aIF9evXZ9CgQTRt2jTgbfJVWFgY48eP57bbbqNKlSrceeedl41SXqyDv//+m2bNmlGrVi0aNGjA7bffTps2bdJc35slS5bQrl27jNw0EQCMtbqZp/ivXr169uLP3iUV1pL08zTi5j4FCef5JE9vWvd7gVJFCzgdWY6zbdu2y2694JSff/6ZsWPH8vHHHzsdSsC2bNnCBx98wNixYy97PyYmhkmTJlGpUiWHIstcXbt2ZfTo0V6311OfM8ZssNbW87iCiBcaWRPJSMYQUuceIoet51SpFtwfO5XT42NYt/oHpyMTh9SuXZuWLVum+01xM1P16tWvSNQAdu3aRcWKFR2IKPPFxcXRuXPnHJOYirOUrIlkhvzXUbzf5xxuO5HrzXGi5ndh7eRHSYqLdToyccD999+f7jfFDQb79++/4vYW2VVERAR9+vRxOgzJIXLGXiUSDIyhWMMeRD6yno2FWtFg7wccfLUBZ35f6XRkIiISxJSsiWSy3IWKUW/YDBbWmYCNO0ee6bdzbNajcOGM06GJiEgQUrIm4gBjDK069ubQPUuYFdKGwls+4OybDWDXYqdDy9b0gyrJLOprkp6UrIk4qHaFG2nx6FT+WWQMB88kwcddSPxqMMQeT3tl8UtkZCRHjx7Vl6hkOGstR48e1c1yJd3o1h0SEN26I33FJyYxZu4mCqwdy8CwbyFvUUJvfw2qdnQ6tGwjPj6effv2XXqCgEhGioyMpGTJkoSHh1/2vm7dIYFQsiYBUbKWMeZs+oupX8zmpZCJVOIPqNIR2o2B/Nc6HZqIpAMlaxIInQYVCSIdat3A6CH3MDTvGF5N6EHC9gXYtxvAxk9A/1iJiORIStZEgsxN1+bni4ea8/tNA7gtdhS7KQmzB8G0bnDiv06HJyIimUzJmkgQKhAZznv31KX7bbdw68knGR85gKQ9P8LbjWDNREhKcjpEERHJJErWRIKUMYZBLcrz4f2NmRzXmtviX+VI4dow/3GY0haO/O50iCIikgmUrIkEuWYVi/Ltw9HkLlaGensGMbf8v7CHf4N3msLy1yAx3ukQRUQkAylZE8kCShTKzcwHG9OrwY0M+bUyQ695l7jyt8Ki52FSSziwyekQRUQkgyhZE78YYzoYYyaePHnS6VBynMjwUEZ3rcm/u9Xg+73Q9I++LKs9Fnv6b5jYEhY+C/G6h5iISHaj+6xJQHSfNWdt2X+Skd/8yoY9x6lZxPJ20S8otedLKFIBOo6H0o2dDlFEPNB91iQQGlkTyYKqlyjIrIGNmdSnHrGh+Yne3p1/FXjRdXf+KW1g7nC4cNrpMEVEJB0oWRPJoowxtK56LQuGxfBK95p8f6EqdY4+z3f5OmPXvQ8TGsPOhU6HKSIiV0nJmkgWFxpiuLNeKZYMb8GwdlGMONubbhdGcjDWuG6k+9VAOHfM6TBFRCRAStZEsonI8FAGxJRn2YiWNGrelltjX2RCYmcSN80kaXx9+HW20yGKiEgAlKyJZDMFc4czok1lvn/8NvbWHk6n+BfZdjY/fH4vCZ/2htMHnQ5RRET8oF+DSkD0a9CsY9fhM7z+3VZKbPuAR8O/wITlIqTNS4TXvQeMcTo8kRxFvwaVQGhkTSSbK18sH+PvbkDbgS/zVPF32BhXgvBvH+LQhLYkHf3D6fBERCQNStZEcoioUoUYO6gb53t/w/g8g8lz6Gfi3mrI79+8ik1McDo8ERHxQqdBJSA6DZq1JSVZFq5eT/6FI2ic9BM7wquQ2PEtqtSo73RoItmaToNKIDSyJpIDhYQYbm1Sn7pPLWRZ9VEUj99LuVlt+Hrco+w6eNzp8ERExI2SNZEcLCI8lJjuQwl/eD1/FmtJp2MfEDchhremfc7Bk3rOqIhIMFCyJiLkLXw9lYbO4lTnDymZ6xyDfh/AnNf689q8jZyMjXc6PBGRHE3JmohcUiCqM/kf28D5aj15IGQOXVf3YNi/32bisl2cj090OjwRkRxJyZqIXC53IfLd+Q70+ZoSBcOZwkhy/2cEt786j5nr9pKQmOR0hCIiOYqSNRHxrFwLIh5aA42GcHfYIj5LeJR5X31ImzeXs3LnEaejExHJMZSsiYh3EXmhzUuYft9TtEgRpka8yoizrzH0/e95+qtfOH1e17OJiGQ0JWsikrZS9TEPLoPmT9LarmRFvqc4tX4GbV5fxvLfDzsdnYhItqZkTUR8E5YLWj6FeXAZeYuXYXz4W7yS+DL/mPwdT325WaNsIiIZRMmaiPjn2mrQbyHc+iJN2MyyvE9gN3zEbWN/4IcdGmUTEUlvStZExH+hYdDkIcygVUSWjOLl8EmMT3yOf06ZwxOzNnNKo2wiIulGyZqIBK5Iebh3DrR/g9phf7Ao95MU2PgubccuYcn2Q05HJyKSLShZE5GrExIC9e7DDF5DeIWW/F/YdN5P+D9GT/2Sxz/fpCcgiIhcJSVrIpI+CpaAXp9Bt8lUjjzK/Mj/o8SmN2k3diGLf/vb6ehERLIsJWsikn6MgRrdMUPWElq9C8PCvmBa4hO8+eFn/GPmJk6e0yibiIi/lKyJSPrLWxS6vQ+9PqNM3ji+yvUsVX75N+3H/oeFWzXKJiLiDyVrIpJxKrXFDFlDSN176R86lxlJjzF52oc8NmMjJ87FOR2diEiWoGRNRDJWZEHo8Abc+y3XF8zNpxGjqL/leTqPnc/3GmUTEUmTkjURyRxlozGDVkGTh+gZtoRZicP4bNp7DPvsZ46f1SibiIg3StZEJPNE5IFbX8T0X0jhYtcxOeI1bvn1Ke4YO4fvfj3odHQiIkFJyZqIZL4SdQkZ8AO0eJr24ev5ImkYc6eP4+FPfuKYRtlERC6jZE1EnBEWAS2ewAxcTv4bbmJcxNt0+e0x7h77JfN/OeB0dCIiQUPJmog4q3gVQvr9B24bTfOI35iVNIzln73K0OnrOXrmgtPRiYg4TsmaiDgvJBQaDyZk8I9ElmnAS+GTuWf7UPqNncHczRplE5GcTcmaiASPwmUJ6fM1dHyLepH7mZH0DzbNeI6h09ZyRKNsIpJDKVkTAIwx5Ywxk40xs5yORXI4Y6BOH0KHriX8plY8Hf4pD+54kCFjP+LbzX85HZ2ISKZTspYNGGM+MMYcMsZsSfF+G2PMdmPMTmPMk6mVYa3dba3tl7GRivihwPWE9PoEuk+hSp6TTE96gl0zn+aZWes5H5/odHQiIplGyVr2MBVo4/6GMSYUeBtoC1QFehljqhpjahhjvk3xKp75IYv4wBio3pWwh9cTUqM7j4R9RZ/NfXjyzcnsPnzG6ehERDKFkrVswFq7DDiW4u0GwM7kEbM44DOgk7X2F2tt+xSvQ758jjFmgDFmvTFm/eHDh9N5K0RSkacwId0mQu9Z3JgvibFnRrBi/APM3bDT6chERDKckrXsqwSw1216X/J7Hhljihhj3gVqG2Oe8rSMtXaitbaetbZesWLF0jdaEV9UbE3kw2uJrXUvfcw8anzdhskfTdFpURHJ1pSsZV/Gw3vW28LW2qPW2oHW2vLW2tEZGJfI1YksQN4ub5LQ51vy5Ymk3+5hLH21J3v268cHIpI9KVkLEsaYkelc5D6glNt0SUDfZpJthJWLpvA/1rGnygO0iltI7olNWLfgI6fDEhFJd0rWgsdIY8y/jTGTjDGDjDHXXGV564CKxpiyxpgIoCfwzdWHKRJEwnNTuscYjvWaz9nwa6i/+iG2vNGF88d1I10RyT6UrAUPC5wHvsM1IrbKGFPLlxWNMZ8CPwKVjDH7jDH9rLUJwNDk8rYBM621v2ZM6CLOKl6pESVHrGZpyYFUPL6MuHH1ObLyQ7Bez/yLiGQZxupgFhSMMb9aa6u5Td8EvGutvdnBsK5gjOkAdKhQocIDv//+u9PhiFxh1eqV5FkwjCh2cOjaGIr3mgCFSqW9okgmMMZssNbWczoOyVo0shY8jhhj6l6csNbuAILuJ5fW2jnW2gEFCxZ0OhQRj5o0akrRhxYzMd9A8h5cw4VxDYhfPRGSkpwOTUQkIErWgsfDwDRjzDRjzBPGmOnAH04HJZIVlSySn77DRjO55qesiS9P+ILHOf9+Gzii+7KJSNajZC1IWGs3AVHAp8lvLQF6ORaQSBYXERbCw91uIbbH5zzDYOL+2kLihMaw4nVITHA6PBERn+maNQlIvXr17Pr1650OQ8Qne4+d45lpC+l5eBxtQ9eRdF0tQjqNh+trOh2a5DC6Zk0CoZE1Ecn2ShXOw8TBt7Om/psMjBvGyb//i53YAhY9D/HnnQ5PRCRVStZEJEfIFRbKsx2r0anXQNonvcY3thksfw3ei4b/rnE6PBERr5SsiV+MMR2MMRNPnjzpdCgiAWlb43o+ebgNk4o8Tp+4Jzh56hT2g9tg3gi4cMbp8ERErqBkTfyiW3dIdlC6SF6+GNSEMg070uTUKOblbo9dOxEmNIadi5wOT0TkMkrWRCRHyhUWyvOdqvPvu5rwxLl7uI/nOJsUCtO6wuzBcO6Y0yGKiABK1kQkh2tf8wbmPNSMQ9fUoc7hkay8/l7sps9gQiPYqsfpiojzlKyJSI5XtmhevhzchDsaVaD3H7cx/Jo3ictdHGbeAzPugdN/Ox2iiORgStZERIDI8FBe7FyDcb1qs+BIMZoeeZpdtYbDju/g7Qaw8RM9GF5EHKFkTUTETcdartOiRQrm45Y1dZhU/WNsscowe5Drerbje5wOUURyGCVr4hfdukNygnLF8jF7SFN6NbiRUWsS6BH3T060fAn2rnX9YnTNe3owvIhkGj1uSgKix01JTjH75/08/dUvhBjDCy0K0HnfGMyuhVCqEXR8C4rd5HSIkoXocVMSCI2siYikonPtEix4JIaaJQvy6HfH6HH2Hxxq9SYc2Q7vNoVlYyAx3ukwRSQbU7ImIpKGG4vkYXr/hrzSrSbbDp6m2YJrmVJ7Bkk3tYPFL8CklnBgk9Nhikg2pWRNRMQHxhjurF+KRY815+ZKxXlu8RE6/N2fPa0mwplDMLElLHxWD4YXkXSnZE1ExA/FC0Ty7j11ead3Hf4+dYGb5+XnjcrTSKjZE1a87jo1uudHp8MUkWxEyZqISADa1rieRY81p1udEryx4jCtd93JtlYfQWIcTGkDc4fDhdNOhyki2YCSNRGRABXME84r3WsxrV9DEpKSaPttGM+VfJ+4eg/Cuvddt/n4faHTYYpIFqdkTUTkKjWrWJTvhsXQv1lZPtxwhJjNt7Hu5s8gPA9M7wZfDdSD4UUkYErWxC+6Ka6IZ3kiwnimfVW+HNyUgrnDuWNeIo9e8xbnGj0Gv3zuemTVr7P1yCoR8ZtuiisB0U1xRbyLS0jinaW7GL/kd/LlCmNMTAg3b38Bc2AjVG4Pt78G+a9zOkxxgG6KK4HQyJqISDqLCAvhkVYVmftwNGWK5qXfggv0C3+ZE83+CTsXukbZfp6mUTYR8YmSNRGRDHLTtfmZNbAJIztUZfWfJ2myrDpfNpyJLV4Nvh4CH3eG4386HaaIBDklayIiGSg0xHBf07J8NyyGuqWv4bFFZ7jj/P9xKGY07Nvg+sXo6ncgKdHpUEUkSClZExHJBKUK5+Gj+xvw2h21+P3wOZotKsuUqE9JKt0UFjwJH7SBQ785HaaIBCElayIimcQYQ7e6JVn4WHNaV7uW55adot3hh/hvizfg6E54Lxp+eFUPhheRyyhZExHJZMXy5+Ltu+owqU89jsfG0+K74rxe6WMSKrWHJS/CxBbw189OhykiQULJmoiIQ1pXvZbvH2tOj/o38ubqk9z8Zx+2NX8Pzh2FSTfD9/+C+FinwxQRhylZExFxUIHIcEZ3rcGnDzQixEDb7/IzstQHxNXsDSvfhHeawp8rnQ5TRBykZE1EJAg0Ll+EBcNieLB5OaZtPEGzrZ1ZGzMVbCJMbQffPgbnTzkdpog4QMma+EWPmxLJOJHhoTzVtgqzBzelSL5c3PmfCB4qNIHjtQbAhikwoRHs+I/TYYpIJtPjpiQgetyUSMaKT0xi4rLdjF+8k/MJiQwsf5xHzr5J5PEdUONOaPMy5C3idJjiJz1uSgKhZE0ComRNJHMcPXOBKSv/5MNVf3L+wnlevfZ7Op7+jJDIgtDuFajWFYxxOkzxkZI1CYSSNQmIkjWRzHUyNp6Pf/yTySv+oHjsLt7O9wEV4ndgK7XF3P46FLje6RDFB0rWJBBK1iQgStZEnHH2QgKfrv0v7//wO7fHfs2I8FmEhEcQ1mYUpk4fjbIFOSVrEgj9wEBEJAvJmyuM/tHlWPpEK8p2eIJ7Il5n3flSmDkPc2RCGxKP7HY6RBFJZ0rWRESyoMjwUO5uVJrpT/Tir04zeS3XECIObSJ+fEM2z3yR+Hg9skoku9BpUAmIToOKBJfEJMuSdRvJ9/0IGiWsZaupyO6mL9O6RUtyhYU6HZ4k02lQCYSSNQmIkjWR4GSTktjynyncuOZZciedZWpYd8Ji/kHPxuXJExHmdHg5npI1CYROg4qIZCMmJIQabfpRYPjPnCjbngGJM2iyqBsDXp7E20t2cuq8To+KZDVK1kREsiGTtyjF+34Ed82kXL4EPkp6mvBF/+KWl+fz2n+2c+xsnNMhioiPlKyJiGRnN91GxMNrCanXlwFhc5kXNoL1S7+m2b8XM2ruVg6dOu90hCKSBl2zJn4xxnQAOlSoUOGB33//3elwRMQff66Abx6CY7tZVagDg/7uRGxoPnrUK8WDzctR8po8TkeY7emaNQmEkjUJiH5gIJJFxZ2DpaPhx/Ek5CnO9KLDeHFnaayFLrVLMKhFecoVy+d0lNmWkjUJhJI1CYiSNZEsbv9P8PVQOPQr5yp1YXxEfyb/fJr4xCTa1bieIS0rUOX6Ak5Hme0oWZNA6Jo1EZGcqEQdGLAUWj5Dnp1zGbHzHtZ1PM6A6HIs+e0Qbd9czsivtzgdpYigZE1EJOcKi4Dmj8ODy6FIeQrMG8STJ55l1ZDKDGtVkdo3XuN0hCKCkjURESleGe7/Dtq8DH8so+DkZgwrtILOta53OjIRQcmaiIgAhIRCo0EwaJXrFOm3j8KHHeDoLqcjE8nxlKyJiMj/FC4Lfb6GjuPh4C/wThNYOQ4SE5yOTCTHUrImIiKXMwbq3AND1kCFVvD9P2FyKzioHxyIOEHJmoiIeFbgeugxDe6YCif3wcTmsHgUJFxwOjKRHEXJmoiIeGcMVOsCQ9ZCjTtg2SvwXgzsXed0ZCI5hpI1ERFJW57C0OVd6D0LLpyBya1hwVMQd9bpyESyPSVrIiLiu4qtYchqqN8fVk+ACY1h91KnoxLJ1pSsiYiIf3Llh9vHwH3zITQcPurkenRV7AmnIxPJlpSsiYhIYEo3gYErodmjsPETeLshbPvW6ahEsh0layIiErjwSGj1LDywGPIVgxm9Yea9cOaQ05GJZBtK1kRE5OrdEAUPLIGb/wnb58HbDWDTZ2Ct05GJZHlK1sQvxpgOxpiJJ0+edDoUEQk2oeEQM9x1arToTfDVgzD9Djix1+nIRLI0JWviF2vtHGvtgIIFCzodiogEq2I3wX0LoO0rsGcVTGgEaydBUpLTkYlkSUrWREQk/YWEQMMHYfCPULI+zBsOU2+HIzudjkwky1GyJiIiGeea0nDPV9BpAhz61fVg+BWv68HwIn5QsiYiIhnLGKjdG4asg5tug4XPwvs3w4HNTkcmkiUoWRMRkcyR/1ro8THc+RGcOgATW8Ci5yH+vNORiQQ1JWsiIpK5qnaCIWugVk9Y/hq8Fw3/XeN0VCJBS8maiIhkvjyFofMEuPtL18jaB7fBvBGuh8SLyGWUrImIiHMq3OL6xWiDAbB2ouvB8DsXOR2VSFBRsiYiIs7KlQ/avQL3L4CwXDCtK8weDOeOOR2ZSFBQsiYiIsHhxkYwcAVE/8P1qKq3G8LWr52OSsRxStZERCR4hEfCLf+CAUsh/3Uwsw/MuAdO/+10ZCKOUbImIiLB5/qa8MBiaPUs7PjO9WD4n6frwfCSIylZExGR4BQaDs0ehUEroXgV+Hqw63q243ucjkwkUylZExGR4Fa0IvSdB+3GwN61rl+MrnlPD4aXHEPJmoiIBL+QEGjwgOs2H6Ubw/wRMKUNHN7udGQiGU7JmoiIZB2FboTes6DLe3BkB7zbDJaNgcR4pyMTyTBK1kREJGsxxvWoqiFrofLtsPgFmNQS/trodGQiGULJmoiIZE35isMdU6HHdDhzCCbdDAufhfhYpyMTSVdK1kREJGur0t71YPiou2DF665To3tWOR2VSLpRsiYiIllf7mug03i4ZzYkxsGUtjD3H3DhtNORiVw1JWsiIpJ9lG8Jg1dDo8GwbjK83Qh+/97pqESuipI1ERHJXiLyQpvR0O8/rofET+8OXz6oB8NLlqVkTfxijOlgjJl48uRJp0MREUldqQbw4DKIGQFbZsH4+rDlSz2ySrIcJWviF2vtHGvtgIIFCzodiohI2sJywc3/BwN+gIIlYdZ9MONuOHXA6chEfKZkTUREsr/rqkP/RdD6edi5EN5uCD99pFE2yRKUrImISM4QGgZNH4FBq+C6GvDNQ/BRJzj2h9ORiaRKyZqIiOQsRcrDvXOg/euw/yd4pwn8OAGSEp2OTMQjJWsiIpLzhIRAvftdN9MtEw3fPQWTb4VD25yOTOQKStZERCTnKlgC7poBXd+HY7vh3Wj44RVIiHM6MpFLlKyJiEjOZgzUvAOGroOqHWHJKJjYAvZvcDoyEUDJmoiIiEveotD9A+j5KcQeg/dbwX/+CXHnnI5McjglayIiIu4qt3Ndy1anD6waB+82hT+WOx2V5GBK1kRERFKKLAgd3nT9atRa+LA9zBkG5/X0Fsl8StZERES8KRvjui9b46Hw04d6MLw4QsmaiIhIaiLywG2joN9CyF0ITu13OiLJYcKcDkBERCRLKFnX9YzR0HCnI5EcRsmaiIiIr8IinI5AciCdBhUREREJYkrWRERERIKYkjURERGRIKZkTURERCSIKVkTERERCWJK1kRERESCmJI1ERERkSBmrLVOxyBZkDHmMLDHw6yCQMqH56V8ryhwJINCS4un+DKjHF+XT2u51OZ7m+dLm4Bz7eJUm/izTnq3i69tpX0l8OWCdV8pba0tFuC6klNZa/XSK91ewMS03gPWB1N8mVGOr8untVxq873N86VNnGwXp9rEyXbxta20r2Rem/jTVk62i14586XToJLe5vj4nlPSKxZ/y/F1+bSWS22+t3lqk6tfJ73bxZ+2cor2Fd8+RyTD6TSoZDpjzHprbT2n45DLqV2Cj9okOKldJLNpZE2cMNHpAMQjtUvwUZsEJ7WLZCqNrImIiIgEMY2siYiIiAQxJWsiIiIiQUzJmoiIiEgQU7ImjjPGlDPGTDbGzHI6FvkfY0xnY8wkY8zXxphbnY5HwBhTxRjzrjFmljFmkNPxiIsxJq8xZoMxpr3TsUj2pGRNMoQx5gNjzCFjzJYU77cxxmw3xuw0xjwJYK3dba3t50ykOYuf7TLbWvsA0Bfo4UC4OYKfbbLNWjsQuBPQrSMyiD9tkuwJYGbmRik5iZI1yShTgTbubxhjQoG3gbZAVaCXMaZq5oeWo03F/3Z5Jnm+ZIyp+NEmxpiOwApgUeaGmaNMxcc2Mca0ArYCf2d2kJJzKFmTDGGtXQYcS/F2A2Bn8khaHPAZ0CnTg8vB/GkX4/JvYL619qfMjjWn8HdfsdZ+Y61tAvTO3EhzDj/bpCXQCLgLeMAYo+9VSXdhTgcgOUoJYK/b9D6goTGmCDAKqG2MecpaO9qR6HIuj+0CPAS0AgoaYypYa991Irgcytu+0gLoCuQC5mV+WDmaxzax1g4FMMb0BY5Ya5MciE2yOSVrkpmMh/estfYoMDCzg5FLvLXLOGBcZgcjgPc2WQoszdxQJJnHNrn0h7VTMy8UyWk0XCuZaR9Qym26JPCXQ7HI/6hdgo/aJPioTcQxStYkM60DKhpjyhpjIoCewDcOxyRql2CkNgk+ahNxjJI1yRDGmE+BH4FKxph9xph+1toEYCjwHbANmGmt/dXJOHMatUvwUZsEH7WJBBs9yF1EREQkiGlkTURERCSIKVkTERERCWJK1kRERESCmJI1ERERkSCmZE1EREQkiClZExEREQliStZEREREgpiSNREREZEgpmRNRLINY8xXxpgXjTHLjTEHjTGtnI5JRORqKVkTkeykOnDCWhsNDAZ6OxyPiMhVU7ImItmCMSYPUBB4PfmtMOCEYwGJiKQTJWsikl1UAzZYaxOTp2sCWxyMR0QkXShZE5Hsojqw0W26JrDZmVBERNKPkjURyS5qcHmyVh2NrIlINmCstU7HICIiIiJeaGRNREREJIgpWRMREREJYkrWRERERIKYkjURERGRIKZkTURERCSIKVkTERERCWJK1kRERESCmJI1ERERkSD2/5/Iks43N8WrAAAAAElFTkSuQmCC\n", + "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", "text/plain": [ "
" ] @@ -380,8 +380,6 @@ } ], "source": [ - "from matplotlib import pyplot as plt\n", - "\n", "a_fit, p_fit = fit_power_law(stddev_data)\n", "\n", "plt.figure()\n", @@ -423,7 +421,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 8, "metadata": { "deletable": false, "nbgrader": { @@ -441,14 +439,22 @@ "source": [ "def number_of_hits_in_d_dimensions(N,d):\n", " \"\"\"Simulates number of hits in case of N trials in the d-dimensional direct-sampling pebble game.\"\"\"\n", - " # YOUR CODE HERE\n", - " raise NotImplementedError()\n", + "\n", + " def random_in_d_dimensional_box():\n", + " \"\"\"Returns a random position in the box [-1,1)^d.\"\"\"\n", + " return rng.uniform(-1,1,d)\n", + " \n", + " number_hits = 0\n", + " for i in range(N):\n", + " position = random_in_d_dimensional_box()\n", + " if is_in_circle(position):\n", + " number_hits += 1\n", " return number_hits" ] }, { "cell_type": "code", - "execution_count": null, + "execution_count": 9, "metadata": { "deletable": false, "editable": false,