chem_descriptors.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import pandas as pd\n",
    "from rdkit import Chem\n",
    "from rdkit.Chem import Descriptors\n",
    "from tqdm import tqdm\n",
    "import matplotlib.pyplot as plt"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "100%|██████████| 26832/26832 [00:11<00:00, 2364.27it/s]\n"
     ]
    }
   ],
   "source": [
    "df = pd.read_csv('mnx_chem_bioreachable.csv')\n",
    "df = df.dropna(subset='SMILES')\n",
    "smi_list = df['SMILES'].values.tolist()\n",
    "smi_list = list(set(smi_list))\n",
    "mol_list = [Chem.MolFromSmiles(smi) for smi in tqdm(smi_list)]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "100%|██████████| 26832/26832 [00:00<00:00, 189149.30it/s]\n",
      "100%|██████████| 26832/26832 [00:27<00:00, 981.60it/s] \n",
      "100%|██████████| 26832/26832 [00:00<00:00, 720754.98it/s]\n",
      "100%|██████████| 26832/26832 [00:02<00:00, 9583.14it/s]\n",
      "100%|██████████| 26832/26832 [00:00<00:00, 72513.57it/s]\n",
      "100%|██████████| 26832/26832 [00:00<00:00, 966022.02it/s]\n",
      "100%|██████████| 26832/26832 [00:00<00:00, 83871.33it/s]\n",
      "100%|██████████| 26832/26832 [01:28<00:00, 302.41it/s]\n"
     ]
    }
   ],
   "source": [
    "descriptors = {\n",
    "    'Molecular Weight': [Descriptors.ExactMolWt(mol) for mol in tqdm(mol_list)],\n",
    "    'LogP': [Descriptors.MolLogP(mol) for mol in tqdm(mol_list)],\n",
    "    'Molecular MR?': [Descriptors.MolMR(mol) for mol in tqdm(mol_list)],\n",
    "    '# H Acceptors': [Descriptors.NumHAcceptors(mol) for mol in tqdm(mol_list)],\n",
    "    '# H Donors': [Descriptors.NumHDonors(mol) for mol in tqdm(mol_list)],\n",
    "    '# Rings': [Descriptors.RingCount(mol) for mol in tqdm(mol_list)],\n",
    "    'TPSA?': [Descriptors.TPSA(mol) for mol in tqdm(mol_list)],\n",
    "    'QED?': [Descriptors.qed(mol) for mol in tqdm(mol_list)],\n",
    "}"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAEICAYAAACktLTqAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAT/0lEQVR4nO3df7TkdX3f8efLXRcVDGpYG1zAhS6hbtommlsUmxhr1CzKSmNtwkYTNFv22Bxy0qQ9CkdPT9pTczQxaWpLpBvFrQZBQoyyZi21trrYorIYk4C4ybqCXEBZpW6QGBF994/5rgzD3N25d2b2zv3s83HOPTvz+c58v++5XF73e9/fH59UFZKktjxmuQuQJE2e4S5JDTLcJalBhrskNchwl6QGGe6S1CDDXTMpyauTfGLK2/j1JH8wzW0sRZIfT7J3xNc+P8n8tGvSymO4a+KS3J7kwSQnDYx/NkklWb9MpU1Nki1JPjcw9pEFxi453Lqq6oaqOmtCde1I8h8msS6tLIa7puWLwJZDT5L8A+Dxy1fOZCVZNTD0ceAZSdZ2y1cDPww8YWDsHGD30axVxybDXdPyHuAX+p5fCLy7/wVJTkzy7iQHktyR5I1Jhv5MJvl73V7vfUn2JvmZvmWPT/Lb3ToOJvlEN/aolkX3V8ULF9jGHyb5creO3Ul+qG/ZjiRvT7IryQPAP+l/b1XdDewHntcNPQu4lV7o9489BtiT5Lgkb03ypSRfSXJ5ksd323pE3UmeleRPk9zf1fi+wb3xJP86yb1J7knymm5sG/BK4HVJvpFk57DPrTYZ7pqWTwLfl+QZ3V7uzwKD/e3/DJwInAH8BL1fBq8ZXFGS44GPAO8FnkrvL4Lf6wvftwI/CjwXeArwOuC7S6j5w8CZ3TY+A1w5sPzngDcBTwSGHQ/YzcNB/jzghu51/WOfrKoHgbcAPwj8CLABWAf828EVJlkD/DGwo/tsVwE/PfCyH6D3fVwHbAUuS/LkqtrefYbfrKoTqmrzkb4Baofhrmk6tPf+IuDzwF2HFvQF/qVVdX9V3Q78NvDzQ9ZzHnB7Vb2rqh6qqs8AfwS8otvT/0XgV6rqrqr6TlX936r61mKLraorulq+Bfw68MNJTux7yQer6v9U1Xer6m+HrKJ/L/3H6YX7DQNjH08S4CLgV6vqvqq6H/gN4IIh63wOsBp4W1V9u6reD3x64DXfBv59t3wX8A1gIj17rVyrl7sANe099PZmT2egJQOcBKwB7ugbu4Pe3uegpwPPTvL1vrHV3fpPAh4HfGGcQrtfNm8C/jmwlof3/E8CDnaP7zzCanYD70zyZHqh/Mqq+kaSk7uxHwN+t1v/E4CbeznfKwEY7OMDPA24qx55h7/BOr5WVQ/1Pf8b4IQj1KrGGe6amqq6I8kXgZfQaxf0+yq9Pc6nA4fOKDmNvr37PncCH6+qFw0u6Pbc/xb4u8CfDSx+gF6IHnrtKnrBOszPAecDLwRup9fm+H/0Qvd7H2mB9/YWVu1PcjewDfhSVX2jW3RjN3YCvXbVt4BvAj9UVcM+b797gHVJ0hfwpzL6LzNv+3qMsi2jadsKvKCqHugfrKrvANcAb0ryxCRPB36NR/flAT4E/GCSn0/y2O7rHyV5RlV9F7gC+J0kT0uyKsk5SY4D/hJ4XJKXJnks8EbguAXqfCK90P0avV8Iv7HEz3tD9zlu6Bv7RDe2p6q+2dX8+8B/TPJUgCTrkvzUkPXdCHwHuDjJ6iTnA2cvop6v0DumoWOM4a6pqqovVNWeBRb/Mr296/30AvC99IJ6cB33Ay+m15O+G/gyvQOSh4L63wB/AdwE3Ncte0xVHQR+CXgHvb8IHgAWuuDn3fTaQnfR+0vik4v5nH0+Tu+AbP8B1xu6sf5TIF8P7AM+meSvgf/JkD55d/D15fR+SX4deBW9X3ajHlN4J7AxydeTfGAxH0QrW5ysQ1pZknwKuLyq3rXctWh2uecuzbgkP5HkB7q2zIXAPwT++3LXpdnmAVVp9p1F7/jECfQOpL6iqu5Z3pI062zLSFKDbMtIUoNmoi1z0kkn1fr165e7DElaUW6++eavVtXQazdmItzXr1/Pnj0LnS0nSRomyR0LLbMtI0kNMtwlqUGGuyQ1yHCXpAZNJdyTHJ/k5iTnTWP9kqTDGynck1zRTeF1y8D4pm7Ks30Dk/6+nt4VdZKkZTDqnvsOYFP/QHdv7MuAc4GNwJYkG7v5KT9H71ajkqRlMNJ57lW1O8n6geGzgX1VtR8gydX0Jjs4ATieXuB/M8mu7v7Vj9BN3rsN4LTTTlvyB5AkPdo4FzGt45HTfc0Dz66qiwGSvBr46rBgB+gm790OMDc3t+Qb3Ky/5E+Gjt/+5pcudZWStOKNE+4ZMva9kK6qHUdcQbIZ2Lxhw4YxyhjO0Jd0LBvnbJl5enM5HnIKvVlyRlZVO6tq24knnnjkF0uSRjZOuN8EnJnk9CRr6E2Bdt1kypIkjWPUUyGvojdR71lJ5pNsraqHgIuB64HbgGuq6tbFbDzJ5iTbDx48uNi6JUmHMerZMlsWGN8F7FrqxqtqJ7Bzbm7uoqWuQ5L0aN5+QJIatKzhbltGkqZjWcPds2UkaTpsy0hSg2zLSFKDbMtIUoNsy0hSgwx3SWqQPXdJapA9d0lqkG0ZSWqQ4S5JDTLcJalBHlCVpAaNM83e2Jbjlr8LTb8HTsEnqR22ZSSpQYa7JDXIcJekBhnuktQgz5aRpAZ5+wFJapBtGUlqkOEuSQ0y3CWpQct6heqsWejqVa9clbTSuOcuSQ0y3CWpQYa7JDXIi5gkqUFexCRJDbItI0kNMtwlqUGGuyQ1yHCXpAYZ7pLUIMNdkhpkuEtSgwx3SWqQ4S5JDTLcJalBEw/3JM9IcnmSa5P8y0mvX5J0ZCOFe5Irktyb5JaB8U1J9ibZl+QSgKq6rapeC/wMMDf5kiVJRzLqnvsOYFP/QJJVwGXAucBGYEuSjd2ylwGfAD46sUolSSMbKdyrajdw38Dw2cC+qtpfVQ8CVwPnd6+/rqqeC7xyoXUm2ZZkT5I9Bw4cWFr1kqShxplDdR1wZ9/zeeDZSZ4PvBw4Dti10JurajuwHWBubq7GqGPqnFtV0kozTrhnyFhV1ceAj420gmQzsHnDhg1jlCFJGjTO2TLzwKl9z08B7l7MCpysQ5KmY5xwvwk4M8npSdYAFwDXTaYsSdI4Rj0V8irgRuCsJPNJtlbVQ8DFwPXAbcA1VXXrYjbuHKqSNB0j9dyrassC47s4zEHTEda7E9g5Nzd30VLXIUl6tGW9/YB77pI0Hcsa7h5QlaTp8MZhktQgw12SGmTPXZIaZM9dkhpkW0aSGmS4S1KD7LlLUoPsuUtSg2zLSFKDDHdJapDhLkkN8oCqJDXIA6qS1CDbMpLUIMNdkho00kxMGm79JX8ydPz2N7/0KFciSY/knrskNcizZSSpQZ4tI0kNsuc+BfbiJS03e+6S1CDDXZIaZLhLUoMMd0lqkOEuSQ0y3CWpQV7EJEkN8iImSWqQbRlJapDhLkkNMtwlqUGGuyQ1yHCXpAYZ7pLUIMNdkhrk/dyPIu/zLulocc9dkho0lXBP8k+T/H6SDyZ58TS2IUla2MjhnuSKJPcmuWVgfFOSvUn2JbkEoKo+UFUXAa8GfnaiFUuSjmgxe+47gE39A0lWAZcB5wIbgS1JNva95I3dcknSUTRyuFfVbuC+geGzgX1Vtb+qHgSuBs5Pz1uAD1fVZ4atL8m2JHuS7Dlw4MBS65ckDTFuz30dcGff8/lu7JeBFwKvSPLaYW+squ1VNVdVc2vXrh2zDElSv3FPhcyQsaqqtwFvG3PdkqQlGnfPfR44te/5KcDdo77ZyTokaTrGDfebgDOTnJ5kDXABcN2ob3ayDkmajsWcCnkVcCNwVpL5JFur6iHgYuB64Dbgmqq6dRHrdM9dkqZg5J57VW1ZYHwXsGspG6+qncDOubm5i5byfknScN5+QJIatKzhbltGkqZjWcPdA6qSNB3e8ncGeCtgSZNmz12SGmTPXZIaZM9dkhpkW0aSGmRbRpIaZFtGkhpkW0aSGmS4S1KDDHdJapAHVCWpQct6+wFv+bs03q5A0pF4b5kZtlCIS9KR2HOXpAYZ7pLUIMNdkhrk2TKS1CBvPyBJDbItI0kNMtwlqUGGuyQ1yHCXpAYZ7pLUIG8/cAzwXjTSscdwb4j3opF0iBcxSVKDvIhJkhrkAVVJapDhLkkNMtwlqUGGuyQ1yHCXpAYZ7pLUIC9iOoYt5aInr2qVVgb33CWpQYa7JDXIcJekBk083JOckeSdSa6d9LolSaMZ6YBqkiuA84B7q+rv941vAv4TsAp4R1W9uar2A1sN9zZ5+2BpZRh1z30HsKl/IMkq4DLgXGAjsCXJxolWJ0lakpHCvap2A/cNDJ8N7Kuq/VX1IHA1cP6oG06yLcmeJHsOHDgwcsGSpCMbp+e+Driz7/k8sC7J9ye5HHhmkksXenNVba+quaqaW7t27RhlSJIGjXMRU4aMVVV9DXjtSCtINgObN2zYMEYZkqRB4+y5zwOn9j0/Bbh7MStwsg5Jmo5xwv0m4MwkpydZA1wAXDeZsiRJ4xgp3JNcBdwInJVkPsnWqnoIuBi4HrgNuKaqbl3Mxp1DVZKmY6See1VtWWB8F7BrqRuvqp3Azrm5uYuWug5J0qN5+wFJatCyhrttGUmajmUNd8+WkaTpsC0jSQ2yLSNJDbItI0kNsi0jSQ0y3CWpQePcOGxs3jisHYudxMNJP6TpsucuSQ2yLSNJDTLcJalBhrskNcgDqmqSB2x1rPOAqiQ1yLaMJDXIcJekBhnuktQgw12SGuTZMlrRFjorRpo1R/sMLs+WkaQG2ZaRpAYZ7pLUIMNdkhpkuEtSgwx3SWqQ4S5JDTLcJalBXsSkqVopFxl5i2C1xouYJKlBtmUkqUGGuyQ1yHCXpAYZ7pLUIMNdkhpkuEtSgwx3SWqQ4S5JDTLcJalBhrskNWji95ZJcjzwe8CDwMeq6spJb0OSdHgj7bknuSLJvUluGRjflGRvkn1JLumGXw5cW1UXAS+bcL2SpBGM2pbZAWzqH0iyCrgMOBfYCGxJshE4Bbize9l3JlOmJGkxRmrLVNXuJOsHhs8G9lXVfoAkVwPnA/P0Av6zHOaXR5JtwDaA0047bbF16xgzqVsHL3Y9h3v9QrcD9vbBR7ZSvkcrpc5hxjmguo6H99ChF+rrgPcD/yzJ24GdC725qrZX1VxVza1du3aMMiRJg8Y5oJohY1VVDwCvGWkFTtYhSVMxzp77PHBq3/NTgLsXswIn65Ck6Rgn3G8CzkxyepI1wAXAdZMpS5I0jlFPhbwKuBE4K8l8kq1V9RBwMXA9cBtwTVXdupiNJ9mcZPvBgwcXW7ck6TBGPVtmywLju4BdS914Ve0Eds7NzV201HVIkh7N2w9IUoOWNdxty0jSdCxruHu2jCRNR6pquWsgyQHgjiW+/STgqxMs52ix7qNnJdYM1n00rcSaAZ5eVUOvAp2JcB9Hkj1VNbfcdSyWdR89K7FmsO6jaSXWfCQeUJWkBhnuktSgFsJ9+3IXsETWffSsxJrBuo+mlVjzYa34nrsk6dFa2HOXJA0w3CWpQSs63BeYw3W5ajk1yf9OcluSW5P8Sjf+lCQfSfJX3b9P7nvPpV3te5P8VN/4jyb5i27Z25IMu3f+JGtfleRPk3xoBdX8pCTXJvl89z0/Z4XU/avdz8ctSa5K8rhZrHvYvMmTrDPJcUne141/ashMb5Os+7e6n5M/T/LHSZ40a3VPRVWtyC9gFfAF4AxgDfBnwMZlrOdk4Fnd4ycCf0lvbtnfBC7pxi8B3tI93tjVfBxwevdZVnXLPg2cQ29ClA8D50659l8D3gt8qHu+Emr+b8C/6B6vAZ4063XTm6nsi8Dju+fXAK+exbqB5wHPAm7pG5tYncAvAZd3jy8A3jfFul8MrO4ev2UW657Kz9tyFzDGf8RzgOv7nl8KXLrcdfXV80HgRcBe4ORu7GRg77B66d06+ZzuNZ/vG98C/Ncp1nkK8FHgBTwc7rNe8/fRC8kMjM963YempnwKvTuyfqgLnpmsG1g/EJITq/PQa7rHq+ldHZpp1D2w7KeBK2ex7kl/reS2zEJzuC677k+1ZwKfAv5OVd0D0P371O5lC9W/rns8OD4tvwu8Dvhu39is13wGcAB4V9dOekeS42e97qq6C3gr8CXgHuBgVf2PWa+7zyTr/N57qjc3xEHg+6dW+cN+kd6e+CNqGKhvFutetJUc7kPncD3qVQxIcgLwR8C/qqq/PtxLh4zVYcYnLsl5wL1VdfOobxkydlRr7qym96f326vqmcAD9NoEC5mJurse9fn0WgBPA45P8qrDvWXI2HJ8v49kKXUe9c+Q5A3AQ8CVR6hhpupeqpUc7mPP4TppSR5LL9ivrKr3d8NfSXJyt/xk4N5ufKH657vHg+PT8I+BlyW5HbgaeEGSP5jxmg/VMV9Vn+qeX0sv7Ge97hcCX6yqA1X1beD9wHNXQN2HTLLO770nyWrgROC+aRWe5ELgPOCV1fVUVkLd41jJ4T5Tc7h2R9PfCdxWVb/Tt+g64MLu8YX0evGHxi/ojr6fDpwJfLr7c/f+JM/p1vkLfe+ZqKq6tKpOqar19L5//6uqXjXLNXd1fxm4M8lZ3dBPAp+b9brptWOek+QJ3fZ+kt4UlbNe9yGTrLN/Xa+g97M3rb9QNwGvB15WVX8z8Hlmtu6xLXfTf5wv4CX0zkr5AvCGZa7lx+j9efbnwGe7r5fQ68d9FPir7t+n9L3nDV3te+k72wGYA27plv0XjsIBG+D5PHxAdeZrBn4E2NN9vz8APHmF1P3vgM9323wPvTM1Zq5u4Cp6xwW+TW9vdesk6wQeB/whsI/emSlnTLHuffT65If+v7x81uqexpe3H5CkBq3ktowkaQGGuyQ1yHCXpAYZ7pLUIMNdkhpkuEtSgwx3SWrQ/weTICjOiHx53wAAAABJRU5ErkJggg==",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAEICAYAAACktLTqAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAARm0lEQVR4nO3dfYxldX3H8ffHRTAFS9VdW9wFF7KUspq01gnaB9OmRbuIKz60lg3xkbKlCY1NYyoGarT2D7APUSLVrBVXW4US68NuXYPWSEkbWlkexMWVutAlDBB2lRStGgj67R/3LF6Gmdk7c+feO/ub9yuZzL2/M/fc754789lzv+fc80tVIUlqy1MmXYAkaekZ7pLUIMNdkhpkuEtSgwx3SWqQ4S5JDTLcJalBhrtWjCT7k5y5hOvbnuTRJP+X5KEkX0ryC0u1fmkYhrs0nPdW1XHAOuAAsH2y5Ug9hrtWtCTHJHlfkvu7r/clOaZv+Z8leaBb9gdJKsmGmeupqh8AnwSeP876pbkY7lrpLgFeDPwS8IvAGcClAEk2AX8KnAlsAH5jrpUkOQ44D7h1tOVKgzHctdKdB/xFVR2oqoPAu4HXd8teB3y0qu7o9szfPcvj35bkf4F9wHHAm0ZfsnR4R026AGnCngPc03f/nm7s0LLdfcvuneXxf11Vl46oNmnR3HPXSnc/8Ny++yd1YwAP0DtQesiJ4ypKGpbhrpXmqUmedugLuBq4NMmaJKuBdwL/2P3stcCbk5ye5Ke6ZdIRwXDXSrML+GHf19PotV5uB74O3AL8JUBVfQG4AvgKvZ76jd06HhlvydLCxck6pMEkOR3YAxxTVY9Nuh5pPu65S/NI8uokRyd5BnA5sNNg15HAcJfm94fAQeAu4EfAH022HGkwtmUkqUHuuUtSg5bFh5hWr15d69evn3QZknREufnmm79dVWtmW7Yswn39+vXs3r378D8oSXpcknvmWmZbRpIaZLhLUoMMd0lqkOEuSQ0aSbgnOTbJzUleMYr1S5LmN1C4J7kqyYEke2aMb0pyZ5J9SS7uW/R2elfUkyRNwKB77tuBTf0DSVYBVwJnARuBLUk2drPLfwN4cAnrlCQtwEDnuVfVDUnWzxg+A9hXVXcDJLkGOIfeVGPH0gv8HybZVVU/nrnOJFuBrQAnnXTSov8BkqQnG+ZDTGt54rRj08CLquoigCRvAr49W7ADVNU2YBvA1NSUF7iRlon1F39+1vH9l5095ko0jGHCPbOMPR7SVbX9sCtINgObN2zYMEQZkqSZhjlbZponzim5jp/MPTmQqtpZVVuPP/74IcqQJM00TLjfBJya5OQkRwPnAjuWpixJ0jAGPRXyanrzR56WZDrJ+d1sNBcB1wF7gWur6o6FPHmSzUm2PfzwwwutW5I0j0HPltkyx/guehMOL0pV7QR2Tk1NXbDYdUiSnszLD0hSgyZ6PXfPlpGOHJ4ieWSZ6J67Z8tI0mjYlpGkBk003D1bRpJGw7aMJDXItowkNchwl6QGeSqktELNdWqj2mDPXZIaZFtGkhpkuEtSgwx3SWqQH2KSpAZ5QFWSGmRbRpIaZLhLUoMMd0lqkOEuSQ3ybBlJapBny0hSg2zLSFKDDHdJapDhLkkNMtwlqUGGuyQ1yHCXpAYZ7pLUID/EJEkNmugE2VW1E9g5NTV1wSTrkLR48020vf+ys8dYifrZlpGkBhnuktQgw12SGmS4S1KDDHdJapDhLkkNMtwlqUGGuyQ1yHCXpAYZ7pLUoCW//ECS04G3AquBL1fVB5f6OSQdGea6NIGXJRi9gcI9yVXAK4ADVfX8vvFNwPuBVcDfV9VlVbUXuDDJU4APj6BmSQsw37Vf1K5B2zLbgU39A0lWAVcCZwEbgS1JNnbLXgn8O/DlJatUkjSwgcK9qm4AHpoxfAawr6rurqpHgWuAc7qf31FVvwqcN9c6k2xNsjvJ7oMHDy6ueknSrIbpua8F7u27Pw28KMlvAq8BjgF2zfXgqtoGbAOYmpqqIeqQJM0wTLhnlrGqquuB6wdaQbIZ2Lxhw4YhypAkzTTMqZDTwIl999cB9y9kBVW1s6q2Hn/88UOUIUmaaZhwvwk4NcnJSY4GzgV2LE1ZkqRhDBTuSa4GbgROSzKd5Pyqegy4CLgO2AtcW1V3LOTJnUNVkkZjoJ57VW2ZY3wX8xw0HWC9zqEqSSPg5QckqUETDXfbMpI0GhMNd8+WkaTRsC0jSQ2yLSNJDbItI0kNsi0jSQ0y3CWpQfbcJalB9twlqUG2ZSSpQYa7JDXIcJekBnlAVZIa5AFVSWrQMHOoSlpG1l/8+UmXoGXEnrskNchwl6QG2ZaRNHZztZD2X3b2mCtpl2fLSFKDPFtGkhpkz12SGmS4S1KDDHdJapDhLkkNMtwlqUGGuyQ1yHCXpAb5ISZJapAfYpKkBtmWkaQGGe6S1CDDXZIaZLhLUoMMd0lqkJN1SEcY50rVINxzl6QGGe6S1CDDXZIaNJJwT/KqJB9O8rkkLxvFc0iS5jZwuCe5KsmBJHtmjG9KcmeSfUkuBqiqz1bVBcCbgN9f0oolSYe1kLNltgMfAD5+aCDJKuBK4KXANHBTkh1V9Y3uRy7tlkvSYc11JtD+y84ecyVHvoH33KvqBuChGcNnAPuq6u6qehS4BjgnPZcDX6iqW5auXEnSIIbtua8F7u27P92N/TFwJvC7SS6c7YFJtibZnWT3wYMHhyxDktRv2A8xZZaxqqorgCvme2BVbQO2AUxNTdWQdUjN8cNKGsawe+7TwIl999cB9w/6YCfrkKTRGDbcbwJOTXJykqOBc4Edgz7YyTokaTQWcirk1cCNwGlJppOcX1WPARcB1wF7gWur6o4FrNM9d0kagYF77lW1ZY7xXcCuxTx5Ve0Edk5NTV2wmMdLkmbn5QckqUETDXfbMpI0GhMNdw+oStJo2JaRpAYZ7pLUIHvuktQge+6S1CDbMpLUIMNdkho07FUhh5JkM7B5w4YNkyxD0jLnJB4LZ89dkhpkW0aSGmS4S1KD7LlLE+aMSxoFe+6S1CDbMpLUIMNdkhpkuEtSgwx3SWqQ4S5JDfKSv5LUIE+FlKQG2ZaRpAYZ7pLUIMNdkhpkuEtSgwx3SWrQRK8KKa0kXv1R4+R57pLUIM9zl6QG2XOXpAbZc5d0xJrrOMb+y84ecyXLj3vuktQgw12SGmS4S1KDDHdJapDhLkkNMtwlqUGGuyQ1yHCXpAYtebgnOSXJR5J8aqnXLUkazEDhnuSqJAeS7JkxvinJnUn2JbkYoKrurqrzR1GsJGkwg+65bwc29Q8kWQVcCZwFbAS2JNm4pNVJkhZloHCvqhuAh2YMnwHs6/bUHwWuAc4Z9ImTbE2yO8nugwcPDlywJOnwhum5rwXu7bs/DaxN8qwkHwJekOQdcz24qrZV1VRVTa1Zs2aIMiRJMw1zVcjMMlZV9R3gwiHWK0ka0jDhPg2c2Hd/HXD/QlaQZDOwecOGDUOUIS0vTqc3eV4KeLi2zE3AqUlOTnI0cC6wYyErcCYmSRqNQU+FvBq4ETgtyXSS86vqMeAi4DpgL3BtVd2xkCd3DlVJGo2B2jJVtWWO8V3ArsU+eVXtBHZOTU1dsNh1SJKezMsPSFKDJjqHqgdUdSTzwGk7WjwAO9E9dw+oStJo2JaRpAZNNNw9W0aSRsO2jCQ1yLaMJDXIcJekBtlzl6QG2XOXpAbZlpGkBhnuktQgw12SGuS1ZSRpDMZ9/RoPqEpSg2zLSFKDDHdJapDhLkkNMtwlqUFH/NkyLc6gouXD2ZZ0pPJsGUlqkG0ZSWqQ4S5JDTLcJalBhrskNchwl6QGGe6S1CDDXZIadMR/iEkrw0I/TOSH2LTS+SEmSWqQbRlJapDhLkkNMtwlqUGGuyQ1yHCXpAYZ7pLUIMNdkhpkuEtSgwx3SWqQ4S5JDVrya8skORb4O+BR4Pqq+sRSP4ckaX4D7bknuSrJgSR7ZoxvSnJnkn1JLu6GXwN8qqouAF65xPVKkgYwaFtmO7CpfyDJKuBK4CxgI7AlyUZgHXBv92M/WpoyJUkLMVBbpqpuSLJ+xvAZwL6quhsgyTXAOcA0vYC/jXn+80iyFdgKcNJJJy207kWb79Kxk7pM7Fw1HUmXrV3ov2FS/+aFXjpYbVno679Uv9eTMMwB1bX8ZA8deqG+Fvg08NokHwR2zvXgqtpWVVNVNbVmzZohypAkzTTMAdXMMlZV9X3gzQOtwMk6JGkkhtlznwZO7Lu/Drh/IStwsg5JGo1hwv0m4NQkJyc5GjgX2LE0ZUmShjHoqZBXAzcCpyWZTnJ+VT0GXARcB+wFrq2qOxby5Ek2J9n28MMPL7RuSdI8Bj1bZssc47uAXYt98qraCeycmpq6YLHrkCQ9mZcfkKQGTTTcbctI0mhMNNw9W0aSRiNVNekaSHIQuGfSdXRWA9+edBFzWM61wfKuz9oWbznXt9Jre25Vzfop0GUR7stJkt1VNTXpOmaznGuD5V2ftS3ecq7P2ubmAVVJapDhLkkNMtyfbNukC5jHcq4Nlnd91rZ4y7k+a5uDPXdJapB77pLUIMNdkhq0osM9ye8luSPJj5NM9Y2/NMnNSb7eff+tvmXXd/PG3tZ9PXuctXXL3tHNW3tnkt/pG39hV/O+JFckme2a+0td5z/1bYv9SW7rxtcn+WHfsg+NupY56ntXkvv66nh537JZt+MYa/urJN9McnuSzyT5mW58uWy72eZInogkJyb5SpK93d/FW7vxOV/fCdS4v/v7uy3J7m7smUm+lORb3fdnjK2gqlqxX8DpwGnA9cBU3/gLgOd0t58P3Ne37Ak/O4HaNgJfA44BTgbuAlZ1y74K/Aq9iVS+AJw15u35N8A7u9vrgT3L4DV+F/C2Wcbn3I5jrO1lwFHd7cuBy5fLtgNWddvkFODoblttnGA9JwC/3N1+OvDf3Ws46+s7oRr3A6tnjL0XuLi7ffGh13gcXyt6z72q9lbVnbOM31pVhyYeuQN4WpJjlkNt9OapvaaqHqmq/wH2AWckOQH46aq6sXq/SR8HXjWuert3Ca8Drh7Xcw5p1u04zgKq6ovVu3Q2wH/Sm/BmuXh8juSqehQ4NEfyRFTVA1V1S3f7e/QuM752UvUswDnAx7rbH2OMf5MrOtwH9Frg1qp6pG/so91brz8fR+tjhrnmrl3b3Z45Pi4vAR6sqm/1jZ2c5NYk/5bkJWOsZaaLutbHVX1vi+fajpPyFnrvtg6Z9LZbbtvncUnW03t3/V/d0Gyv7yQU8MWulbu1G/vZqnoAev9BASNp485mmDlUjwhJ/hX4uVkWXVJVnzvMY59H7+3yy/qGz6uq+5I8Hfhn4PX09pLHVdusc9fOMz60AevcwhP32h8ATqqq7yR5IfDZJM+rqu8uRU2D1gd8EHgPvW3xHnqto7cwwu01aG2Htl2SS4DHgE90y8a27eYxlu2zUEmOo/d39ydV9d0kc72+k/BrVXV/dxzuS0m+OaE6gBUQ7lV15mIel2Qd8BngDVV1V9/67uu+fy/JJ+m9fV1UuC+ytrnmrp3miW/rFzyn7VwOV2eSo4DXAC/se8wjwCPd7ZuT3AX8PLB7KWpaSH19dX4Y+Jfu7tBzAA9igG33RuAVwG937bSxbrt5jGX7LESSp9IL9k9U1acBqurBvuX9r+/YHWrlVtWBJJ+hlw0PJjmhqh7oWqcHxlWPbZlZdGctfB54R1X9R9/4UUlWd7efSu+Pcs+Yy9sBnJvkmCQnA6cCX+3e8n0vyYu7VtEbgHnfmSyhM4FvVtXjbaEka5Ks6m6f0tV595jqeVz3B3XIq/nJ6zXrdhxzbZuAtwOvrKof9I0vh223rOZI7n6nPwLsraq/7Ruf6/UdqyTHdu/mSXIsvXf7e+htszd2P/ZGxvc3ueLPlnk1vT2UR4AHgeu68UuB7wO39X09GzgWuBm4nd6B1vczojMs5qqtW3YJvTMZ7qTvjBhgit4v1F3AB+g+gTyG7bgduHDG2Gu7bfQ14BZg84Re438Avt69ZjuAEw63HcdY2z56fe1Dv2MfWmbb7uX0zkq5i14baew19NXy6/RaL7f3ba+Xz/f6jrm+U7rX62vda3dJN/4s4MvAt7rvzxxXTV5+QJIaZFtGkhpkuEtSgwx3SWqQ4S5JDTLcJalBhrskNchwl6QG/T8Qs8jpklVfzwAAAABJRU5ErkJggg==",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAEICAYAAACktLTqAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAASM0lEQVR4nO3df7BcZ13H8ffHSCuWWoUEB9OWBFIrcRwF1vJDR9FBTAuxKAoNKqCZRtA6Oo6jQfmj/+CIin+gBSdKrUVorfxsJExlHEoHqNqUKdASK2lJp5ciTS3EWn+Uwtc/9kSW7e7N7t17u3efvF8zd+7d55w9+zzd9HPP/Z6zz5OqQpLUlm+YdwckSavPcJekBhnuktQgw12SGmS4S1KDDHdJapDhroWV5FVJPrLGr3Fpkr9ey9eQ1oLhrrlIciTJQ0k2DrXfkqSSbJlT19ZUkuu78X3vUPt7u/bndY8vTfLlJP+Z5EtJPpbkOQP7Py/JR5McS/KZJBc8uiPReme4a54+C+w6/iDJ9wCPnV93VleSDWM2/SvwioH9ngA8Gzg6tN/fVNXjgI3Ah4C/Hdh2NvA64And96uSfOMqdV0NMNw1T29jIOSAVwJXDu6Q5IwkVyY5muSuJK9LMvLfbZLvSvLBJPcnuT3JSwe2PTbJG7tjHEvyka7teUmWho5zJMnzx7zG3yb5t+4YNyT57oFtVyR5S5IDSR4EfmTMuN8OvGwg/HcB7wEeGrVzVT3cPWdzkk1d25VV9aFu2w3A6cA3j3k9nYQMd83TPwLfkuRpXdC9DBiub/8JcAbwFOCH6f8y+IXhAyU5Dfgg8A7gifQD880D4ftHwDOB5wKPB34L+OoK+vwB4JzuNT5OP3QHvRx4Pf2wHXc94B7g08ALusevYOiX2qAkp3T7/DvwxaFtG4A3A++pqv+YZiBqm3/Gad6On71/GPgX4HPHNwwE/tOr6gHggSRvBH4eeOvQcV4EHKmqv+wefzzJu4CfTnII+EXg2VV1/Pgf615jqs5W1eUD/bsU+GKSM6rqWNf8vqr6aPfz/yxzqCuBVyS5E/jWqrpxRF9emuRF9H9RfAl4SXemPuhP6P+i+fGpBqLmGe6at7fRLyts5ZFnrxuBU4C7BtruAjaPOM6TgWcl+dJA2zd2x98IfBNwxywd7X7ZvB74GWATXzvz3wgcD/e7Jzzcu4E30j8bf9uYfa6pqp/rLjq/i/5fHtcP9OexwC8BW6vqPycfiU4GhrvmqqruSvJZ4AJg99Dm+4Av0w/uT3dtZzNwdj/gbuDDVfVjwxu6Gv3/AE8FPjG0+UEGatVdgG8a092XAxcCzweO0C8XfREYPOWeaJrVqvqvJB8AXtP1a7l970vyS8BNSd5RVZ/vNj2Rfmn1nkleUycXa+5aD3YDP1pVDw42VtVXgGuA1yc5PcmTgd/gkXV5gL8DvjPJzyd5TPf1/UmeVlVfBS4H/jjJdyTZkOQ5SU6lf+fKNyV5YZLH0L/z5NQx/Twd+F/6Z9vfDPzejOP+HeCHq+rIiXasqn8BrqN/reC4JeBJI0o1kuGu+auqO6rq4JjNv0r/7PpO+hco30E/qIeP8QD9C5QX0T+T/TfgDXwtqH8T+BRwE3B/t+0bulr5LwN/Qf8vggfph+YoV9IvC32O/l8S/zjNOEf0+Z6qmuZDWH8I7EnyxO7xZuDwMrdc6iQWF+uQpPZ45i5JDTLcJalBhrskNchwl6QGrYv73Ddu3FhbtmyZdzckaaHcfPPN91XVyM9lrItw37JlCwcPjrsTTpI0SpK7xm2zLCNJDTLcJalBhrskNWhNwj3JaUlu7qYrlSQ9yiYK9ySXJ7k3ya1D7Tu6FW8OJ9k7sOm36U/4JEmag0nP3K8Adgw2dJMVXQacD2wHdiXZ3i1P9mngC6vYT0nSFCa6FbKqbhixGv15wOGquhMgydX057p+HHAa/cD/7yQHuilXv06SPcAegLPPPnvFA5AkPdIs97lv5utXnVkCnlVVlwAkeRVw36hgB6iqfcA+gF6v59SUkrSKZgn3UYtP/n9IV9UVJzxAshPYuW3bthm6MdqWve8f2X7k91+46q8lSevNLHfLLAFnDTw+kymX+6qq/VW154wzzpihG5KkYbOE+03AOUm2JjmF/go4165OtyRJs5j0VsirgBuBc5MsJdndrdt4Cf11HQ/RX6n9tmlePMnOJPuOHTt24p0lSRNbF8vs9Xq9WunEYeNq6+NYc5fUiiQ3V1Vv1DanH5CkBs013C3LSNLamGu4e7eMJK0NyzKS1CDLMpLUIMsyktQgyzKS1CDDXZIaZM1dkhpkzV2SGmRZRpIaZLhLUoNmWaxjIS030ZiTiklqhRdUJalBXlCVpAZZc5ekBhnuktQgw12SGmS4S1KDvFtGkho01/vcq2o/sL/X6108z34cN+4eeO9/l7RoLMtIUoMMd0lqkOEuSQ0y3CWpQYa7JDXIcJekBhnuktQgP8QkSQ1yyl9JapBlGUlqkOEuSQ0y3CWpQYa7JDXIcJekBs11yt9F4VTAkhaNZ+6S1CDDXZIaZLhLUoMMd0lq0KqHe5KnJfmzJO9M8prVPr4k6cQmCvcklye5N8mtQ+07ktye5HCSvQBVdaiqXg28FOitfpclSScy6Zn7FcCOwYYkG4DLgPOB7cCuJNu7bT8BfAT4h1XrqSRpYhOFe1XdANw/1HwecLiq7qyqh4CrgQu7/a+tqucCPzvumEn2JDmY5ODRo0dX1ntJ0kizfIhpM3D3wOMl4FlJngf8FHAqcGDck6tqH7APoNfr1Qz9kCQNmSXcM6Ktqup64PqJDpDsBHZu27Zthm5IkobNcrfMEnDWwOMzgXumOYCLdUjS2pgl3G8CzkmyNckpwEXAtavTLUnSLCa9FfIq4Ebg3CRLSXZX1cPAJcB1wCHgmqq6bZoXdw1VSVobE9Xcq2rXmPYDLHPRdILj7gf293q9i1d6DEnSIzn9gCQ1aK7hbllGktbGXMPdu2UkaW1YlpGkBlmWkaQGWZaRpAZZlpGkBhnuktSgWSYOm9miTxy2Ze/7R7Yf+f0XPso9kaSvZ81dkhpkWUaSGmS4S1KDDHdJapAfYpKkBnlBVZIaZFlGkhpkuEtSgwx3SWqQ4S5JDXL6gTXgtASS5s27ZSSpQZZlJKlBhrskNchwl6QGzfWC6snGC62SHi2euUtSgwx3SWqQ4S5JDXLKX0lqkB9ikqQGWZaRpAYZ7pLUIO9zXwe8/13SavPMXZIaZLhLUoMMd0lqkOEuSQ0y3CWpQYa7JDXIcJekBq1JuCd5cZI/T/K+JC9Yi9eQJI03cbgnuTzJvUluHWrfkeT2JIeT7AWoqvdW1cXAq4CXrWqPJUknNM2Z+xXAjsGGJBuAy4Dzge3AriTbB3Z5XbddkvQomjjcq+oG4P6h5vOAw1V1Z1U9BFwNXJi+NwAfqKqPjzpekj1JDiY5ePTo0ZX2X5I0wqw1983A3QOPl7q2XwWeD/x0klePemJV7auqXlX1Nm3aNGM3JEmDZp04LCPaqqreBLzphE9OdgI7t23bNmM3JEmDZj1zXwLOGnh8JnDPpE92sQ5JWhuzhvtNwDlJtiY5BbgIuHb2bkmSZjHNrZBXATcC5yZZSrK7qh4GLgGuAw4B11TVbVMc0zVUJWkNTFxzr6pdY9oPAAdW8uJVtR/Y3+v1Ll7J8yVJozn9gCQ1aK7hbllGktbGXNdQtSyzPNdWlbRSlmUkqUGWZSSpQXMNdz/EJElrw7KMJDXIcJekBllzl6QGWXOXpAZZlpGkBs31Q0xaGT/cJOlEPHOXpAZ5QVWSGuQFVUlqkDX3hliLl3ScNXdJapDhLkkNMtwlqUHeLSNJDfJuGUlqkGUZSWqQt0LqEcbdUgneViktCs/cJalBhrskNchwl6QGWXM/CTgtgXTy8cxdkho01zP3JDuBndu2bZtnN05ay90VM+1z/CtAWl/8EJMkNciyjCQ1yHCXpAYZ7pLUIMNdkhpkuEtSgwx3SWqQ4S5JDTLcJalBhrskNciJw7QqnJZAWl9W/cw9yVOSvDXJO1f72JKkyUwU7kkuT3JvkluH2nckuT3J4SR7AarqzqravRadlSRNZtIz9yuAHYMNSTYAlwHnA9uBXUm2r2rvJEkrMlG4V9UNwP1DzecBh7sz9YeAq4ELJ33hJHuSHExy8OjRoxN3WJJ0YrPU3DcDdw88XgI2J3lCkj8Dnp7kteOeXFX7qqpXVb1NmzbN0A1J0rBZ7pbJiLaqqn8HXj3RAVysQ5LWxCxn7kvAWQOPzwTumeYALtYhSWtjlnC/CTgnydYkpwAXAdeuTrckSbOY9FbIq4AbgXOTLCXZXVUPA5cA1wGHgGuq6rZpXjzJziT7jh07Nm2/JUnLmKjmXlW7xrQfAA6s9MWraj+wv9frXbzSY0iSHsm5ZSSpQXMNd8sykrQ25hru3i0jSWvDsowkNciyjCQ1yLKMJDXIsowkNchwl6QGzXWZPScOO3mNW5ZvHJfrk6ZjzV2SGmRZRpIaZLhLUoMMd0lqkB9ikqQGeUFVkhpkWUaSGmS4S1KDDHdJapDhLkkNcvoBLYRx0xVMOy3Bah1Hmtaj/W/Pu2UkqUGWZSSpQYa7JDXIcJekBhnuktQgw12SGmS4S1KDDHdJapAfYtKamnat1HlZzX76gSitB36ISZIaZFlGkhpkuEtSgwx3SWqQ4S5JDTLcJalBhrskNchwl6QGGe6S1CDDXZIaZLhLUoNWfW6ZJKcBbwYeAq6vqrev9mtIkpY30Zl7ksuT3Jvk1qH2HUluT3I4yd6u+aeAd1bVxcBPrHJ/JUkTmLQscwWwY7AhyQbgMuB8YDuwK8l24Ezg7m63r6xONyVJ05ioLFNVNyTZMtR8HnC4qu4ESHI1cCGwRD/gb2GZXx5J9gB7AM4+++xp+y2tW9NOHzxuiuBxx2lhSuFFGdui9HOUWS6obuZrZ+jQD/XNwLuBlyR5C7B/3JOral9V9aqqt2nTphm6IUkaNssF1Yxoq6p6EPiFiQ7gYh2StCZmOXNfAs4aeHwmcM80B3CxDklaG7OE+03AOUm2JjkFuAi4dnW6JUmaxaS3Ql4F3Aicm2Qpye6qehi4BLgOOARcU1W3TfPiSXYm2Xfs2LFp+y1JWsakd8vsGtN+ADiw0hevqv3A/l6vd/FKjyFJeiSnH5CkBs013C3LSNLamGu4e7eMJK2NVNW8+0CSo8BdK3z6RuC+VezOeuLYFpNjW0yLOLYnV9XIT4Gui3CfRZKDVdWbdz/WgmNbTI5tMbU2Ni+oSlKDDHdJalAL4b5v3h1YQ45tMTm2xdTU2Ba+5i5JeqQWztwlSUMMd0lq0EKH+5g1XBdGkiNJPpXkliQHu7bHJ/lgks90379tYP/XdmO9PcmPz6/no41aa3cl40nyzO6/y+Ekb0oyau2AR9WYsV2a5HPd+3dLkgsGti3E2JKcleRDSQ4luS3Jr3XtC/++LTO2hX/fJlJVC/kFbADuAJ4CnAJ8Atg+735NOYYjwMahtj8A9nY/7wXe0P28vRvjqcDWbuwb5j2Gob7/EPAM4NZZxgP8M/Ac+gvCfAA4f52O7VLgN0fsuzBjA54EPKP7+XTgX7v+L/z7tszYFv59m+Rrkc/c/38N16p6CDi+huuiuxD4q+7nvwJePNB+dVX9b1V9FjhM/7/BulFVNwD3DzVPNZ4kTwK+papurP7/VVcOPGduxoxtnIUZW1V9vqo+3v38AP3puzfTwPu2zNjGWZixTWKRw33cGq6LpIC/T3Jzt2A4wLdX1eeh/48TeGLXvqjjnXY8m7ufh9vXq0uSfLIr2xwvXSzk2JJsAZ4O/BONvW9DY4OG3rdxFjncR67h+qj3YjY/UFXPAM4HfiXJDy2zbwvjHTRuPIs0zrcATwW+D/g88MaufeHGluRxwLuAX6+q/1hu1xFtiza2Zt635SxyuM+8huu8VdU93fd7gffQL7N8ofszkO77vd3uizreacez1P083L7uVNUXquorVfVV4M/5WplsocaW5DH0w+/tVfXurrmJ923U2Fp5305kkcN9oddwTXJaktOP/wy8ALiV/hhe2e32SuB93c/XAhclOTXJVuAc+hd51rupxtOVAB5I8uzujoRXDDxnXTkefp2fpP/+wQKNrevHW4FDVfXHA5sW/n0bN7YW3reJzPuK7ixfwAX0r4DfAfzuvPszZd+fQv/K/CeA2473H3gC8A/AZ7rvjx94zu92Y72ddXi1HriK/p+5X6Z/trN7JeMBevT/h7sD+FO6T1Kvw7G9DfgU8En6wfCkRRsb8IP0SwyfBG7pvi5o4X1bZmwL/75N8uX0A5LUoEUuy0iSxjDcJalBhrskNchwl6QGGe6S1CDDXZIaZLhLUoP+D1zqGNCYjzrTAAAAAElFTkSuQmCC",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAEICAYAAABcVE8dAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAStklEQVR4nO3dfbBcdX3H8ffH8KTgYIHUakJMEIYxYxX1Dj6OxY6jQUWstZZUR1QG6rRYbWs1VKfVTh/QqY7VUp2IDBYslPGhEo2ljpZSlakEBQUpNSKaGDVBJNbWisC3f+xB1+u9N3vvbnb38nu/Znbu7m93z/neX3I/9+z3nHtOqgpJ0n3f/SZdgCRpPAx8SWqEgS9JjTDwJakRBr4kNcLAl6RGGPiS1AgDXxOX5HNJjktyTJLP7+O1leTYWWNvTHLxPt53aJIfJNk6ipqXKsna7ns4YJJ1qE0GviYqyYHAw4DtwOOABQN/CC8AfgQ8I8lD9tM69jt/UWgYBr4m7ZHAl6v3J98z7L/APx14N/BF4EX9TyR5SpLPJrkjyY4kL+3G75/krUm+nmRvkk8nuX/33BP63nN9kpP6lndlkr/qPrnsTfKRJEd0T1/Vfb2j+8TxxCT3S/KGbj27k/x9ksO7Zd37ieCMJN8APpXkkCQXJ/lut/5rkjx4P82b7kMMfE1EkpcluQP4DPDE7v4fAm/uQmzdCNe1BjgJeH93e8ms5z4OvBNYCZwAXNc9/df0PnU8CTgCeC1wT5JVwMeAP+/GXwN8MMnKvtW+BHg58FDgLuAd3fhTu68PqqrDqupq4KXd7WnAMcBhwN/O+jZ+BXgE8Ex6v7wOB44GjgReAfxwUZOiNlWVN28TuwH/Ti9k19AL2uzj9QV8H7ij7/Z/wMULvOcNwHXd/YcCdwOP6R6fA3x4jvfcj16IPnqO514HXDRr7Arg9O7+lcC5fc+tB+4EVgBru+/hgL7nPwn8Tt/j44EfAwf0vf6YvudfDnwWeNSk//28La+bW/gauyRHdFvxe+ltPV8J3Ewv6L6X5NX7WMRjq+pB996Ac/fx+pfQ27KnqnYB/0ZvKxl6W8lfneM9RwGHzPPcw4Df6L6HO7pPJ08B+vcN7Oi7/3XgwG6Zc3lo95r+1x8A9Ldp+pd3Eb1fMJcm2ZXkLd2+EGlBBr7Grqpu74L6t4Hzu/v/DJzShfjbR7WuJE8CjgPOSfLtJN8GHg9s7HaA7gAePsdbb6P3yWGu53bQ28J/UN/t0Krq/8VzdN/9NfS22G+jt7U+2y56v0T6X38X8J2+sZ+8r6p+XFVvqqr19H5hPoe+NpU0HwNfk9R/VM5jgGv3wzpOBz5Br61yQnd7JPAA4GR6W/5PT/LCJAckOTLJCVV1D3AB8LYkD02yotvBejBwMXBKkmd244ckOSnJ6r71vjjJ+iQPAP4M+EBV3Q3sAe6h16u/1yXA7ydZl+Qw4C+Bf6yqu+b6hpI8LckvJ1lBr731Y3ptKmlBBr4m6XHA55McCdxdVd8b5cKTHAK8EHhnVX277/Y1em2R06vqG8Cz6O0wvp3efoRHd4t4DfAl4JruuTcD96uqHcCpwB/TC/AdwB/xsz9PFwEXAt+m1xr6PYCq+l/gL4DPdO2gJ9D7xXIRvSN4vkbvk8UrF/jWfgn4AL2wv4lei2rBv0OQoNtBJml0klxJbyfy+ZOuRernFr4kNcLAl6RG2NKRpEa4hS9JjZiKEzEdddRRtXbt2kmXIUnLyrXXXntbVa3c9yt7piLw165dy7Zt2yZdhiQtK0m+vu9X/ZQtHUlqhIEvSY0w8CWpEQa+JDViooGf5JQkm/fu3TvJMiSpCRMN/KraUlVnHX744ZMsQ5KaYEtHkhph4EtSI6biD6+GsXbTx+Ycv/XcZ4+5Ekmabm7hS1IjDHxJaoSBL0mNMPAlqREGviQ1wsCXpEbsl8BPcmiSa5M8Z38sX5K0eAMFfpILkuxOcsOs8Q1Jbk6yPcmmvqdeB1w2ykIlScMZdAv/QmBD/0CSFcB5wMnAemBjkvVJng58GfjOCOuUJA1poL+0raqrkqydNXwisL2qbgFIcilwKnAYcCi9XwI/TLK1qu4ZXcmSpKUY5tQKq4AdfY93Ao+vqrMBkrwUuG2+sE9yFnAWwJo1a4YoQ5I0iGF22maOsfrJnaoLq+qj8725qjZX1UxVzaxcOfBF1yVJSzRM4O8Eju57vBrYtZgFeAEUSRqfYQL/GuC4JOuSHAScBly+mAV4ARRJGp9BD8u8BLgaOD7JziRnVNVdwNnAFcBNwGVVdeNiVu4WviSNz6BH6WycZ3wrsHWpK6+qLcCWmZmZM5e6DEnSYDy1giQ1YqKBb0tHksZnooHvTltJGh9bOpLUCFs6ktQIWzqS1IhhzqUz1dZu+tic47ee++wxVyJJ08EeviQ1wh6+JDXCHr4kNcKWjiQ1wsCXpEbYw5ekRtjDl6RG2NKRpEYY+JLUCANfkhph4EtSIzxKR5Ia4VE6ktSI++zZMucz31k0wTNpSrpvs4cvSY0w8CWpEQa+JDXCwJekRhj4ktQIj8OXpEZ4HL4kNcKWjiQ1wsCXpEYY+JLUCANfkhph4EtSIwx8SWqEgS9JjTDwJakRBr4kNWLkF0BJ8gjgVcBRwCer6l2jXsf+Mt/FUbwwiqT7goG28JNckGR3khtmjW9IcnOS7Uk2AVTVTVX1CuCFwMzoS5YkLcWgLZ0LgQ39A0lWAOcBJwPrgY1J1nfPPRf4NPDJkVUqSRrKQIFfVVcBt88aPhHYXlW3VNWdwKXAqd3rL6+qJwEvmm+ZSc5Ksi3Jtj179iyteknSwIbp4a8CdvQ93gk8PslJwPOBg4Gt8725qjYDmwFmZmZqiDokSQMYJvAzx1hV1ZXAlUMsV5K0HwxzWOZO4Oi+x6uBXYtZgBdAkaTxGSbwrwGOS7IuyUHAacDli1mAF0CRpPEZ9LDMS4CrgeOT7ExyRlXdBZwNXAHcBFxWVTcuZuVu4UvS+AzUw6+qjfOMb2WBHbMDLHcLsGVmZubMpS5DkjQYT60gSY2YaODb0pGk8Zlo4LvTVpLGx5aOJDXClo4kNcKWjiQ1wpaOJDXCwJekRtjDl6RG2MOXpEbY0pGkRhj4ktQIe/iS1Ah7+JLUiGEucdiMtZs+Nuf4rec+e8yVSNLS2cOXpEYY+JLUCANfkhrhUTqS1AiP0pGkRtjSkaRGGPiS1AgDX5IaYeBLUiMMfElqhIEvSY3wOHxJaoTH4UtSI2zpSFIjDHxJaoSBL0mNMPAlqRFe8WoIXglL0nLiFr4kNcLAl6RGGPiS1AgDX5IasV8CP8nzkrwnyUeSPGN/rEOStDgDB36SC5LsTnLDrPENSW5Osj3JJoCq+qeqOhN4KfCbI61YkrQki9nCvxDY0D+QZAVwHnAysB7YmGR930ve0D0vSZqwgQO/qq4Cbp81fCKwvapuqao7gUuBU9PzZuDjVfX5uZaX5Kwk25Js27Nnz1LrlyQNaNge/ipgR9/jnd3YK4GnAy9I8oq53lhVm6tqpqpmVq5cOWQZkqR9GfYvbTPHWFXVO4B37PPNySnAKccee+yQZUiS9mXYLfydwNF9j1cDuwZ9s+fDl6TxGXYL/xrguCTrgG8CpwG/NXRVy5zn2JE0jRZzWOYlwNXA8Ul2Jjmjqu4CzgauAG4CLquqGxexTC9xKEljMvAWflVtnGd8K7B1KSuvqi3AlpmZmTOX8n5J0uA8tYIkNWKigW9LR5LGZ6KB71E6kjQ+tnQkqRG2dCSpEbZ0JKkRtnQkqRHD/qXtUFo7l45/gStpkmzpSFIjbOlIUiMMfElqhIEvSY3wOHxJaoQ7bSWpEbZ0JKkRBr4kNcLAl6RGGPiS1AhPrTAFPOWCpHHwKB1JaoQtHUlqhIEvSY0w8CWpEQa+JDXCwJekRkz0sEwtzMM1JY2SZ8uUpEZ4HL4kNcIeviQ1wh7+MmRvX9JSuIUvSY1wC/8+xC1/SQtxC1+SGmHgS1IjDHxJaoSBL0mNMPAlqREjD/wkxyR5b5IPjHrZkqSlG+iwzCQXAM8BdlfVI/vGNwB/A6wAzq+qc6vqFuAMA396eLimJBh8C/9CYEP/QJIVwHnAycB6YGOS9SOtTpI0MgMFflVdBdw+a/hEYHtV3VJVdwKXAqcOuuIkZyXZlmTbnj17Bi5YkrQ0w/TwVwE7+h7vBFYlOTLJu4HHJDlnvjdX1eaqmqmqmZUrVw5RhiRpEMOcWiFzjFVVfRd4xUALSE4BTjn22GOHKEOSNIhhtvB3Akf3PV4N7FrMAjwfviSNzzCBfw1wXJJ1SQ4CTgMuH01ZkqRRG/SwzEuAk4CjkuwE/rSq3pvkbOAKeodlXlBVNy5m5bZ0JsvDNaW2DBT4VbVxnvGtwNalrryqtgBbZmZmzlzqMiRJg/Ei5pLUCC9iLkmN8ORpktQIWzqS1AhbOpLUCFs6ktQIA1+SGmEPX5IaYQ9fkhphS0eSGmHgS1IjDHxJaoQ7bSWpEe60laRG2NKRpEYY+JLUCANfkhph4EtSIwa6xOH+4jVtJbViGq4h7VE6ktQIWzqS1AgDX5IaYeBLUiMMfElqhIEvSY0w8CWpER6Hr58z3/HCCxnVscRLWfdcxnlss7RceBy+JDXClo4kNcLAl6RGGPiS1AgDX5IaYeBLUiMMfElqhIEvSY0w8CWpEQa+JDVi5KdWSHIo8HfAncCVVfX+Ua9DkrR4A23hJ7kgye4kN8wa35Dk5iTbk2zqhp8PfKCqzgSeO+J6JUlLNGhL50JgQ/9AkhXAecDJwHpgY5L1wGpgR/eyu0dTpiRpWAO1dKrqqiRrZw2fCGyvqlsAklwKnArspBf617HAL5QkZwFnAaxZs2axdWvKzHeWy0mdtdKzbk6Haft/sVjLvf7Zhtlpu4qfbslDL+hXAR8Cfj3Ju4At8725qjZX1UxVzaxcuXKIMiRJgxhmp23mGKuq+h/gZQMtwPPhS9LYDLOFvxM4uu/xamDXYhbg+fAlaXyGCfxrgOOSrEtyEHAacPliFpDklCSb9+7dO0QZkqRBDHpY5iXA1cDxSXYmOaOq7gLOBq4AbgIuq6obF7Nyt/AlaXwGPUpn4zzjW4GtI61IkrRfTPTUCrZ0JGl8vIi5JDXCk6dJUiNSVZOugSR7gK8v8e1HAbeNsJxRm+b6prk2mO76rG3pprm+aa4Nfr6+h1XVwH+5OhWBP4wk26pqZtJ1zGea65vm2mC667O2pZvm+qa5Nhi+Pls6ktQIA1+SGnFfCPzNky5gH6a5vmmuDaa7Pmtbummub5prgyHrW/Y9fEnSYO4LW/iSpAEY+JLUiGUd+PNcU3eS9dya5EtJrkuyrRs7Isknknyl+/oLY6zn565FvFA9Sc7p5vLmJM+cQG1vTPLNbv6uS/KsCdV2dJJ/TXJTkhuTvKobn5a5m6++ic9fkkOSfC7J9V1tb+rGJz53C9Q28XmbVeeKJF9I8tHu8ejmrqqW5Q1YAXwVOAY4CLgeWD/hmm4Fjpo19hZgU3d/E/DmMdbzVOCxwA37qofedYmvBw4G1nVzu2LMtb0ReM0crx13bQ8BHtvdfyDwX10N0zJ389U38fmjd2Gkw7r7BwL/ATxhGuZugdomPm+z1vsHwD8AH+0ej2zulvMW/k+uqVtVdwL3XlN32pwKvK+7/z7geeNacVVdBdw+YD2nApdW1Y+q6mvAdnpzPM7a5jPu2r5VVZ/v7v83vdN/r2J65m6++uYztvqq5wfdwwO7WzEFc7dAbfMZ678rQJLVwLOB82fVMZK5W86BP981dSepgH9Jcm16F2kHeHBVfQt6P6jAL06suoXrmZb5PDvJF7uWz70fXSdWW5K1wGPobQ1O3dzNqg+mYP66lsR1wG7gE1U1NXM3T20wBfPWeTvwWuCevrGRzd1yDvw5r6k79ip+1pOr6rHAycDvJnnqhOtZjGmYz3cBDwdOAL4FvLUbn0htSQ4DPgi8uqq+v9BL5xibRH1TMX9VdXdVnUDvsqcnJnnkAi+fhtqmYt6SPAfYXVXXDvqWOcYWrG85B/7Q19Qdtara1X3dDXyY3ser7yR5CED3dffkKoQF6pn4fFbVd7ofyHuA9/DTj6djry3JgfTC9P1V9aFueGrmbq76pmn+unruAK4ENjBFcze7timatycDz01yK70W9a8muZgRzt1yDvyhr6k7SkkOTfLAe+8DzwBu6Go6vXvZ6cBHJlPhT8xXz+XAaUkOTrIOOA743DgLu/c/defX6M3f2GtLEuC9wE1V9ba+p6Zi7uarbxrmL8nKJA/q7t8feDrwn0zB3M1X2zTMG0BVnVNVq6tqLb08+1RVvZhRzt3+3uO8P2/As+gdofBV4PUTruUYenvMrwduvLce4Ejgk8BXuq9HjLGmS+h9RP0xva2BMxaqB3h9N5c3AydPoLaLgC8BX+z+Mz9kQrU9hd5H4y8C13W3Z03R3M1X38TnD3gU8IWuhhuAP9nXz8EU1DbxeZuj1pP46VE6I5s7T60gSY1Yzi0dSdIiGPiS1AgDX5IaYeBLUiMMfElqhIEvSY0w8CWpEf8P77Zjl7klLogAAAAASUVORK5CYII=",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAEICAYAAACktLTqAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAQ0klEQVR4nO3df6zddX3H8edrRcwCDOWHxhWwxcuYjdkUb3DDxPGHula8oG5OOuPEECqbNTPZ4uqyRPaHSV0ys5mgps6uZmMQ4nBS24lmGyEzRCkEZysyO6zjWqRFpFNcJIX3/jinena4pz33nnt67v3c5yO5ued8zvl+z/t+8u3rfPs5n/P5pqqQJLXl5yZdgCRp8RnuktQgw12SGmS4S1KDDHdJapDhLkkNMtwlqUGGu5aFJF9NclGSC5Pcd4LnVpKpvrYbkvz9gOev6W7zo+7Po0k+n+R1i/k3SCeT4a4lL8lzgBcD+4FXAscN9xE8r6pOB34V+BLw2STXjOm15pQO/11qZB5EWg5eBnyjOl+nnmZ84Q5AVX2vqv4auAH48LGwTfLSJHcmeSLJviRXHtsmyY4kNybZleSHSb6S5CU9j1+W5J4kR7q/L+t57M4kH0ryZeDHwIVJrknyUHdf307y9nH+zWqP4a4lK8m7kjwBfBn49e7tP6ITuE8kWTvmEm4DXgBc3P3fw07gi9229wI3Jbm45/kbgT8Hnk/nfxkf6v4dZwG7gI8CZwMfAXYlObtn23cAm4AzgMPd526oqjOAy4D7x/MnqlWGu5asqvrbqnoecC/wa8CvAHuBX6iq51XVt4+z+X3dN4Anum8KWxZQwsHu77O6r386sLWqnqqqfwU+TyfQj7mtqr5aVUeBm4CXd9uvAL5VVX9XVUer6mbgm8BMz7Y7qmpfd9ujwDPAy5L8fFU9UlX7FlC/VjDDXUtSkrO6wXyEzpnrncCDwMXAD5K87wS7uKT7BvC87hvE1gWUsbr7+3HgF4GHq+qZnse/0/McgO/13P4xnTcDutt+p2/f/ds+fOxGVT0JvA24HnikO9TzywuoXyuY4a4lqaoe74byu4G/6d7+AjDTDey/OgllvBk4ROdN5SBwft+HnRcA3x1iPwfpfCDcq3/b/7c8a1XdUVWvA15E5yz/k/MrXSud4a6lrnd2zCvoDNGMVZIXJtkMfBD4QPds/SvAk8D7kzwnyeV0hlVuGWKXu4FfSvK7SU5J8jZgHZ1hnUGvf2WS04CfAD8Cnh7179LKYrhrqXslnfHzs4Gnq+oHY3ytJ5I8CXwdeAPw1qraDlBVTwFXAhuAx4CPAb9XVd880U6r6vvAG+l8GPx94P3AG6vqsQGb/Fz3uQfpDAn9BvAHI/xdWoHixTokqT2euUtSgwx3SWqQ4S5JDTLcJalBp0y6AIBzzjmn1qxZM+kyJGlZuffeex+rqnPnemyi4Z5kBpiZmppiz549kyxFkpadJP3ffP6piQ7LVNXOqtp05plnTrIMSWqOY+6S1CDDXZIaZLhLUoMMd0lqkOEuSQ0y3CWpQRMN9yQzSbYdOXJkkmVIUnMm+iWmqtoJ7Jyenr5uoftYs2XXnO0Htl6x0F1K0rLnsIwkNchwl6QGGe6S1CDDXZIaZLhLUoMMd0lqkOEuSQ0y3CWpQYa7JDVoLOGe5LQk9yZ54zj2L0k6vqHCPcn2JIeS7O1rX5/kwST7k2zpeehPgFsXs1BJ0vCGPXPfAazvbUiyCrgR2ACsAzYmWZfktcA3gEcXsU5J0jwMtXBYVd2VZE1f86XA/qp6CCDJLcBVwOnAaXQC/3+T7K6qZ/r3mWQTsAngggsuWPAfIEl6tlFWhVwNPNxzfxZ4VVVtBkhyDfDYXMEOUFXbgG0A09PTNUIdkqQ+o4R75mj7aUhX1Y4T7iCZAWampqZGKEOS1G+U2TKzwPk9988DDs5nB1W1s6o2nXnmmSOUIUnqN0q43wNclGRtklOBq4Hb57MDr8QkSeMx7FTIm4G7gYuTzCa5tqqOApuBO4AHgFurat98Xtwzd0kaj2Fny2wc0L4b2L2oFUmSRuYFsiWpQRMNd4dlJGk8XDhMkhrksIwkNchhGUlqkMMyktSgUZYfWNLWbNk1Z/uBrVec5Eok6eRzzF2SGuSYuyQ1yDF3SWqQ4S5JDXLMXZIa5Ji7JDXIYRlJapDhLkkNMtwlqUGGuyQ1yNkyktQgZ8tIUoMclpGkBhnuktQgw12SGmS4S1KDDHdJapDhLkkNcp67JDXIee6S1CCHZSSpQYa7JDXolEkXcLKt2bJr4GMHtl5xEiuRpPHxzF2SGmS4S1KDDHdJapDhLkkNMtwlqUGGuyQ1aNHDPclLk3wiyWeS/P5i71+SdGJDhXuS7UkOJdnb174+yYNJ9ifZAlBVD1TV9cDvANOLX7Ik6USGPXPfAazvbUiyCrgR2ACsAzYmWdd97Erg34F/WbRKJUlDGyrcq+ou4PG+5kuB/VX1UFU9BdwCXNV9/u1VdRnw9kH7TLIpyZ4kew4fPryw6iVJcxpl+YHVwMM992eBVyW5HHgL8Fxg96CNq2obsA1genq6RqhDktRnlHDPHG1VVXcCdw61g2QGmJmamhqhDElSv1Fmy8wC5/fcPw84OJ8duJ67JI3HKOF+D3BRkrVJTgWuBm5fnLIkSaMYdirkzcDdwMVJZpNcW1VHgc3AHcADwK1VtW8+L+5l9iRpPIYac6+qjQPad3OcD02H2O9OYOf09PR1C92HJOnZvEC2JDXIC2RLUoNcOEySGuSwjCQ1yGEZSWqQwzKS1CDDXZIaNMraMiNbamvLrNmya872A1uvOMmVSNJoHHOXpAY5LCNJDTLcJalBhrskNcgvMUlSg/xAVZIa5LCMJDXIcJekBhnuktQgP1CVpAZNdPmB5XKZPZclkLTcOCwjSQ0y3CWpQYa7JDXIcJekBhnuktQgw12SGuQ8d0lqkAuHSVKDHJaRpAYZ7pLUIMNdkho00bVlljvXnJG0VHnmLkkNMtwlqUGGuyQ1yHCXpAYZ7pLUoLGEe5I3Jflkks8lef04XkOSNNjQ4Z5ke5JDSfb2ta9P8mCS/Um2AFTVP1XVdcA1wNsWtWJJ0gnN58x9B7C+tyHJKuBGYAOwDtiYZF3PU/6s+7gk6SQaOtyr6i7g8b7mS4H9VfVQVT0F3AJclY4PA/9cVffNtb8km5LsSbLn8OHDC61fkjSHUcfcVwMP99yf7ba9F3gt8NtJrp9rw6raVlXTVTV97rnnjliGJKnXqMsPZI62qqqPAh894cbJDDAzNTU1YhmSpF6jnrnPAuf33D8PODjsxq7nLknjMWq43wNclGRtklOBq4HbRy9LkjSKoYdlktwMXA6ck2QW+GBVfSrJZuAOYBWwvar2zWOfTQ7LuFqkpEkbOtyrauOA9t3A7oW8eFXtBHZOT09ft5DtJUlzc/kBSWrQRMM9yUySbUeOHJlkGZLUnImGu7NlJGk8HJaRpAY5LCNJDXJYRpIaNOryA5oH579LOlkclpGkBjksI0kNcraMJDXIcJekBhnuktQgP1CVpAb5gaokNchhGUlqkOEuSQ3yG6pLgN9clbTYDPdlyDcDSSfibBlJapCzZSSpQX6gKkkNMtwlqUGGuyQ1yNkyS9igWTGSdCKeuUtSgwx3SWqQ89wlqUHOc5ekBjksI0kNMtwlqUGGuyQ1yHCXpAYZ7pLUIMNdkhpkuEtSgwx3SWrQood7kguTfCrJZxZ735Kk4QwV7km2JzmUZG9f+/okDybZn2QLQFU9VFXXjqNYSdJwhj1z3wGs721Isgq4EdgArAM2Jlm3qNVJkhZkqHCvqruAx/uaLwX2d8/UnwJuAa5a5PokSQswypj7auDhnvuzwOokZyf5BPCKJB8YtHGSTUn2JNlz+PDhEcqQJPUb5UpMmaOtqur7wPUn2riqtgHbAKanp2uEOiRJfUY5c58Fzu+5fx5wcD47cD13SRqPUcL9HuCiJGuTnApcDdw+nx24nrskjcewUyFvBu4GLk4ym+TaqjoKbAbuAB4Abq2qffN5cc/cJWk8hhpzr6qNA9p3A7sX+uJVtRPYOT09fd1C9yFJejaXH5CkBo0yW2ZkSWaAmampqUmW0Yw1W3bN6/kHtl4xpkokTZoXyJakBjksI0kNmmi4O1tGksbDYRlJapDDMpLUIMNdkhrkVMgVbNDUSadISsufY+6S1CCHZSSpQYa7JDXIee6S1CDH3CWpQQ7LSFKDDHdJapDhLkkNMtwlqUF+Q1Xz4rdapYU52f92nC0jSQ1yWEaSGmS4S1KDDHdJapDhLkkNMtwlqUGGuyQ1yHnuWhSLNYd33HOBB+1/EOfva7lynrskNchhGUlqkOEuSQ0y3CWpQYa7JDXIcJekBhnuktQgw12SGmS4S1KDDHdJatCiLz+Q5DTgY8BTwJ1VddNiv4Yk6fiGOnNPsj3JoSR7+9rXJ3kwyf4kW7rNbwE+U1XXAVcucr2SpCEMOyyzA1jf25BkFXAjsAFYB2xMsg44D3i4+7SnF6dMSdJ8DDUsU1V3JVnT13wpsL+qHgJIcgtwFTBLJ+Dv5zhvHkk2AZsALrjggvnWrTGa78qJJ2tfk9i/xmPcq38uluVS51xG+UB1NT87Q4dOqK8GbgN+K8nHgZ2DNq6qbVU1XVXT55577ghlSJL6jfKBauZoq6p6EnjXUDtwPXdJGotRztxngfN77p8HHJzPDlzPXZLGY5Rwvwe4KMnaJKcCVwO3z2cHSWaSbDty5MgIZUiS+g07FfJm4G7g4iSzSa6tqqPAZuAO4AHg1qraN58X98xdksZj2NkyGwe07wZ2L2pFkqSRTXT5AYdlJGk8vEC2JDXIhcMkqUGpqknXQJLDwHcWuPk5wGOLWE6r7Kfh2VfDsZ+GM85+enFVzfkt0CUR7qNIsqeqpiddx1JnPw3PvhqO/TScSfWTwzKS1CDDXZIa1EK4b5t0AcuE/TQ8+2o49tNwJtJPy37MXZL0bC2cuUuS+hjuktSgZR3uA67hKiDJgSRfT3J/kj3dtrOSfCnJt7q/nz/pOk+2ua4HfLx+SfKB7vH1YJLfnEzVkzGgr25I8t3ucXV/kjf0PLbi+irJ+Un+LckDSfYl+cNu++SPqapalj/AKuC/gAuBU4GvAesmXddS+QEOAOf0tf0FsKV7ewvw4UnXOYF+eQ1wCbD3RP1C59rAXwOeC6ztHm+rJv03TLivbgD+eI7nrsi+Al4EXNK9fQbwn92+mPgxtZzP3H96Ddeqego4dg1XDXYV8Onu7U8Db5pcKZNRVXcBj/c1D+qXq4BbquonVfVtYD+d425FGNBXg6zIvqqqR6rqvu7tH9JZ/nw1S+CYWs7hPugaruoo4ItJ7u1ejBzghVX1CHQOSuAFE6tuaRnULx5jc9uc5D+6wzbHhhtWfF8lWQO8AvgKS+CYWs7hPuc1XE96FUvXq6vqEmAD8J4kr5l0QcuQx9izfRx4CfBy4BHgL7vtK7qvkpwO/CPwvqr6n+M9dY62sfTTcg73ka/h2rKqOtj9fQj4LJ3/+j2a5EUA3d+HJlfhkjKoXzzG+lTVo1X1dFU9A3ySnw0prNi+SvIcOsF+U1Xd1m2e+DG1nMN95Gu4tirJaUnOOHYbeD2wl07/vLP7tHcCn5tMhUvOoH65Hbg6yXOTrAUuAr46gfqWjGOB1fVmOscVrNC+ShLgU8ADVfWRnocmfkwNdZm9paiqjiY5dg3XVcD2muc1XBv2QuCzneOOU4B/qKovJLkHuDXJtcB/A2+dYI0T0b0e8OXAOUlmgQ8CW5mjX6pqX5JbgW8AR4H3VNXTEyl8Agb01eVJXk5nKOEA8G5Y0X31auAdwNeT3N9t+1OWwDHl8gOS1KDlPCwjSRrAcJekBhnuktQgw12SGmS4S1KDDHdJapDhLkkN+j+u63ef1a65EQAAAABJRU5ErkJggg==",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAEICAYAAACktLTqAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAQYElEQVR4nO3dfaykZXnH8e/Pxa0VIlbY+rIvnrVL0I2toiOIfbPVNEthxdgX2WiLFtliSquNtV2bxpekTdCYxpqgZpUVUyy42RLdLau0taXUhigLNhFYSTe4yCkoq+hq1BTBq3/MUKeHOcucM2d2ztzn+/mHmXvmeeY69+757cP13PM8qSokSW153KQLkCQtPcNdkhpkuEtSgwx3SWqQ4S5JDTLcJalBhrt0DEk+neTCSdchLVRc567WJPkC8BrgYWBPVb3gGO8t4PtAAUeBTwBvraqHj0et0rh45K6mJHk88EzgEPBC4NYhNnteVZ0E/DLwauD3xlehdHwY7mrNc4E7qvu/pB2GC3cAquoQ8B/A8x8ZS3JDkjf0Hr8uyeeSvDfJt5J8Jck5fe/dmOTGJN9N8s9JLk9yVe+1JyS5Ksk3k3w7yc1Jnro0P7L0aIa7mpDk9Um+TTecz+49fgvw7l6YbhxiH88GfpHuUf98zgLuBE4F3gNckSS91/4O+AJwCvBO4Hf6trsQOBlY33v9EuAHQ/540oIZ7mpCVX20qp4M3AK8GPg54DbgSVX15Kr6yjE2vzXJ94CDwA3AB47x3rur6sO9nvzHgKcDT02yAXgR8PaqerCqPgfs7dvuh3RDfVNVPVxVt1TVdxb1w0pDMNw19ZI8pXd0fhR4Cd2AvhM4HfhWkjc/xi5eAJxEt99+FnDiMd77tUceVNX3ew9PAp4BPNA3BnBP3+O/Ba4Hrklyb5L39M4PSGNhuGvqVdUDvaP23wc+0nv8GWBr76j9fUPso6pqN3AT8PZFlHEf8JQkT+wbW9+3/x9W1buqajPdf4DOA353EZ8jDcVwV0v6V8ecQbdFs1CXAduTPG0hG1XV3cAB4J1JVic5G9j6yOtJfiXJzyZZBXyHbpvG5ZYaG8NdLXkh3f75KcDDVfWthe6gqr4E/Bvw1kV8/muAs4FvAn9Jd838//Reexqwh26wH+x9xlWL+AxpKH6JSRqTJJ8AvlxV75h0LVp5PHKXlkiSFyX5mSSPS7IFOB/45ITL0gp1wqQLkBryNOBaukseZ4E3VtUXJ1uSVirbMpLUINsyktSgZdGWOfXUU2tmZmbSZUjSVLnlllu+UVVrBr22LMJ9ZmaGAwcOTLoMSZoqSe6e77WJtmWSbE2y8+jRo5MsQ5KaM9Fwr6p9VbX95JNPnmQZktQcT6hKUoMMd0lqkOEuSQ0y3CWpQYa7JDXIcJekBi2LLzGNYmbHdQPHD1927nGuRJKWD4/cJalBhrskNWgs4Z7kxCS3JDlvHPuXJB3bUOGeZFeS+5PcNmd8S5I7kxxKsqPvpT8Ddi9loZKk4Q175H4lsKV/oHcX98uBc4DNwLYkm5O8HLgD+PoS1ilJWoChVstU1Y1JZuYMnwkcqqq7AJJcQ/eekScBJ9IN/B8k2V9VP5q7zyTbge0AGzZsWPQPIEl6tFGWQq4F7ul7PgucVVWXAiR5HfCNQcEOUFU7gZ0AnU7He/1J0hIaJdwzYOz/QrqqrnzMHSRbga2bNm0aoQxJ0lyjrJaZBdb3PV8H3LuQHXg9d0kaj1HC/WbgtCQbk6wGLgD2Lk1ZkqRRDNWWSXI18FLg1CSzwDuq6ooklwLXA6uAXVV1+0I+fJxtGS9LIGklG3a1zLZ5xvcD+xf74VW1D9jX6XQuXuw+JEmP5g2yJalB3iBbkhrkhcMkqUGGuyQ1yJ67JDXInrskNci2jCQ1yLaMJDXItowkNci2jCQ1yHCXpAaNcj33qTTfBcXAi4pJaocnVCWpQZ5QlaQG2XOXpAYZ7pLUIMNdkhpkuEtSg1wtI0kNcrWMJDXItowkNchwl6QGGe6S1CDDXZIaZLhLUoMMd0lqkOEuSQ3yS0yS1CC/xCRJDbItI0kNMtwlqUGGuyQ1yHCXpAYZ7pLUIMNdkhpkuEtSgwx3SWqQ4S5JDVrycE/ynCQfSrInyRuXev+SpMc2VLgn2ZXk/iS3zRnfkuTOJIeS7ACoqoNVdQnw20Bn6UuWJD2WYY/crwS29A8kWQVcDpwDbAa2Jdnce+0VwOeAzy5ZpZKkoQ0V7lV1I/DAnOEzgUNVdVdVPQhcA5zfe//eqnoJ8Jr59plke5IDSQ4cOXJkcdVLkgY6YYRt1wL39D2fBc5K8lLgVcBPAPvn27iqdgI7ATqdTo1QhyRpjlHCPQPGqqpuAG4YagfJVmDrpk2bRihDkjTXKKtlZoH1fc/XAfcuZAdez12SxmOUcL8ZOC3JxiSrgQuAvUtTliRpFMMuhbwauAk4Pclskouq6iHgUuB64CCwu6puX8iHe5s9SRqPoXruVbVtnvH9HOOk6RD73Qfs63Q6Fy92H5KkRxvlhOrIltsJ1Zkd1w0cP3zZuce5EkkajTfIlqQGTfTIfVp4RC9p2nhVSElq0ETD3dUykjQe9twlqUG2ZSSpQbZlJKlBtmUkqUG2ZSSpQYa7JDXIcJekBnlCVZIa5AlVSWqQbRlJapDhLkkNMtwlqUGGuyQ1yNUyktQgV8tIUoNsy0hSgwx3SWqQ4S5JDTLcJalBhrskNeiESRcwzWZ2XDdw/PBl5x7nSiTp//PIXZIa5JeYJKlBfolJkhpkW0aSGmS4S1KDDHdJapDhLkkNMtwlqUGGuyQ1yHCXpAYZ7pLUIMNdkho0lnBP8sokH07yqSS/No7PkCTNb+hwT7Iryf1JbpszviXJnUkOJdkBUFWfrKqLgdcBr17SiiVJj2khR+5XAlv6B5KsAi4HzgE2A9uSbO57y1/0XpckHUdDh3tV3Qg8MGf4TOBQVd1VVQ8C1wDnp+vdwKer6tZB+0uyPcmBJAeOHDmy2PolSQOM2nNfC9zT93y2N/aHwMuB30xyyaANq2pnVXWqqrNmzZoRy5Ak9Rv1TkwZMFZV9X7g/Y+5cbIV2Lpp06YRy5Ak9Rv1yH0WWN/3fB1w77Abez13SRqPUcP9ZuC0JBuTrAYuAPaOXpYkaRQLWQp5NXATcHqS2SQXVdVDwKXA9cBBYHdV3b6AfXqbPUkag6F77lW1bZ7x/cD+xXx4Ve0D9nU6nYsXs70kabBRT6iOpNUTqjM7rhs4fviyc49zJZJWKm+QLUkN8sJhktQgw12SGjTRcHe1jCSNhz13SWqQbRlJapBLIZcBl05KWmq2ZSSpQbZlJKlBhrskNchwl6QGuc5dkhrkCVVJapBtGUlqkOEuSQ0y3CWpQRP9hupKM983USVpqblaRpIa5GoZSWqQPXdJapDhLkkNMtwlqUGGuyQ1yHCXpAYZ7pLUINe5S1KDXOcuSQ2yLSNJDTLcJalBhrskNchwl6QGGe6S1CDDXZIaZLhLUoMMd0lqkOEuSQ0y3CWpQUse7kmeleSKJHuWet+SpOEMFe5JdiW5P8ltc8a3JLkzyaEkOwCq6q6qumgcxUqShjPskfuVwJb+gSSrgMuBc4DNwLYkm5e0OknSopwwzJuq6sYkM3OGzwQOVdVdAEmuAc4H7hhmn0m2A9sBNmzYMGy9AmZ2XDdw/PBl5x7nSiQtV6P03NcC9/Q9nwXWJjklyYeAM5K8bb6Nq2pnVXWqqrNmzZoRypAkzTXUkfs8MmCsquqbwCVD7SDZCmzdtGnTCGVIkuYa5ch9Fljf93wdcO9CduDNOiRpPEYJ95uB05JsTLIauADYuzRlSZJGMexSyKuBm4DTk8wmuaiqHgIuBa4HDgK7q+r2hXy491CVpPEYdrXMtnnG9wP7F/vhVbUP2NfpdC5e7D4kSY820csPeOQuSeMx0XD3hKokjYcXDpOkBhnuktQge+6S1CB77pLUINsyktQg2zKS1CDbMpLUINsyktQgw12SGmS4S1KDRrlZx8i8WcfxsdDb8s33/mNtI2l58YSqJDXItowkNchwl6QGGe6S1CDDXZIa5GoZjdWxVt4M4moctWqhq9ZG5WoZSWqQbRlJapDhLkkNMtwlqUGGuyQ1yHCXpAYZ7pLUINe5L2MLXSO+0PdLapfr3CWpQbZlJKlBhrskNchwl6QGGe6S1CDDXZIaZLhLUoMMd0lqkOEuSQ0y3CWpQYa7JDVoya8tk+RE4APAg8ANVfXxpf4MSdKxDXXknmRXkvuT3DZnfEuSO5McSrKjN/wqYE9VXQy8YonrlSQNYdi2zJXAlv6BJKuAy4FzgM3AtiSbgXXAPb23Pbw0ZUqSFmKotkxV3ZhkZs7wmcChqroLIMk1wPnALN2A/0+O8Y9Hku3AdoANGzYstG4tgeV4ieCF1nT4snPHVImGdaw/s2n585nvZ5iW+gcZ5YTqWn58hA7dUF8LXAv8RpIPAvvm27iqdlZVp6o6a9asGaEMSdJco5xQzYCxqqrvAa8fagferEOSxmKUI/dZYH3f83XAvQvZgTfrkKTxGCXcbwZOS7IxyWrgAmDv0pQlSRrFsEshrwZuAk5PMpvkoqp6CLgUuB44COyuqtsX8uFJtibZefTo0YXWLUk6hmFXy2ybZ3w/sH+xH15V+4B9nU7n4sXuQ5L0aBO9/IBH7pI0HhMNd0+oStJ4eOEwSWpQqmrSNZDkCHD3Ijc/FfjGEpbTCudlMOdlMOdlsOU+L8+sqoHfAl0W4T6KJAeqqjPpOpYb52Uw52Uw52WwaZ4X2zKS1CDDXZIa1EK475x0AcuU8zKY8zKY8zLY1M7L1PfcJUmP1sKRuyRpDsNdkho01eE+zz1cV5wk65P8a5KDSW5P8qbe+FOS/FOS/+r996cmXevxlmRVki8m+Yfe8xU/JwBJnpxkT5Iv9/7enL3S5ybJH/d+f25LcnWSJ0zznExtuB/jHq4r0UPAW6rqOcCLgT/ozcUO4LNVdRrw2d7zleZNdK9a+gjnpOtvgM9U1bOB59GdoxU7N0nWAn8EdKrqucAqupcxn9o5mdpwp+8erlX1IPDIPVxXnKq6r6pu7T3+Lt1f1LV05+Njvbd9DHjlRAqckCTrgHOBj/QNr+g5AUjyJOCXgCsAqurBqvo2zs0JwE8mOQF4It2bD03tnExzuM93D9cVrXcj8zOAzwNPrar7oPsPAPDTEyxtEt4H/Cnwo76xlT4nAM8CjgAf7bWsPpLkRFbw3FTVfwPvBb4K3Accrap/ZIrnZJrDfeA9XI97FctIkpOAvwfeXFXfmXQ9k5TkPOD+qrpl0rUsQycALwA+WFVnAN9jitoN49DrpZ8PbASeAZyY5LWTrWo00xzuI9/DtSVJHk832D9eVdf2hr+e5Om9158O3D+p+ibg54FXJDlMt2X3q0muYmXPySNmgdmq+nzv+R66Yb+S5+blwFeq6khV/RC4FngJUzwn0xzu3sO1J0no9k8PVtVf9720F7iw9/hC4FPHu7ZJqaq3VdW6qpqh+3fjX6rqtazgOXlEVX0NuCfJ6b2hlwF3sLLn5qvAi5M8sff79DK6566mdk6m+huqSX6dbl91FbCrqv5qshVNRpJfAP4d+BI/7i//Od2++25gA92/vL9VVQ9MpMgJSvJS4E+q6rwkp+CckOT5dE80rwbuAl5P92Bvxc5NkncBr6a7+uyLwBuAk5jSOZnqcJckDTbNbRlJ0jwMd0lqkOEuSQ0y3CWpQYa7JDXIcJekBhnuktSg/wVmpmErXqOqjgAAAABJRU5ErkJggg==",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXoAAAEICAYAAABRSj9aAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAQvUlEQVR4nO3df5Dcd13H8efLlFRoK78asKaNKZPaMYNI4UypKCK/TIFQBMRE0KK1GZQ6MPwh6eA4+JfgDOpUIhilFsbSEis/EghTEO1Uxmqb0qJJQySEQI9Qkootip2Bwts/7ltYj7vrXnb39vaT52Pm5nY/u/vdzztz98r33t/Pfr+pKiRJ7fqBcU9AkjRaBr0kNc6gl6TGGfSS1DiDXpIaZ9BLUuMMeklqnEGvJiX5n56v7yR5oOf+q5K8Jcm3uvv3JfnnJBd1r12Z5O1JprvHv5DkT+Z4j7ckqSQbZo2/IskdSb6e5M7Zj0tLzaBXk6rq9Ie+gC8Bm3rGru2e9v7u8VXAp4APJAlwJTAFbADOAH4euKN3+93zfhX4GnDprLc/C7gceCzwPuCaEZQo9c2g10mvqr4FvAf4YeDxwE8BH6yqozXjSFW9d9bLfhb4EeD1wOYkK3u292dVtbeqvg38E/DEJSlEmodBr5NeklOB1wDTVXUv8C/AG5P8dpKf6PbeZ7sU2A28v7v/4jm2exrwx7hHrzEz6HUye2WS+4C7gacDL+3G/xB4G/AqYC/w5STfbc8keRTwS8D7ur8GbuD72zcAO4GjwJtGNH+pL6eMewLSGO2sqlfPHuxaLtuB7UkeCfwGcHWSW6vqAPCLwIPAnu4l1wJ/n2RVVR0HSPJjwM8BT6iqB5egFmle7tFLC6iqB6pqO/BfwPpu+FLgdOBLSe4B/hZ4BLCl56VnAfdV1f8u5XyluRj00ixJ3pDk2UkemeSUrm1zBnBHktXAc5npyT+1+/pJZlo9ve2bW4ALlnTi0jwMeun7PQC8HbgHuBd4HfDyqjrMzJLKO6vq41V1z0NfwFXAU5I8udvGhcws2ZTGLl54RJLa5h69JDXOoJekxhn0ktQ4g16SGrcsPjB15pln1tq1a8c9DUmaKLfffvu9VbXq4Z63LIJ+7dq17N27d9zTkKSJkuSL/TzP1o0kNc6gl6TGjTXok2xKsuP+++8f5zQkqWljDfqq2l1VWx/96EePcxqS1DRbN5LUOINekhpn0EtS4wx6SWrcsvjA1CDWbvvonONH3vqiJZ6JJC1P7tFLUuMMeklqnEEvSY0z6CWpcSMJ+iSnJbk9yYtHsX1JUv/6CvokVyc5lmTfrPGNSQ4mOZRkW89DbwJ2DnOikqQT0+8e/TXAxt6BJCuA7cDFwHpgS5L1SZ4H3AV8dYjzlCSdoL7W0VfVzUnWzhreAByqqsMASa4HLgFOB05jJvwfSLKnqr4ze5tJtgJbAdasWXPCBUiSFjbIB6ZWA3f33J8GLqyqKwCSvAa4d66QB6iqHcAOgKmpqRpgHpKkBQwS9Jlj7LuBXVXXPOwGkk3ApnXr1g0wDUnSQgZZdTMNnNNz/2zg6GI24PnoJWn0Bgn624DzkpybZCWwGdi1mA14hSlJGr1+l1deB9wCnJ9kOsllVfUgcAVwI3AA2FlV+xfz5u7RS9Lo9bvqZss843uAPUOdkSRpqLw4uCQ1bqzno6+q3cDuqampy4e9bc9TL0kzPKmZJDXO1o0kNW6sQe+qG0kaPVs3ktQ4g16SGmePXpIaZ49ekhpn60aSGmfQS1LjDHpJapwHYyWpcR6MlaTG2bqRpMYZ9JLUOINekhpn0EtS41x1I0mNc9WNJDXO1o0kNc6gl6TGGfSS1LhTxj2B5WTtto/OOX7krS9a4plI0vCcdEE/X5hLUqts3UhS41xHL0mNcx29JDXO1o0kNc6gl6TGGfSS1DiDXpIaZ9BLUuMMeklqnEEvSY0z6CWpcQa9JDVu6EGf5MeTvCvJDUl+a9jblyQtTl9Bn+TqJMeS7Js1vjHJwSSHkmwDqKoDVfVa4JXA1PCnLElajH736K8BNvYOJFkBbAcuBtYDW5Ks7x57CfAp4JNDm6kk6YT0FfRVdTPwtVnDG4BDVXW4qr4JXA9c0j1/V1X9NPCq+baZZGuSvUn2Hj9+/MRmL0l6WINceGQ1cHfP/WngwiTPBl4GnArsme/FVbUD2AEwNTVVA8xDkrSAQYI+c4xVVd0E3DTAdiVJQzTIqptp4Jye+2cDRxezAS88IkmjN0jQ3wacl+TcJCuBzcCuxWzAC49I0uj1u7zyOuAW4Pwk00kuq6oHgSuAG4EDwM6q2r+YN3ePXpJGr68efVVtmWd8DwsccO1ju7uB3VNTU5ef6DYkSQvzFAiS1LixBr2tG0kavbEGvQdjJWn0bN1IUuNs3UhS42zdSFLjbN1IUuMMeklqnD16SWrcIGevHNikfDJ27baPzjl+5K0vWuKZSNLi2bqRpMYZ9JLUOHv0ktQ419FLUuNs3UhS4wx6SWqcQS9JjTPoJalxrrqRpMa56kaSGmfrRpIaN9Zz3Uw6z4EjaRK4Ry9JjTPoJalxBr0kNc6gl6TGGfSS1Dg/MCVJjfMDU5LUOFs3ktQ4g16SGmfQS1LjDHpJapxBL0mNM+glqXEGvSQ1zqCXpMYZ9JLUuJEEfZKXJvnLJB9O8oJRvIckqT99B32Sq5McS7Jv1vjGJAeTHEqyDaCqPlRVlwOvAX55qDOWJC3KYi4leA3wDuC9Dw0kWQFsB54PTAO3JdlVVXd1T/m97vGTipcYlLSc9L1HX1U3A1+bNbwBOFRVh6vqm8D1wCWZ8TbgY1X16bm2l2Rrkr1J9h4/fvxE5y9JehiD9uhXA3f33J/uxn4HeB7wiiSvneuFVbWjqqaqamrVqlUDTkOSNJ/FtG7mkjnGqqquAq562Bcnm4BN69atG3AakqT5DLpHPw2c03P/bOBovy/2fPSSNHqD7tHfBpyX5Fzgy8Bm4Ff6ffHJtkfvQVpJ47CY5ZXXAbcA5yeZTnJZVT0IXAHcCBwAdlbV/n636R69JI1e33v0VbVlnvE9wJ6hzUiSNFReHFySGufFwSWpcZ7UTJIaZ+tGkhpn60aSGmfrRpIaZ9BLUuPs0UtS4+zRS1LjbN1IUuMMeklq3KBnr9QQeFZLSaPkwVhJapwHYyWpcfboJalxBr0kNc6gl6TGGfSS1DhX3UhS41x1I0mNs3UjSY3zk7HLmJ+YlTQM7tFLUuPco59A7ulLWgyDviH+ByBpLrZuJKlxrqOXpMa5jl6SGmfrRpIaZ9BLUuMMeklqnEEvSY0z6CWpcQa9JDXOoJekxnkKhJOYp0yQTg7u0UtS44Ye9EmelOTdSW4Y9rYlSYvXV9AnuTrJsST7Zo1vTHIwyaEk2wCq6nBVXTaKyUqSFq/fPfprgI29A0lWANuBi4H1wJYk64c6O0nSwPo6GFtVNydZO2t4A3Coqg4DJLkeuAS4q59tJtkKbAVYs2ZNv/PVCZjvoKukk8MgPfrVwN0996eB1Uken+RdwAVJrpzvxVW1o6qmqmpq1apVA0xDkrSQQZZXZo6xqqr/BF47wHYlSUM0yB79NHBOz/2zgaOL2YAXHpGk0Rsk6G8DzktybpKVwGZg12I24IVHJGn0+l1eeR1wC3B+kukkl1XVg8AVwI3AAWBnVe1fzJu7Ry9Jo9fvqpst84zvAfac6JtX1W5g99TU1OUnug1J0sI8BYIkNW6sQW/rRpJGb6xB78FYSRo9WzeS1Lixno8+ySZg07p168Y5Dc2y0CkTPFe9NHls3UhS42zdSFLjDHpJapzLKyWpcfboJalxtm4kqXEGvSQ1znX0Gor51t6Pa939cpuPNE726CWpcbZuJKlxBr0kNc6gl6TGGfSS1DhX3WhRFjqz5WKev9jVL66i0aQb58+wq24kqXG2biSpcQa9JDXOoJekxhn0ktQ4g16SGmfQS1LjDHpJapwfmNJY+AEoaen4gSlJapytG0lqnEEvSY0z6CWpcQa9JDXOoJekxhn0ktQ4g16SGmfQS1LjDHpJatzQT4GQ5DTgz4FvAjdV1bXDfg9JUv/62qNPcnWSY0n2zRrfmORgkkNJtnXDLwNuqKrLgZcMeb6SpEXqt3VzDbCxdyDJCmA7cDGwHtiSZD1wNnB397RvD2eakqQT1VfrpqpuTrJ21vAG4FBVHQZIcj1wCTDNTNjfyQL/kSTZCmwFWLNmzWLnrUbNd1bLUT9/vrNmLrR9z7T58Cb9LKWTPv+HDHIwdjXf23OHmYBfDXwAeHmSdwK753txVe2oqqmqmlq1atUA05AkLWSQg7GZY6yq6hvAr/e1Ac9HL0kjN8ge/TRwTs/9s4Gji9mA56OXpNEbJOhvA85Lcm6SlcBmYNdwpiVJGpZ+l1deB9wCnJ9kOsllVfUgcAVwI3AA2FlV+xfz5kk2Jdlx//33L3bekqQ+9bvqZss843uAPSf65lW1G9g9NTV1+YluQ5K0sLGeAsE9ekkaPS8OLkmN86RmktS4VNW450CS48AXT/DlZwL3DnE642Y9y1dLtYD1LGf91vKjVfWwnzhdFkE/iCR7q2pq3PMYFutZvlqqBaxnORt2LbZuJKlxBr0kNa6FoN8x7gkMmfUsXy3VAtaznA21lonv0UuSFtbCHr0kaQEGvSQ1bqKDfp5r1i4rc11vN8njknwiyee674/teezKrp6DSX6hZ/zpSf69e+yqJHNdD2DkkpyT5B+THEiyP8nrJ7WmJD+Y5NYkn+lq+YNJraVXkhVJ7kjyke7+xNaT5Eg3jzuT7J3kepI8JskNST7b/f5ctGS1VNVEfgErgM8DTwJWAp8B1o97XnPM81nA04B9PWN/BGzrbm8D3tbdXt/VcSpwblffiu6xW4GLmLngy8eAi8dUz1nA07rbZwD/0c174mrq3vf07vYjgH8FnjGJtcyq643A+4CPNPDzdgQ4c9bYRNYDvAf4ze72SuAxS1XLWH4Qh/SPdhFwY8/9K4Erxz2veea6lv8f9AeBs7rbZwEH56qBmVNAX9Q957M941uAvxh3Xd1cPgw8f9JrAh4FfBq4cJJrYeYCQJ8EnsP3gn6S6znC9wf9xNUD/BDwBboFMEtdyyS3bua7Zu0keGJVfQWg+/6Ebny+mlZ3t2ePj1VmLhh/ATN7whNZU9fmuBM4Bnyiqia2ls6fAr8LfKdnbJLrKeDjSW5PsrUbm8R6ngQcB/66a6v9VZLTWKJaJjno57xm7ZLPYrjmq2nZ1ZrkdODvgDdU1dcXeuocY8umpqr6dlU9lZk94Q1JnrzA05d1LUleDByrqtv7fckcY8umns4zq+ppwMXA65I8a4HnLud6TmGmhfvOqroA+AYzrZr5DLWWSQ76ga9ZO0ZfTXIWQPf9WDc+X03T3e3Z42OR5BHMhPy1VfWBbniia6qq+4CbgI1Mbi3PBF6S5AhwPfCcJH/D5NZDVR3tvh8DPghsYDLrmQamu78YAW5gJviXpJZJDvpJvmbtLuDS7valzPS5HxrfnOTUJOcC5wG3dn/S/XeSZ3RH2H+t5zVLqnv/dwMHquqPex6auJqSrErymO72I4HnAZ9lAmsBqKorq+rsqlrLzO/DP1TVq5nQepKcluSMh24DLwD2MYH1VNU9wN1Jzu+GngvcxVLVMo4DLEM8wPFCZlZ9fB5487jnM88crwO+AnyLmf+NLwMez8wBs8913x/X8/w3d/UcpOdoOjDFzA/554F3MOugzhLW8zPM/Kn4b8Cd3dcLJ7Em4CnAHV0t+4Df78YnrpY5ans23zsYO5H1MNPX/kz3tf+h3/EJruepwN7u5+1DwGOXqhZPgSBJjZvk1o0kqQ8GvSQ1zqCXpMYZ9JLUOINekhpn0EtS4wx6SWrc/wGIwbpC+392cwAAAABJRU5ErkJggg==",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXYAAAEICAYAAABLdt/UAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAOfUlEQVR4nO3df4xl5V3H8fdHFmoVHNqyGrJAd3EXLH9ItCMlVBPUWheadWtDCltsU0rYoIGoaUxXY6ym/4DxD4PQ4NoS0qisiKTspotoqhRNwbJryq8SzLpW2FLDLrRjoImw7dc/5hKvw8zumb137p373PcrmWTuc8898z0nM5/7zHOe+5xUFZKkdnzfuAuQJA2XwS5JjTHYJakxBrskNcZgl6TGGOyS1BiDXZIaY7BLUmMMdk2NJB9N8kSS7yT5rySfTjLTe+73k7yW5OW+r2/3vbaSvNJrfzHJF5Nc2ff8miS3JXkuybeT3JXkzWM4TMlg13RI8nHgZuC3gBngYmA98HdJTu5t9ldVdWrf1+kLdnNhVZ0KnA/cCdya5JO9504GvgW8E3g7sAG4ceWOSFpaXFJArUvyQ8DzwMeq6u6+9lOBg8yH/QZgY1X9yhL7KGBTVR3oa7sC+HNgXVW9uGD7PwVerqqPD/t4pOOxx65pcAnw/cC9/Y1V9TJwP/DeE9zvfcAa4KL+xiTvBrYBf3mC+5UGYrBrGpwBHKmqo4s8901gbe/7D/bGx1//+sdj7bSqXgOOAG99vS3JJmA38/8d7B9O+dLyrBl3AdIIHAHOSLJmkXA/Ezjc+/7upYZiFtMbm18LvNTXfA1wX1XdM0jB0iDssWsaPAz8D/CB/sYkPwhcBnzpBPe7FTgKfKWv7Uzmx/OlsTHY1byqmgP+APiTJJuTnJxkPfDXzPfm/2I5+0vy1iRXA7cBNy+4cPobwE1DKVw6QQ7FaCpU1R8meRH4I2Aj8Cbme+rvqapXkgBcmeT9C156blW90Pv+sd7smFeBx4DfrKqFF0hvBr4BfGpljkQ6Pqc7aiol+Rjzvfh3V9Wz465HGiaDXVMryYeB16pq17hrkYbJYJekxnjxVJIaM9aLp0m2AFtOO+20684777xxliJJE2f//v1HqmrtwvZVMRQzOztb+/btG3cZkjRRkuyvqtmF7Q7FSFJjDHZJaozBLkmNMdglqTEGuyQ1xmCXpMYY7JLUGINdkhrT7LK963d8YdH2r9/0vhFXIkmjZY9dkhpjsEtSYwx2SWqMwS5JjRlrsCfZkmTn3NzcOMuQpKaMNdirak9VbZ+ZmRlnGZLUFIdiJKkxBrskNcZgl6TGGOyS1BiDXZIaY7BLUmMMdklqjMEuSY0x2CWpMQa7JDWm2RttLMUbcEhqnT12SWqMwS5JjTHYJakxBrskNcZgl6TGGOyS1BiDXZIaM/Hz2Jealy5J02roPfYk70hye5J7kvzqsPcvSTq2TsGe5I4kLyR5ckH75iTPJDmQZAdAVT1dVdcDHwRmh1+yJOlYuvbY7wQ29zckOQm4DbgMuADYluSC3nO/BPwz8MWhVSpJ6qRTsFfVQ8BLC5ovAg5U1cGqehXYBWztbb+7qi4Brh5msZKk4xvk4uk64Lm+x4eAdyW5FPgA8CZg71IvTrId2A5wzjnnDFCGJKnfIMGeRdqqqh4EHjzei6tqJ7ATYHZ2tgaoQ5LUZ5BZMYeAs/senwU8P1g5kqRBDRLsjwKbkmxIcgpwFbB7OTtIsiXJzrm5uQHKkCT16zrd8S7gYeD8JIeSXFtVR4EbgAeAp4G7q+qp5fzwqtpTVdtnZmaWW7ckaQmdxtiratsS7Xs5xgVSabXxDlqaBhO/pMCwHGtpAv/oJU0Sg11Ncg0hTbOxBnuSLcCWjRs3jrMMTQCDWupurMFeVXuAPbOzs9eNsw6tHuMKcMfe1RLXY5ekxhjsktQYx9ilY3CIRpPIMXatqFaD0emxWs2c7qjODLNuWn0z0+Qw2DvwD1XSJDHYV4BvBFqMvxcaFYNdY+EHjqSV46yYARhO/8dzIa0eY53H7rK9kjR8DsWMkGOskkbBYJ9iDp9IbTLYp4ABvrr5n5yGzWCXVqnlviH7RqDXjfXiqTezlqThc62YCeS/7pKOxaGYVWBYQe1YuiRwPXZJao49dqkRDtHpdQb7KubQiqQT4VCMJDXGHrvUOIdopo/z2CWpMa7uKEmNcYxdkhrjGLuk/8cx+clnj12SGmOPXZpSfk6iXfbYJakxBrskNcZgl6TGOMYuqRNny0wOe+yS1BiXFJCkxrikgCQ1xqEYSWqMwS5JjXFWjKQVcaxPtjqTZmXZY5ekxhjsktQYh2IkDWSYi4n5IajhsMcuSY0x2CWpMQa7JDXGMXZJI+dNPlaWPXZJaozBLkmNcXVHSWqMqztKUmMcipGkxhjsktQYg12SGmOwS1JjDHZJaozBLkmNMdglqTGuFSNp1XOd9uUx2CVNLAN/cQ7FSFJjDHZJaozBLkmNMdglqTEGuyQ1xmCXpMYY7JLUGINdkhpjsEtSYwx2SWqMwS5JjRl6sCd5f5I/S3JfkvcOe/+SpGPrFOxJ7kjyQpInF7RvTvJMkgNJdgBU1eer6jrgo8CVQ69YknRMXVd3vBO4Ffjc6w1JTgJuA34BOAQ8mmR3VX2tt8nv9p6XpJGa9lUfO/XYq+oh4KUFzRcBB6rqYFW9CuwCtmbezcD9VfWvS+0zyfYk+5LsO3z48InWL0laYJAx9nXAc32PD/XabgTeA1yR5PqlXlxVO6tqtqpm165dO0AZkqR+g9xoI4u0VVXdAtwywH4lSQMYpMd+CDi77/FZwPODlSNJGtQgwf4osCnJhiSnAFcBu5ezgyRbkuycm5sboAxJUr9OQzFJ7gIuBc5Icgj4ZFV9NskNwAPAScAdVfXUcn54Ve0B9szOzl63vLIlafmmZbZMp2Cvqm1LtO8F9g61IknSQFxSQJIaM9Zgd4xdkoZvrMFeVXuqavvMzMw4y5CkpjgUI0mNMdglqTEGuyQ1xounktQYL55KUmMcipGkxhjsktQYg12SGuPFU0lqjBdPJakxg9xBSZKa0Npyvga7JC1hUgPfi6eS1BiDXZIaY7BLUmOc7ihJjXG6oyQ1xqEYSWqMwS5JjTHYJakxBrskNcZgl6TGGOyS1BjnsUtSY5zHLkmNcShGkhpjsEtSYwx2SWqMwS5JjTHYJakxBrskNcZgl6TGGOyS1Bg/eSpJjfGTp5LUGIdiJKkxBrskNcZgl6TGGOyS1BiDXZIaY7BLUmMMdklqjMEuSY0x2CWpMQa7JDXGYJekxhjsktQYV3eUpMa4uqMkNWbNuAuQpEmzfscXFm3/+k3vG3Eli3OMXZIaY7BLUmMMdklqjMEuSY0x2CWpMQa7JDXGYJekxhjsktQYg12SGmOwS1JjDHZJaozBLkmNMdglqTEGuyQ1xmCXpMYY7JLUGG+0IUlDslpuwGGPXZIaY7BLUmOGHuxJzk3y2ST3DHvfkqTj6xTsSe5I8kKSJxe0b07yTJIDSXYAVNXBqrp2JYqVJB1f1x77ncDm/oYkJwG3AZcBFwDbklww1OokScvWKdir6iHgpQXNFwEHej30V4FdwNauPzjJ9iT7kuw7fPhw54IlScc2yBj7OuC5vseHgHVJ3pbkduAnkvz2Ui+uqp1VNVtVs2vXrh2gDElSv0HmsWeRtqqqF4HrB9ivJGkAg/TYDwFn9z0+C3h+sHIkSYMapMf+KLApyQbgG8BVwIeWs4MkW4AtGzduHKAMSVrdlvpEKqzMp1K7Tne8C3gYOD/JoSTXVtVR4AbgAeBp4O6qemo5P7yq9lTV9pmZmeXWLUlaQqcee1VtW6J9L7B3qBVJkgbikgKS1JixBnuSLUl2zs3NjbMMSWrKWIPdMXZJGj6HYiSpMQa7JDXGYJekxnjxVJIak6oadw0kOQz85wm89AzgyJDLmTTTfg6m/fjBczDNx//2qnrDKoqrIthPVJJ9VTU77jrGadrPwbQfP3gOpv34F+MYuyQ1xmCXpMZMerDvHHcBq8C0n4NpP37wHEz78b/BRI+xS5LeaNJ77JKkBQx2SWrMRAR7ks1JnklyIMmORZ5Pklt6zz+e5CfHUedK6XD8V/eO+/EkX05y4TjqXEnHOwd92/1Uku8muWKU9a20Lsef5NIkX03yVJIvjbrGldbh72AmyZ4kj/XOwTXjqHNVqKpV/QWcBPw7cC5wCvAYcMGCbS4H7mf+BtsXA/8y7rpHfPyXAG/pfX9ZS8ff9Rz0bfcPzN/85Ypx1z3i34HTga8B5/Qe//C46x7DOfgd4Obe92uBl4BTxl37OL4mocd+EXCgqg5W1avALmDrgm22Ap+reY8Apyc5c9SFrpDjHn9VfbmqvtV7+AjzNxZvSZffAYAbgb8BXhhlcSPQ5fg/BNxbVc8CVNU0noMCTksS4FTmg/3oaMtcHSYh2NcBz/U9PtRrW+42k2q5x3Yt8/+9tOS45yDJOuCXgdtHWNeodPkdOA94S5IHk+xP8pGRVTcaXc7BrcA7gOeBJ4Bfr6rvjaa81aXTPU/HLIu0LZyj2WWbSdX52JL8LPPB/tMrWtHodTkHfwx8oqq+O99ha0qX418DvBP4eeDNwMNJHqmqf1vp4kakyzn4ReCrwM8BPwr8fZJ/qqr/XuHaVp1JCPZDwNl9j89i/h15udtMqk7HluTHgc8Al1XViyOqbVS6nINZYFcv1M8ALk9ytKo+P5IKV1bXv4EjVfUK8EqSh4ALgVaCvcs5uAa4qeYH2Q8k+Q/gx4CvjKbE1WMShmIeBTYl2ZDkFOAqYPeCbXYDH+nNjrkYmKuqb4660BVy3ONPcg5wL/Dhhnpo/Y57DqpqQ1Wtr6r1wD3ArzUS6tDtb+A+4GeSrEnyA8C7gKdHXOdK6nIOnmX+PxaS/AhwPnBwpFWuEqu+x15VR5PcADzA/JXxO6rqqSTX956/nflZEJcDB4DvMP/O3YSOx/97wNuAT/d6rEerodXuOp6DZnU5/qp6OsnfAo8D3wM+U1VPjq/q4er4O/Ap4M4kTzA/dPOJqprK5XxdUkCSGjMJQzGSpGUw2CWpMQa7JDXGYJekxhjsktQYg12SGmOwS1Jj/hev4ETbeFcfjwAAAABJRU5ErkJggg==",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "for key in descriptors.keys():\n",
    "    plt.hist(descriptors[key], bins=50)\n",
    "    plt.yscale('log')\n",
    "    plt.title(key)\n",
    "    plt.show()"
   ]
  }
 ],
 "metadata": {
  "interpreter": {
   "hash": "c45e1f4cff16ded6653d5122b243e74746b1533e5e96f730958c4db3332ffa40"
  },
  "kernelspec": {
   "display_name": "Python 3.10.0 ('AIHacks')",
   "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.10.0"
  },
  "orig_nbformat": 4
 },
 "nbformat": 4,
 "nbformat_minor": 2
}