From 066d393b325250aeb1e54cd716dcd67b83407e3f Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Ga=C3=A9tan=20Facchinetti?= Date: Fri, 29 Mar 2024 18:22:28 +0100 Subject: [PATCH] cleaning --- example.ipynb | 74 ++------------------------------------------------- 1 file changed, 2 insertions(+), 72 deletions(-) diff --git a/example.ipynb b/example.ipynb index 4ce8fbc..0a829d0 100644 --- a/example.ipynb +++ b/example.ipynb @@ -38,27 +38,6 @@ "z_pbh, xe_pbh, Tm_pbh = pyhy.call_run_hyrec(cosmo(), inj_pbh())" ] }, - { - "cell_type": "code", - "execution_count": 4, - "metadata": {}, - "outputs": [ - { - "data": { - "text/plain": [ - "array([9.39846976e-01, 1.74273336e+00, 2.71873230e+00, ...,\n", - " 2.18012745e+04, 2.18040000e+04, 2.18067255e+04])" - ] - }, - "execution_count": 4, - "metadata": {}, - "output_type": "execute_result" - } - ], - "source": [ - "Tm_pbh" - ] - }, { "cell_type": "code", "execution_count": 5, @@ -75,51 +54,9 @@ }, { "cell_type": "code", - "execution_count": 6, - "metadata": {}, - "outputs": [ - { - "data": { - "text/plain": [ - "array([1.000e+00, 2.000e+00, 3.000e+00, ..., 7.998e+03, 7.999e+03,\n", - " 8.000e+03])" - ] - }, - "execution_count": 6, - "metadata": {}, - "output_type": "execute_result" - } - ], - "source": [ - "z_c" - ] - }, - { - "cell_type": "code", - "execution_count": 7, + "execution_count": null, "metadata": {}, - "outputs": [ - { - "data": { - "text/plain": [ - "" - ] - }, - "execution_count": 7, - "metadata": {}, - "output_type": "execute_result" - }, - { - "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAkMAAAG1CAYAAAD6GvACAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8g+/7EAAAACXBIWXMAAA9hAAAPYQGoP6dpAABa0klEQVR4nO3deXxU5f3+/9fMZN8J2RcI+yJLEAiCC8FiAS0qxUKtC2jVLujXCqhYrUrbn/ipS6k1ra2KuLUureKOCyWgKLIGRDaBAAGyEBKykm3m/P44ZCBCQgJJTjJzPXmcx2TOnJnznpwkc3Gf+9y3zTAMAxEREREvZbe6ABERERErKQyJiIiIV1MYEhEREa+mMCQiIiJeTWFIREREvJrCkIiIiHg1hSERERHxagpDIiIi4tV8rC6go3O5XBw6dIjQ0FBsNpvV5YiIiEgzGIZBWVkZCQkJ2O1Nt/0oDJ3BoUOHSE5OtroMEREROQs5OTkkJSU1uY3C0BmEhoYC5jczLCzM4mpERESkOUpLS0lOTnZ/jjdFYegM6k+NhYWFKQyJiIh0Ms3p4qIO1CIiIuLVFIZERETEq+k0WStwuVzU1NRYXYaIpXx9fXE4HFaXISLSYgpD56impobs7GxcLpfVpYhYLiIigri4OA1DISKdisLQOTAMg9zcXBwOB8nJyWccx0DEUxmGQWVlJQUFBQDEx8dbXJGISPMpDJ2Duro6KisrSUhIICgoyOpyRCwVGBgIQEFBATExMTplJiKdhpoyzoHT6QTAz8/P4kpEOob6/xTU1tZaXImISPMpDLUC9Y8QMel3QUQ6I4UhERER8WoKQyIiIuLVFIakSYsXLyYiIuKM29lsNpYsWdLm9bSWzMxMbDYbR48ebXSb5r53ERHp3BSGpEnTp09n586d7vsPP/wwqampltXT2UKXiIh0fLq0XpoUGBjovmRaRESkMdnr9/G3j77mQF0ZUz/MIjGvAEfNMZal9eaDsedhnHR9hWEDA3PF43/6L4NWv0p49whrCsdLWobef/99+vXrR58+fXjuuefabD+GARUV1iyG0fzvRUREhHtYgKysLGw2G/PmzXNvc8stt3D99dcDDU8VLV68mPnz57Np0yZsNhs2m43Fixe7n1dYWMiUKVMICgqiT58+vPvuuw32vWLFCtLS0vD39yc+Pp558+ZRV1fnfjwlJYWFCxc2eE5qaioPP/yw+3GAKVOmYLPZ3Pe/b+/evdhsNl577TXGjBlDQEAAgwYNYsWKFadsu2rVKoYMGUJAQAAXXHABW7ZsOWWbjz/+mAEDBhASEsLEiRPJzc097X5FRLyN4TLY/JflrEn+MdFj+rNosD+vpfeie+mnjM55g7T897AH57JqeC++PP/E8tWwXqwe1pPVw3oyrHgZtZXWDsfh8S1DdXV1zJ49m+XLlxMeHs7w4cOZMmUKXbt2bfV9VVZCSEirv2yzlJdDcPCZt7v44ospKytj48aNjBgxghUrVhAVFUVmZqZ7mxUrVnDvvfee8tzp06ezZcsWli5dymeffQZAeHi4+/H58+fzpz/9iccee4y//vWvXHfddezbt4/IyEgOHjzI5ZdfzsyZM3nppZfYvn07t956KwEBAe6wcyZr164lJiaGF154gYkTJ55xUL+7776bhQsXMnDgQJ588kkmT55MdnZ2g2N/991385e//IW4uDh++9vfMnnyZHbu3Imvry8AlZWVPP7447z88svY7Xauv/565s6dy6uvvtqsmkVEOiTDAJcLw2XgchrmbZ2TGqcLl8vA5XThdBo4XSfuu1wGhmFQW+di/aYclm/dzi7fIpb+7i7sx/9Hftu/VvLdwAHsGXwrFYP8sAcFEBvlz7xVBWCzYf4Dw2Zzt8asuelZLogPtexbAV4QhtasWcN5551HYmIiAJMmTeKTTz7h2muvtbgya4SHh5OamkpmZiYjRowgMzOTu+66i/nz51NeXk5JSQm7du1i7Nixpzw3MDCQkJAQfHx8iIuLO+XxmTNnur+vjzzyCE899RRr1qxh4sSJ/O1vfyM5OZmnn34am81G//79OXToEPfeey8PPvhgs6YyiY6OBk7Mf3Umt99+O1OnTgXg73//O0uXLuX555/nnnvucW/z0EMPcdlllwHw4osvkpSUxNtvv820adMAc/DAZ555hl69erlf8/e///0Z9y0iciaVZVVs3JhD9r7DHCos5nB5GUV1VVz/7iZ8KkrxqSzlXxMHsnJET254bzk//Xg5DsPJ9p6JTP+/+3DZbbhsdlx2G4atfjGDhst20jrMdd/OvIk+hw4CcN8tt/Dodddx15tv8uTf/gZATnQ03d54o3nFxwAxyUAyHw4bTcixIcT+/nYWXHPeKZte0jrfrjbV4cPQypUreeyxx1i/fj25ubm8/fbbXH311Q22ycjI4LHHHiMvL4+hQ4fy17/+lbS0NAAOHTrkDkIAiYmJHDx4sE1qDQoyW2is0JLZQMaOHUtmZiZz5szh888/Z8GCBbzxxht88cUXFBUVkZCQQJ8+fVpcw5AhQ9xfBwcHExYW5p6ratu2bYwePbrBoHwXXngh5eXlHDhwgG7durV4f2cyevRo99c+Pj6MGDGCbdu2NbpNZGQk/fr1a7BNUFCQOwiBOedW/XsSEWkOZ52TzE+/ZfmGHXzrLOFQuIPcmBAOxUTidDggGUgOBszm/cefvIMuxz9MXgy7i839u1P5lUGScz8AeQRyODKi5YWcZkzUk3tY2Jrb3wLoerSEXvuLuLAkkMFvvkP3nlEtr6cD6fBhqKKigqFDh3LzzTfz4x//+JTHX3/9dWbPns0zzzzDqFGjWLhwIRMmTGDHjh3ExMS0a602W/NOVVktPT2dRYsWsWnTJnx9fenfvz/p6elkZmZSXFx82lah5qg/tVTPZrPhcrma/Xy73Y7xvV9Gq6d1ON17+n6NIiL1ao/VUrgln5x3N7L92y08c3EM23rFcDQsFC6MBqIbbO9fU0PMkRK6lBwjrLyW8CqDlWm/IcInEMLDGenflR4bK+g3cibfps/CJ8AHh8PG+yVgdziw2+3Y7DYcDjt2uw27zYbDbsPmsGO32cx1djvYbAS8s4Y8Hwc2u41bDBfX2yH4joc5ctcfsdltBNhsZOM68Xr1r+mwub/2OX4bEtTh40OLdPh3M2nSJCZNmtTo408++SS33norN910EwDPPPMMH3zwAYsWLWLevHkkJCQ0aAk6ePCgu9XodKqrq6murnbfLy0tbYV30bHU9xv685//7A4+6enpPProoxQXFzNnzpxGn+vn5+fufN0SAwYM4L///S+GYbhbh1atWkVoaChJSUmAeRrs5M7JpaWlZGdnN3gdX1/fZu9/9erVXHKJ2UBbV1fH+vXruf3220/Zpr5Vqri4mJ07dzJgwIAWvz8R8R51tU6e/8cy+j37HqE52wirzOeVK0fzlxuu5PpPP+Xpvz5FPJDUtSsz/t9/AAiorqb/7jz6HHbS2yeUPvGxDEvtwaAhifj4fK+bwAMnvmzsv6ZDGlnfXLGNrO/c7Ttnr8OHoabU1NSwfv167rvvPvc6u93O+PHj+eqrrwBIS0tjy5YtHDx4kPDwcD766CN+97vfNfqaCxYsYP78+W1eu5W6dOnCkCFDePXVV3n66acBuOSSS5g2bRq1tbVNtgylpKSQnZ1NVlYWSUlJhIaG4u/vf8Z9/vrXv2bhwoXccccd3H777ezYsYOHHnqI2bNnu/sLXXrppSxevJjJkycTERHBgw8+eEon6ZSUFJYtW8aFF16Iv78/Xbp0aXSfGRkZ9OnThwEDBvDnP/+Z4uJibr755gbb/P73v6dr167ExsZy//33ExUVdcppWBGReqtX7eaXm79g06DurK5dxvBi87R6OAMoCQ3hSHgYTuxk+/cnNyKNezJzGXFeH664fAxBEzSpd0fVqcNQYWEhTqeT2NiGGTc2Npbt27cDZl+RJ554gnHjxuFyubjnnnuavJLsvvvuY/bs2e77paWlJCcnt80bsNDYsWPJysoiPT0dMPvLDBw4kPz8fPr169fo86ZOncpbb73FuHHjOHr0KC+88AIzZ8484/4SExP58MMPufvuuxk6dCiRkZH8/Oc/54EHTvwX6L777iM7O5sf/ehHhIeH84c//OGUlqEnnniC2bNn8+yzz5KYmMjevXsb3eejjz7Ko48+SlZWFr179+bdd98lKirqlG3uvPNOvvvuO1JTU3nvvffw89MfLBE5Vda6/VxzYCMHB3QnpLKSj4ZfT+2YBAJ7JzIiLoRPywNJ+ckvMR54mN7+DnoDF1tdtDSLzehEHSBsNluDDtT1naO//PLLBh1h77nnHlasWMHXX399zvssLS0lPDyckpISwsLCGjxWVVVFdnY2PXr0ICAg4Jz3Ja1j79699OjRg40bN1o6WrY30u+EeKq6WicX/W0xXw/tRY+cPF7rlUbamJ5WlyVNaOrz+/s69aCLUVFROBwO8vPzG6zPz89v1qXXIiIizfH/PfkuXw/tRdCxKv4ZPVhByMN06jDk5+fH8OHDWbZsmXudy+Vi2bJlDVqKzkZGRgYDBw5k5MiR51qmiIh0YlWVNbyYaF5Yc+3X+YyfeOpYOtK5dfg+Q+Xl5ezatct9v77zbmRkJN26dWP27NnMmDGDESNGkJaWxsKFC6moqHBfXXa2Zs2axaxZs9zNbNJ5pKSk6PJ3EWk1T2R8RPbIOLoeLeWPvzx1iBfp/Dp8GFq3bh3jxo1z36/v3DxjxgwWL17M9OnTOXz4MA8++CB5eXmkpqaydOnSUzpVi4iInI33fQuBcCZlHSHuamunjZC20ak6UFtBHahFmk+/E+JpPv9kC5f4FWJzufgioCdjxqRYXZI0k9d0oBYREWlLL6xaC8CozfsUhDyYwlAj1IFaRMTLGQaru5m9ScYeU99RT6Yw1IhZs2axdetW1q5da3UpIiJigc8/ymJbr2QcTie3TEu3uhxpQwpD0qTFixcTERFxxu1sNhtLlixp9f2f6XX37t2LzWYjKyur1fctIt7t8KJPmPP660xYuZ3evSKtLkfakMKQNGn69Ons3LnTff/hhx/WqM4i4hX6fvoOjz/zDHN2nXn+RencOvyl9WKtwMBAAgMDrS5DRKRdle8vYkCpOaVTz19OsLgaaWtqGWqEp3agfv/994mIiMDpdAKQlZWFzWZj3rx57m1uueUWrr/+eqDhabLFixczf/58Nm3ahM1mw2azsXjxYvfzCgsLmTJlCkFBQfTp04d33323yVpSUlL4wx/+wLXXXktwcDCJiYlkZGScsl1ubi6TJk0iMDCQnj178p///OeUbfbs2cO4ceMICgpi6NChfPXVVy391oiIuH380nKWn5/KpuBUul/keZN1S0MKQ404qw7UhgEVFdYszRwu6uKLL6asrIyNGzcCsGLFCqKiosjMzHRvs2LFCvds9iebPn06c+bM4bzzziM3N5fc3FymT5/ufnz+/PlMmzaNzZs3c/nll3PddddRVFTUZD2PPfYYQ4cOZePGjcybN48777yTTz/9tME2v/vd75g6dSqbNm3iuuuu46c//Snbtm1rsM3999/P3LlzycrKom/fvlx77bXU1dU163siIvJ9/6nJ57InnuC3836DzWZ1NdLWFIZaU2UlhIRYs1RWNqvE8PBwUlNT3eEnMzOTu+66i40bN1JeXs7BgwfZtWsXY8eOPeW5gYGBhISE4OPjQ1xcHHFxcQ1Ooc2cOZNrr72W3r1788gjj1BeXs6aNWuarOfCCy9k3rx59O3blzvuuINrrrmGP//5zw22+clPfsItt9xC3759+cMf/sCIESP461//2mCbuXPncsUVV9C3b1/mz5/Pvn37GkzjIiLSbIZBwsEdJBw+TFpYotXVSDtQGPJCY8eOJTMzE8Mw+Pzzz/nxj3/MgAED+OKLL1ixYgUJCQn06dOnxa87ZMgQ99fBwcGEhYVRUFDQ5HO+P6Hu6NGjT2n1ac42J+87Pj4e4Iz7FhE5ncPr9vHEoqfInvYzfnnVuU36LZ2DOlC3pqAgKC+3bt/NlJ6ezqJFi9i0aRO+vr7079+f9PR0MjMzKS4uPm2rUHP4+vo2uG+z2XC5XGf1Wueyb9vxNu322reIeJZ9L68kGtgeOJwh3YOtLkfagcJQa7LZILjj/+LU9xv685//7A4+6enpPProoxQXFzNnzpxGn+vn5+fufN0aVq9efcr9AQMGnLLuxhtvbHB/2LBhrVaDiMjJCteswWWzkdf3EoaceXPxADpN1ghPvZoMoEuXLgwZMoRXX33V3VH6kksuYcOGDezcubPJlqGUlBSys7PJysqisLCQ6urqc6pl1apV/OlPf2Lnzp1kZGTw5ptvcueddzbY5s0332TRokXs3LmThx56iDVr1nD77bef035FRBrz2xnnE/3223yjUae9hsJQIzx9Oo6xY8fidDrdYSgyMpKBAwcSFxdHv379Gn3e1KlTmThxIuPGjSM6Opp///vf51THnDlzWLduHcOGDeOPf/wjTz75JBMmNBzTY/78+bz22msMGTKEl156iX//+98MHDjwnPYrInI6RXty2dI7maLwcNJ+eL7V5Ug7sRlGM6/J9lKlpaWEh4dTUlJCWFhYg8eqqqrIzs6mR48eBAQEWFRh55WSksJvfvMbfvOb31hdirQS/U5IZ/fKU+9ww5Bwoo6UkD/lSux2XVffWTX1+f19ahkSERE57uu8QwD0zy5VEPIiCkMiIiLHfRthBqABNc2/Qlc6P11NJpbZu3ev1SWIiJzgdLK9ZzQAo3v1tLgYaU9qGRIREQG2fL6V3Kiu2FwurvjheVaXI+1IYagRnnxpvYiInCpz9VYAUg4UEtPFz+JqpD0pDDXC0y+tFxGRhraUHgYgJbfG4kqkvSkMiYiIALvDzdteNH0ZtngehSEREZGaGranxAAwoncvi4uR9qYwJCIiXu+7L77hQKwZhq4Y19fiaqS9KQxJkxYvXkxERMQZt7PZbCxZsqTN62ktmZmZ2Gw2jh492ug2zX3vp/O73/2O2267rclt0tPTWzz69pIlS+jduzcOh6NVR+5uzvGbN28ed9xxR6vtU6Qj+fzr7wDodqCQpEhfi6uR9qYwJE2aPn06O3fudN9/+OGHSU1NtayezhC68vLy+Mtf/sL999/f6q/9i1/8gmuuuYacnBz+8Ic/tPrrgzn+k81mIysrq8H6uXPn8uKLL7Jnz5422a+IlQZ9sIaiyZOZ9/wmq0sRCygMSZMCAwOJiYmxuoxO5bnnnmPMmDF07969VV+3vLycgoICJkyYQEJCAqGhoa36+mcSFRXFhAkT+Pvf/96u+xVpD0G7NtOlvJzuXXtYXYpYQGGoEZ46ztD7779PREQETqcTgKysLGw2G/PmzXNvc8stt3D99dcDDU8VLV68mPnz57Np0yZsNhs2m43Fixe7n1dYWMiUKVMICgqiT58+vPvuuw32vWLFCtLS0vD39yc+Pp558+ZRV1fnfjwlJYWFCxc2eE5qaioPP/yw+3GAKVOmYLPZ3Pe/r75l47XXXmPMmDEEBAQwaNAgVqxYccq2q1atYsiQIQQEBHDBBRewZcuWU7b5+OOPGTBgACEhIUycOJHc3NzT7rfea6+9xuTJkxusq6io4MYbbyQkJIT4+HieeOKJU55XXV3N3LlzSUxMJDg4mFGjRpGZmQmYp/Xqw8+ll16KzWYjMzOTI0eOcO2115KYmEhQUBCDBw/m3//+d4PXPdP39ft69DA/DIYNG4bNZiM9Pd392OTJk3nttdeafP8inY5hkFBotgiFXTTE4mLECgpDjTiXcYYqnM4WL3Uul/v5dS4XFU4nx44HljO9bktcfPHFlJWVsXHjRsAMKFFRUe4P3fp1J38A1ps+fTpz5szhvPPOIzc3l9zcXKZPn+5+fP78+UybNo3Nmzdz+eWXc91111FUVATAwYMHufzyyxk5ciSbNm3i73//O88//zx//OMfm117/bF44YUXyM3NPeOxufvuu5kzZw4bN25k9OjRTJ48mSNHjpyyzRNPPMHatWuJjo5m8uTJ1NbWuh+vrKzk8ccf5+WXX2blypXs37+fuXPnNrrPoqIitm7dyogRI07Zz4oVK3jnnXf45JNPyMzMZMOGDQ22uf322/nqq6947bXX2Lx5Mz/5yU+YOHEi3333HWPGjGHHjh0A/Pe//yU3N5cxY8ZQVVXF8OHD+eCDD9iyZQu33XYbN9xwA2vWrDnzN7QR9c/97LPPyM3N5a233nI/lpaWxoEDBzSViniUypw8rv/jHO761a+JmTjA6nLEApqbrA2EfP55i5/zxsCB/OT46ai3CwuZtnUrY8PDyRw2zL1NyurVFJ70QV3POE1waUx4eDipqalkZmYyYsQIMjMzueuuu5g/fz7l5eWUlJSwa9cuxo4de8pzAwMDCQkJwcfHh7i4uFMenzlzJtdeey0AjzzyCE899RRr1qxh4sSJ/O1vfyM5OZmnn34am81G//79OXToEPfeey8PPvggdvuZc3l0tDlnUERExGn3/3233347U6dOBeDvf/87S5cu5fnnn+eee+5xb/PQQw9x2WWXAfDiiy+SlJTE22+/zbRp0wCora3lmWeeoVevXu7X/P3vf9/oPvfv349hGCQkJLjXlZeX8/zzz/PKK6/wgx/8oMG+Tn7eCy+8wP79+93PnTt3LkuXLuWFF17gkUcecZ+ujIyMdL//xMTEBuHsjjvu4OOPP+aNN94gLS3tjN+j06n/Pnft2vWU73N9bfv27Wu0ZU6ks/ly+Td8dMEFrBiSyhN9NUGrN1IY8kJjx44lMzOTOXPm8Pnnn7NgwQLeeOMNvvjiC4qKikhISKBPnz4tft0hQ040LwcHBxMWFkZBQQEA27ZtY/To0dhsNvc2F154IeXl5Rw4cIBu3bqd+xv7ntGjR7u/9vHxYcSIEWzbtq3RbSIjI+nXr1+DbYKCgtxBCCA+Pt79nk7n2LFjAAQEBLjX7d69m5qaGkaNGnXKvup98803OJ1O+vZteElvdXU1Xbt2bXR/TqeTRx55hDfeeIODBw9SU1NDdXU1QUFt8wc9MDAQMFvMRDyFz8pveTbvM9ZFj8F++USryxELKAy1gfKLL27xc/xPCglToqIov/jiU85h7r3ggnOszJSens6iRYvYtGkTvr6+9O/fn/T0dDIzMykuLj5tq1Bz+Po2vBzVZrPhOun035nY7XYMw2iwrvY0LWHt6XTv6fs1niwqKgqA4uJidwtLc5SXl+NwOFi/fj0Oh6PBYyEhIY0+77HHHuMvf/kLCxcuZPDgwQQHB/Ob3/yGmpoT0wm05ve1/rRnS96bSEcXsWYDt2z5kK7Dx1hdilhEfYbaQLDD0eLF56TTRD52O8EOB4Hf+1Bs7LktVd9v6M9//rM7+NSHoczMzNP2F6rn5+fn7nzdEgMGDOCrr75q8KG8atUqQkND3aeLoqOjG3ROLi0tJTs7u8Hr+Pr6Nnv/q1evdn9dV1fH+vXrGTBgQKPbFBcXs3PnzlO2aYlevXoRFhbG1q1bG6zz9fXl66+/PmVf9YYNG4bT6aSgoIDevXs3WJo6Jbhq1Squuuoqrr/+eoYOHUrPnj0bvC407/t6Mj8/c4LK032ft2zZgq+vL+edpxm9xXOE798MgGOYOk97K4UhL9SlSxeGDBnCq6++6g4+l1xyCRs2bGDnzp1NtgylpKSQnZ1NVlYWhYWFVFdXN2ufv/71r8nJyeGOO+5g+/btvPPOOzz00EPMnj3b3V/o0ksv5eWXX+bzzz/nm2++YcaMGae0kqSkpLBs2TLy8vIoLi5ucp8ZGRm8/fbbbN++nVmzZlFcXMzNN9/cYJvf//73LFu2jC1btjBz5kyioqK4+uqrm/WeTsdutzN+/Hi++OIL97qQkBB+/vOfc/fdd/O///3Pva+T+0n17duX6667jhtvvJG33nqL7Oxs1qxZw4IFC/jggw8a3V+fPn349NNP+fLLL9m2bRu/+MUvyM/Pb7BNc76vJ4uJiSEwMJClS5eSn59PSUmJ+7HPP/+ciy++2H26TKTTq63l85FxfHneeYSlD7K6GrGIwpCXGjt2LE6n0x2GIiMjGThwIHFxcQ36snzf1KlTmThxIuPGjSM6OvqUy7gbk5iYyIcffsiaNWsYOnQov/zlL/n5z3/OAw884N7mvvvuY+zYsfzoRz/iiiuu4Oqrr27QXwfgiSee4NNPPyU5OZlhJ3UuP51HH32URx99lKFDh/LFF1/w7rvvuk9jnbzNnXfeyfDhw8nLy+O9995zt4ycrVtuuYXXXnutwSnCxx57jIsvvpjJkyczfvx4LrroIoYPH97geS+88AI33ngjc+bMoV+/flx99dWsXbu2yf5UDzzwAOeffz4TJkwgPT2duLi4U8Jcc76vJ/Px8eGpp57iH//4BwkJCVx11VXux1577TVuvfXWFn5HRDquI1t2M/O393Lh008TdUm81eWIRWxGUx0ghNLSUsLDwykpKSEsrOFMxlVVVWRnZ9OjR48GHWbFWnv37qVHjx5s3LjRktGyDcNg1KhR3HXXXe6r6zzBRx99xJw5c9i8eTM+PqfvbqjfCelslvz9PaYMCKXL0XKOXHVFg4s8pHNr6vP7+9QyJNLKbDYb//znPxsMKOkJKioqeOGFFxoNQiKd0eZD5mnlbgdLFYS8mP6qibSB1NRUS+dwawvXXHON1SWItLq9zgoA4ksVhLyZwlAjMjIyyMjIOKsrp8RaKSkpTV7+LiJS72CoeYIkydH0aRTxbDpN1ohzmY5DREQ6AcMgJ84MQX3j1XnamykMiYiIV6rZf5A9iWYIGjmst8XViJUUhlqBTsmImFoy4riI1Tau2k61nx9+NbWM6R9qdTliIfUZOge+vr7YbDYOHz5MdHS0rkQQr2UYBjU1NRw+fBi73X7OYzWJtIesXQcgIYWE3BL8fPT325spDJ0Dh8NBUlISBw4cYO/evVaXI2K5oKAgunXr1mB0bZGOandZMZBC3BHPGgZDWk5h6ByFhITQp08fyycUFbGaw+HAx8dHLaTSaRzwM0NQXI2ml/F2CkOtwOFwNDnXk4iIdDwHoswQ1CM06gxbiqdTW7aIiHifigrCqsoIqaxkUJ/G5/8T76CWIRER8TqVW/bw/m9/yxG6YBQesbocsZhahkRExOvkfr4LgH0+vYnqqn5u3k5hSEREvE55lhmGCiM02KLoNJmIiHihl2IdLH3hBUZtrOWHVhcjllPLkIiIeJ19ITa2pqRQGR9jdSnSASgMiYiI17n/+Rf5ZO5cfmDosnrRaTIREfE2VVUMPbgZ+0GDNfP7Wl2NdABqGWpERkYGAwcOZOTIkVaXIiIirahmRzZ2DEoIo9v5ahkSsBmacr1JpaWlhIeHU1JSQlhYmNXliIjIOVr/3Pu898V/CMhxcO9nz6MZZDxTSz6/dZpMRES8yvoDecyfOZOUfUeYpyAk6DSZiIh4mX3lpQBEFTstrkQ6CoUhERHxKod8agGIqfG3uBLpKBSGRETEq+SG+wGQEBBhbSHSYSgMiYiI96ir40BsBAC9ExKsrUU6DIUhERHxGq79OWTHxwGQ2j/e4mqko1AYEhERr7Fj0z4qAwOxO11c2DvI6nKkg1AYEhERr/HNd4cAiCksJSRAH4Fi0k+CiIh4jV1FRwGILqy2thDpUBSGRETEaxxyVQEQVaHRFuUEhSEREfEaBQEuAKKNQIsrkY5EYUhERLxGQYQ50GJ8aKTFlUhHornJRETEO9TUMHD/HorDg+ibeKnV1UgHojAkIiJewcg5wN/+spBjBJC7q9LqcqQD0WkyERHxCkc3ZgOwlxSSktWBWk5QGBIREa9wJGsvBpDr3wM/P6urkY5Ep8lERMQrvFdbyYMffMD56wtQjyE5mVqGRETEKxyqraQ8KIjaoFCrS5EOxivC0JQpU+jSpQvXXHON1aWIiIhFbvrvR+y44QamfmN1JdLReEUYuvPOO3nppZesLkNERCwUU7CLvgcOkBLXzepSpIPxijCUnp5OaKiaRUVEvFZ1NZFV5iStEakp1tYiHY7lYWjlypVMnjyZhIQEbDYbS5YsOWWbjIwMUlJSCAgIYNSoUaxZs6b9CxURkc5r/35u/82dPHDjzwk9L8LqaqSDsfxqsoqKCoYOHcrNN9/Mj3/841Mef/3115k9ezbPPPMMo0aNYuHChUyYMIEdO3YQExMDQGpqKnV1dac895NPPiEhIaHN34OIiHRsB7fs5e9XXQXA/+tueTuAdDCWh6FJkyYxadKkRh9/8sknufXWW7npppsAeOaZZ/jggw9YtGgR8+bNAyArK6vV6qmurqa6utp9v7S0tNVeW0RErPHtrjwYmUxYSSUxYQ6ry5EOpkPH45qaGtavX8/48ePd6+x2O+PHj+err75qk30uWLCA8PBw95KcnNwm+xERkfazu7AYgK5HqiyuRDqiDh2GCgsLcTqdxMbGNlgfGxtLXl5es19n/Pjx/OQnP+HDDz8kKSmpySB13333UVJS4l5ycnLOun4REekYDtQcA6BrqWFxJdIRWX6arD189tlnzd7W398ff3//NqxGRETaW56fGYKianwtrkQ6og7dMhQVFYXD4SA/P7/B+vz8fOLi4iyqSkREOpuCUPP//nF+YRZXIh1Rhw5Dfn5+DB8+nGXLlrnXuVwuli1bxujRo9t03xkZGQwcOJCRI0e26X5ERKSNGQaHuwQBkBjZ1eJipCOy/DRZeXk5u3btct/Pzs4mKyuLyMhIunXrxuzZs5kxYwYjRowgLS2NhQsXUlFR4b66rK3MmjWLWbNmUVpaSnh4eJvuS0RE2lBREbldIwHo0T3e4mKkI7I8DK1bt45x48a578+ePRuAGTNmsHjxYqZPn87hw4d58MEHycvLIzU1laVLl57SqVpEROR06rL3cSgqCoBB3UIsrkY6IpthGOpa34T6lqGSkhLCwnSuWUSks9n5ynv0SwrF5nJRfmE6Qf4duoeItJKWfH7rJ0JERDzad/sLAIg4ekxBSE5LPxWNUAdqERHPcPCIOeBiRHGNxZVIR6Uw1IhZs2axdetW1q5da3UpIiJyDkJyDvGjL7+k916FITk9hSEREfFoY1av5b3772fa1zarS5EOSmFIREQ8WkixOa2STw/NNSmnpzAkIiKey+kkoLYAAwjq383qaqSDUhhqhDpQi4h4gLw8Lvh7BqEffkjOwECrq5EOSmGoEepALSLiAXJyONS1KxWBgfSO0yTccnqWj0AtIiLSVqq/20/ObTN4t+sPGb35bavLkQ5KLUMiIuKxyrbmEFxVhSMvmOgu+siT09NPhoiIeKyq78wryY6GJWPTlfXSCJ0mExERj/Wlw8nHd9/NsfzB3GZ1MdJhqWWoEbqaTESk8/smwp9Fl1/O9uGxVpciHZjCUCN0NZmISOd3JND8mOtiC7a4EunIFIZERMQzVVdzOCIIgKiQSIuLkY5MYUhERDzTgQMcjIoCILmrwpA0TmFIREQ8U06OOwz1i9KAi9I4hSEREfFIdftyOHQ8DA3W6NPSBIUhERHxSPv35FLn44PNZTAk0c/qcqQDUxhqhC6tFxHp3LILigEIKakhyF8fd9I4/XQ0QpfWi4h0bocqygEIKTEsrkQ6OoUhERHxSAVGLQChlb4WVyIdncKQiIh4pMLjAy6GGSEWVyIdncKQiIh4nvJyCiLMEBQZ0MXiYqSjUxgSERHPk5PDsO++Y9Kqr+kfrAEXpWmatV5ERDzPgQPcvmQJ6Uu+49DH91pdjXRwahkSERHPk5Nj3pBMYqLFtUiHpzDUCI0zJCLSeVXtyaE4JIT9CkPSDDbDMDQAQxNKS0sJDw+npKSEsLAwq8sREZFm2HjLbzj/+qsJPlpD2VU/xGazuiJpby35/FbLkIiIeJzcUnPARb9Km4KQnJHCkIiIeJwxX35FxcSJXPOMkpCcma4mExERjxNYmIN/dTXhEclWlyKdgFqGRETEs5SW4l9dBoB/rySLi5HOQC1DIiLiWXJyeOS66/g2sQcxfeqsrkY6AbUMiYiIZ8nJ4YMLLuBfk36AI7ba6mqkE1AYEhERz3LgAAejogDo19Xf4mKkM1AYEhERj+Lcn8Ohrl0BGBSnMCRnpjAkIiIe5UDOYWp9fcFlMDTRz+pypBNQGGqEpuMQEemc9peUAhBU4iLIXx9zcmb6KWnErFmz2Lp1K2vXrrW6FBERaYHcWrPTdFCpPuKkefSTIiIinsMwKPQ1p9wMqgq0uBjpLBSGRETEcxw9Sl6EOSlnmE2Ta0vzKAyJiIjnOOmy+lj/IIuLkc5CYUhERDxHzonL6rsF6bJ6aR6FIRER8Rw5Oe6WoT4RCkPSPApDIiLiOU46TTYgRmMMSfMoDImIiMeo2ZvDtf/7H32+KGJYUoDV5UgnoTAkIiIeo253Dk8/9RQX/6GSbpE+VpcjnYTCkIiIeI4DBwCoik7GZrO4Fuk0FJtFRMQzGAZVpYVUh4TgSkqyuhrpRNQyJCIinqGggEUTxhP53nusuanK6mqkE1EYEhERz7BvH0Vh5qjTXX3UeVqaT2FIREQ8w759PPLcc3w24WGmVna3uhrpRBSGRETEM+zbB0BeTSK94tUlVppPYagRGRkZDBw4kJEjR1pdioiINMfevQDsozuJidaWIp2LzTAMw+oiOrLS0lLCw8MpKSkhLEwzIIuIdFTOH03mRxdfxJ6i/rz1qys4L0WtQ96sJZ/f+kkRERGPkF9whKWjRgGQEq8TH9J8+mkREZHOzzA4VF0JgO9RO8H++niT5tNPi4iIdH5Hj5IXFAiAX7kuq5eWURgSEZHOb98+DkZHAxBSHWhxMdLZKAyJiEjnt28fB6OiAOhq+FlcjHQ2ZxWGdu/ezQMPPMC1115LQUEBAB999BHffvttqxYnIiLSLCeFoThff4uLkc6mxWFoxYoVDB48mK+//pq33nqL8vJyADZt2sRDDz3U6gWKiIic0d697jDUPUhhSFqmxWFo3rx5/PGPf+TTTz/Fz+9EU+Sll17K6tWrW7U4ERGRZjmpZah3hMKQtEyLw9A333zDlClTTlkfExNDYWFhqxQlIiLSInv2uMPQwBj1GZKWaXEYioiIIDc395T1GzduJFHjn4uISHszDCpzcig+Psrw0AS1DEnLtDgM/fSnP+Xee+8lLy8Pm82Gy+Vi1apVzJ07lxtvvLEtahQREWnc4cPk+h8PQFV2unfV5ArSMi0OQ4888gj9+/cnOTmZ8vJyBg4cyCWXXMKYMWN44IEH2qJGERGRxu3e7T5F5lvij91us7gg6WxaHIZWrVrFs88+y+7du3n//fd55ZVX2L59Oy+//DLPPfdcW9QoIiLSuF27iCop4aq31xC/I9rqaqQTanEYmjhxInfffTfx8fFcfvnlTJs2jS5dunDllVcyb968tqhRRESkcbt3M3DfPi5/ahcX7ehpdTXSCbU4DC1fvpy3336bkSNHsnXrVj744AMGDRrE0aNHycrKaoMSRUREmrBrFwC76UVSksW1SKfU4jA0ZswYsrKyGDRoEOeffz5TpkzhrrvuYsWKFXTv3r0tahQREWnc7t3sj4nh27A+JCQaVlcjndBZTcexc+dO1q1bR1JSEj4+PuzYsYPKysrWrk1EROTMdu3ip7/7HR+8E8mBHhrvTlquxWHo0UcfZfTo0Vx22WVs2bKFNWvWsHHjRoYMGcJXX33VFjWKiIicXkkJFBZy7Pil9f2jNMaQtJzNMIwWtSnGx8ezaNEiJk2a5F5XW1vLb3/7W5566imqq6tbvUgrlZaWEh4eTklJCWHHB/QSEZEOYt06GDmSPGKJ9z3I3j02uied1UkP8TAt+fw+q+k4Tg5CAL6+vjz22GN88sknLX25NpeTk0N6ejoDBw5kyJAhvPnmm1aXJCIirWXrVgC2MQCHy0FSvIKQtFyLh+mMOj6w1emMHTv2nIppCz4+PixcuJDU1FTy8vIYPnw4l19+OcHBwVaXJiIi52rbNgC2MpC4OHA4LK5HOiWPH7M8Pj6e+Ph4AOLi4oiKiqKoqEhhSETEE2zdyrLzz+fRaydhK9oLpFhckHRGlrcnrly5ksmTJ5OQkIDNZmPJkiWnbJORkUFKSgoBAQGMGjWKNWvWnNW+1q9fj9PpJDk5+RyrFhGRDmHbNrZ168aBESG4epZbXY10UpaHoYqKCoYOHUpGRsZpH3/99deZPXs2Dz30EBs2bGDo0KFMmDCBgoIC9zapqakMGjTolOXQoUPubYqKirjxxhv55z//2WQ91dXVlJaWNlhERKQDqq5uMC9ZNLqSTM6O5afJJk2adEqH7JM9+eST3Hrrrdx0000APPPMM3zwwQcsWrTIPf3HmUa+rq6u5uqrr2bevHmMGTOmyW0XLFjA/PnzW/YmRESk/e3cCS4Xe2MTAUjwUxiSs2N5y1BTampqWL9+PePHj3evs9vtjB8/vtljGhmGwcyZM7n00ku54YYbzrj9fffdR0lJiXvJyck56/pFRKQNHe88vTu2GwA9QhWG5Ox06DBUWFiI0+kkNja2wfrY2Fjy8vKa9RqrVq3i9ddfZ8mSJaSmppKamso333zT6Pb+/v6EhYU1WEREpAPavBmAA5HmabK+XRSG5OxYfpqsrV100UW4XC6ryxARkda2cSMGcLhrCACD4xWG5Ox06JahqKgoHA4H+fn5Ddbn5+cTFxfXpvvOyMhg4MCBjBw5sk33IyIiZ2njRo6GhFAXYH6UDUvys7gg6aw6dBjy8/Nj+PDhLFu2zL3O5XKxbNkyRo8e3ab7njVrFlu3bmXt2rVtuh8RETkLeXmQm0tOTAwAtlIfugRrxEU5O5afJisvL2fXrl3u+9nZ2WRlZREZGUm3bt2YPXs2M2bMYMSIEaSlpbFw4UIqKircV5eJiIgX2rgRgN39hwHgX6ZTZHL2LA9D69atY9y4ce77s2fPBmDGjBksXryY6dOnc/jwYR588EHy8vJITU1l6dKlp3SqFhERL3I8DG3rMQSA0CqFITl7loeh9PR0DMNocpvbb7+d22+/vZ0qEhGRDm/1agC2xvYBoKuhMCRnr0P3GbKSOlCLiHRQLhesWgXA9ogkAOJ9FIbk7CkMNUIdqEVEOqjt26GoCAIDce3sDW8lcr6/xoSTs2f5aTIREZEWOd4qxKhRHPssHrbDpKusLUk6N7UMiYhI57JyJQDGhRexf7+5qls3C+uRTk9hSEREOg+XCz7+GICC4eOojC+D0FqSky2uSzo1nSZrREZGBhkZGTidTqtLERGRehs2wOHDEBLCsrhUeG499mI/AgLGWF2ZdGJqGWqEOlCLiHRAH31k3o4fz+5CGxT5EqgBF+UcKQyJiEjn8Z//mLdXXEHYni4w9UImfjDM2pqk01MYEhGRzuGbb2DzZvDzg6lT2bfPXJ2SrI8yOTf6CRIRkc7hpZfM2yuugC5d3FeSde9uXUniGdSBWkREOr7SUvjnP82vj0/U/b+Lv4XRtdSl9AJCratNOj2FoUboajIR8QrV1bBnDxw4AHl5J5bCQqisNJeKCjh2DGw2cDjMxc8PwsOhSxeIiDBvExPNZppu3SA52dymtWRkmIGof3+zZQg42v0oRNSSENJ6uxHvZDPONEuqlystLSU8PJySkhLCwjTcu4h0UnV1sHUrrFtn9r3ZsQN27oTsbHPsntZms0GPHjBkCAwebC7DhkGvXuZjLbF/PwwcaIayl1+G66+nqMJJ17WfA7Br4IX0ivFt/fcgnVpLPr/VMiQi4omOHIEVK8zRmteuhY0bzdad0wkNNVt04uIgPt68jYqCkBAICjKXgABzW6fTXKqroaQEiovh6FFzfwcOmMFl3z6oqjJbnPbsgSVLTuwrJgYuushcLrzQDEi+TQSZsjKYOtUMQhdeCNddB8CavdXm45UOekTpo0zOjX6CREQ8QVkZLF8O//ufebt586nbhIXB8OGQmmqeburXD/r2NcNPS1trmmIYUFAA335rtkLVXwW2aZO5/q23zAUgMBDS0mDMGHMZMcIMTDU15vu4+27zdbp2hVdecdeZlVsFPuB/1B+7vRVrF6+kMCQi0lnt3QvvvQfvvw+ZmWaAONnAgZCefiJk9OkD9na4iNhmg9hYc7n00hPrq6pg/Xr44gtzWbXKbFlascJc6tntDU/dxcaa7zElxb1qW1EVxEDYsYC2fz/i8RSGREQ6C8OArCx48014912zxeRkPXvCZZfBuHFmCIqNtaLKxgUEmKe6LrwQ7r3XDDw7dsCXX5rLqlXm/fogFBsL06fD/febrUUn2VluhqE4Q2FIzp3CkIhIR2YYsGULvP46vPEGfPfdicccDrPvzY9+ZC79+rXu6a62ZrfDgAHm8vOfm+tqaswr2QICzCvUGnk/B1xm/6eeAYHtVa14MIWhRujSehGx1Pbt8NprZgjavv3E+oAA89LyqVNh4kQzMHgSPz9ISDjjZkX+Zhg6L1JhSM6dLq0/A11aLyLtprTUDD+LFsHq1SfW+/nBpEnmKaPJk82rvLyc/f0vMELqeCtiBFNS9f2QU+nSehGRzsIwzMvfFy0y+wLVX/7ucMCECfDTn8KVV5oDHAoAuwtqMULqALi4l1qG5NwpDImIWOHAAXjxRXjhBdi9+8T6/v3h5pvhhhvMS97lFKuyzcBoL/YjKtRhcTXiCRSGRETaS3W1eRXYokXwyScnrpoKDTVbgG66CS64oHN1grbA+twqiIDgUl1JJq1DYUhEpK1t2mQGoFdegaKiE+vHjjVbgaZOheBg6+rrZOzZwbC6J0MGteLcZ+LVFIZERNpCURH8+99mCNqw4cT6xESYOdNceve2qrpOrXRLMLwRzA/Ps7oS8RQKQyIircXpNKfDWLQI3n7bPC0G5txbV19ttgJddpnZOVrO2p495m3PntbWIZ5DYagRGmdIRJotOxsWLzaX/ftPrB8yxBxM8Gc/Myc+lVaxJbQQegTQvWcQ0A7Ti4jH0zhDZ6BxhkTktCorzclGFy0yJxStFxFhzqx+883mjOzqDN2qCkqdxG74HICdAy6kT2wTM96LV9M4QyIibcEwYO1aMwD9+9/mIIlgBp7LLjMD0FVXmaNES5tY/10tfBuGPbKWPukKQtI6FIZERM6koMC8EmzRooaTo/boYV4OP2MGdOtmXX1e5OjOALj9fMZcBFxrdTXiKRSGREROp7YWPvzQHBTxgw+gzhzxmIAAuOYasxVo7FhzslFpN/XTtPXvb20d4lkUhkRETrZlixmAXnnFbBGql5ZmBqCf/lRTY1ho23YDsCkMSatSGBIRKS42+wC98AKsW3difWysOS3GTTfBwIHW1Sdu71+xHn7kxB41ANBFLdI6FIZExDs5nfDZZ2YAWrKk4ZhAkyebAWjCBPO+dAh1ToNjsZXg72JInD6+pPXop0lEvMv27fDyy/DSS+ZkqfWGDjUD0M9+BtHR1tUnjVq9pwr8XVBj46LeumJPWo/CkIh4vrw8eO01sx/Q+vUn1kdGmmMC3XSTOSaQdGiZ2ZXgB/6Fgfj7qOO6tB6FIRHxTOXl5pQYr7xing6rnyHexwcmTjQvh588Gfz9ra1Tmu2r/ApIhuhKTWorrUthqBGajkOkE6qrg08/NQPQkiXmKNH1Ro+G66+HadM0NUYnta26AoD+viEWVyKeRmGoEbNmzWLWrFnu4bxFpINyuWD1avM02GuvweHDJx7r08cMQNddB716WVejtIrcwHIALohWy5C0LoUhEel8DAPWrIHXX4c332zYETo6Gq691gxBI0ZobjAPUV7loirGbOmb2FctQ9K6FIZEpHMwDLPzc30A2rfvxGOhoXDllWYL0PjxuhzeAy3dWgm+BlQ4GN1T/bykdSkMiUjHZRiwcSO88Ya5ZGefeCwkxAxA06aZ4wFpclSPtmxvOURC2OEQ7Ha19knrUhgSkY7F6YSvvjKvBFuyBPbsOfFYUJB5Bdi0aTBpEgQGWlamtK91xRUQCd3r1F9IWp/CkIhYr6oKli0zA9C77zbsBB0YCFdcYQagK64wA5F4nT02s/P0kFD1F5LWpzAkItY4etScFX7JEvjoI3NcoHoREWYL0NVXww9/aJ4SE6/lchkUR5cBcFl3/SxI61MYEpH2s3MnfPCBGYJWrIDa2hOPJSaa4WfKFLjkEnWCFrfdu8H41fk4Bpcy9W2FIWl9CkMi0naqqmDlyhMBaNeuho8PHHgiAA0frsvg5bTWrrXBwSBGJAURon7y0gYUhkSkdR04YAafDz4wp8E4eRRoX18YO9bs+3P55dC3r3V1SqexZo15m5ZmbR3iuRSGROTclJaap7w+/dQMP9u2NXw8IcEMPldcAT/4gTkmkEgL/DcsG37iQ78LYgE/q8sRD6QwJCItU1tr/lf9s8/MAPT11+acYPXsdrjgghMBaOhQnf6Ss1ZZ7eLAhTlwqYveMV1RGJK2oDAkIk1zuWDrVli+3AxAy5dDWVnDbXr3Nkd+vuwyGDcOunSxplbxOBu2uOCl7vj2r+AHv9W4UtI2FIZEpCGnE7KyzI7PK1bA559DUVHDbbp2NU951QeglBQrKhUvsO5zH/hXd8ZPAh+H1dWIp1IYEvF2NTXmnF8rVpgB6IsvTm35CQqCMWPMAHTZZTBsmHk6TKSNrVhh3o4da20d4tkUhhqRkZFBRkYGTqfT6lJEWldhodnPZ/Vq+PJLc+qLY8cabhMWBhdfbI73c8kl5mXvGvdH2lmd0+BTVwFEhTN2rK6pl7ZjMwzDsLqIjqy0tJTw8HBKSkoICwuzuhyRlqmthc2bzeBTv3x/rB8wT3vVB5+xY2HIEHDonIRY662N5UwtWQdVdirGXUSQv1ojpfla8vmtliERT2EYsG+fecqrPvisW2cOfPh9/fubV3yNGmW2AA0YoNNe0uG8tq0EEiAyN1xBSNqUwpBIZ+RymS08GzY0XIqLT922Sxcz9FxwgbmkpelqL+kUvqw8CsD59ghL6xDPpzAk0tHV1sL27bBx44nQs3Fjw4lN6/n6wqBBDcNPnz5q9ZFOp6rWxaE4M9z/pHeEtcWIx1MYEukoiovN0PP9Zc+ehoMa1gsIgNRUOP/8E8t554GfBqWTzu+F1WUYIXXYynyYOUn9NaVtKQyJtKeKCjPc7N5tLjt2mMv27VBQ0PjzQkPNy9lPDj79+oGPfoXFM7266wj0gOT8SPx8NIK5tC39JRVpTYYBR46YQWfXrhOhp37Jy2v6+UlJZufm+qVfP/M2MVFTWohX2eBrDvQ5oUukxZWIN1AYEmkJwzBPZ+3bB/v3N7ytDzylpU2/Rpcu0KuXufTteyL49O0LISHt8z5EOrDNB6o5lmT2ibtztMKQtD2FIZGTVVTAoUPmcvCgGXK+H3xO13H5+xITzfm66kPPyYuu5BJp0uOrjkAsBOWEcl66+sBJ21MYEu9QU2Oeojp48ETYqQ88J399pladejEx0L07dOt24rZnTzPs9OgBgZpQUuRsfVhxGICLXFEWVyLeQmFIOieXy5w89PDh0y8FBSe+zsszb5srONhs2UlIOBF2Tg4+yckKOyJt5Lv8Wo50Ny+pn5sWbXE14i0UhsR6hgGVlXD0qNkfp7GAc3LIOXLEDEQt4etrBpyEhBNh5/tfJySY83KJiCUeWVkI0RBwMJjL0oOsLke8hMKQnDuXy+xHU1Z2YiktNYNNfcBp7Ov629ras9t3RARERzdcYmIa3o+NNQNPZKQGHxTp4LYtDYLYWC7trYsJpP0oDHkTl8ucp6qqypyl/OSvKyoahpkzLSeHn4qK1qnP4TA7F0dFnRpwThd2oqI0k7qIB8nLg3UvhoMznCe3W12NeBOFIascPXpqIDnbr5u7XU1N274nHx9zcMCQEPNUU5cu5hIR0fC2sXXBwRpLR8SLvfIKOJ3mLDL9+lldjXgThSGr9OpldgC2io+POZ1D/RIcbAaZc1n8/RVmROSsuFwGfzq4H3p05aabdIpM2pfCkFUCAszgEBhofl1/ezZfn81zNI2DiHQgi78u4/BV2TBhH1ekjUEfT9Ke9NNmlb17zUCilhQRET542w4hUfSK9SFxoj6apH3pJ84q6vgrIgKYPQY+ejoEjg3iueWG1eWIF9J1xiIiYqlnnzWv9Rg6FMaOVWu5tD+PD0NHjx5lxIgRpKamMmjQIJ599lmrSxIRkeOqal08kpcN8ce48071HBBr2AzD8Og2SafTSXV1NUFBQVRUVDBo0CDWrVtH165dm/X80tJSwsPDKSkpIUwjE4uItKo7lxTwVMRW7MV+lEy8gJBAj/8/urSTlnx+e/xPncPhICjIHNK9uroawzDw8PwnItIp1DkN/lmxD4CLS+IVhMQylv/krVy5ksmTJ5OQkIDNZmPJkiWnbJORkUFKSgoBAQGMGjWKNWvWtGgfR48eZejQoSQlJXH33XcTFaWZkEVErPbbDwupSqyACgcvXJFkdTnixSwPQxUVFQwdOpSMjIzTPv76668ze/ZsHnroITZs2MDQoUOZMGECBQUF7m3q+wN9fzl06BAAERERbNq0iezsbP71r3+Rn5/faD3V1dWUlpY2WEREpHW5XAZ/Ld0LwEUFSfSI1hW2Yp0O1WfIZrPx9ttvc/XVV7vXjRo1ipEjR/L0008D4HK5SE5O5o477mDevHkt3sevf/1rLr30Uq655prTPv7www8zf/78U9arz5CISOu5573DPBb6LVQ62Dn8AvrEKgxJ6/KYPkM1NTWsX7+e8ePHu9fZ7XbGjx/PV1991azXyM/Pp6ysDDADzcqVK+nXxKQ39913HyUlJe4lJyfn3N6EiIg0UF7lYmHlHgAuyk9SEBLLdehBFwsLC3E6ncTGxjZYHxsby/btzZvSeN++fdx2223ujtN33HEHgwcPbnR7f39//P39z6luERFp3A3/PUht4jHsR315fUqy1eWIdOww1BrS0tLIysqyugwREQF2F9TyTph5Bdm11T1IiPD4jyHpBDr0abKoqCgcDscpHZ7z8/OJi4tr031nZGQwcOBARo4c2ab7ERHxJj96dw9GaB3+h4J57sdt+3dcpLk6dBjy8/Nj+PDhLFu2zL3O5XKxbNkyRo8e3ab7njVrFlu3bmXt2rVtuh8REW/xt8+Psr13LgCPJvQmwLdDfwSJF7G8fbK8vJxdu3a572dnZ5OVlUVkZCTdunVj9uzZzJgxgxEjRpCWlsbChQupqKjgpptusrBqERFpicoqg7tydkIC9N0Vx29u6WJ1SSJuloehdevWMW7cOPf92bNnAzBjxgwWL17M9OnTOXz4MA8++CB5eXmkpqaydOnSUzpVi4hIx/XoIzZq3u2L723ZfPSTXlaXI9JAhxpnqCPS3GQiIudm1Sq45BJwueD112HaNKsrEm/gMeMMWUkdqEVEzt3+ojqmza3A5YIbblAQko5JLUNnoJYhEZGz43IZ9Fy0jX1JhUS/2I9d/4hFf0alvahlSERELPfk0y72FdaBr4s/zApQEJIOy/IO1CIi4nmWL4d5sx1gDOaOZ0r5xa3hVpck0ii1DImISKvauLuWn0wzcDrh+p/Z+MstCkLSsSkMNUIdqEVEWm5HXg0XrN3AkZt3MGyEi3/+E2w2q6sSaZo6UJ+BOlCLiDRPQamTnkuyqOhWhqPQn3Ujh5Pa3c/qssRLqQO1iIi0q/IqF4Pf/JaKbmXYynxY0m+IgpB0GgpDIiJyTkqPuej96rcU9CqCajt/7zKYHw0OtroskWbT1WQiInLWSo+56POvbynodQSq7SxwDOIXF6nDtHQuCkMiInJW8krqGPzfrRQebxH6P59B3PODSKvLEmkxnSZrhK4mExFp3LeHauj17iYKexZBlZ0/+QxWEJJOS1eTnYGuJhMRaejTbZVc/s1m6mKqsJX68mz0YH4+Wn8fpWNpyee3TpOJiEizLfikmPtrvsWIqcPncADvnzeECQODrC5L5JzoNJmIiJyRYcCCBfDbO3ww/JwE7w9l/ehhCkLiEdQyJCIiTSoocvGrW+289RZAKFdmpvLy70IJC9T/p8UzKAyJiEijnlxezD2Hd+DcOhA/vzCefhpu1aSr4mEUhkRE5BQ1NfDgg/B/frlwaRVBv97LitFDGDHC6spEWp/CUCMyMjLIyMjA6XRaXYqISLta+YXBr35hY+tWIKQPQ+L8+WRGd2J1wZh4KF1afwa6tF5EvMX+ojouf2sP3+ZXwwODiI628Y9/wJQpVlcm0nK6tF5ERJqtzmlw+5J8nrXtwdW7BnrD5LvLWDwvjEiNoyheQGFIRMSL/e3zo9y9dzeVyWUA+OYH8mhMX2b/SS3h4j0UhkREvNCHWyq4bfVeDvY+DMlApYMJhd147ZokIoIcVpcn0q4UhkREvMjHWyv51Zf7yO6RD70BJ/TPjufNST0YlOhndXkillAYEhHxAh9vreTXX+5jT30IAuJ2d+Wvw3pwzQ9CrC1OxGIaPlRExEMZBqxYAROm1zDx0Fr29M4HB8Ts7srLQcPJ/flgrjlfQUhELUON0DhDItJZVVa7ePT9Uj5cEMH69QB+0D+K2G4uFgzszk0/V+dokZNpnKEz0DhDItJZ7N4Nf3vBycLzv8YVUQM3pBFQFMSMGXDHb1yc118nA8R7aJwhEREvUV7l4smlpXyeEcFnnwE4YEEI9v5l/GzeMZ78aRDR0aBeESKNUxgSEelk6pwGGZ+X8MyufHbEHsYIrYNNF2CzBfDDH8JP+vXj2km+BPkrAIk0h8KQiEgn4HIZ/GdjBU9symd9eAHOrtXuq8LsxX5ce+8x/vDjAHr0APC3slSRTkdhSESkg6qpM3j+q1IW7SwkK/QwdTFV0PP4gxUO+uRGc2uPWO6cHIGfj83SWkU6M4UhEZEO5Ngx+HNmMS8dLGBnVCFGRC30Ov5gjY2E/V35WWws9/8gUiNFi7QShSEREQu5XAbvbK5g17JgPv3YxsqVUP3/8uHyPABs5T6k5HVlamwUcy7pQtwP9WdbpLXpt0pEpJ0VFMDy5bD0Y4NXxq2jLrkCnh8G28IBiNoUS3xfO9f3iOL2yyLUEVqkjSkMiYi0sbV7q1i8sYTlR46yz1VJ5S9SAZu5dAuC6GMMvvwYN90azoQJMGBAF2y2LhZXLeI9FIZERFpRndPg0+3H+O+2ElYcLWFv2FGz43MXzAUgoYrBXQP54Q9hxIW9mXhhfyImqv+PiFUUhhqh6ThEpDl25NXw76wy/pdbylajlKKYMoyQOojCXACcEHQohAE1EfwwLoJfbfIjuf4xXQYvYjlNx3EGmo5DROrtK6zlyy115GcFsmYNfLmlln1PrTp1w2o7oXkhnOeMYFJiODenhZPURf/3FGlPmo5DROQcuFwG6/ZVs2uLg10bfcnKgi/88zn8i22wIQLmpB7f0hcOBOLnsJFcHsqIkDAm9wpjyuhgdXoW6UQUhkTEqx0ormPptgq+OFDJN2UVZNvLOdq1HCOsDp7rA+8mmhv2CQLAL7KOSVfB8OEwahQMGT6SuK4KPiKdmcKQiHiF+tCz6mAl35RWsNdWwdHwSnNaC4CY40s9J0T3q2HiDZCaCkOGBdNz4IX0TPeFmSe/soKQSGenMCQiHqOuDnZkO3lrz1G2FtQQuTqe776D7dsh557NMKgUojGXk9iL/YgoDqabEcTgkGDGp4Rw5eBgIn5w8hVedhR8RDyTwpCIdCoul0HWgRpW7qlkXf4xtpdXUrcrmKq349mzB2qDnfD2N5AIPBcDNccDzd5g7IlVZuhxBXFeaDAXJgQzsX8QPaJ9LX1PImIthSER6XD2F9Xx1Z5jbMirYtvRKvZWHSPPVkVJUBU1kVXg7zI3jD2+5HWFHfEA+Ff54sgOI8Lpx9R5ToalOOjbF/r070NM136WvScR6bgUhkSk3VVUwP79Bm/lFPNt0THiN8Wzf5ed7GzYMn4n1RMPmRuGHF++zwm+hYFElAeSSCDDuofxs0+hb19ISrJht59/mifpFJeInJ7CkIi0qtJjLjbkVLPpUDXbiqrILq/mQE01hbYqnAcCcf2lD8XFxzf+4FuIc8LdXeCAebUW5/kBYCvxJagkgMiaABLtgfQJDmBQVAAjkgK5IMVfl66LSKtRGBKRZquogNxcOHQIvjhcyleVJRyoqaLAqOaofzVVoVW4ImpPPCHy+FLPCIbjQSgszIbxXQRBIXDZjQbDo6FnT4jqnkSP7knER+jPk4i0D/21EREKSp1sOljN9oIaviuupuyInajt0Rw6ZIafr65fT1VsJfz6fNgfbD7pxiK4ae/pX7DGhm9xACGV/kQ6/Yl3BNA9yJ/zkgO5cgskJ4M5IOzg0zxZf5ZEpH3pr46Ih6qshIICyMo7xvqiCnLKa8k9VsPhulqKjBpKHDVUBlZTE1YDQSfNwRcBFAbD4yddfz7DCcFO6FpDUGEwCQng5xtC+a4YYm3+dAv0p09YAINi/Tk/yZ9+sb7Y7bb2fssiImdFYUikk6iqdbGroI4dBTXsKaplX2kNFUftxH8XTUGBGXxWXrGZishKfB8czLHtx1twrs+Hn++FoDPs4Jgd31J/go75EVUTxJV3QUICxMdDRcIAkqLsjPhfANHhYLNBw5lIRUQ6L4UhkXbidEJJibmsLSrnUFkt8SVhVB11UFICG2zFbAk6QplRRwVOjtnrqPapo8a/jrrAOoyQ2hMXRDmALkBxMCw4qQXnx1UQXUVdYA0QjL8/hFQHciwnlOAaX8JdfnS1+xLt50e3ED/6dPGjX7Q/QxP9SGiyj05om31fRESspjAk0oQ6p0FeqRNXhYPKchtlZbC1rJIdxyoprq7jaLWTklonZU4n5a46Kgwnx3BSZXOaQca3jjp/J+wNwjl36IkXfnMzRNXA3OGw63jQ+Fkp3Hqg6YJcYCvzxa/Cl8BqP7oeC+TyOyAmxlwKo/oS5Wsj7fVgesVBSAjYbPWD8YiIyOkoDDUiIyODjIwMnE7nmTcWS7lcBpU1BofLnPhU+3Cs0kZlJWwrr2RfVRUl1U5Ka1yU1Tgpr3NRXuek0uWk0uXi2PHwcszmpNpRh9/+ECJe6UNZGZSVQcUbX5j9aa4bBYcCzR3emgs/y4HgFhRZcuJXLTAQ6g4HYqtz0HcIJPSEiAio6BFGwZ5kwh0+RPg6iArwISrAh9gQH+JDfegf60efGF/8fJrqixNxFt9BERHvZjMMw7C6iI6stLSU8PBwSkpKCDMvf5EzcLkMyqpdlBwzqK0C3zofqqqg7JiLzZXllNe66F8dQVUVVFVBlnGUHCo5VufimNNFpdPlDitVOKnCSY3NRY3dSa3DRZ3Dif/WCAIX9aGyEsqPGbg+WWHu/KoxUGqOU8NvdsJVh1pW/OZwuHPYifv/XQWRtdh/MZywglBCQ6FuwiFKL87Fv84Hf5eDQBwE4yDY7kOow0GYr4MIPwdd/H2ICfIhNsRBcrgvg7oGEh4Ofn6t830WEZHGteTzWy1DnZDLBbW1UFVjUFRVR3UthBq+1Naa6/dUHaOkxklVnYuqOoPqOsO8dbqorjOodppLjct1/NYgrjyEuOIwqqrgiLOGL5L346q1cf66XlRVQXU1ZJ2/hyMxZdTZXTgdLpzHb10+LgwfFy5fF/i6wO+kfP1ePDx5fAqEECe8t8H8euIl4DreAeaBQ/CDghZ9D6r2BlCSV3/PBjU28DOwB7sIAYKDobrOn8qDwfjUOfB12fFzOvDHQYBhJ8DmIMjuIMhuJ9THQbifgwg/H5L6+nHBSggNNZfawDRiwu1EbLNjd4/xl3B8ERERT6AwZJGf/hQyL9pKVUg1LpthLnYDl92Fy25gOAwMu4HhcJlfOwxwGNie74HxZrL5IgNK4W8bITcAfnbBiRf/x7fQt7zhDh3Hl8ZaJZYnwz+OJ+c4J/z7AByzs+GXvU5sc3459Cw+/fMb4+ciIAD8/cE/1M6RQn8cdXb6DjMIckBAAOQZoZTscuKLHV/s+GMnyOYg0G4n2O4g2MdBqK+dEF8ztIT520ka68+gjWboCQqCOv8xdA1xEHLo5FGJux9fzoUm8BQR8XQKQxZZvhwKriyFhKoWPc/wdZ244zz+we/jwuEAX1/zFMyxCl+cR32xOe3YXDbsThs2lx27y2Yuhg2Hy27eGjYc2EiOCaL/DWY4sYX4snF3MgE2OxP+aK4LCIDdEUkcK4ohyMdOsO/xxc9OqJ+dED87oQHmEnZ8iQi0E3yJDfuC+oIdwGjzy5+e/K6Sjy/nQqFFRETOjvoMnUFb9Rl64w1Y7SjE8HHh57Dh77Dh77AT4GPD38dGgI+NAB/zfoCvjUBf836Uvw+RAT74+oLDx8DmYxDoZ8OhAe5ERETc1GeoE5g2Daad84B1tuOLiIiInC1N+ywiIiJeTWFIREREvJrCkIiIiHg1hSERERHxagpDIiIi4tUUhkRERMSrKQyJiIiIV1MYEhEREa+mMCQiIiJeTWFIREREvJrCkIiIiHg1hSERERHxagpDIiIi4tU0a/0ZGIYBQGlpqcWViIiISHPVf27Xf443RWHoDMrKygBITk62uBIRERFpqbKyMsLDw5vcxmY0JzJ5MZfLxaFDhwgNDcVmswEwcuRI1q5de9rtG3vs++tLS0tJTk4mJyeHsLCwtim+BZp6T+39ei15bnO2PdM2zT1mTa3vSMeztY/lubxmax/LM22n3822fb2O8rvZ2GP63Wyb53XW303DMCgrKyMhIQG7veleQWoZOgO73U5SUlKDdQ6Ho9ED0dhjja0PCwuz/BcUmn5P7f16LXluc7Y90zYtPWZNvV5HOJ6tfSzP5TVb+1ieaTv9brbt63WU383GHtPvZts8rzP/bp6pRaieOlCfhVmzZrX4saae0xG0dn3n8noteW5ztj3TNi09Zt52LM/lNVv7WJ5pO/1utu3rdZTfzcYe87ZjeS6vqd/NhnSazCKlpaWEh4dTUlJi+f9W5NzpeHoOHUvPouPpOdryWKplyCL+/v489NBD+Pv7W12KtAIdT8+hY+lZdDw9R1seS7UMiYiIiFdTy5CIiIh4NYUhERER8WoKQyIiIuLVFIZERETEqykMiYiIiFdTGOqA3n//ffr160efPn147rnnrC5HztGUKVPo0qUL11xzjdWlyDnKyckhPT2dgQMHMmTIEN58802rS5KzdPToUUaMGEFqaiqDBg3i2WeftbokaQWVlZV0796duXPntuh5urS+g6mrq2PgwIEsX76c8PBwhg8fzpdffknXrl2tLk3OUmZmJmVlZbz44ov85z//sbocOQe5ubnk5+eTmppKXl4ew4cPZ+fOnQQHB1tdmrSQ0+mkurqaoKAgKioqGDRoEOvWrdPf2k7u/vvvZ9euXSQnJ/P44483+3lqGepg1qxZw3nnnUdiYiIhISFMmjSJTz75xOqy5Bykp6cTGhpqdRnSCuLj40lNTQUgLi6OqKgoioqKrC1KzorD4SAoKAiA6upqDMNAbQOd23fffcf27duZNGlSi5+rMNTKVq5cyeTJk0lISMBms7FkyZJTtsnIyCAlJYWAgABGjRrFmjVr3I8dOnSIxMRE9/3ExEQOHjzYHqXLaZzr8ZSOpTWP5/r163E6nSQnJ7dx1XI6rXEsjx49ytChQ0lKSuLuu+8mKiqqnaqX72uN4zl37lwWLFhwVvtXGGplFRUVDB06lIyMjNM+/vrrrzN79mweeughNmzYwNChQ5kwYQIFBQXtXKk0h46nZ2mt41lUVMSNN97IP//5z/YoW06jNY5lREQEmzZtIjs7m3/961/k5+e3V/nyPed6PN955x369u1L3759z64AQ9oMYLz99tsN1qWlpRmzZs1y33c6nUZCQoKxYMECwzAMY9WqVcbVV1/tfvzOO+80Xn311XapV5p2Nsez3vLly42pU6e2R5nSTGd7PKuqqoyLL77YeOmll9qrVDmDc/ndrPerX/3KePPNN9uyTGmmszme8+bNM5KSkozu3bsbXbt2NcLCwoz58+c3e59qGWpHNTU1rF+/nvHjx7vX2e12xo8fz1dffQVAWloaW7Zs4eDBg5SXl/PRRx8xYcIEq0qWJjTneErn0ZzjaRgGM2fO5NJLL+WGG26wqlQ5g+Ycy/z8fMrKygAoKSlh5cqV9OvXz5J6pWnNOZ4LFiwgJyeHvXv38vjjj3Prrbfy4IMPNnsfPq1etTSqsLAQp9NJbGxsg/WxsbFs374dAB8fH5544gnGjRuHy+Xinnvu0dUNHVRzjifA+PHj2bRpExUVFSQlJfHmm28yevTo9i5XzqA5x3PVqlW8/vrrDBkyxN2n4eWXX2bw4MHtXa40oTnHct++fdx2223ujtN33HGHjmMH1dy/tedCYagDuvLKK7nyyiutLkNayWeffWZ1CdJKLrroIlwul9VlSCtIS0sjKyvL6jKkDcycObPFz9FpsnYUFRWFw+E4pZNefn4+cXFxFlUlZ0vH07PoeHoOHUvP0h7HU2GoHfn5+TF8+HCWLVvmXudyuVi2bJlOm3RCOp6eRcfTc+hYepb2OJ46TdbKysvL2bVrl/t+dnY2WVlZREZG0q1bN2bPns2MGTMYMWIEaWlpLFy4kIqKCm666SYLq5bG6Hh6Fh1Pz6Fj6VksP57nfhGcnGz58uUGcMoyY8YM9zZ//etfjW7duhl+fn5GWlqasXr1ausKlibpeHoWHU/PoWPpWaw+npqbTERERLya+gyJiIiIV1MYEhEREa+mMCQiIiJeTWFIREREvJrCkIiIiHg1hSERERHxagpDIiIi4tUUhkRERMSrKQyJiIiIV1MYEhEREa+mMCQiIiJeTWFIRLzO3r17sdlspyzp6elWlyYiFvCxugARkfaWnJxMbm6u+35eXh7jx4/nkksusbAqEbGKZq0XEa9WVVVFeno60dHRvPPOO9jtajAX8TZqGRIRr3bzzTdTVlbGp59+qiAk4qUUhkTEa/3xj3/k448/Zs2aNYSGhlpdjohYRKfJRMQr/fe//+Xaa6/lo48+4gc/+IHV5YiIhRSGRMTrbNmyhVGjRjF79mxmzZrlXu/n50dkZKSFlYmIFRSGRMTrLF68mJtuuumU9WPHjiUzM7P9CxIRSykMiYiIiFfTpRMiIiLi1RSGRERExKspDImIiIhXUxgSERERr6YwJCIiIl5NYUhERES8msKQiIiIeDWFIREREfFqCkMiIiLi1RSGRERExKspDImIiIhXUxgSERERr/b/A8s9sfmbzxH8AAAAAElFTkSuQmCC", - "text/plain": [ - "
" - ] - }, - "metadata": {}, - "output_type": "display_data" - } - ], + "outputs": [], "source": [ "# plotting xe\n", "fig = plt.figure()\n", @@ -136,13 +73,6 @@ "ax.legend()" ] }, - { - "cell_type": "code", - "execution_count": null, - "metadata": {}, - "outputs": [], - "source": [] - }, { "cell_type": "code", "execution_count": null,