{ "cells": [ { "cell_type": "markdown", "id": "95a11fb0-6c93-4c91-8434-69f6c2a76edf", "metadata": {}, "source": [ "# Simple Sample" ] }, { "cell_type": "code", "execution_count": 1, "id": "356a31e0-59f8-4cb1-91bb-2a5b50ec3723", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAkAAAAHgCAYAAABNbtJFAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/GU6VOAAAACXBIWXMAAA9hAAAPYQGoP6dpAAB7Y0lEQVR4nO3dd3wT5eMH8E+6d5ltKZQWCpS9oTK+glCpbPQryJA9HIwvgnwVGUUEcQJuUBl+RaSiKCiKDJmCgAwFQZbMMspsaYGWNs/vj/wuuaRJm2vTJHf5vF+vvJpeLpfnkhufe57n7nRCCAEiIiIiD+Ll6gIQERERORsDEBEREXkcBiAiIiLyOAxARERE5HEYgIiIiMjjMAARERGRx2EAIiIiIo/DAEREREQehwGIiIiIPA4DEHmsLVu2QKfT4euvv3Z1UYiK7c0330T16tXh7e2Nxo0bAwDi4uIwZMgQl5aLyN0xAJHdli5dCp1OB51Ohx07dhR4XQiBmJgY6HQ6dOvWzQUlJMn333+Pdu3aISIiAkFBQahevTr69OmDdevWubpoOHLkCGbMmIEzZ864uiiqt379evz3v/9FmzZtsGTJErz66qul8jnLly/H/PnzS2XazsRlj+QYgEixgIAALF++vMDwrVu34sKFC/D393dBqUjy1ltvoUePHtDpdJg8eTLmzZuHf//73zhx4gRWrFjh6uLhyJEjePnll7kTcoBffvkFXl5eWLRoEQYNGoQuXboAAI4dO4ZPPvnEYZ+jpQDEZY8kPq4uAKlPly5dsHLlSrz77rvw8TEtQsuXL0ezZs1w7do1F5ZO+4QQuHfvHgIDAwu8lpeXh1deeQUPP/ww1q9fX+D19PR0ZxTRYQqbVzL8noGBgfDz8zMbbs9BSHZ2NoKDg0uraKrHZU/7WANEivXr1w/Xr1/Hhg0bjMNyc3Px9ddfo3///lbfk52djYkTJyImJgb+/v5ISEjAW2+9BSFEgXGXLVuGZs2aITAwEOXKlUPfvn1x/vx5s3FOnDiBf//734iKikJAQACqVKmCvn37IiMjwzjOhg0b0LZtW5QpUwYhISFISEjASy+9VODz8vPz8dJLLyEqKgrBwcHo0aNHgc/bvn07evfujapVq8Lf3x8xMTF47rnncPfuXbPx2rdvj/bt2xf4jCFDhiAuLs5smF6vx/z581GvXj0EBAQgMjISTz31FG7evGk2XlxcHLp164aff/4ZzZs3R2BgIBYuXGj1e7527RoyMzPRpk0bq69HREQYn0t9oFJTU4ucfwDYvXs3HnnkEYSHhyMoKAjt2rXDr7/+WmC8tLQ0DB8+HNHR0fD390e1atXwzDPPIDc3F0uXLkXv3r0BAA899JCxSXXLli2FzuuZM2eg0+mwdOnSAp+n0+kwY8YM4/8zZsyATqfD8ePH8eSTTyI8PBwVK1bEtGnTIITA+fPn0bNnT4SFhSEqKgpvv/12gWmeO3cOf//9t9XvUE76Dr/66iu8/PLLqFy5MkJDQ/H4448jIyMDOTk5GD9+PCIiIhASEoKhQ4ciJyfH+P527dqhUaNGVqedkJCA5ORkm5+t0+mwZMkSZGdnG79H6fux7AMkNV9v3boVzz77LCIiIlClShUAwO3btzF+/HjExcXB398fERERePjhh7F//34AhmV67dq1OHv2rPFzLJdlS0Wte1z2bC975DysASLF4uLi0KpVK3z55Zfo3LkzAOCnn35CRkYG+vbti3fffddsfCEEevTogc2bN2P48OFo3Lgxfv75Z0yaNAlpaWmYN2+ecdzZs2dj2rRp6NOnD0aMGIGrV6/ivffew4MPPogDBw6gTJkyyM3NRXJyMnJycjB27FhERUUhLS0NP/zwA27duoXw8HD89ddf6NatGxo2bIiZM2fC398fJ0+etLrRnD17NnQ6HV544QWkp6dj/vz5SEpKwsGDB41HfytXrsSdO3fwzDPPoHz58tizZw/ee+89XLhwAStXrizW9/jUU09h6dKlGDp0KMaNG4fTp0/j/fffx4EDB/Drr7/C19fXOO6xY8fQr18/PPXUUxg5ciQSEhKsTjMiIgKBgYH4/vvvMXbsWJQrV67Ictgz/7/88gs6d+6MZs2aISUlBV5eXliyZAk6dOiA7du3o2XLlgCAixcvomXLlrh16xZGjRqF2rVrIy0tDV9//TXu3LmDBx98EOPGjcO7776Ll156CXXq1AEA418l81qUJ554AnXq1MFrr72GtWvXYtasWShXrhwWLlyIDh064PXXX8cXX3yB559/Hi1atMCDDz5ofO+gQYOwdetWqwHdmjlz5iAwMBAvvvgiTp48iffeew++vr7w8vLCzZs3MWPGDPz2229YunQpqlWrhunTpwMABg4ciJEjR+Lw4cOoX7++cXp79+7F8ePHMXXqVJuf+fnnn+Pjjz/Gnj178OmnnwIAWrduXWg5n332WVSsWBHTp09HdnY2AODpp5/G119/jTFjxqBu3bq4fv06duzYgaNHj6Jp06aYMmUKMjIycOHCBeO6GhISYvMzHL3uedqyR04kiOy0ZMkSAUDs3btXvP/++yI0NFTcuXNHCCFE7969xUMPPSSEECI2NlZ07drV+L7vvvtOABCzZs0ym97jjz8udDqdOHnypBBCiDNnzghvb28xe/Zss/EOHTokfHx8jMMPHDggAIiVK1faLOu8efMEAHH16lWb42zevFkAEJUrVxaZmZnG4V999ZUAIN555x3jMGk+5ebMmSN0Op04e/ascVi7du1Eu3btCow7ePBgERsba/x/+/btAoD44osvzMZbt25dgeGxsbECgFi3bp3NeZGbPn26ACCCg4NF586dxezZs8W+ffuKPf96vV7UrFlTJCcnC71ebxzvzp07olq1auLhhx82Dhs0aJDw8vISe/fuLfB50ntXrlwpAIjNmzcXGMfWvJ4+fVoAEEuWLCnwHgAiJSXF+H9KSooAIEaNGmUclpeXJ6pUqSJ0Op147bXXjMNv3rwpAgMDxeDBg82m2a5dO2HP5lH6DuvXry9yc3ONw/v16yd0Op3o3Lmz2fitWrUyWw5u3bolAgICxAsvvGA23rhx40RwcLDIysoq9PMHDx4sgoODCwyPjY01mydp3W3btq3Iy8szGzc8PFyMHj260M/p2rWrWbkL48h1zxOXPXIeNoFRsfTp0wd3797FDz/8gNu3b+OHH36w2fz1448/wtvbG+PGjTMbPnHiRAgh8NNPPwEAVq1aBb1ejz59+uDatWvGR1RUFGrWrInNmzcDAMLDwwEAP//8M+7cuWP1M8uUKQMAWL16NfR6faHzMmjQIISGhhr/f/zxx1GpUiX8+OOPxmHyfgDZ2dm4du0aWrduDSEEDhw4UOj0rVm5ciXCw8Px8MMPm81rs2bNEBISYpxXSbVq1QptDpF7+eWXsXz5cjRp0gQ///wzpkyZgmbNmqFp06Y4evRogfGLmv+DBw/ixIkT6N+/P65fv24sa3Z2Njp27Iht27ZBr9dDr9fju+++Q/fu3dG8efMCn6PT6ewqv5J5LcyIESOMz729vdG8eXMIITB8+HDj8DJlyiAhIQH//POP2Xu3bNlid+0PYPgO5TV2iYmJEEJg2LBhZuMlJibi/PnzyMvLA2BYlnv27Ikvv/zS+Hn5+flITU1Fr169HN5HZ+TIkfD29jYbVqZMGezevRsXL150yGc4ct3zxGWPnIcBiIqlYsWKSEpKwvLly7Fq1Srk5+fj8ccftzru2bNnER0dbbahA0xVz2fPngVg6NcjhEDNmjVRsWJFs8fRo0eNHXirVauGCRMm4NNPP0WFChWQnJyMDz74wKz/zxNPPIE2bdpgxIgRiIyMRN++ffHVV19Z3SDXrFnT7H+dTocaNWqYnSly7tw5DBkyBOXKlUNISAgqVqyIdu3aAYDZ59rrxIkTyMjIQERERIF5zcrKKtBZuVq1aoqm369fP2zfvh03b97E+vXr0b9/fxw4cADdu3fHvXv3zMYtav5PnDgBABg8eHCBsn766afIyclBRkYGrl69iszMTLOmnOJQOq+2VK1a1ez/8PBwBAQEoEKFCgWGW/a7csRnAUBMTEyB4Xq93myZGTRoEM6dO4ft27cDADZu3IgrV65g4MCBJSqTNda+2zfeeAOHDx9GTEwMWrZsiRkzZpRop+zIdY/LHpUm9gGiYuvfvz9GjhyJy5cvo3PnzsYjv+LS6/XQ6XT46aefChylAub9Dt5++20MGTIEq1evxvr16zFu3DjMmTMHv/32G6pUqYLAwEBs27YNmzdvxtq1a7Fu3TqkpqaiQ4cOWL9+vdXp25Kfn4+HH34YN27cwAsvvIDatWsjODgYaWlpGDJkiNmGXafTWa05yM/PLzCvERER+OKLL6x+ZsWKFc3+L+6ZKGFhYXj44Yfx8MMPw9fXF5999hl2795tDG/2kObvzTffNF5oz1JISAhu3LhRrDJasjavto7gLb9XOWu/sa3fXUltj72fZe/nJScnIzIyEsuWLcODDz6IZcuWISoqCklJSSUqkzXWvts+ffrgX//6F7799lusX78eb775Jl5//XWsWrXK2MdP6Wc4at3jskeliQGIiu3RRx/FU089hd9++w2pqak2x4uNjcXGjRtx+/Zts1og6Syb2NhYAEB8fDyEEKhWrRpq1apV5Oc3aNAADRo0wNSpU7Fz5060adMGCxYswKxZswAAXl5e6NixIzp27Ii5c+fi1VdfxZQpU7B582aznYt0lCkRQuDkyZNo2LAhAODQoUM4fvw4PvvsMwwaNMg4nvwsOEnZsmWtHj1LtVyS+Ph4bNy4EW3atHHaabbNmzfHZ599hkuXLpkNL2r+4+PjARjCVGE75YoVKyIsLAyHDx8utBz2NkfIlS1bFgBw69Yts+GW36saeXt7o3///li6dClef/11fPfdd1abqkpTpUqV8Oyzz+LZZ59Feno6mjZtitmzZxsDkNLfzFHrHpc9Kk1sAqNiCwkJwUcffYQZM2age/fuNsfr0qUL8vPz8f7775sNnzdvHnQ6nXEj+9hjj8Hb2xsvv/xygaMiIQSuX78OAMjMzDT2oZA0aNAAXl5exlOMrR0RSkeQ8tOQAeB///sfbt++bfz/66+/xqVLl4zlknZE8jIJIfDOO+8U+Iz4+Hj8/fffuHr1qnHYH3/8UeAMmD59+iA/Px+vvPJKgWnk5eUV2Nja686dO9i1a5fV16S+VpZnthQ1/82aNUN8fDzeeustZGVlFZiuNK9eXl7o1asXvv/+e/z+++8FxpO+P6lfi5J5DAsLQ4UKFbBt2zaz4R9++KHd01DC3tPgHWXgwIG4efMmnnrqKWRlZeHJJ590yufm5+cXaMKNiIhAdHS02XoSHBxsd1OvI9c9T1z2yHlYA0QlMnjw4CLH6d69Ox566CFMmTIFZ86cQaNGjbB+/XqsXr0a48ePNx7lxcfHY9asWZg8eTLOnDmDXr16ITQ0FKdPn8a3336LUaNG4fnnn8cvv/yCMWPGoHfv3qhVqxby8vLw+eefw9vbG//+978BADNnzsS2bdvQtWtXxMbGIj09HR9++CGqVKmCtm3bmpWvXLlyaNu2LYYOHYorV65g/vz5qFGjBkaOHAkAqF27NuLj4/H8888jLS0NYWFh+Oabb6y23Q8bNgxz585FcnIyhg8fjvT0dCxYsAD16tVDZmamcbx27drhqaeewpw5c3Dw4EF06tQJvr6+OHHiBFauXIl33nnHZp+qwty5cwetW7fGAw88gEceeQQxMTG4desWvvvuO2zfvh29evVCkyZNFM2/l5cXPv30U3Tu3Bn16tXD0KFDUblyZaSlpWHz5s0ICwvD999/DwB49dVXsX79erRr1w6jRo1CnTp1cOnSJaxcuRI7duxAmTJl0LhxY3h7e+P1119HRkYG/P390aFDB7NrFFkzYsQIvPbaaxgxYgSaN2+Obdu24fjx44q/I3soPQ2+pJo0aYL69etj5cqVqFOnDpo2beqUz719+zaqVKmCxx9/HI0aNUJISAg2btyIvXv3ml2jplmzZkhNTcWECRPQokULhISE2DzoceS654nLHjmRE884I5WTnwZfGMvT4IUQ4vbt2+K5554T0dHRwtfXV9SsWVO8+eabZqe2Sr755hvRtm1bERwcLIKDg0Xt2rXF6NGjxbFjx4QQQvzzzz9i2LBhIj4+XgQEBIhy5cqJhx56SGzcuNE4jU2bNomePXuK6Oho4efnJ6Kjo0W/fv3E8ePHjeNIp+J++eWXYvLkySIiIkIEBgaKrl27mp3aLoQQR44cEUlJSSIkJERUqFBBjBw5Uvzxxx9WT49dtmyZqF69uvDz8xONGzcWP//8c4HT4CUff/yxaNasmQgMDBShoaGiQYMG4r///a+4ePFiod+nLffv3xeffPKJ6NWrl4iNjRX+/v4iKChINGnSRLz55psiJyenWPMvhOHyA4899pgoX7688Pf3F7GxsaJPnz5i06ZNZuOdPXtWDBo0SFSsWFH4+/uL6tWri9GjR5t99ieffCKqV68uvL29zU5LLmxe79y5I4YPHy7Cw8NFaGio6NOnj0hPT7d5KrLladi2Thlv166dqFevXoFh9mwepe/Q8pIMttYVW2UTQog33nhDABCvvvpqkZ8rUXoavGV5cnJyxKRJk0SjRo1EaGioCA4OFo0aNRIffvih2XhZWVmif//+okyZMgJAoafEO3rdE8Kzlj1yHp0Q7IFF5Im2bNmChx56CCtXrixWbRM51jvvvIPnnnsOZ86cKXAWkdZw2SN3wD5AREQuJoTAokWL0K5dO82HHyJ3wT5AREQukp2djTVr1mDz5s04dOgQVq9e7eoiEXkMBiAiIhe5evUq+vfvjzJlyuCll15Cjx49XF0kIo/BPkBERETkcdgHiIiIiDwOAxARERF5HAYgIiIi8jgMQERERORxGICIiIjI4zAAERERkcdhACIiIiKPwwBEREREHocBiIiIiDwOAxARERF5HAYgIiIi8jgMQERERORxGICIiIjI4zAAERERkcdhACIiIiKPwwBEREREHocBiIiIiDwOAxARERF5HAYgIiIi8jgMQERERORxGICIiIjI4zAAERERkcdhACIiIiKPwwBEREREHocBiIiIiDwOAxARERF5HAYgIiIi8jg+ri6AO9Lr9bh48SJCQ0Oh0+lcXRwiIiKygxACt2/fRnR0NLy8Cq/jYQCy4uLFi4iJiXF1MYiIiKgYzp8/jypVqhQ6DgOQFaGhoQAMX2BYWJiLS0NERET2yMzMRExMjHE/XhgGICukZq+wsDAGICIiIpWxp/sKO0ETERGRx2EAIiIiIo/DAEREREQehwGIiIiIPA4DEBEREXkcBiAiIiLyOAxARERE5HEYgIiIiMjjMAARERGRx2EAIiIiIo/DAEREREQehwGIiIiIPA5vhupM33wDHD8O6HT2P7y8lI3viIeXF+DjA3h72/e3qNf8/Q3/ExFRQfn5wP37QF6e6W9hz6X/8/MBIQC93vAoyfPCXgMM/0sP+f8lff7ii6X//drAAORM//sfsGaNq0vhGn5+QFCQ4REYaP43MhKIiQEaNgSaNgVq12ZgIiJ1u3UL2LcP2LkTOHMGuHIFSE8HsrOBu3cNjzt3DH/v33d1aV2HAchDSInXE+XmGh63bhU9bmAgUL8+UKtWwUdISKkXlYjILvn5wD//AIcPA0eOGGr4T5wwPK5dc3XpqAg6ITx5r2xdZmYmwsPDkZGRgbCwMMdNePdu4PJl86rEoh5SVaQzH3q9YcXOy7P9t7DX5OPk5QH37pmOdqQjHumvEjodUK8e0KqV4dGxI1C1quN+H7KPEMC5c4aN/LlzwIULQGYmcPs2kJNjavoMDQUiIgyPqlWB6tWBKlVYuycRwrCTPHMGOHvW8F3euGE4SMjONjVJBwUB5csD5coBcXGGA4Hq1Q21quRct28DmzYBO3YYHn/8Ydi+2Uv6PaUacPnD39+07vj4AL6+RT/39jZ0WZC6SjjiubXXdDpT+a39X5LnPXuW7DexoGT/zQBkRakFIDInhGFDf/kycOoUcOAAsH+/4fHPP/bVmDVvDjz+ODBoEFCpUumX2VOdPAl89x2weTOwdy9w9WrxphMcDHTpAgwYAPToYdoYlpbMTGDtWuDXX4E//wQuXgRu3jS8FhAAVK4M1KgBtGkDdOoE1KxZuuUBDLUFn35q6BN44ULxpuHrC7RsCbRrByQnA23bGnZW5Hj37gFffw189RWwfr0h5BdFWq7q1DEsWw0aGLZP5cvzAKCUKdp/Cxd7//33RWxsrPD39xctW7YUu3fvLnT8efPmiVq1aomAgABRpUoVMX78eHH37l2zcS5cuCAGDBggypUrJwICAkT9+vXF3r177S5TRkaGACAyMjKKNU/kAHfvCnH4sBCrVgnx2mtCDBsmRNOmQnh7W6+38vUVYuBAIU6dcnXJrUtPF2LlSiEmTBCiRw8h6tUTonx5IUJDhfD3F6JMGSFq1xaifXshnn1WiMWLhTh/3rVlvn9fiC++EKJ169Kpa2zVSggF66UiZ84IMWSIEH5+ysr00ENC/PCDEHq948t09aoQI0YI4eXl+O+yalUhXnlFiBs3HF9uJbKzhVi3TogXXxSiWzch6tQRIipKiOBgw29RvrwQ1asL0a6dEKNHC7FggRDHjrm2zLZkZRm+04gI2997rVpC/PvfQkyfLkRqqhB//GF4H7mMkv23S2uAUlNTMWjQICxYsACJiYmYP38+Vq5ciWPHjiEiIqLA+MuXL8ewYcOwePFitG7dGsePH8eQIUPQt29fzJ07FwBw8+ZNNGnSBA899BCeeeYZVKxYESdOnEB8fDzi4+PtKhdrgNxYdjbw++/Atm3At98aao3k/PyASZOA6dNd30SQlwekpgJLlhhqTqSzKZRo1Ajo1QsYPtzQUdxZtmwBxo0DDh0q+Fr58obahwYNgNhYQ7nKljU0efn7G+bz/n1DU87Vq8ClS4Zmnn/+AbZvB65fN03LxwdYsMAwf44gBPDmm8C0aYY+Z5bKlTM8dDogK8vUJG3p0UeBRYsM8+UIf/4JdOsGnD9vGubra6i5qVXL8D3GxhqaC8PDDTVlgOG7vHPH0DSWnm6oKT1+3NC59sSJgp8TGgpMmQJMnGj4bp3h7l1DbdayZYbl3Nr3XpSaNYFhwwzLQcWKji+jUt9+a1j+LWvoKlUyLBudOhlqdipUcE35yCbV1AC1bNlSjB492vh/fn6+iI6OFnPmzLE6/ujRo0WHDh3Mhk2YMEG0adPG+P8LL7wg2rZtq6gc9+7dExkZGcbH+fPnWQOkFqdOCfHSS0KUK2d+ZNa6tRAXLrimTHq9EMuXG450bR05+vkJUa2a4Qi5YUMh4uOFCAmxPb6/v6EmLD+/9Ms+Y0bBz69fX4jXXxfi6NGS1Y7cvy/Ed98Zarvk01+0qORlz8831LDIp1u2rBBjxwqxYYP12pGcHCF+/93w3daqZf7eunWFuHy55OU6etRQ8yFNNyzM8F3evFmy6aalCfH550J06VKwVqljR+fUBm3aJERsrO3lNjBQiJgYw3fZqJFhmS9btvDxp04V4vbt0i+7Nbm5Qowfb14mb28h+vYVYtu20l//qMSU1AC5LADl5OQIb29v8e2335oNHzRokOjRo4fV93zxxRciPDzc2Ex26tQpUbt2bTF79mzjOHXq1BHjx48Xjz/+uKhYsaJo3Lix+PjjjwstS0pKigBQ4MEApCIZGUK88IKhKUzacEVHC/H3384tx/nzQjzySMENe3y8EJMmGZpXzp+3vSG9eVOIzZsNIaRFi4LT6dfPECJKg14vxMSJ5p/XrJkQv/zi+Cah3Fwhxo0zb8IsaXPYtGnmZZ8wQVnIyM8XYsUK8zDdsqUhJBVXdrZh5y+fXlpa8adny6lTQowcKYROZ/7b3brl+M+SLF1asEm6alUhRo0S4quvhDh3zvZyfu2aEFu2CPHGG0I8+KB5uQEhKlc2hFZnys01NE/Ly9G5s/O3IVQiqghAaWlpAoDYuXOn2fBJkyaJli1b2nzfO++8I3x9fYWPj48AIJ5++mmz1/39/YW/v7+YPHmy2L9/v1i4cKEICAgQS5cutTlN1gBpyJ49ho2wtAGrUkWIixed89l//ilEpUrmG9CkJMOGvrgB4sIFIZ57znwHMXKkY8st+egj02fodIZaitI84tXrDbUz0mc2aCBEXl7xprVtm+k78vIyBJniOnnSUGshleuFF4o/rcmTzeevNAOJEIbvoWJF02d27Vo6v+G335ovkx06CLF9e/GX85MnDcuCj49pml5eQnzyiUOLbVN+vqEPoTyQf/hh6fQFo1Kl2QC0efNmERkZKT755BPx559/ilWrVomYmBgxc+ZM4zi+vr6iVatWZu8bO3aseOCBB+wuGztBq1x6uqG6XdqYtWpVerUmkrNnzcNP5cpCrF7tuOmvXm1eu7VkieOmLYShmUbeYbiIWlOHyc0VokkT0+d++qnyaeTlCdG4sWkar71W8nLt22dodpR2hsePK5/G2bOmafj5Gb5jZzh82LwW6913HTv9U6cMzXjS9MeMKX5wtXT8uOGgQR7Ev/7aMdMuzKxZ5s3TGzeW/mdSqVBFACpOE1jbtm3F888/bzbs888/F4GBgSL//49yqlatKoYPH242zocffiiio6PtLhsDkAZcvmzeN2HWrNL7rPv3hWjTxvRZLVoYzvhxtGXLTJ8RHi7ElSuOma5ebzj7TJr2uHGOma69tm0zfXadOsqPuletMr2/SRPH7YynTjVNd/Bg5e//739N7580yTFlstf69abPDg52XLObXi9EcrJp2k884fhakvx8Q62n9BmhoaV7due+faamPC8vw/JEqqVk/+2yC0f4+fmhWbNm2LRpk3GYXq/Hpk2b0KpVK6vvuXPnDrwsrnXh/f/XVBBCAADatGmDY8eOmY1z/PhxxMbGOrL45O4iI4EVK0zX3Jg923ChudLw/vuG68wAhgvV/fhj6ZwdMmCA4XpHAJCRAaSkOGa6GzYYzvoCDBfYmzPHMdO117/+ZTgbCgCOHjWVxV7vv296Pnu2466z8sILQJkyhucrVii7su+9e4azyADD2V4TJzqmTPZ6+GHg6acNz7OzHfebrlsH/Pyz4XlMDPDJJ46/lpOXF/D224blHTBcfHDECEMccjQhgDFjDBdvBYCpUw1neZFnKP08ZtuKFSuEv7+/WLp0qThy5IgYNWqUKFOmjLj8/2deDBw4ULz44ovG8VNSUkRoaKj48ssvxT///CPWr18v4uPjRZ8+fYzj7NmzR/j4+IjZs2eLEydOiC+++EIEBQWJZcuW2V0u1gBpiPxIcsgQx0//2jXDNXykz9i+3fGfIXf5suGIWKqqd0T/pgcfNJW/JH1nSuLLL01lGDHC/vedOmV6X82aju/vMmGCafrz59v/vtWrTe/r39+xZbJXerqh9kdaVhxRY/ivf5nmKzW15NMrTGamEHFxps+zaC1wiG++MU2/dm1DkyypmiqawCTvvfeeqFq1qvDz8xMtW7YUv/32m/G1du3aicGyquf79++LGTNmiPj4eBEQECBiYmLEs88+K25anOnx/fffi/r16wt/f39Ru3btIs8Cs8QApCE3b5oCio+PoV+GI02fbtqADh3q2GnbMmmS6TOnTSvZtI4eNU0rIcF1p/lmZwsREGAoR6VK9pdj/nxT+V991fHlOnzYNP2OHe1/38iRpvetWeP4ctnr+edN5XjjjZJNa88e07Tq1nXOsmLZvOnI5ja93vxMy++/d9y0yWVUFYDcEQOQxsj7cpTkjB5L2dmmzqalEa5suXDBdN2XuLiS7RTkO8i5cx1XxuLo1s1UFntPiZd3mD182PFl0utN13Py8bHvLK78fFOH+IAAw3LiKidOmL6fWrVKtqzIz9hzVid5vV6I5s1Nn+vIU+N37jRNt1EjnvGlEaroA0TkNGPHmq6K+/nnpvb+klq1ynCFXgDo29d5N2atXNnQxwMwXGFZ6n+klBCGK/gChu9H6l/kKt26mZ5v3lz0+HfvAlu3Gp7HxQF16zq+TDqdqVx5eYYbYRbl+HHD1a8BoEMHw40vXaVGDcP9wqRyHTlSvOncv2/oBwUY7qHWp49jylcUnQ74739N/3/yieOmvWSJ6fn48aV/XzpyOwxApH0REYYbcAKGm2HasxOzx//+Z3r+1FOOmaa9pA6iALB6dfGmceQIcPq04fmDDxpuceFKrVubnu/dW/T4Bw4YdsyAIWiU1g6sQwfT8z17ih5fXvY2bRxfHqV69TI9//774k3jt99MN8Dt1s1wuw5n6dXLdHuM1asNt1gpqXv3DLepAYCQEKB375JPk1SHAYg8w8CBpufFDQxy6emmIFWtmvN3dJ07m3b4GzYUbxo//GB63r17yctUUnXqmGpL7AlA8nFatCidMllO253KZS/5byv/zZWQL2Pymjpn8PUF+vUzPM/JMdynq6S2bgUyMw3PH33UdO818igMQOQZkpMNG1IAWLu25KfUrl9vurlp797Orz6vUAFo2tTw/I8/gCtXlE9j2zbT886dHVOukvDxMc3TmTOmGgdb5EGjZctSKxaiow0PwHAj3qJuaisvV/PmpVcue8XHG264ChhqsO7dUz4NeQCSml+d6YknTM9/+qnk01u71vS8Z8+ST49UiQGIPENoqKkvxNmzxe8LIZFvhKXmNWeT74ikvjD2EsLQrAEYwpS0g3Q1eY3J/v2FjysFDX9/w53pS5NUrsxM4ORJ2+Pp9cDBg4bn8fGOu5t8SUk1lPfvG+4kr0R2tqnpr25dUxh0ppYtTd/lhg2G/lgl8eOPhr++vq4JdOQWGIDIczzyiOn59u3Fn44QhhogAAgLM++74kzyz1W6Uzt+3NSB+4EH3KcDqLwjc2FB4/594NQp03uk2r3S0rix6fnff9seLy3NVMNSGp2yi0t+cdldu5S99+BBU62Xq/o0+fiYgsqtW/Y1RdqSlmZadlq1MqzD5JEYgMhzyDfeO3cWfzonTpiuCvzgg6W/87WlWTPTc6UBSL4TdFWAsyY+3vRc2klZc/as6Wy+GjVKt0yA/eWSv+aMctmrJAFIvmzJlzlnS0oyPS/J+iuff+kK5OSRGIDIczRpYmguAUq2Ad292/T8gQdKVqaSqFTJcMsPwNBcpKRf06FDpufu0E9FIg8NhQUNee2QPJyUFvlnFFYzJX/NnQJQ3bqmDuby394e8gAk9dFyhcRE03P5OqiUfN13p/BPTscARJ7D39+0sz91ynAmV3FIfWcA1wYgnc50RH7zpqHjsL3kfaDcqammcmVTSC0saDi7psVdg5m9vLyAhATD81OnlHWElgKQj0/p97UqTL16prO15OugUvIaIFeuv+RyDEDkWeSdbJUeCUukDqE6netPc27UyPT86FH73yeNGxbmmk6ttnh5GW7ICgD//GP7jCtnB42KFQ3Xi7H8bEvu2gQGmIKuXm9oxrVHXp6hvxhgCFABAaVTNnt4e5vWt/PnTRebVEKvB/780/C8enXXX/uKXIoBiDxL/fqm54cPK39/fr7pfTVquL4DZc2apueF7ZjlsrIMfWgAw7V33KUDtEQKDvfuGS5caY2zg4ZOZ/qcs2dNF2C0JP0G3t7OuzK4verUMT239yzI8+dN8+oOZwrK+yAVZ/09cwa4c8fw3JW1WeQWGIDIs8gD0F9/KX//2bOm5oN69RxTppKQ7/ztDUDHjpmeu1Pzl0QeHGwFIKm5z8/PeTVYUs1UXp7tcknDo6Nd1zneFvlvbW9tobv1aZLPQ3EuZSFf5+XbAvJIDEDkWeQb0OIcQcp3HPIjalcpTgCSmjQAoHZtx5bHEaSO3YDtCzxKwyMjDc1mzhAVVfDz5fLzTWcHyufBXciXV/kyUBh3DkBKmnwl8nXeHQ5gyKUYgMizhIYCsbGG54cPK78itLt1Ho6KMnUMtTcAnT9vei59F+6kqAAkDxoREc4pE2BeLmsd6K9dM/VZcscAJK9Zky8DhZEvU/LmVlcpTjOenDwAsQbI4zEAkeeRaj1u3y76dguW3K0GSN435fRp+66QK9/5xcSUTrlKoqgAdP26a4KGPGxZK5d8mDsGoKAgU6ffc+fse4+71QCFh5uaPItTAyTNj07nHn2ayKUYgMjzVKtmeq7k1HHAvOlAOq3Y1eR9U9LSih5f7QHIVUHDXculhPR7p6WZLiRZGHlfq8qVS61YikgHHteuma5mbi9pfqKjTZdbII/FAESeJy7O9FxpAJKOnOWnRbuavBOwPacGSwHIy8twMUV3465Bw13LpYQUgPLzgcuXix5f6tQdFeW8vlZFka+/9tZkAYazv6SmS/k0yGO5yRJN5ETFDUD375tqWNzpFGd5iLFnpyYFoOhow8Xt3E1RTU3y/jfu1AdITQEIKLof0P37pr5W7hSU5euekgAkXfoBYAAiAAxA5ImKG4AuXjT1PXHXAFRUDdC9e6Z+T+7Y/AUYOnVLHbvdqaZFSzVAQNEBSD4/DECkQQxA5HnkZz4pCUDyja07nT2lJABduGB67q4BCDAFCHeqaQkOBgIDC5ZBIi+rFgKQfFnSQgCSr+vutP6SyzAAkeeJjDRd0r+4AUitNUDy192lU6s1UoC4caPgVZfl4cOZTWA6nalcaq0Bkv/mSpYVdw1A8lqdosjXddYAERiAyBPpdKYjYXmNSFG0EICuXzc9r1ixdMrjCOXKmZ5nZJi/5sp5kD7vxo2C9ym7dcv0vEwZZ5VIGfn3Jf8erXHXAFSlium5khogtdR+ktMwAJFnko7Qb98G7t617z3uGoAqVDDcewooOgBJnVoB974RpPwea5YBSP6/s4NGeLjhrxCGe6rJSeUKCTH9Hu5G/pvLlwVr3DUABQSY1l8lAUheQye/qjd5LAYg8kxFndFjjfwMK3e6g7q3t6kpSEkNUIUKpVemkpKCBgBkZpq/Jv3v6+v8a7nYUy75OO5GHoCU1AC5W2CQ1r8rVwrWxNkired+fu79G5HTMACRZ7LnflOW5FeNdrfwIO2g0tMLv72HFmqApKARFub8O9nLy2UZgKRyysdxN76+pvIVVQMkPzBwtwAkNeXl55s3PRZGfv84Zy835JYYgMgzyTvP2lsDJAWg8HD3u4qs1GcmP79g04ycWmqACgsa8gDkbLbKJf/e3TkAAabfvagaoJs3Tc/lfbLcgXz9ted2Nnq9aTxndpwnt8YARJ6pODVAUlByx87DZcuanst3XJbUUgNkq6lJCFNNiyuaMeSfKa+ZkodOd29ekX73GzcKvx2GtBz5+5tO/3cX8nXQngDkqvvHkVtjACLPVNTVhi3dv2+qanfHACTvDFxYk4BaApCtJrCcHNNp8e5UAyQvo1pqgIQofFmRXnPHM9qUBiBXXTqB3BoDEHkmpZ2g5cHBHQOQvTVAUrNHWJihP4i7shU05M/dKQC5ulxK2HsmmLQcyZctdyFfB+1Zf9VwkUpyOgYg8kxKm8DkR5nueAQpP0q3pwnMnfv/ALabmuRBwxVNTfbUALl7E5j8t7fVDygvz3CJCMA9a4CU9gFSw0UqyekYgMgzKe0ELR/H3WuAbDVr5OebwpG7ByB3bWqyJ5i5ew2Q/Le3VQMknzd3rwGyJwC56ga65NYYgMgzhYaazuSyZwMqH8fdA5CtGqCbN02nyLtz/x/AdidoVwcNe5rA3L0GyJ5rAcmXIS0EIHdvwiaXYAAiz6TTmXZUlteZscbdA5A9naDVtJO21Qna1fPgrjVTShR2MUeJu9/WQ2kAcvf5IZdgACLPJW0I7QlA7n4EaU8NkNSnAzDUgLkzd+1s7K7lUkL+28uXCTl3rwGSd+JnAKJiYgAizyUdCWdmFn05fVfef8oe9nSCVlMACggw3LIAMP/uXV3TYqsPkJo6QSsNQO64vOt0pnLZqsWSU9PvQ07DAESeS9qACmF7RyBx9Y63KPZ0gpZfrC8kpFSL4xDS9+xOfW3k35taa4Dk82DrquHyZcgda4AA0/dsTw2ufH4YgOj/MQCR55JvCIu6n5Crd7xF0VoNEFB0AHJF0PD2NgUIW+Vy9+9WC01ggHkNbmH3vwNMIcnf31C7SAQGIPJk8tBQ1FGku9cA+fkBQUGG51oJQNL83LljGuYOtVjS58q/T3kZg4OdWx6l7AlAaugzI62HeXnA3buFj+vOV7Uml2EAIs9l7+0jANMRvpeX++7gpB2CrZ2aPDyoIQBJ95+6e9d0hC8PGlJAcjbpc+U7Xflzd7tvliV5cLS1rLi6ps0e9pzNJpEOYBiASIYBiDyXrQ6t1kivh4UZOmC6IymYyUOCnHxnp4Y+QPIgkZNj+OsOQUMezCTuUC57yQO8rQCkhhotW5dKsKTXmwKSOzZfk8swAJHnKk4NkDtvQKUdVXa29dfV2gQGmHbI7hA0iqqZcvcA5OVlCsC2OkHLlyF3DUD2HsDcvm36nVgDRDIMQOS5ilsD5K6kwHDvnuG2F5bU2gQGmIKPOwUgvd50Z3qpXD4+7n2TWYm1fkxy8gDkqqbGoti6JpMlngFGNjAAkeeytwYoN9cQKgD33oDKj9StdQpVWw2QuwYgeSCwLJe71/5IpN9fzU1g9h7AuPs1vMhlGIDIcympQpeooQYIsN4MprY+QNaChrRj9vJyXU2LPORYNs25a22JJXkAsnYKubT86HTue9o4a4CohBiAyHPZWwOklqvIyo/UrXWEVnMNkLWg4arO6NZqpqTyqaUGSArA+fmmDuZy0vy48nsuCmuAqIQYgMhz2bsBVcMpwYB5ALJWA+QO19BRorAmMFcGDXctlxJFXQtIWn7ctfkLYA0QlRgDEHku+U7A1plTgPtfBFFibxNYcLChCcnduWvQcNdyKWFvAHLnJj17T4OXrwtqqPkkp1HBVpColBQVGCTufhsMib1NYGqo/QHct7OxZbnu3zeddefOgUGuqAAkLT/uXANk74UQ1XBGG7kEAxB5Ln9/U/8GWxcPBLRTAyQ1ganlKNhaHyB36GtjWS41XQNIUtgNUYUw7wPkruytAXKHq4eTW2IAIs+l01m/35QltVw/p6gaIGk+1FIDZNnUlJ9vuu6OK3dkluVyh1PzlSqsBujePdOZYe5cAySfB1sXdAQYgMgmBiDybEXdPsLyNXfegBbWCTo/33A9I8C950HOsqnJXYJGYQFILd+tfB6ka1xJ1HAVaMB6E6k1arimEbkEAxB5NmkjWlgfIHfZ8RalsCYw+anO7npdF0uWTU3u8jtY7njV2AQmXwYKC0DuHOjk81DYAYxa5oecjgGIPJs9TWBqrAGynB/5Tk6NAcidgoa7BjMlCgtAaqkx0ems35jWklrWX3I6BiDybPIAZO2KuIB6dnCF1QBpIQC5y+/gruVSQgs1QIDp+9bCAQw5nVsEoA8++ABxcXEICAhAYmIi9uzZU+j48+fPR0JCAgIDAxETE4PnnnsO9yxX4v/32muvQafTYfz48aVQclI96QhXrzf1kbGklg1oYX2A1BiACusD5K6doN15+ZCzNwC5cw0QoLwG193nh5zK5QEoNTUVEyZMQEpKCvbv349GjRohOTkZ6enpVsdfvnw5XnzxRaSkpODo0aNYtGgRUlNT8dJLLxUYd+/evVi4cCEaNmxY2rNBamXPtYDUcoSv5SYwd2pqctemOSW00AQG2NcEJl+v1fL7kFO4PADNnTsXI0eOxNChQ1G3bl0sWLAAQUFBWLx4sdXxd+7ciTZt2qB///6Ii4tDp06d0K9fvwK1RllZWRgwYAA++eQTlC1bttAy5OTkIDMz0+xBHkIegGwdRaqlBohNYM7hrmenKaGVJjAlNUB+foCPT+mXiVTDpQEoNzcX+/btQ1JSknGYl5cXkpKSsGvXLqvvad26Nfbt22cMPP/88w9+/PFHdOnSxWy80aNHo2vXrmbTtmXOnDkIDw83PmJiYkowV6Qq9gQgtezgtFYD5K5nWxVWM+XugUGitRqg3FzT1bgtqeGijuQSLo3D165dQ35+PiIjI82GR0ZG4u+//7b6nv79++PatWto27YthBDIy8vD008/bdYEtmLFCuzfvx979+61qxyTJ0/GhAkTjP9nZmYyBHkK1gC5LzaBlR6t1QABht/C2kU+GYDIBpc3gSm1ZcsWvPrqq/jwww+xf/9+rFq1CmvXrsUrr7wCADh//jz+85//4IsvvkCAnRt6f39/hIWFmT3IQxR1B3XAfMfrzuGhsDCn9gDkTk1N7louJeytAXL30GDPxRDVcGd7cgmX1gBVqFAB3t7euHLlitnwK1euICoqyup7pk2bhoEDB2LEiBEAgAYNGiA7OxujRo3ClClTsG/fPqSnp6Np06bG9+Tn52Pbtm14//33kZOTA29v79KbKVIXJTVAAQHufRd1f3/Tc8sz2tQYgPz8DN+3Xu9eZ1tpvQ+QfNmRL1PuyNr94iyxBohscOnW3M/PD82aNcOmTZuMw/R6PTZt2oRWrVpZfc+dO3fgZbETkgKNEAIdO3bEoUOHcPDgQeOjefPmGDBgAA4ePMjwQ+aUBCB334B6eZk6ecqv/AyoMwDpdKay3rvnPjVxljtd+XerhQAkX3bcPQAVVQOUn2+aH3dff8npXN4lfsKECRg8eDCaN2+Oli1bYv78+cjOzsbQoUMBAIMGDULlypUxZ84cAED37t0xd+5cNGnSBImJiTh58iSmTZuG7t27w9vbG6Ghoahfv77ZZwQHB6N8+fIFhhMp6gSthp2bnx+Ql6eNAAQYdsB37hjmx13mQR4KcnLUFRgk9gYgPz/nlKe4iqoBcpdaQ3JLLg9ATzzxBK5evYrp06fj8uXLaNy4MdatW2fsGH3u3DmzGp+pU6dCp9Nh6tSpSEtLQ8WKFdG9e3fMnj3bVbNAamZPHyC11AABpsCghSYwwBQo7t1zn/uZeXkBvr6GO9NbBjMtBCA1NYEVVQOkpos6ktO5PAABwJgxYzBmzBirr23ZssXsfx8fH6SkpCAlJcXu6VtOg8hIazVA0g6rsBogNcyHRNpRu1vQ8Pc3BCB3CmZKFHY3eDXVaBVVA6SmDt3kdG7co5PICYoKQHq9aQehhg2oPQFILTtpwHx+3CloyIOZmgKDRN60peYAVFQNEAMQFYIBiDxbUQFIbTUn0o5N601grt4xy4OZO9VM2cuyg7mcfNlRex8gBiAqBAMQebai7gWmtg2o1mqA3LkJDFBvExhgOwC5U9AsCmuAqAQYgMizyXdYlqEBUN81XrQWgOQ74Nu3Tc9dPQ9qbwIDtBGAiqoBUutyT07BAESezfKUZktqO4KUmizy883vjeQu19BRSv77ZGRYH+4Kam8CA0zLgWXNiZpOgy+qBkhNYY6cjgGIPJt8A2/ZbwZQbw0QYD4/aj0Slpc1M9P03NU7M+nzc3PVH4DUfBp8UTVAaurPRE7HAESeTWs1QLbmR60ByFYNkKvnwVow8/MzdC5WCy00gbEGiEqAAYg8m/yosKgApIYaIFvzo7UA5OqdmbVyubpMSknLQW6u4XIPEvly4+MWl4qzragaIDU155HTMQCRZyvsBqKA+s7w8ZQmMF9f55dFTv49S+VS0/cK2D4BQFpu/P3dv0arqJMY1NScR07HAESeragmMLX1IdByE5h0mYKAANfvmOXfoVQute1gbd0OQ1pu1DA/RfXhYxMYFYIBiDxbURtQrQUgb2/3b9aQs7bTcocdmbuWSwlPCEBqW3/JqRiAyLPJm1K0UANka4cg7eDUVPsDWC+vO+yYrZVBzd+tPABJy42al3cJa4CoEAxA5Nl0Otu3j7AcpoYdQlE1QGrbCbhr0HDXYKYEa4DIwzEAEdm6ejKgvg2orQCkpqN6OXdtanLXcikhr/2UL+daCkCsAaJCMAARSRtRLQQgW6fB379f8HU1cNeaFnetmVJCvixIywdgWm7UsKwoCUBqmB9yKgYgIvlVfS2pLQDZOg1e2sG5+vRxpdw1aLhrMFNCvixIy4cQ5qfBuzslTWBqmB9yKgYgInubwNSwAbXVBKalAOQOv4O7lksJawFIXhOkhvmx1YwnYRMYFYIBiMgTOkGrNQC5a02Lu9ZMKWEtAKktMGjtJAZyKgYgIi11grbVJKDWAOSuQcNdg5kS1gKQ2pZ3oPAApLZAR07FAESkpSNIazVAQmgrALnDjsxdy6WEteYjNQYGLa2/5FQMQETShl6vB/LyzF9T2wbUWgCSz5PaApC71rS4a82UEtbOAtNaAFLj/JDTMAARFXYmidoCkLV5kXdsVcM8yLlr0HDXYKaEJzSBqXF+yGkYgIgKuyGq2jag1uZFHoDUVgPkrk1N7louJbTQCRpgDRAVGwMQUWEBSG0XUtNaALJW0+KuNUDuUC4ltBaArJ3EIA3z8lLXTYDJKRiAiLTaBKaFAOSuNS3uWi4ligpAaljeAfuawNQyL+RUDEBEWm0Ck8ounwcGIMdw13IpYe0sMLVd+BMwD0BCmL+mpvuakdMxABFpqQaITWDOoYVO0Fo7CwywfRanGtZdcjoGICKt1gBpIQAFBxcc5g47ZmsBKCTE+eUoCa31AQIKHsCwBogKwQBEpKUaoKLu76SGeZCzFirCwpxfDkuhofYNc2dFLStqCcv2rL9qW+7JKRiAiOypAdLpAG9v55WpuLSyU5N4eRWsBXKHoKHVAKTGi2ayBoiKiQGIyJ4A5OdnCEHuTn6qr7QzU3MAAgrWArlD0AgMNIQzOS00gckDkFpOG2cAomJiACKypwpdLRtQ+U5N2pmp+SwwoGDgcYcApNO5Z7mUsHYWmBoDkLUzHwHDGWFsAqNCMAAR2VsDpAbynZYWmsAA9w0a7loue1k7C0yNAcjWAYx8uVfLAQw5FQMQkZY6UVqrAVJ7ALJsWnKHTtBAwcCjhSYw+bKi9gCkthMYyOkYgIhYA+Te3LWmRYsBSEs1QGo8pZ+cigGIiDVA7s0dO0EDBWui1HCWoJzWzwJjDRAVgQGIyJNqgNQyH3KWgccdrgQNuE8QKy7WAJGHYwAi0lINkJeX6XR9rZ4F5i6XI9BSAFLzWWD21ACpcbmnUscARGTt2jmA4TRa6chYLQEIMG3stdIHyF371qg9AGn9LDA1NueRUzEAEVnrNwOotw+BtOPSSh8gdw0a7tIUV1xaPAtM3uzFAERFYAAislUDpNYApLUaoKAgV5fAOndpiisurfcBUuO8kFMxABFpLQBprQZI7UHDXTEAkYdjACKyduYUoN4AJO3YtBKAhHB1CaxTezDT+mnwamzOI6diACLSag2QVk6DZwAqHfLrFrEGiDwQAxCR1gKQZQ2Q2k8HrlPH9Dwx0XXlsNSkiel5166uK0dx6XSm5UGLp8GrcV7IqbhUENk6C0ytTUeF1QCpaT4kHTsCzzwDHDoELF3q6tKY9OkD/PQTcO4csHChq0tTPH5+huVDK2eB8TR4UkAlSzhRKbJVA6TWDajW+gDpdMCHH7q6FAV5eQGffebqUpSM5RmDaqw1kZczP9/0XI3zQk7FJjAiewKQmu7zpLUaICo9WghA8nXT1vqrlnkhpyrWUnHr1i3s2bMH6enp0Ov1Zq8NGjTIIQUjchpbZ4GpdQOqtdPgqfQUFoDUsqzYOoBRY3MeOZXipeL777/HgAEDkJWVhbCwMOhkZ0LodDoGIFIfWxtQeXW6mjag8iYw+e085K8RAdqoAbKnBlct80JOpbgJbOLEiRg2bBiysrJw69Yt3Lx50/i4ceNGaZSRqHTZ6gSt1g2oZZ8ItZ7NRqVPC2eB2dMHiMGfrFAcgNLS0jBu3DgEuevl6YmU0toRpOUF7lgDRLZIgVjNZ4Fpbf0lp1EcgJKTk/H777+XRlmIXENrG1DL+WEAIlvYBEYeTPFS0bVrV0yaNAlHjhxBgwYN4GuxQe3Ro4fDCkfkFFrbgFo26TEAkS0MQOTBFNcAjRw5EufPn8fMmTPRu3dv9OrVy/h49NFHi1WIDz74AHFxcQgICEBiYiL27NlT6Pjz589HQkICAgMDERMTg+eeew737t0zvj5nzhy0aNECoaGhiIiIQK9evXDs2LFilY08gD1nganxNHiATWBUOGl50OsNDzX2m+FZYFRMigOQXq+3+ciXd0CzU2pqKiZMmICUlBTs378fjRo1QnJyMtLT062Ov3z5crz44otISUnB0aNHsWjRIqSmpuKll14yjrN161aMHj0av/32GzZs2ID79++jU6dOyM7OVlw+8gBaO4JkDRDZy7K/mBqXeV4HiIrJ5UvF3LlzMXLkSAwdOhQAsGDBAqxduxaLFy/Giy++WGD8nTt3ok2bNujfvz8AIC4uDv369cPu3buN46xbt87sPUuXLkVERAT27duHBx98sBTnhlTJ1llgaj0NnjVAZC/58pCba1r+dTrDla7VQGtXcienKdYSvnXrVnTv3h01atRAjRo10KNHD2zfvl3xdHJzc7Fv3z4kJSWZCuTlhaSkJOzatcvqe1q3bo19+/YZm8n++ecf/Pjjj+jSpYvNz8nIyAAAlCtXzurrOTk5yMzMNHuQB9HaEaRloOOOgGyRLw/5+aawrKblXWs1uOQ0igPQsmXLkJSUhKCgIIwbNw7jxo1DYGAgOnbsiOXLlyua1rVr15Cfn4/IyEiz4ZGRkbh8+bLV9/Tv3x8zZ85E27Zt4evri/j4eLRv396sCUxOr9dj/PjxaNOmDerXr291nDlz5iA8PNz4iImJUTQfpHJeXqajXS1sQC1rgNTal4lKn2X4l5YVtS7vvBcYKaA4AM2ePRtvvPEGUlNTjQEoNTUVr732Gl555ZXSKKOZLVu24NVXX8WHH36I/fv3Y9WqVVi7dq3Nzx49ejQOHz6MFStW2Jzm5MmTkZGRYXycP3++tIpP7sry9hGWz9W0AbU8IlZjswY5h61lRc3Lu7XnapofchrFS8U///yD7t27Fxjeo0cPm7UwtlSoUAHe3t64cuWK2fArV64gKirK6numTZuGgQMHYsSIEQCABg0aIDs7G6NGjcKUKVPgJdvAjxkzBj/88AO2bduGKlWq2CyHv78//P39FZWdNMbHx9AHQgv3ArPs2CodFatpHsg5LGtPpGVeTU2lPAuMiknx4WBMTAw2bdpUYPjGjRsVNx35+fmhWbNmZtPT6/XYtGkTWrVqZfU9d+7cMQs5AOD9/9W4Qgjj3zFjxuDbb7/FL7/8gmrVqikqF3mgomqA1NR0pIWjenIOLSwrrAGiYlK8VEycOBHjxo3DwYMH0bp1awDAr7/+iqVLl+Kdd95RXIAJEyZg8ODBaN68OVq2bIn58+cjOzvbeFbYoEGDULlyZcyZMwcA0L17d8ydOxdNmjRBYmIiTp48iWnTpqF79+7GIDR69GgsX74cq1evRmhoqLE/UXh4OAIDAxWXkTyA/AaiErVuQG2d2qymeSDn8JQApKYaLXIaxUv5M888g6ioKLz99tv46quvAAB16tRBamoqevbsqbgATzzxBK5evYrp06fj8uXLaNy4MdatW2fsGH3u3DmzGp+pU6dCp9Nh6tSpSEtLQ8WKFdG9e3fMnj3bOM5HH30EAGjfvr3ZZy1ZsgRDhgxRXEbyANZqgLRwGrxad2rkHLZum6KmZUVrZ3GS0xRrqXj00UeLfdVna8aMGYMxY8ZYfW3Lli1m//v4+CAlJQUpKSk2pyc1hRHZTUudoG2dBq+meSDn0EJYZhMYFRNPCSECTBtILXSCtnUavJr6MZFzaO00eHaCJgXsWirKlSuH48ePo0KFCihbtix0Op3NcW/cuOGwwhE5DWuAyBPZqgFSU58Z+UkxvA4QKWDXUjFv3jyEhoYanxcWgIhUSUudoG3VAKlpHsg5tNAEptMZymt51XN2gqYi2LWUDx482PicnYhJk7R0GrxlDRCvA0S22LoOkNqWlaICkNrmh5xCcR+g/fv349ChQ8b/V69ejV69euGll15Cbm6uQwtH5DRaagJjDRDZSyvLipbWX3IaxQHoqaeewvHjxwEYrgr9xBNPICgoCCtXrsR///tfhxeQyCl4Gjx5IvkykZNjfbgaMABRMSgOQMePH0fjxo0BACtXrkS7du2wfPlyLF26FN98842jy0fkHNIGMj8fkC6joNYNKC+ESPaSLxP37lkfrgZSEzXPAiMFFAcgIQT0ej0Aw+0vunTpAsBwi4xr1645tnREzmLtVFq1BiDWAJG9bAUgtXUaZg0QFYPiANS8eXPMmjULn3/+ObZu3YquXbsCAE6fPm28ejOR6lh2HJb/BdS1AbV1M1Q1deQm55AvE2quASoqAKkt0JFTKA5A8+fPx/79+zFmzBhMmTIFNWrUAAB8/fXXxnuDEamOVmuA1LxTo9KnlWVF3oQtUev6S06jeKlo2LCh2VlgkjfffNN4M1Ii1SkqAKlp2ZYf7aq5YyuVPnaCJg+muAbo/PnzuHDhgvH/PXv2YPz48fjf//4HX1YzklpptQbo7l3rw4kA7SwrRQUgL971iQpSvFT0798fmzdvBgBcvnwZDz/8MPbs2YMpU6Zg5syZDi8gkVNYXg8FUO9p8PIDETU3a1Dp01oTmLWzwHx8DFeLJrKgOAAdPnwYLVu2BAB89dVXqF+/Pnbu3IkvvvgCS5cudXT5iJyDNUDkiTzhLDAu92SD4gB0//59+Pv7AzCcBt+jRw8AQO3atXHp0iXHlo7IWbR0FphWjuqp9GllWbF2HSA13tiVnEpxAKpXrx4WLFiA7du3Y8OGDXjkkUcAABcvXkT58uUdXkAip9BqDZB8p6amjtzkHJ5wGrza5oWcRnEAev3117Fw4UK0b98e/fr1Q6NGjQAAa9asMTaNEamOVgOQms/sodKnlRogeXn//0K9DEBUFMVLRvv27XHt2jVkZmaibNmyxuGjRo1CUFCQQwtH5DRaOg1eKzs1Kn1aWVYs118/PwYgKlKxlgxvb2+z8AMAcXFxjigPkWtYOwuMNUCkdVoOQPKzwIissGvJaNq0KTZt2oSyZcuiSZMm0BVySuH+/fsdVjgip7FWA6TW0+C1slOj0qe1s8CAgjW4XO7JBruWjJ49exrP/OrVq1dplofINXgWGHkirSwrhQUgtYU5chq7lvKUlBSrz4k0o6g+QGq6kqxWdmpU+rSyrLAGiIqhREtGVlYW9FKP+/8XFhZWogIRuURRG1A1XUlWKzs1Kn1aWVbkJykwAJGdFB/Wnj59Gl27dkVwcDDCw8NRtmxZlC1bFmXKlCnQMZpINQrrBK22DSivA0T2ki8Tar5quLUDGHaCpiIoXjKefPJJCCGwePFiREZGFtohmkg1CqsBUltw0MrF7aj0aaUGSF7e/HxACPUewJDTKF4y/vjjD+zbtw8JCQmlUR4i1yisE7TaNqA8DZ7sZWtZUVvHYcsDGHnXDLXNCzmN4iawFi1a4Pz586VRFiLXsdaHQDoNXm3BQStH9VT6tLKsWAYgtV7ElJxK8VL+6aef4umnn0ZaWhrq168PX4t03bBhQ4cVjshp5BtJKfhooQYoN9f6cCJAuwFIfg0vBiCyQfFSfvXqVZw6dQpDhw41DtPpdBBCQKfTIV++4BGphVYDkD3DyXMxAJEHU7yUDxs2DE2aNMGXX37JTtCkHQxA5Inky4QQ1oerAQMQFYPipfzs2bNYs2YNatSoURrlIXINTwhA3BGQJVvLhNqWecs+fAxAZAfFnaA7dOiAP/74ozTKQuQ6hQUgtW1AtbJTo9KnldpC1gBRMSheyrt3747nnnsOhw4dQoMGDQp0gu7Ro4fDCkfkNFqqAdLpDPNj2R9PbfNBpc/WMqG2U8ctrwPEAER2ULxFfPrppwEAM2fOLPAaO0GTasnv9SVdQ0Stp8EDhjIzAFFRtFoDJL8OkJru40dOpXgpt7z3F5EmaKkGCDCUWX5hO2kYkZxWAxBrgMgODonGt27dcsRkiFxHawHI2kZfjfNBpYsBiDyY4gD0+uuvIzU11fh/7969Ua5cOVSuXJmdo0m9LAOQXm86LVhtOwPAepnVOB9UuhiAyIMpDkALFixATEwMAGDDhg3YuHEj1q1bh86dO2PSpEkOLyCRU1gGIPml9NW2MwCsl5k7ArKklTMGGYCoGBQv5ZcvXzYGoB9++AF9+vRBp06dEBcXh8TERIcXkMgpCgtAatyAsgaI7KGVs8B4HSAqBsU1QGXLljXeDHXdunVISkoCAAgheAYYqZdlAFL7BpQBiOyhlSYwra2/5BSKl/LHHnsM/fv3R82aNXH9+nV07twZAHDgwAFeHZrUS2sbUAYgsoeXl+G6UfLbYADqW1a0tv6SUyheyufNm4e4uDicP38eb7zxBkJCQgAAly5dwrPPPuvwAhI5hbVO0BI1XkeEAYjs5eMD3L9fcJiaFLb+MgCRDYqXcl9fXzz//PMFhj/33HMOKRCRS1heCFHtR5AMQGQvLQYg+fqrxgMYcopiLeUnTpzA5s2bkZ6eXuDCiNOnT3dIwYicSmtV6AxAZC8tLCtaW3/JKRQv5Z988gmeeeYZVKhQAVFRUdDpdMbXdDodAxCpk9Y2oFrYqZFzWFsu1HwWmBbWX3IKxVvEWbNmYfbs2XjhhRdKozxErqG1PgTWyqzG+aDSp4WrhjMAUTEobhy9efMmevfuXRplIXIdrfUhYA0Q2UsLywoDEBWD4i177969sX79+tIoC5HraG0DqoWdGjmHFpYV+UGKFtZfcgrFS3mNGjUwbdo0/Pbbb2jQoAF8LdqKx40b57DCETkNAxB5Ki0sK1pbf8kpFC/lH3/8MUJCQrB161Zs3brV7DWdTscAROqktT5AvBcY2UtrAUgLl7Egp1C8lJ8+fbo0ykHkWp7QB0htZ/aQc2hhWdHaAQw5hQq37ESlwBMuhOjn5/xykPuztlyouQZICwcw5BTFWsovXLiANWvW4Ny5c8jNzTV7be7cuQ4pGJFTaa0PgRaO6sk5rC0Xalvmtbb+klMoDkCbNm1Cjx49UL16dfz999+oX78+zpw5AyEEmjZtWhplJCp9WtuAMgCRvSyXC29vww1S1URr6y85heK6wcmTJ+P555/HoUOHEBAQgG+++Qbnz59Hu3btin19oA8++ABxcXEICAhAYmIi9uzZU+j48+fPR0JCAgIDAxETE4PnnnsO9+7dK9E0ycN5ws1QGYDIGsvlQm3NXwADEBWL4i370aNHMWjQIACAj48P7t69i5CQEMycOROvv/664gKkpqZiwoQJSElJwf79+9GoUSMkJycjPT3d6vjLly/Hiy++iJSUFBw9ehSLFi1CamoqXnrppWJPk0hzG1D2ASJ7WS4XDEDkIRQHoODgYGO/n0qVKuHUqVPG165du6a4AHPnzsXIkSMxdOhQ1K1bFwsWLEBQUBAWL15sdfydO3eiTZs26N+/P+Li4tCpUyf069fPrIZH6TSJNLcBtVZm1gCRNZbLhRqXE62tv+QUigPQAw88gB07dgAAunTpgokTJ2L27NkYNmwYHnjgAUXTys3Nxb59+5CUlGQqkJcXkpKSsGvXLqvvad26Nfbt22cMPP/88w9+/PFHdOnSpdjTzMnJQWZmptmDPIzWNqBauLYLOQebwMhDKV7S586di6ysLADAyy+/jKysLKSmpqJmzZqKzwC7du0a8vPzERkZaTY8MjISf//9t9X39O/fH9euXUPbtm0hhEBeXh6efvppYxNYcaY5Z84cvPzyy4rKThqjteuIWO7EfH3V17GVnIMBiDyUohqg/Px8XLhwAVWrVgVgaA5bsGAB/vzzT3zzzTeIjY0tlULKbdmyBa+++io+/PBD7N+/H6tWrcLatWvxyiuvFHuakydPRkZGhvFx/vx5B5aYVEFr1xGxFoCIrNFiHyC1H8CQUyha0r29vdGpUyccPXoUZcqUKfGHV6hQAd7e3rhy5YrZ8CtXriAqKsrqe6ZNm4aBAwdixIgRAIAGDRogOzsbo0aNwpQpU4o1TX9/f/j7+5d4fkjFtH4hRAYgskXrNUBqPIAhp1C8ZNSvXx///POPQz7cz88PzZo1w6ZNm4zD9Ho9Nm3ahFatWll9z507d+BlsUB7///CL4Qo1jSJNFeFzgBE9tJ6AFLj+ktOoXhJnzVrFp5//nm88soraNasGYKDg81eDwsLUzS9CRMmYPDgwWjevDlatmyJ+fPnIzs7G0OHDgUADBo0CJUrV8acOXMAAN27d8fcuXPRpEkTJCYm4uTJk5g2bRq6d+9uDEJFTZOoAK1tQC13YjwFnmzhWWDkoRQHIOlsqx49ekAn61QphIBOp0O+fMGzwxNPPIGrV69i+vTpuHz5Mho3box169YZOzGfO3fOrMZn6tSp0Ol0mDp1KtLS0lCxYkV0794ds2fPtnuaRAVo/UKIatypkXNosQ8QAxDZQfGSvmTJEsTExBhrWyR6vR7nzp0rViHGjBmDMWPGWH1ty5YtZv/7+PggJSUFKSkpxZ4mUQFa24AyAJG9tNAEJj9I0cL6S06heEkfNmwYLl26hIiICLPh169fR1JSEgYPHuywwhE5DQMQeSotBCCtrb/kFIrr9qWmLktZWVkICAhwSKGInE5rR5DsA0T2YgAiD2X3kj5hwgQAgE6nw7Rp0xAUFGR8LT8/H7t370bjxo0dXkAip9DpDCFIOgVe7dcRsSwza4DIFi32AVL7+ktOYfeSfuDAAQCGGqBDhw7BT7bS+Pn5oVGjRnj++ecdX0IiZ/H2NgUgtV9HhE1gZC+tnQWmhet4kVPYHYA2b94MABg6dCjeeecdxae7E7k9KehoYQPKJjCyl9abwNR4AENOUayzwIg0SdqIaqEPAWuAyF5aD0BqXH/JKRiNiSTyAKT2PgQMQGQvLfYBYgAiOzAAEUls1QCpsQpddpICAAYgss1y2bBcdtSAAYiKQYVbdqJSoqUmMItb1LAPENlkGYAslx01YACiYmAAIpJoKQCxBojsxRog8lAMQEQSLQUgy6N4BiCyxbJ2kDVA5CEYgIgktjpBq7EPEJvAyF5aaAKzvJK72k9iIKdQ4ZadqJRoqQaITWBkLy00gUlXcge0sf6SUzAAEUm0dCFENoGRvbRQAwRo6yxOcgouGUQSLdUAMQCRvbTQBwjQ1vpLTsEARCTR0oUQLZsx2AeIbNFCExjAJjBSjAGISKKlKnQt3OCSnEPrTWAMQGSDCrfsRKVEyxtQBiCyhQGIPBQDEJFEyxtQBiCyxbJ5VK1NYFpef6lUMAARSbTUB8gS+wCRLVqsAdLa+kulggGISKKlPkCWWANEtmgxALEGiOyggS07kYPIg05enum5FjagOp2rS0DuSitngUnrqeV1vLRwAEOlgksGkUQedHJzrQ9Xq/v3XV0Ccle8DhB5KAYgIgkDEHkiyxogtfYXYwAihRiAiCTyDaU8MGihCl0e6IjkLAOQWptLGYBIIQ1s2YkcxFYA0sIGlDVAZItWOsgzAJFCDEBEEi03gbEGiGzx8XF1CRyDAYgUYgAikmgtANWta3peqZLrykHuTa1NXpYYgEghBiAiidYC0MqVQGQk0KgRMGqUq0tD7uz554GAAOCjj1xdkuKzdSFELfTho1KhkbpPIgfQWifounWBCxcM86WVo3wqHW++Cbz6qrr7A1mrAVLrwQs5BQMQkUQedLRQAwRop38HlT41hx/A+oUQ1XrwQk7BpYNIorUmMCJPYq0Gl+suFYIBiEjCAESkXtZqcLnuUiEYgIgkWusDRORJrB3AMABRIbhlJ5Jo+UKIRFrHAEQKMQARSdgERqReDECkEAMQkYQBiEi9GIBIIQYgIgkDEJF68SwwUogBiEjCTtBE6sUaIFKIW3YiiRYvhEjkKeTrKS+ESHbg0kEkYRMYkXpZW0+57lIhGICIJNaOIC2HE5F7YgAihRiAiCS2Npa8kSiR+2MAIoUYgIgk3IASqRfXX1KIAYhIwg0okXpx/SWFGICIJNyAEqkX119SiAGISGJtY8nTaInUgQGIFOLWnUjCDSiRenH9JYUYgIgk1mp7uAElUgdr6y9rcKkQXDqIJDyCJFIvrr+kEAMQkYQbUCL14vpLCjEAEUnYCZpIvRiASCFu3Ykk3IASqRfXX1KIAYhIwg0okXpx/SWF3CIAffDBB4iLi0NAQAASExOxZ88em+O2b98eOp2uwKNr167GcbKysjBmzBhUqVIFgYGBqFu3LhYsWOCMWSE14waUSL24/pJCLg9AqampmDBhAlJSUrB//340atQIycnJSE9Ptzr+qlWrcOnSJePj8OHD8Pb2Ru/evY3jTJgwAevWrcOyZctw9OhRjB8/HmPGjMGaNWucNVukRuwDRKReDECkkMu37nPnzsXIkSMxdOhQY01NUFAQFi9ebHX8cuXKISoqyvjYsGEDgoKCzALQzp07MXjwYLRv3x5xcXEYNWoUGjVqVGjNEpHVjaWPj/PLQUTKMQCRQi4NQLm5udi3bx+SkpKMw7y8vJCUlIRdu3bZNY1Fixahb9++CA4ONg5r3bo11qxZg7S0NAghsHnzZhw/fhydOnWyOo2cnBxkZmaaPcgD8UKIROrFGlxSyKVLx7Vr15Cfn4/IyEiz4ZGRkbh8+XKR79+zZw8OHz6MESNGmA1/7733ULduXVSpUgV+fn545JFH8MEHH+DBBx+0Op05c+YgPDzc+IiJiSn+TJF68QiSSL24/pJCqo7HixYtQoMGDdCyZUuz4e+99x5+++03rFmzBvv27cPbb7+N0aNHY+PGjVanM3nyZGRkZBgf58+fd0bxyd1wA0qkXlx/SSGXdnCoUKECvL29ceXKFbPhV65cQVRUVKHvzc7OxooVKzBz5kyz4Xfv3sVLL72Eb7/91nhmWMOGDXHw4EG89dZbZs1tEn9/f/j7+5dwbkj1uAElUi+uv6SQS2uA/Pz80KxZM2zatMk4TK/XY9OmTWjVqlWh7125ciVycnLw5JNPmg2/f/8+7t+/Dy+Ltl9vb2/o9XrHFZ60hxtQIvXi+ksKufwUlwkTJmDw4MFo3rw5WrZsifnz5yM7OxtDhw4FAAwaNAiVK1fGnDlzzN63aNEi9OrVC+XLlzcbHhYWhnbt2mHSpEkIDAxEbGwstm7div/973+YO3eu0+aLVIhngRGpFwMQKeTyrfsTTzyBq1evYvr06bh8+TIaN26MdevWGTtGnzt3rkBtzrFjx7Bjxw6sX7/e6jRXrFiByZMnY8CAAbhx4wZiY2Mxe/ZsPP3006U+P6Ri3IASqRfXX1LI5QEIAMaMGYMxY8ZYfW3Lli0FhiUkJEAIYXN6UVFRWLJkiaOKR56CG1Ai9eL6Swqp+iwwIofiBpRIvXgdL1KIAYhIwg0okXrxQoikEJcOIglrgIjUi+svKcQARCThBpRIvbj+kkIMQEQSngZPpF4MQKQQAxCRhBtQIvXi+ksKMQARSbgBJVIvrr+kEAMQkYQbUCL14vpLCjEAEUm4ASVSL66/pBADEJGEG1Ai9eL6SwoxABFJeCFEIvXihRBJIS4dRBKeBk+kXqwBIoUYgIgk3IASqRfXX1KIAYhIwg0okXpx/SWFGICIJNyAEqkX119SiAGISMINKJF6cf0lhRiAiCTcgBKpF8/iJIUYgIgkDEBE6sX1lxRiACKSWDuC5GnwROrA6wCRQlw6iCQ8giRSL66/pBADEJGEG1Ai9eL6SwoxABFJuAElUi+uv6QQAxCRhBtQIvXi+ksKMQARSbgBJVIvrr+kEAMQkYQbUCL14vpLCjEAEUl4GjyRejEAkUIMQERylhtMbkCJ1IEBiBRiACKSs6wF4gaUSB14IURSiEsHkRxrgIjUiTVApBADEJEcAxCROjEAkUIMQERyDEBE6sQARAoxABHJMQARqRMDECnEAEQkZ7nB5GnwROrAAEQKMQARybEGiEidrJ3xxfWXCsEARCTHAESkTjqd4SHH9ZcKwQBEJMcARKReXH9JAQYgIjleCJFIvSzXV14IkQrBpYNIjkeQROrF9ZcUYAAikuMGlEi9uP6SAgxARHI8DZ5IvRiASAEGICI5bkCJ1IvrLynAAEQkxw0okXpx/SUFGICI5LgBJVIvrr+kAAMQkRw3oETqxfWXFGAAIpLjBpRIvbj+kgIMQERyvBAikXrxQoikAM/xLSYhBPLy8pCfn+/qomier68vvJ0VRHgaPJF6sQaIFODWvRhyc3Nx6dIl3Llzx9VF8Qg6nQ5VqlRBSEhI6X8YN6BE6sX1lxRgAFJIr9fj9OnT8Pb2RnR0NPz8/KCzvAMxOYwQAlevXsWFCxdQs2bN0q8J4gaUSL24/pICDEAK5ebmQq/XIyYmBkFBQa4ujkeoWLEizpw5g/v37zMAEZFt7MNHCrCHWDF5sXOd0zi1ho0BiEi9uP6SAtyLE8lxA0qkXlx/SQEGICI5bkCJ1MtyfWX/TCoEAxA5XPv27TF+/Hi7x1+6dCnKlClTauVRJCzM/H+eBk+kHvIAFBrqunKQKrhFAPrggw8QFxeHgIAAJCYmYs+ePTbHbd++PXQ6XYFH165dzcY7evQoevTogfDwcAQHB6NFixY4d+5cac8KqV2lSub/swaISD3k62tEhOvKQarg8gCUmpqKCRMmICUlBfv370ejRo2QnJyM9PR0q+OvWrUKly5dMj4OHz4Mb29v9O7d2zjOqVOn0LZtW9SuXRtbtmzBn3/+iWnTpiEgIMBZs0VqxQBEpF5375qeMwBREVwegObOnYuRI0di6NChqFu3LhYsWICgoCAsXrzY6vjlypVDVFSU8bFhwwYEBQWZBaApU6agS5cueOONN9CkSRPEx8ejR48eiPDwFaJ9+/YYO3Ysxo8fj7JlyyIyMhKffPIJsrOzMXToUISGhqJGjRr46aefjO/ZunUrWrZsCX9/f1SqVAkvvvgi8vLyjK9nZ2dj0KBBCAkJQaVKlfD2228X+NycnBw8//zzqFy5MoKDg5GYmIgtW7Y4Y5aVYwAiUi/5gbOHb++paC4NQLm5udi3bx+SkpKMw7y8vJCUlIRdu3bZNY1Fixahb9++CA4OBmC4UOHatWtRq1YtJCcnIyIiAomJifjuu+9sTiMnJweZmZlmD8WaNweqVHHuo3lzxcX87LPPUKFCBezZswdjx47FM888g969e6N169bYv38/OnXqhIEDB+LOnTtIS0tDly5d0KJFC/zxxx/46KOPsGjRIsyaNcs4vUmTJmHr1q1YvXo11q9fjy1btmD//v1mnzlmzBjs2rULK1aswJ9//onevXvjkUcewYkTJ5R/z6WNAYhIvRiASAnhQmlpaQKA2Llzp9nwSZMmiZYtWxb5/t27dwsAYvfu3cZhly5dEgBEUFCQmDt3rjhw4ICYM2eO0Ol0YsuWLVank5KSIgAUeGRkZBQY9+7du+LIkSPi7t275i9UriwE4NxH5cp2fMsm7dq1E23btjX+n5eXJ4KDg8XAgQMLfH+7du0SL730kkhISBB6vd74+gcffCBCQkJEfn6+uH37tvDz8xNfffWV8fXr16+LwMBA8Z///EcIIcTZs2eFt7e3SEtLMytLx44dxeTJk4UQQixZskSEh4fbLLfN77w0HDtm/h3n5pb+ZxKRY8jX3SlTXF0acoGMjAyb+29Lqj7FZdGiRWjQoAFatmxpHKbX6wEAPXv2xHPPPQcAaNy4MXbu3IkFCxagXbt2BaYzefJkTJgwwfh/ZmYmYmJilBUmKqoYc1BCxfjMhg0bGp97e3ujfPnyaNCggXFYZGQkACA9PR1Hjx5Fq1atzC5E2KZNG2RlZeHChQu4efMmcnNzkZiYaHy9XLlySEhIMP5/6NAh5Ofno1atWmblyMnJQfny5RWXv9SxBohIG1gDREVwaQCqUKECvL29ceXKFbPhV65cQVQRO/fs7GysWLECM2fOLDBNHx8f1K1b12x4nTp1sGPHDqvT8vf3h7+/fzHmQOb330v2fifx9fU1+1+n05kNk8KOFCRLKisrC97e3ti3b1+B21g45eamSlmeOssrfhOpEwMQFcGlW3c/Pz80a9YMmzZtMg7T6/XYtGkTWrVqVeh7V65ciZycHDz55JMFptmiRQscO3bMbPjx48cRGxvruMJ7gDp16mDXrl0QQhiH/frrrwgNDUWVKlUQHx8PX19f7N692/j6zZs3cfz4ceP/TZo0QX5+PtLT01GjRg2zR1Ehl4io2Lh9oSK4/PB2woQJ+OSTT/DZZ5/h6NGjeOaZZ4xnJQHAoEGDMHny5ALvW7RoEXr16mW1GWXSpElITU3FJ598gpMnT+L999/H999/j2effbbU50dLnn32WZw/fx5jx47F33//jdWrVyMlJQUTJkyAl5cXQkJCMHz4cEyaNAm//PILDh8+jCFDhpjdJ61WrVoYMGAABg0ahFWrVuH06dPYs2cP5syZg7Vr17pw7goxdqzhb8eOri0HESnzzjuGvzVqAG3auLYs5PZc3gfoiSeewNWrVzF9+nRcvnwZjRs3xrp164x9Uc6dO1fgxqPHjh3Djh07sH79eqvTfPTRR7FgwQLMmTMH48aNQ0JCAr755hu0bdu21OdHSypXrowff/wRkyZNQqNGjVCuXDkMHz4cU6dONY7z5ptvIisrC927d0doaCgmTpyIjIwMs+ksWbIEs2bNwsSJE5GWloYKFSrggQceQLdu3Zw9S/aZNw948kmgUSNXl4SIlBg7FvjXv4D4eMCiuZ/Ikk7I2zcIgKETdHh4ODIyMhBmcWuEe/fu4fTp06hWrRovrOgk/M6JiMgehe2/Lbm8CYyIiIjI2RiAiIiIyOMwABEREZHHYQAiIiIij8MAVEzsO+48/K6JiMjRGIAUkq6afOfOHReXxHPk5uYCQIErSRMRERWXy68DpDbe3t4oU6YM0v//rsNBQUFm98oix9Lr9bh69SqCgoLg48PFlYiIHIN7lGKQbuEghSAqXV5eXqhatSqDJhEROQwDUDHodDpUqlQJERERuH//vquLo3l+fn4FrgZORERUEgxAJeDt7c1+KURERCrEw2oiIiLyOAxARERE5HEYgIiIiMjjsA+QFdKF9zIzM11cEiIiIrKXtN+25wK6DEBW3L59GwAQExPj4pIQERGRUrdv30Z4eHih4+gE7zNQgF6vx8WLFxEaGurwa89kZmYiJiYG58+fR1hYmEOn7Q44f+qn9XnU+vwB2p9Hzp/6ldY8CiFw+/ZtREdHF3n5FNYAWeHl5YUqVaqU6meEhYVpdsEGOH9aoPV51Pr8AdqfR86f+pXGPBZV8yNhJ2giIiLyOAxARERE5HEYgJzM398fKSkp8Pf3d3VRSgXnT/20Po9anz9A+/PI+VM/d5hHdoImIiIij8MaICIiIvI4DEBERETkcRiAiIiIyOMwABEREZHHYQBysNmzZ6N169YICgpCmTJlrI5z7tw5dO3aFUFBQYiIiMCkSZOQl5dX6HRv3LiBAQMGICwsDGXKlMHw4cORlZVVCnOgzJYtW6DT6aw+9u7da/N97du3LzD+008/7cSS2y8uLq5AWV977bVC33Pv3j2MHj0a5cuXR0hICP7973/jypUrTiqxMmfOnMHw4cNRrVo1BAYGIj4+HikpKcjNzS30fe78G37wwQeIi4tDQEAAEhMTsWfPnkLHX7lyJWrXro2AgAA0aNAAP/74o5NKqtycOXPQokULhIaGIiIiAr169cKxY8cKfc/SpUsL/FYBAQFOKrEyM2bMKFDW2rVrF/oeNf1+1rYnOp0Oo0ePtjq+Gn67bdu2oXv37oiOjoZOp8N3331n9roQAtOnT0elSpUQGBiIpKQknDhxosjpKl2PlWIAcrDc3Fz07t0bzzzzjNXX8/Pz0bVrV+Tm5mLnzp347LPPsHTpUkyfPr3Q6Q4YMAB//fUXNmzYgB9++AHbtm3DqFGjSmMWFGndujUuXbpk9hgxYgSqVauG5s2bF/rekSNHmr3vjTfecFKplZs5c6ZZWceOHVvo+M899xy+//57rFy5Elu3bsXFixfx2GOPOam0yvz999/Q6/VYuHAh/vrrL8ybNw8LFizASy+9VOR73fE3TE1NxYQJE5CSkoL9+/ejUaNGSE5ORnp6utXxd+7ciX79+mH48OE4cOAAevXqhV69euHw4cNOLrl9tm7ditGjR+O3337Dhg0bcP/+fXTq1AnZ2dmFvi8sLMzstzp79qyTSqxcvXr1zMq6Y8cOm+Oq7ffbu3ev2bxt2LABANC7d2+b73H33y47OxuNGjXCBx98YPX1N954A++++y4WLFiA3bt3Izg4GMnJybh3757NaSpdj4tFUKlYsmSJCA8PLzD8xx9/FF5eXuLy5cvGYR999JEICwsTOTk5Vqd15MgRAUDs3bvXOOynn34SOp1OpKWlObzsJZGbmysqVqwoZs6cWeh47dq1E//5z3+cU6gSio2NFfPmzbN7/Fu3bglfX1+xcuVK47CjR48KAGLXrl2lUELHe+ONN0S1atUKHcddf8OWLVuK0aNHG//Pz88X0dHRYs6cOVbH79Onj+jatavZsMTERPHUU0+VajkdJT09XQAQW7dutTmOre2RO0pJSRGNGjWye3y1/37/+c9/RHx8vNDr9VZfV9NvJ4QQAMS3335r/F+v14uoqCjx5ptvGofdunVL+Pv7iy+//NLmdJSux8XBGiAn27VrFxo0aIDIyEjjsOTkZGRmZuKvv/6y+Z4yZcqY1agkJSXBy8sLu3fvLvUyK7FmzRpcv34dQ4cOLXLcL774AhUqVED9+vUxefJk3LlzxwklLJ7XXnsN5cuXR5MmTfDmm28W2mS5b98+3L9/H0lJScZhtWvXRtWqVbFr1y5nFLfEMjIyUK5cuSLHc7ffMDc3F/v27TP77r28vJCUlGTzu9+1a5fZ+IBhnVTTbwWgyN8rKysLsbGxiImJQc+ePW1ub9zBiRMnEB0djerVq2PAgAE4d+6czXHV/Pvl5uZi2bJlGDZsWKE33lbTb2fp9OnTuHz5stlvFB4ejsTERJu/UXHW4+LgzVCd7PLly2bhB4Dx/8uXL9t8T0REhNkwHx8flCtXzuZ7XGXRokVITk4u8may/fv3R2xsLKKjo/Hnn3/ihRdewLFjx7Bq1SonldR+48aNQ9OmTVGuXDns3LkTkydPxqVLlzB37lyr41++fBl+fn4F+oBFRka63e9lzcmTJ/Hee+/hrbfeKnQ8d/wNr127hvz8fKvr2N9//231PbbWSTX8Vnq9HuPHj0ebNm1Qv359m+MlJCRg8eLFaNiwITIyMvDWW2+hdevW+Ouvv0r9xs9KJSYmYunSpUhISMClS5fw8ssv41//+hcOHz6M0NDQAuOr+ff77rvvcOvWLQwZMsTmOGr67ayRfgclv1Fx1uPiYACyw4svvojXX3+90HGOHj1aZEc9NSnOPF+4cAE///wzvvrqqyKnL++/1KBBA1SqVAkdO3bEqVOnEB8fX/yC20nJ/E2YMME4rGHDhvDz88NTTz2FOXPmuPWl6ovzG6alpeGRRx5B7969MXLkyELf6+rfkIDRo0fj8OHDhfaRAYBWrVqhVatWxv9bt26NOnXqYOHChXjllVdKu5iKdO7c2fi8YcOGSExMRGxsLL766isMHz7chSVzvEWLFqFz586Ijo62OY6afju1YQCyw8SJEwtN6ABQvXp1u6YVFRVVoCe7dHZQVFSUzfdYdvzKy8vDjRs3bL6npIozz0uWLEH58uXRo0cPxZ+XmJgIwFD74IydZ0l+08TEROTl5eHMmTNISEgo8HpUVBRyc3Nx69Yts1qgK1eulNrvZY3Sebx48SIeeughtG7dGh9//LHiz3P2b2hNhQoV4O3tXeCMu8K++6ioKEXju4sxY8YYT4hQWhPg6+uLJk2a4OTJk6VUOscpU6YMatWqZbOsav39zp49i40bNyquMVXTbweY9mtXrlxBpUqVjMOvXLmCxo0bW31PcdbjYnFYbyIyU1Qn6CtXrhiHLVy4UISFhYl79+5ZnZbUCfr33383Dvv555/dqhO0Xq8X1apVExMnTizW+3fs2CEAiD/++MPBJXO8ZcuWCS8vL3Hjxg2rr0udoL/++mvjsL///tutO0FfuHBB1KxZU/Tt21fk5eUVaxru8hu2bNlSjBkzxvh/fn6+qFy5cqGdoLt162Y2rFWrVm7biVav14vRo0eL6Ohocfz48WJNIy8vTyQkJIjnnnvOwaVzvNu3b4uyZcuKd955x+rravv9JCkpKSIqKkrcv39f0fvc/beDjU7Qb731lnFYRkaGXZ2glazHxSqrw6ZEQgghzp49Kw4cOCBefvllERISIg4cOCAOHDggbt++LYQwLLz169cXnTp1EgcPHhTr1q0TFStWFJMnTzZOY/fu3SIhIUFcuHDBOOyRRx4RTZo0Ebt37xY7duwQNWvWFP369XP6/NmyceNGAUAcPXq0wGsXLlwQCQkJYvfu3UIIIU6ePClmzpwpfv/9d3H69GmxevVqUb16dfHggw86u9hF2rlzp5g3b544ePCgOHXqlFi2bJmoWLGiGDRokHEcy/kTQoinn35aVK1aVfzyyy/i999/F61atRKtWrVyxSwU6cKFC6JGjRqiY8eO4sKFC+LSpUvGh3wctfyGK1asEP7+/mLp0qXiyJEjYtSoUaJMmTLGMy8HDhwoXnzxReP4v/76q/Dx8RFvvfWWOHr0qEhJSRG+vr7i0KFDrpqFQj3zzDMiPDxcbNmyxey3unPnjnEcy3l8+eWXxc8//yxOnTol9u3bJ/r27SsCAgLEX3/95YpZKNTEiRPFli1bxOnTp8Wvv/4qkpKSRIUKFUR6eroQQv2/nxCGnXnVqlXFCy+8UOA1Nf52t2/fNu7rAIi5c+eKAwcOiLNnzwohhHjttddEmTJlxOrVq8Wff/4pevbsKapVqybu3r1rnEaHDh3Ee++9Z/y/qPXYERiAHGzw4MECQIHH5s2bjeOcOXNGdO7cWQQGBooKFSqIiRMnmh0FbN68WQAQp0+fNg67fv266NevnwgJCRFhYWFi6NChxlDlDvr16ydat25t9bXTp0+bfQfnzp0TDz74oChXrpzw9/cXNWrUEJMmTRIZGRlOLLF99u3bJxITE0V4eLgICAgQderUEa+++qpZbZ3l/AkhxN27d8Wzzz4rypYtK4KCgsSjjz5qFijcyZIlS6wus/IKYrX9hu+9956oWrWq8PPzEy1bthS//fab8bV27dqJwYMHm43/1VdfiVq1agk/Pz9Rr149sXbtWieX2H62fqslS5YYx7Gcx/Hjxxu/j8jISNGlSxexf/9+5xfeDk888YSoVKmS8PPzE5UrVxZPPPGEOHnypPF1tf9+Qhhq8AGIY8eOFXhNjb+dtM+yfEjzodfrxbRp00RkZKTw9/cXHTt2LDDvsbGxIiUlxWxYYeuxI+iEEMJxDWpERERE7o/XASIiIiKPwwBEREREHocBiIiIiDwOAxARERF5HAYgIiIi8jgMQERERORxGICIiIjI4zAAERERkcdhACKiUqfT6fDdd9+57fTkNm3ahDp16iA/P79Upu9Iffv2xdtvv+3qYhCpEgMQEanOpUuX0LlzZwDAmTNnoNPpcPDgQYdM+7///S+mTp0Kb29vh0yvOD777DO0bdu2yPGmTp2K2bNnIyMjwwmlItIWBiAiUp2oqCj4+/s7fLo7duzAqVOn8O9//9vh01Zi9erV6NGjR5Hj1a9fH/Hx8Vi2bJkTSkWkLQxARGTTxx9/jOjoaOj1erPhPXv2xLBhw4z/r169Gk2bNkVAQACqV6+Ol19+GXl5eTane+jQIXTo0AGBgYEoX748Ro0ahaysLLNxFi9ejHr16sHf3x+VKlXCmDFjjK/Jm8CqVasGAGjSpAl0Oh3at2+Pbdu2wdfXF5cvXzab5vjx4/Gvf/3LZrlWrFiBhx9+GAEBAcZhM2bMQOPGjbF48WJUrVoVISEhePbZZ5Gfn4833ngDUVFRiIiIwOzZs82mpdPpsHDhQnTr1g1BQUGoU6cOdu3ahZMnT6J9+/YIDg5G69atcerUKbP33bt3D+vXrzcGoA8//BA1a9ZEQEAAIiMj8fjjj5uN3717d6xYscLmPBGRDQ69tSoRacqNGzeEn5+f2Lhxo3HY9evXzYZt27ZNhIWFiaVLl4pTp06J9evXi7i4ODFjxgzjewCIb7/9VgghRFZWlqhUqZJ47LHHxKFDh8SmTZtEtWrVzO6A/eGHH4qAgAAxf/58cezYMbFnzx4xb948q9Pbs2ePACA2btwoLl26JK5fvy6EEKJWrVrijTfeML4nNzdXVKhQQSxevNjm/DZs2FC89tprZsNSUlJESEiIePzxx8Vff/0l1qxZI/z8/ERycrIYO3as+Pvvv8XixYsFALO7VQMQlStXFqmpqeLYsWOiV69eIi4uTnTo0EGsW7dOHDlyRDzwwAPikUceMfu8H374QdSqVUsIIcTevXuFt7e3WL58uThz5ozYv3+/eOedd8zG/+mnn4Sfn5+4d++ezfkiooIYgIioUD179hTDhg0z/r9w4UIRHR0t8vPzhRBCdOzYUbz66qtm7/n8889FpUqVjP/LA8vHH38sypYtK7Kysoyvr127Vnh5eYnLly8LIYSIjo4WU6ZMsVkm+fROnz4tAIgDBw6YjfP666+LOnXqGP//5ptvREhIiNnnWgoPDxf/+9//zIalpKSIoKAgkZmZaRyWnJws4uLijN+BEEIkJCSIOXPmmJVx6tSpxv937dolAIhFixYZh3355ZciICDA7PNGjhwpnn/+eWOZw8LCzD7b0h9//CEAiDNnztgch4gKYhMYERVqwIAB+Oabb5CTkwMA+OKLL9C3b194eRk2H3/88QdmzpyJkJAQ42PkyJG4dOkS7ty5U2B6R48eRaNGjRAcHGwc1qZNG+j1ehw7dgzp6em4ePEiOnbsWKJyDxkyBCdPnsRvv/0GAFi6dCn69Olj9rmW7t69a9b8JYmLi0NoaKjx/8jISNStW9f4HUjD0tPTzd7XsGFDs9cBoEGDBmbD7t27h8zMTACAEALff/+9sfnr4YcfRmxsLKpXr46BAwfiiy++KPCdBgYGAoDV75qIbGMAIqJCde/eHUIIrF27FufPn8f27dsxYMAA4+tZWVl4+eWXcfDgQePj0KFDOHHihNUwURRph15SERER6N69O5YsWYIrV67gp59+Muu3ZE2FChVw8+bNAsN9fX3N/tfpdFaHWfaVko+j0+lsDpPet2fPHuTl5aF169YAgNDQUOzfvx9ffvklKlWqhOnTp6NRo0a4deuWcRo3btwAAFSsWLHQeSMicz6uLgARubeAgAA89thj+OKLL3Dy5EkkJCSgadOmxtebNm2KY8eOoUaNGnZNr06dOli6dCmys7ONtTG//vorvLy8kJCQgNDQUMTFxWHTpk146KGHipyen58fAFi9bs+IESPQr18/VKlSBfHx8WjTpk2h02rSpAmOHDli13yUhtWrV6Nr165mp+D7+PggKSkJSUlJSElJQZkyZfDLL7/gscceAwAcPnwYVapUQYUKFVxVbCJVYgAioiINGDAA3bp1w19//YUnn3zS7LXp06ejW7duqFq1Kh5//HF4eXnhjz/+wOHDhzFr1iyr00pJScHgwYMxY8YMXL16FWPHjsXAgQONzUQzZszA008/jYiICHTu3Bm3b9/Gr7/+irFjxxaYXkREBAIDA7Fu3TpUqVIFAQEBCA8PBwAkJycjLCwMs2bNwsyZM4ucz+TkZHz22WfF+YocYs2aNWbl/OGHH/DPP//gwQcfRNmyZfHjjz9Cr9cjISHBOM727dvRqVMnVxSXSNXYBEZERerQoQPKlSuHY8eOoX///mavJScn44cffsD69evRokULPPDAA5g3bx5iY2OtTisoKAg///wzbty4gRYtWuDxxx9Hx44d8f777xvHGTx4MObPn48PP/wQ9erVQ7du3XDixAmr0/Px8cG7776LhQsXIjo6Gj179jS+5uXlhSFDhiA/Px+DBg0qcj4HDBiAv/76C8eOHbPna3GoU6dO4eTJk0hOTjYOK1OmDFatWoUOHTqgTp06WLBgAb788kvUq1cPgOGU+e+++w4jR450enmJ1E4nhBCuLgQRUWkZPnw4rl69ijVr1tg1/qRJk5CZmYmFCxeWcsnMzZ07Fxs3bsSPP/5o93s++ugjfPvtt1i/fn0ploxIm1gDRESalJGRgR07dmD58uVWm85smTJlCmJjYwt0aC5tVapUweTJkxW9x9fXF++9914plYhI21gDRESa1L59e+zZswdPPfUU5s2b5+riEJGbYQAiIiIij8MmMCIiIvI4DEBERETkcRiAiIiIyOMwABEREZHHYQAiIiIij8MARERERB6HAYiIiIg8DgMQEREReZz/A6totRPphvpEAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "import nexus as nx\n", "import numpy as np\n", "import matplotlib.pyplot as plt\n", "\n", "site = nx.Hyperfine(magnetic_field = 33,\n", " isotropic = True)\n", "\n", "sample = nx.SimpleSample(thickness=3000,\n", " composition = [[\"Fe\", 1]],\n", " density = 7.874,\n", " isotope = nx.lib.moessbauer.Fe57,\n", " abundance = 0.02119,\n", " lamb_moessbauer = 0.796,\n", " hyperfine_sites = [site])\n", "\n", "beam = nx.Beam(polarization = 0)\n", "\n", "exp = nx.Experiment(beam = beam,\n", " objects = [sample],\n", " isotope = nx.lib.moessbauer.Fe57)\n", "\n", "velocities = np.linspace(-10, 10, 512)\n", "\n", "moessbauer_spectrum = nx.MoessbauerSpectrum(experiment = exp,\n", " velocity = velocities,\n", " id = \"my first spectrum\")\n", "\n", "moessbauer_spectrum.Calculate()\n", "\n", "moessbauer_spectrum.Plot(velocity = True)" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.11.9" }, "widgets": { "application/vnd.jupyter.widget-state+json": { "state": {}, "version_major": 2, "version_minor": 0 } } }, "nbformat": 4, "nbformat_minor": 5 }