{ "cells": [ { "cell_type": "code", "execution_count": 1, "id": "604ace42-84fa-4391-bd92-0209e50a4dda", "metadata": { "tags": [] }, "outputs": [], "source": [ "import numpy as np\n", "import matplotlib.pyplot as plt\n", "from scipy.optimize import curve_fit\n", "import h5py" ] }, { "cell_type": "markdown", "id": "214268c5-b46c-48bb-8abd-d82131611df3", "metadata": {}, "source": [ "## load HDF5 file with context manager" ] }, { "cell_type": "code", "execution_count": 5, "id": "82b0c7a2-c365-49de-88c3-32a92024ac9c", "metadata": { "tags": [] }, "outputs": [], "source": [ "with h5py.File(\n", " r\"kaust_2024.nxs\", \"r\"\n", ") as f:\n", " # get the data\n", " data = f[list(f.keys())[0]][\"data\"][\"Simple_Sweep\"]\n", " SMU_mesI1 = data[\"SMU_mesI1\"][:]\n", " SMU_mesV1 = data[\"SMU_mesV1\"][:]\n", " SMU_setV1 = data[\"SMU_setV1\"][:]" ] }, { "cell_type": "markdown", "id": "95dcbf9c-53f4-40b6-923f-c5c0f8f9aa8b", "metadata": { "tags": [] }, "source": [ "## Fit the data: a * x + b" ] }, { "cell_type": "code", "execution_count": 18, "id": "6dd2f7e6-58dc-4c31-b73a-783460c8cee5", "metadata": { "tags": [] }, "outputs": [], "source": [ "fit_result = curve_fit(lambda x, a, b: a * x + b, SMU_mesV1, SMU_mesI1)\n", "a, b = fit_result[0]" ] }, { "cell_type": "markdown", "id": "07926076-9d78-4afa-9849-6997f2a6768b", "metadata": {}, "source": [ "## plot the data" ] }, { "cell_type": "code", "execution_count": 19, "id": "42375a8d-9aa9-453a-8d18-566dd5221364", "metadata": { "tags": [] }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAnYAAAHWCAYAAAD6oMSKAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAACBCklEQVR4nO3dd3gUZdfH8e+mJ5CEEkhRCMVQIihNKUoNoQkIFqoREFFUCBB4BFEEVJoFlAcVUYoFBRXBhmiogkREpHcwFIHQSagpu/f7B2/2cUkCCQQ25fe5rlwyM2dmzmQTOdwz9xmLMcYgIiIiIvmei7MTEBEREZHcocJOREREpIBQYSciIiJSQKiwExERESkgVNiJiIiIFBAq7EREREQKCBV2IiIiIgWECjsRERGRAkKFnYiIiEgBocJOpJCbNWsWFovF/uXm5kZwcDBdunRh9+7dN+28o0aNwmKx5GifCxcuMGrUKJYvX35zksrnmjRp4vBZpn+1atUqQ+yePXuIioqibNmyeHt7U7FiRWJiYjh58qRDXLly5TI9psViwcvL65o5TZ48mXr16hEQEICnpydly5alS5cubN26NdvXdf78ecaPH0/NmjUpWrQoRYoUoUaNGowdO5bz589niLdYLPTr1y/bxxcpSNycnYCI5A0zZ86kSpUqXLp0id9++40xY8awbNkyduzYQfHixXP9fE8++WSmBcfVXLhwgdGjRwOXixjJqEKFCsyePdthXbFixRyWjx8/Tr169fDz8+PVV1+lbNmyrF+/npEjR7Js2TLWrVuHi8vlf/fPnz+f5ORkh/0PHDhA586d6dix4zXzOXnyJK1bt+buu++mePHi/P3334wfP566deuybt06KleufNX9jx49SvPmzdm7dy/R0dG8/vrrACxdupTXXnuNL774gsWLFxMYGHjNXEQKAxV2IgJAtWrVqFOnDnC5aLJarYwcOZIFCxbQq1evXD/f7bffzu23357rxy3svL29qVev3lVjvv32W06ePMncuXOJiIgAoGnTpiQnJzN8+HA2btxIzZo1Aez//beff/4ZuFycX0t6IZ6ucePG1KtXj/DwcGbPns0rr7xy1f0ff/xxduzYwbJly7j//vvt6yMjI3nggQdo2rQpPXr0YNGiRdfMRaQw0K1YEclUepF39OhRh/V//vkn7du3p0SJEnh5eVGzZk2+/PJLh5gLFy4wZMgQypcvj5eXFyVKlKBOnTp88cUX9pjMbsUuXbqUJk2aULJkSby9vSlbtiwPP/wwFy5cYN++fZQqVQq4XCyk3w7s2bMncPnWYq9evQgLC8PHx4fbbruNdu3asXnzZodzLF++HIvFwhdffMGLL75ISEgIfn5+NG/enJ07d2b4PixatIiIiAj8/f3x8fGhatWqjBs3Lsffk7zE3d0dAH9/f4f16SN7V7vFaoxh5syZVKhQgWbNml3X+dM/Rze3q48t/Pnnn/zyyy/07t3boahLd//99/PEE0/w888/s27dugzbP/30U6pWrYqPjw933303P/zwg8P29J/BTZs28eijj+Lv70+JEiWIiYkhLS2NnTt30qpVK3x9fSlXrpx9tFAkL1NhJyKZio+PB6BSpUr2dcuWLeO+++7jzJkzTJ06lW+//ZYaNWrQuXNnZs2aZY+LiYnh/fffJzo6mkWLFvHpp5/y6KOPZnh+69/27dvHAw88gIeHBzNmzGDRokWMHz+eIkWKkJKSQnBwsH1Upnfv3sTFxREXF8eIESMAOHz4MCVLlmT8+PEsWrSId999Fzc3N+rWrZtpwTZ8+HD279/PRx99xLRp09i9ezft2rXDarXaY6ZPn06bNm2w2WxMnTqV77//nujoaP75558cf0+yYrVaSUtLu+aXzWa75rEA9u7dS4kSJXBzc6NixYq8+OKLXLx40SGmQ4cOlC1blsGDB7N161bOnTvHr7/+yvjx42nXrh1Vq1bN8viLFy9m//79PPHEEzl6RtJqtZKcnMyOHTt48sknKV269DVHgmNjY+35ZiV9W3psuh9//JEpU6bwyiuvMG/ePEqUKEHHjh35+++/MxyjU6dO3H333cybN48+ffowadIkBg0aRIcOHXjggQeYP38+zZo1Y+jQoXzzzTfZvmYRpzAiUqjNnDnTAOb33383qamp5uzZs2bRokUmKCjINGrUyKSmptpjq1SpYmrWrOmwzhhj2rZta4KDg43VajXGGFOtWjXToUOHq5535MiR5t//C/r6668NYDZs2JDlPsePHzeAGTly5DWvKy0tzaSkpJiwsDAzaNAg+/ply5YZwLRp08Yh/ssvvzSAiYuLM8YYc/bsWePn52fuv/9+Y7PZsjxPdr8nWQkNDTXANb+yc80vvviiee+998zSpUvNjz/+aPr162fc3NxMo0aNMuRx+PBhU79+fYdzPProo+bSpUtXPUfnzp2Nq6ur+eeff66Zz795enraz1OpUiWzbdu2a+7Tt29fA5gdO3ZkGbN9+3YDmGeeeca+DjCBgYEmKSnJvi4hIcG4uLiYcePG2del/wy+9dZbDsesUaOGAcw333xjX5eammpKlSplHnrooWxdr4iz6Bk7EQHI8FxW1apV+fbbb+23y/bs2cOOHTt48803AUhLS7PHtmnThh9++IGdO3dStWpV7r33XmbPns2wYcNo1aoVdevWxdvb+6rnr1GjBh4eHjz11FM8++yzNGzYkAoVKmQ7/7S0NF5//XU+++wz9uzZQ2pqqn3b9u3bM8S3b9/eYfmuu+4CYP/+/dSrV4/Vq1eTlJTEs88+m+XIVE6+J1n5/vvvM0xOyExISMg1Y1577TWH5TZt2lCuXDmGDBnCt99+a5/scPr0aR588EEuXLjA7NmzKVOmDFu2bOHVV1+lffv2/Pjjj5neJj116hQLFiygVatW3HbbbdfM599Wr15NSkoKe/fuZdKkSTRt2pQlS5Zw55135ug4VzLGAGT4jJo2bYqvr699OTAwkNKlS7N///4Mx2jbtq3DctWqVdm4cSOtW7e2r3Nzc+OOO+7IdH+RvESFnYgA8Mknn1C1alXOnj3L3Llz+eCDD+jatSs//fQT8L9n7YYMGcKQIUMyPcaJEyeAyy0ubr/9dubOncuECRPw8vKiZcuWvPHGG4SFhWW6b8WKFVm8eDGvv/46zz33HOfPn6dChQpER0czYMCAa+YfExPDu+++y9ChQ2ncuDHFixfHxcWFJ598MsOtSICSJUs6LHt6egLYY48fPw5w1QkeOfmeZCU8PNxenFxN+izVnHrssccYMmQIv//+u72wmzBhAhs2bGD//v0EBwcD0LBhQ6pUqUKzZs2YPXs2PXr0yHCszz77jOTk5GxNmrhSrVq1gMv/gGjfvj133HEHw4cP59tvv81yn7JlywKXHwvIavbsvn37AChTpozD+is/X7j8GWf2s1CiRAmHZQ8PD3x8fDI8a+jh4UFSUlKW+YrkBSrsRAS4PEqRPmGiadOmWK1WPvroI77++mseeeQRAgICAHjhhRd46KGHMj1G+l++RYoUYfTo0YwePZqjR4/y008/MWzYMNq1a8eOHTuyzKFhw4Y0bNgQq9XKn3/+yX//+18GDhxIYGAgXbp0uWr+n332GY8//jhjx451WH/ixIkM7T6yI/0B/38/T3elnHxPslKxYsVsjQKNHDmSUaNGXTMuK/8uDDds2MBtt91mL+rS3XPPPQBs2bIl02NMnz6dwMDADCNcOeXr60uVKlXYtWvXVeMiIyMZPny4fZQwMwsWLLDHiogKOxHJwuuvv868efN4+eWXeeihh6hcuTJhYWFs3LgxQ/F0NYGBgfTs2ZONGzfy9ttvc+HCBXx8fK66j6urK3Xr1qVKlSrMnj2bv/76iy5dumQYVfs3i8Vi357uxx9/5NChQ9xxxx3ZzjddgwYN8Pf3Z+rUqXTp0iXT27HX+z35t9y8FZuZjz/+GHC81R4SEsKSJUs4dOiQwy3VuLg4IPNRyj///JNNmzbx/PPPX3M267WcOHGCzZs3c9999101rk6dOrRo0YLp06cTFRWVIX7VqlXMmDGDVq1aUbt27RvKSaSgUGEnIpkqXrw4L7zwAs8//zyff/45jz32GB988AGtW7emZcuW9OzZk9tuu41Tp06xfft2/vrrL7766isA6tatS9u2bbnrrrsoXrw427dv59NPP6V+/fpZFnVTp05l6dKlPPDAA5QtW5ZLly4xY8YMAJo3bw5cHukJDQ3l22+/JSIighIlShAQEEC5cuVo27Yts2bNokqVKtx1112sW7eON95447p75RUtWpS33nqLJ598kubNm9OnTx8CAwPZs2cPGzduZMqUKQDZ/p5kpXr16teV35VWrlzJmDFj6NixIxUqVODSpUv89NNPTJs2jWbNmtGuXTt77HPPPcfs2bOJjIxk2LBh9mfsXnvtNQIDA+nevXuG40+fPh24PCM5K+kF9J49ewBITEwkMjKSbt26ERYWhre3N7t27eKdd94hOTmZkSNHXvO6PvnkE5o3b06LFi2Ijo62991bunQp77zzDlWqVMnW7GORQsPZszdExLnSZ8WuXbs2w7aLFy+asmXLmrCwMJOWlmaMMWbjxo2mU6dOpnTp0sbd3d0EBQWZZs2amalTp9r3GzZsmKlTp44pXry48fT0NBUqVDCDBg0yJ06csMdcOSs2Li7OdOzY0YSGhhpPT09TsmRJ07hxY/Pdd9855LR48WJTs2ZN+yzLHj16GGOMOX36tOndu7cpXbq08fHxMffff79ZuXKlady4sWncuLF9//RZsV999ZXDcePj4w1gZs6c6bB+4cKFpnHjxqZIkSLGx8fHhIeHmwkTJjjEZOd7crPt3r3btGnTxtx2223G09PTeHl5merVq5sxY8ZkOtP1r7/+Mh07djS33367/TN68sknzYEDBzLEXrhwwfj7+5tGjRpdNYfQ0FATGhpqX7506ZJ58sknTdWqVU3RokWNm5ubuf32281jjz1mtm7dmu1rO3funBk7dqypUaOG8fHxMT4+Puauu+4yr732mjl37lyGeMA899xzmeaX/vNizP9+Bo8fP+4Q16NHD1OkSJEM+zdu3Njceeed2c5bxBksxmTjqV0RERERyfPUoFhERESkgFBhJyIiIlJAqLATERERKSBU2ImIiIgUECrsRERERAoIFXYiIiIiBYQaFOchNpuNw4cP4+vrm+VLx0VERKRwMcZw9uxZQkJCrvneaBV2ecjhw4czvMhaREREBODgwYPXfJuOCrs8xNfXF7j8wfn5+Tk5GxEREckLkpKSKFOmjL1OuBoVdnlI+u1XPz8/FXYiIiLiIDuPaWnyhIiIiEgBocJOREREpIBQYSciIiJSQOgZu3zGZrORkpLi7DQkl7i7u+Pq6ursNEREpIBQYZePpKSkEB8fj81mc3YqkouKFStGUFCQeheKiMgNU2GXTxhjOHLkCK6urpQpU+aaDQol7zPGcOHCBY4dOwZAcHCwkzMSEZH8ToVdPpGWlsaFCxcICQnBx8fH2elILvH29gbg2LFjlC5dWrdlRUTkhmjYJ5+wWq0AeHh4ODkTyW3phXpqaqqTMxERkfxOhV0+o+ewCh59piIikltU2ImIiIgUECrsRG6i5cuXY7FYOHPmjLNTERGRQkCFnYiIiEgOTYrdxeQluzPdNnnJbibF7rrFGV2mwq6QyKs/gPmBMYa0tDRnpyEiInmIq4uFiZn83Tp5yW4mxu7C1cU5z0+rsCsknPUD2KRJE/r378/AgQMpXrw4gYGBTJs2jfPnz9OrVy98fX2pWLEiP/30k32fbdu20aZNG4oWLUpgYCBRUVGcOHHCvn3RokXcf//9FCtWjJIlS9K2bVv27t1r356SkkK/fv0IDg7Gy8uLcuXKMW7cOAD27duHxWJhw4YN9vgzZ85gsVhYvnw58L/bpz///DN16tTB09OTlStXYozh9ddfp0KFCnh7e3P33Xfz9ddfO1zvwoULqVSpEt7e3jRt2pR9+/bl/jdVREScLjoijJjISv/7uzX1kv3v1JjISkRHhDklLxV2hUSGH0C4ZT+AH3/8MQEBAfzxxx/079+fZ555hkcffZQGDRrw119/0bJlS6Kiorhw4QJHjhyhcePG1KhRgz///JNFixZx9OhROnXqZD/e+fPniYmJYe3atSxZsgQXFxc6duxofyPH5MmT+e677/jyyy/ZuXMnn332GeXKlctx3s8//zzjxo1j+/bt3HXXXbz00kvMnDmT999/n61btzJo0CAee+wxVqxYAcDBgwd56KGHaNOmDRs2bODJJ59k2LBhufI9FBGRvCc6IoyY5mH8s/QDjr5WhS8Xr3JqUQeAkTwjMTHRACYxMTHDtosXL5pt27aZixcv3tA53lm8y4QO/cGEDV9oQof+YN5ZvOuGjnctjRs3Nvfff799OS0tzRQpUsRERUXZ1x05csQAJi4uzowYMcK0aNHC4RgHDx40gNm5c2em5zh27JgBzObNm40xxvTv3980a9bM2Gy2DLHx8fEGMOvXr7evO336tAHMsmXLjDHGLFu2zABmwYIF9phz584ZLy8vs3r1aofj9e7d23Tt2tUYY8wLL7xgqlat6nDeoUOHGsCcPn06y+9Rbn22IiJyi11MNOarJ4wZ6WfMSD/z0Utdb8pprlYfXEkjdoVMdEQYHq4upFhteLi63JJ/Vdx11132P7u6ulKyZEmqV69uXxcYGAhcfvvCunXrWLZsGUWLFrV/ValSBcB+u3Xv3r1069aNChUq4OfnR/ny5QE4cOAAAD179mTDhg1UrlyZ6Ohofvnll+vKu06dOvY/b9u2jUuXLhEZGemQ2yeffGLPa/v27dSrV8+hL139+vWv69wiIpLHHd4A0xrDlq9JMy68ae3Ka6lds3ye/VbRK8UKmclLdtuLuhSrjclLdt/04s7d3d1h2WKxOKxLL4RsNhs2m4127doxYcKEDMdJf5dqu3btKFOmDB9++CEhISHYbDaqVatGSkoKALVq1SI+Pp6ffvqJxYsX06lTJ5o3b87XX39tf8euMcZ+3Kze+FCkSBH7n9Nv8/7444/cdtttDnGenp4ZjikiIgWUMbDmA4gdAdYU/jEBrK4xgSEdH8Hj/x9xApx2O1aFXSFy5TN1k/PAD+CVatWqxbx58yhXrhxubhl/PE+ePMn27dv54IMPaNiwIQCrVq3KEOfn50fnzp3p3LkzjzzyCK1ateLUqVOUKlUKgCNHjlCzZk0Ah4kUWQkPD8fT05MDBw7QuHHjLGMWLFjgsO7333+/5rFFRCRvmfT/kwoz/N144RR7p/ei4snlAPxsrcO++ybwdKvLd3jS4535d6sKu0Iis4kSeeEH8ErPPfccH374IV27duU///kPAQEB7Nmzhzlz5vDhhx9SvHhxSpYsybRp0wgODubAgQMZJihMmjSJ4OBgatSogYuLC1999RVBQUEUK1YMFxcX6tWrx/jx4ylXrhwnTpzgpZdeumZevr6+DBkyhEGDBmGz2bj//vtJSkpi9erVFC1alB49etC3b1/eeustYmJiePrpp1m3bh2zZs26Sd8pERG5WdI7ScC//m48sIak2Y9TMTkBq8WdX8sNYOdtnYluXslh3/R4q805d3FU2BUSVpvJdKaOs38ArxQSEsJvv/3G0KFDadmyJcnJyYSGhtKqVStcXFywWCzMmTOH6OhoqlWrRuXKlZk8eTJNmjSxH6No0aJMmDCB3bt34+rqyj333MPChQvtt2FnzJjBE088QZ06dahcuTKvv/46LVq0uGZur776KqVLl2bcuHH8/fffFCtWjFq1ajF8+HAAypYty7x58xg0aBDvvfce9957L2PHjuWJJ564Kd8rERG5ORwGPoyNaM8fsS15DT+snPEqQ7HHP6NpSA2aXmN/Z7AYPRiUZyQlJeHv709iYiJ+fn4O2y5dukR8fDzly5fHy8vLSRnKzaDPVkQkb/rwp9+pvHoIjVw3A7CjVCuqPPkRePre0jyuVh9cSbNiRURERK709wr6bH2cRq6buWg8eCHtaao8O+eWF3U5pcJOREREJJ01DZaOgU8ehHNH2Wm7nYfTxvBFWmMmL93j7OyuSc/YiYiIiAAkHoJ5T8KB1QB8ntaUM41eZWGL6nmyk0RmVNiJiIiI7PoZ5veFi6dIcS3C4Iu9CIvomac7SWRGhZ2IiIgUXmkpsGQ0xE25vBxcg8+CXyasSNk830kiMyrsREREpFDI0Hj4VDx8/QQc/guAv4K7Uqv3Ozzh5pnlMfLqSF06FXYiIiJSKDg0Hg7cDN8PgOQkLrn50e9CH+66oyu1rlLU5Qcq7ERERKRQiI4Iw9V6iWLLhoLbEgAO+93NI8d60yWyQZ4fjcsOFXYiIiJSOBzfxXN7nga3rdiMhQ9sD/LmsYcZEFm1QBR1oD52cpMZY3jqqacoUaIEFouFYsWKMXDgQGenJSIihc2Gz2FaYzi2FYqU4gnrC0xI7YSrq3uBKeogHxZ27733nv3VS7Vr12blypVXjV+xYgW1a9fGy8uLChUqMHXq1Awx8+bNIzw8HE9PT8LDw5k/f77D9l9//ZV27doREhKCxWJhwYIFGY7Rs2dPLBaLw1e9evVu6FoLgkWLFjFr1ix++OEHjhw5wq5du3j11Vft28uVK8fbb7/tvARFRKRAmBS7i8lLdmfckHyO7e91hQXPQOoFKN+Yj6p9yvK0ani4upBitWW+Xz6Vrwq7uXPnMnDgQF588UXWr19Pw4YNad26NQcOHMg0Pj4+njZt2tCwYUPWr1/P8OHDiY6OZt68efaYuLg4OnfuTFRUFBs3biQqKopOnTqxZs0ae8z58+e5++67mTJlylXza9WqFUeOHLF/LVy4MHcuPB/bu3cvwcHBNGjQgKCgIEqXLo2vb95+HYuIiOQ/6RMjHIq0I5s4/XZ9qh5biA0XaPYS/73tdV5bcYqYyErsGtOamMhKGffLxyzGmLzbjOUKdevWpVatWrz//vv2dVWrVqVDhw6MGzcuQ/zQoUP57rvv2L59u31d37592bhxI3FxcQB07tyZpKQkfvrpJ3tMq1atKF68OF988UWGY1osFubPn0+HDh0c1vfs2ZMzZ85kOpqXXVd7yW9+fFF8z549+fjjj+3LoaGhlCtXjho1avD222/TpEkTVqxY4bBPPvpxzDX58bMVEcmL0t8OEdM8jGi/FaT9NBw3k8JZj0B8u3/M5D0Bl7dHVnK4/Wrf74r1ecXV6oMr5ZvJEykpKaxbt45hw4Y5rG/RogWrV6/OdJ+4uDhatGjhsK5ly5ZMnz6d1NRU3N3diYuLY9CgQRliruf24PLlyyldujTFihWjcePGjBkzhtKlS2cZn5ycTHJysn05KSkp+ycz5vKQsjO4+4DFcs2wd955h4oVKzJt2jTWrl2Lq6srjz76qH37N998w913381TTz1Fnz59bmbGIiJSCERHhOGRdpayK54D1z9wA+JLNKT8k5+ATwmsuzIv3vJD4+HsyjeF3YkTJ7BarQQGBjqsDwwMJCEhIdN9EhISMo1PS0vjxIkTBAcHZxmT1TGz0rp1ax599FFCQ0OJj49nxIgRNGvWjHXr1uHpmXlPnHHjxjF69Ogcnccu9QKMDbm+fW/U8MPgUeSaYf7+/vj6+uLq6kpQUFCG7SVKlMDV1RVfX99Mt4uIiOTIwbX03fYEuB4gxbjyhq07L/afbB+MGBRZKctd8+JI3fXIN4VdOssVI0XGmAzrrhV/5fqcHjMznTt3tv+5WrVq1KlTh9DQUH788UceeuihTPd54YUXiImJsS8nJSVRpkyZHJ1XRESk0LPZIO6/sOQVsKWx31aaGNsA1qWVx3fpngJTtGVHvinsAgICcHV1zTCSduzYsQwjbumCgoIyjXdzc6NkyZJXjcnqmNkVHBxMaGgou3dn/TCmp6dnlqN51+Tuc3nkzBncfZxzXhERKZQyvArsXz786Q+a7xpJ+dOXH8v6wVqPf+4fz7yWNe3PzkHBGZG7lnxT2Hl4eFC7dm1iY2Pp2LGjfX1sbCwPPvhgpvvUr1+f77//3mHdL7/8Qp06dXB3d7fHxMbGOjxn98svv9CgQYMbyvfkyZMcPHiQ4ODgGzpOliyWbN0Ozes8PDywWq3OTkNERPIwh1eB/atAmzfvC9pteoEgy2nSXDwZkfwYwU37Et28kkNsYSru8k1hBxATE0NUVBR16tShfv36TJs2jQMHDtC3b1/g8q3NQ4cO8cknnwCXZ8BOmTKFmJgY+vTpQ1xcHNOnT3eY7TpgwAAaNWrEhAkTePDBB/n2229ZvHgxq1atssecO3eOPXv22Jfj4+PZsGEDJUqUoGzZspw7d45Ro0bx8MMPExwczL59+xg+fDgBAQEORahkVK5cOX799Ve6dOmCp6cnAQEBzk5JRETymAwFWtMK/D5rGB32f4irxUBAJT6/fTTBvmEFemJEduSrwq5z586cPHmSV155hSNHjlCtWjUWLlxIaGgoAEeOHHHoaVe+fHkWLlzIoEGDePfddwkJCWHy5Mk8/PDD9pgGDRowZ84cXnrpJUaMGEHFihWZO3cudevWtcf8+eefNG3a1L6c/lxcjx49mDVrFq6urmzevJlPPvmEM2fOEBwcTNOmTZk7d656tl3DK6+8wtNPP03FihVJTk4ulO1ORETk2tILtM9i13Dvisep57INLECNx6DN6zx+lbtYhWGkLl2+6mNX0BW0PnaSPfpsRUSyafdiTn3WkxKWs5w3nhR56L9wd+dr75fP5aSPXb5684SIiIgUQtZU+GUEzH6YEpazbDPlaJsylsknajk7szwnX92KFRERkULm9H6Y1xv+WQvAxuBHufuJKXT89WChmhSRXSrsREREJG/a9h181w8uJZJkfFgZPooHOj8NFM4Zr9mhwk5EREScIsv+dKmX2DijH3cf+QqAw77V+bnKGHo90NghrLDNeM0OFXYiIiLiFJn2pzuxh+Mzu3L3+cvruW8AIc1G0MvVPdNjaKTOkQq7fEaTmAsem83m7BRERJwiw+3UgHWkfDuQUraLXHAvjk+njyCsuTNTzHdU2OUT7u7uWCwWjh8/TqlSpXL8LlvJe4wxpKSkcPz4cVxcXPDw8HB2SiIit1x0RBhu1ouUWhYDbr/iAfzjX5vbe88Gv5v09qYCTIVdPuHq6srtt9/OP//8w759+5ydjuQiHx8fypYti4uLug+JSCGUsIVnd/UGt11YjYUptkcYMOADcHF1dmb5kgq7fKRo0aKEhYWRmprq7FQkl7i6uuLm5qYRWBEpfIyBdTNh0QuQdokEU5zB1v78llYFy7K/9ezcdVJhl8+4urri6qp/xYiISD52KRG+HwBb5wOwzHo3exq8yezW9zJ5yW61MLkBKuxERETk1jm0Dr5+Ak7vw2pxZXxKZ3ybDiS6eWVA/elulAo7ERERyXUZetQZA7+/B7EjwZZKomcICyuPwde/eobiTf3prp8KOxEREcl1Dj3q6peEBc/ArkUALLTey4FaE+jbMut3vWqk7vqosBMREZFcl16YrVz8Hb3WfIBvylHSLB6MSulO6abPEt28kpMzLJjUX0FERERyn81KtNsC5nq+hm/KUf42wbS/NIrSzZ5TUXcTacROREREctfZo/BNH4hfgQuwwNaQ4Sm9SHP10S3Wm0wjdiIiIpJ79i6FqfdB/Apw9+GXSqMYmPIMaa4+pFhtTF6y29kZFmgasRMREZEbZ02DZWNg1STAQGA1Pr19FCN+SyEmshLREWHqUXcLqLATERGRbMvQxgTgzEGY1xsOrrm8XOcJ3vV4gjeWHrAXdaAedbeCCjsRERHJNoc2JhFhsGPh5VYml86QZLxZVXUkbdo+Q0rsLoeiLp161N1cFmOMvrN5RFJSEv7+/iQmJuLn5+fsdERERDI1eclupsRuZW6FRdQ8/AUAG2wV2HDvRHq2berk7AqenNQHGrETERGRHImu4UKX9eMofXgHANPSHiCl8Uv0iwx3cmaiwk5ERESALJ6f+3+Tl+zGajMMCtoE3w+kdMpZThtfYlL78pulNrtU1OUJanciIiIiwP+en7uyJcnkJbt5L3YzkXvHXJ4kkXKWQ361aJ08lt8stdXGJA/RiJ2IiIgAmc9anbxkN98vXspvJaZSMuFvwMKask/SbVcjBkRWVRuTPEaFnYiIiNj9u7ibsnQ3HVjKQu9PcL+QDEUD+ab8KGLW+quNSR6lwk5EREQcREeEMXPpJka5fMSDrqvBBlRsBh2nsT/uNDGRGZ/DUxuTvEHtTvIQtTsREZG84Itvv6feusGUdzlKmnFhTflnue/xV8FFj+Y7Q07qA31CIiIicpkxrPjkVR76qyflXY6Cfxm+qfEh3Xc0YPKyvc7OTrJBt2JFREQELpxi7/ReND65HCxAlbbQ/r908ilBgp8mR+QXKuxEREQKuwNrYF5vKiYexGpxx7XVGLj3KbBYAD0/l5+osBMRESngsmw8bLPx2ycjqL/vfVywQokKuD4yE0JqZDiGRuryBxV2IiIiBVx642H4V4F27hj7P4rivjO/X16u/ii0nQSevk7KUnKDCjsREZECLkOfufKHOf9FL0JTT5Lq4ol727eg5mP2W6+Sf6mwExERKQSiI8KwmDRsy8Zg+3UBRSyGEz4VCej5OZSu4uz0JJeosBMRESkMEg/R/0AMuK0GYK6tGZ0HzgYPHycnJrlJfexEREQKup2LYOr9cGA1Z403g9KiGZryJJNXHnJ2ZpLLNGInIiJSUKWlwJLREDcFgE228qy/dyKT2jaj/BL1piuIVNiJiIjkY1m2MjkVT8KMbgSd2wbA9LTWXGo8guci7wQymVCh4q5AUGEnIiKSj2XaymTLNyTP70eQ9TyX3Pz4+Y6XOV+qSYbiTY2HC55894zde++9R/ny5fHy8qJ27dqsXLnyqvErVqygdu3aeHl5UaFCBaZOnZohZt68eYSHh+Pp6Ul4eDjz58932P7rr7/Srl07QkJCsFgsLFiwIMMxjDGMGjWKkJAQvL29adKkCVu3br2haxUREbmW6IgwYiIrMTF2F+/+shm+Hwhf98LTep7Dfnfj1W81D3bpk+WIXHREGIMiK93apOWmyVeF3dy5cxk4cCAvvvgi69evp2HDhrRu3ZoDBw5kGh8fH0+bNm1o2LAh69evZ/jw4URHRzNv3jx7TFxcHJ07dyYqKoqNGzcSFRVFp06dWLNmjT3m/Pnz3H333UyZMiXL3F5//XUmTpzIlClTWLt2LUFBQURGRnL27Nnc+waIiIhkIjoijFfvc6fZqq6wbiY2Y2FtmV6EDFgKxco4Oz25hSzGmHwz/lq3bl1q1arF+++/b19XtWpVOnTowLhx4zLEDx06lO+++47t27fb1/Xt25eNGzcSFxcHQOfOnUlKSuKnn36yx7Rq1YrixYvzxRdfZDimxWJh/vz5dOjQwb7OGENISAgDBw5k6NChACQnJxMYGMiECRN4+umns3V9SUlJ+Pv7k5iYiJ+fX7b2ERGRQs4Y2PA5LBwCqRc4bvz4j7Ufs14b6uzMJJfkpD7INyN2KSkprFu3jhYtWjisb9GiBatXr850n7i4uAzxLVu25M8//yQ1NfWqMVkdMzPx8fEkJCQ4HMfT05PGjRvn6DgiIiI5knwW5j8N3z4LqRdYZb2TB9NeZ3laNSYv2e3s7MQJ8s3kiRMnTmC1WgkMDHRYHxgYSEJCQqb7JCQkZBqflpbGiRMnCA4OzjImq2NmdZ70/a48zv79+7PcLzk5meTkZPtyUlJSts8pIiKF3JFN8HUvOLkHGy68lfoIXk0Hs7p5FSarlUmhlW8Ku3SWK95jZ4zJsO5a8Veuz+kxcyu3cePGMXr06ByfR0RECr4s25gYw/LPxnL/32/jZlI46xFIr7NP06h5e/r/f6xamRRe+aawCwgIwNXVNcNI2rFjxzKMlKULCgrKNN7NzY2SJUteNSarY2Z1Hrg8chccHJzt47zwwgvExMTYl5OSkihTRg+5iohIFm1MLp5hz0c9aXJy2eXlSq35rFgMjTyLq5WJAPmosPPw8KB27drExsbSsWNH+/rY2FgefPDBTPepX78+33//vcO6X375hTp16uDu7m6PiY2NZdCgQQ4xDRo0yHZu5cuXJygoiNjYWGrWrAlcfiZwxYoVTJgwIcv9PD098fT0zPZ5RESk8Mgw6lbpDEmfRXFH8hGsFjdcW7wK9Z7hmavcGdJIXeGTbwo7gJiYGKKioqhTpw7169dn2rRpHDhwgL59+wKXR8AOHTrEJ598AlyeATtlyhRiYmLo06cPcXFxTJ8+3WG264ABA2jUqBETJkzgwQcf5Ntvv2Xx4sWsWrXKHnPu3Dn27NljX46Pj2fDhg2UKFGCsmXLYrFYGDhwIGPHjiUsLIywsDDGjh2Lj48P3bp1u0XfHRERKWiiI8LA2Di39C1Sf/0SP4uVRK/b8I/6DG6r5ez0JC8y+cy7775rQkNDjYeHh6lVq5ZZsWKFfVuPHj1M48aNHeKXL19uatasaTw8PEy5cuXM+++/n+GYX331lalcubJxd3c3VapUMfPmzXPYvmzZMgNk+OrRo4c9xmazmZEjR5qgoCDj6elpGjVqZDZv3pyja0tMTDSASUxMzNF+IiJSQJ07bsynDxsz0s+YkX7mhxEtjbl4xtlZyS2Wk/ogX/WxK+jUx05EROziV8I3feDsES4Zd8ZYe/BpWlNiIivrFmshk5P6IF/dihURESnwbFZY8Tr8+joYG7ttt/FHnbd4tX1rSqmNiVyDCjsREZG8IukwzOsD+y8/5z03rQknG73Ksy3uAtTGRK5NhZ2IiMgtlmmPut2xl98iceEkKS4+LKk4jKNBbdTGRHJEhZ2IiMgt5tCjrkk5WPIKrJ4MwFZbKOtqT+Txds1pncX+GqmTrKiwExERucXSC7MvF/9Gpw09CTq7BYBZaS0433gkz0VWc2Z6ko+psBMREXGC6OBtPO3zIp5nz5FoivB8ah/ujHhMo3FyQ1TYiYiI3Eqpl+CXF2HtR3gC621h9E/txzGXQD5QUSc3yMXZCYiIiBQaJ/bAR81h7UcA/Hl7Dx5NGcExl0BSrDYmL9nt5AQlv9OInYiISC7JdLbr//vl87dpsmc8HraL4BPAt+VfZsC6AGIiKxEdEcZk9aiTXKDCTkREJJc4zHZNL85SzrNt+tO0OPr95eVyDZleejiv/nraXtT9O17FndwIFXYiIiK5JENxdmcypz7uTvjFfdhwwaXJMGg0hKQle4mJLKUedZLr9K7YPETvihURKRgmL95FwrKpvOz2CV6WVM55lKJot1lQ7n5npyb5kN4VKyIi4iyXEok+PRbc5wOw3FaTJgO+hiIBTk5MCgPNihUREckth9bBB41g63xSjSvjrY/RK2Uwk38/7ezMpJDQiJ2IiMiNMgbi3oXFo8CWygFbKeJqvcGwDh3x0WxXuYVU2ImIiGRDlq1Mzp8kfvrjlD+1CoAfrfdy8P4J9G1ZC9BsV7m1VNiJiIhkQ6atTPav5tzsHpRPOUaaxYMV5WPYG/Iw0c0rOeyr2a5yq6iwExERyYZ/j7xZjJX+7t9jWzaWotg45R1KiR6ziQiqTsQ19he5mVTYiYiIZFN0RBg+KSeouuIJcN2KC7CtdFvCe38AnkWdnZ6IZsWKiIhk257FPLn5Me5z3cp548nzac8S/uxsFXWSZ6iwExERuRZr6uUZr589DBdOsN1WlofTxvJl2v1MXrLb2dmJ2OlWrIiIyNWcOQBf94Z//gDg07TmJDUaxaIW1ZmsViaSx6iwExGRQi3LNibAD3M/pPmu0XhZz5LsWoSBF5+kakSUPVatTCSvUWEnIiKFWqZtTNKS2TCjP20Pz728fFttPgsaQVWf2zMUb2plInmJCjsRESnUMoy61XDh2Mxu1Di343JAg/7Q7GV6u3lc8xgizqbCTkRECr30wmzPklmc+/UjSlsucdHNH+9OH0GlFk7OTiT7VNiJiIikXCD63Dvg8SkAf9iqcm/0N+AX4uTERHJG7U5ERKRwO7YdPmwK6z/FZixMsT5M15ThTF573tmZieSYRuxERKRwMgb++gR+GgppFzlqivHbXePp93BXbGpjIvmUCjsRESmwsmxlcimJndN7U/n4LwCssN7FrgZv0qd1XUBtTCT/UmEnIiIFVqatTA6v58wnj1H50j/YcOW30GfZWCaK6OaVHfZVGxPJj1TYiYhIgeUw8mYM0UWXYP35JYqZNJI8g/B77FMalrmXhtfYXyS/UGEnIiIFWnREGJ6piZRf0Rdc1+EK7CnZlDuenAnexZ2dnkiuUmEnIiIF24HfeXpbb3D9h2TjxgTbY7zc722wWJydmUiuU7sTEREpmGw2WPkWzGwDSf8QbwukU9qrzEhtweSle5ydnchNoRE7EREpeM4dg2+egr+XATDfeh8J94/l25Y1mKxWJlKAqbATEZF8KctWJnuXcX5Ob4qkniTVxYvhyY9TpmkfoptXAtTKRAo2FXYiIpIvZWhlYk2D5eMwK9+iCIYTPhX5sfJYyhStmKF4UysTKahU2ImISL7075G3oslHeSLhNTgQhwXYHNSR6r3fp4e79zX3FylIVNiJiEi+FR0RRrmTK2m45imwnOOs8WZV1RG07vKcs1MTcYp8Nyv2vffeo3z58nh5eVG7dm1Wrlx51fgVK1ZQu3ZtvLy8qFChAlOnTs0QM2/ePMLDw/H09CQ8PJz58+fn+Lw9e/bEYrE4fNWrV+/GLlZERLKWlgKLhtN+2yCKW86xyVaeDmnjVNRJoZavCru5c+cycOBAXnzxRdavX0/Dhg1p3bo1Bw4cyDQ+Pj6eNm3a0LBhQ9avX8/w4cOJjo5m3rx59pi4uDg6d+5MVFQUGzduJCoqik6dOrFmzZocn7dVq1YcOXLE/rVw4cKb840QESnsTv0NM1rA7+8CMD2tNV2tr7DXWprJS3Y7OTkR57EYY/LNk6N169alVq1avP/++/Z1VatWpUOHDowbNy5D/NChQ/nuu+/Yvn27fV3fvn3ZuHEjcXFxAHTu3JmkpCR++ukne0yrVq0oXrw4X3zxRbbP27NnT86cOcOCBQuu+/qSkpLw9/cnMTERPz+/6z6OiEiBtmUefD8QkpO46OZPvwtPcndEV6IjwuytTGIiK+kZOikwclIf5JsRu5SUFNatW0eLFi0c1rdo0YLVq1dnuk9cXFyG+JYtW/Lnn3+Smpp61Zj0Y+bkvMuXL6d06dJUqlSJPn36cOzYsZxfqIhIITcpdlfmo26pF9k8tSd8/QQkJ3HY726anXvNXtTB5WfuYiIrMTGrY4gUcPlm8sSJEyewWq0EBgY6rA8MDCQhISHTfRISEjKNT0tL48SJEwQHB2cZk37M7J63devWPProo4SGhhIfH8+IESNo1qwZ69atw9PTM9P8kpOTSU5Oti8nJSVd47sgIlLwZWhjAnBsBydmdaP6hb0YLFgaDuZL68N0dXVXKxORf8k3hV06yxXv9jPGZFh3rfgr12fnmNeK6dy5s/3P1apVo06dOoSGhvLjjz/y0EMPZZrbuHHjGD16dJa5i4gURg4NhI0huuQfpH4/mADbJc67l6RIl+lQsSkDs3EMkcIm3xR2AQEBuLq6ZhidO3bsWIbRtHRBQUGZxru5uVGyZMmrxqQf83rOCxAcHExoaCi7d2d9K+CFF14gJibGvpyUlESZMmWyjBcRKSyiI8JwTztP4IpB4LoKd+BAsXsp2/sz8M36/70ihV2+ecbOw8OD2rVrExsb67A+NjaWBg0aZLpP/fr1M8T/8ssv1KlTB3d396vGpB/zes4LcPLkSQ4ePEhwcHCWMZ6envj5+Tl8iYgIcGQTz+x8godcV5FmXJho7UzZ6J9V1IlcQ74ZsQOIiYkhKiqKOnXqUL9+faZNm8aBAwfo27cvcHkE7NChQ3zyySfA5RmwU6ZMISYmhj59+hAXF8f06dPts10BBgwYQKNGjZgwYQIPPvgg3377LYsXL2bVqlXZPu+5c+cYNWoUDz/8MMHBwezbt4/hw4cTEBBAx44db+F3SEQknzMG1n4EPw8HawqHTEmGWPsTl1YJt2V7dYtV5BryVWHXuXNnTp48ySuvvMKRI0eoVq0aCxcuJDQ0FIAjR4449JYrX748CxcuZNCgQbz77ruEhIQwefJkHn74YXtMgwYNmDNnDi+99BIjRoygYsWKzJ07l7p162b7vK6urmzevJlPPvmEM2fOEBwcTNOmTZk7dy6+vr636LsjIpLPXTwN3/aDHT8AEGutzd/3vc4XrerY25iAnp8TuZp81ceuoFMfOxEp6CbF7sLVxZKxODu4lqTPovBLPoLV4sZrKV0p3jSa6OaV7CHqUSeFVU7qg3w1YiciIvlbhlYmNhusnox1ySv4GStnvG5nYaUxFPe/U21MRK6DCjsREbll/t3KxDvlNH1OToA9i3EFdpZqQeXe0+nmlfWIhEbqRK5OhZ2IiNxS0RFh3HbmT+7//VmwnOGi8WB1peeJ6DYErtKXVESuTYWdiIjcOjYrrJjAw5tfB4thl+02BlkH8mP3p5ydmUiBkG/62ImISD6XdBg+bg8rJgCGOWlNeMQ6lq3W2/ReV5FcohE7ERG5+Xb9Agv6woWTpLj6MORiL+6I6MWmiDC1MhHJRSrsRETkhmXZxiQthXUzY6h96FMAjhWpTKdTT/FQ88b2WId3w6LiTuRGqLATEZEblqGNCcDpfSTM6Ebts1svL9/7NHNcH+chN0+1MhG5SVTYiYjIDcsw6ha0leRvniPIeo5Lbr54PTwVqrYlOhvHEJHrl6PCLjExkfnz57Ny5Ur27dvHhQsXKFWqFDVr1qRly5Y0aNDgZuUpIiJ5XHREGK62ZPyWDQO3xXgCR3zvIrj3bChW1tnpiRQK2ZoVe+TIEfr06UNwcDCvvPIK58+fp0aNGkRERHD77bezbNkyIiMjCQ8PZ+7cuTc7ZxERyYuO7+K5PU8T5bYYgA+sDxI8cKmKOpFbKFsjdnfffTePP/44f/zxB9WqVcs05uLFiyxYsICJEydy8OBBhgwZkquJiohIHrbhC/hxMKSe54Tx43nrcyxNq07y8n26xSpyC2WrsNu6dSulSpW6aoy3tzddu3ala9euHD9+PFeSExGRPC75HCwcAhu/AOA3651sb/AWM1rXVxsTESfIVmF3raIOwGq18v3339OhQ4dsxYuISP6QZSuThC2c+rgbJS7ux4YLk1IfwqPpf+jfvAqgNiYiznDDs2J37NjBjBkz+Pjjjzl9+jQpKSm5kZeIiOQRGVqZGAN/ziBt4TBKmBTOeZRmYaVXcS9em/5qYyLiVNdV2J0/f565c+cyffp0fv/9d5o2bcqYMWPo0KFDLqcnIiLO9u+RN4+0s/RNfBu2fYsbEF/ifsr3/oRORUpec38RuflyVNjFxcXx0Ucf8eWXXxIWFkb37t1Zs2YNkydPJjw8/GblKCIiThYdEUbg2S3U/60TuBwnxbiypkJ/Gj4+CiwWZ6cnIv8v24VdeHg4Fy5coFu3bqxZs8ZeyA0bNuymJSciInmAzQa/v0vnTaPAJY0DtlLE2AbwdY/+zs5MRK6QrT52AHv27KFRo0Y0bdqUqlWr3sycREQkrzh/Er7oAr+8BLY0frDWpaN1PH+mVWDykt3Ozk5ErpDtEbv4+HhmzZrFM888w8WLF+natSvdu3fHoiF4EZGCad9vMO9JOHuYNIsHL6dEEdS0L+uaV1IrE5E8Ktsjdrfddhsvvvgie/bs4dNPPyUhIYH77ruPtLQ0Zs2axa5du25mniIikssmxe7KfNTNZiVu5lBss9rC2cOc8i5H20uvENTsGaKbVwIuF3MxkZWYmNUxRMQprmtWbLNmzWjWrBmJiYnMnj2bGTNm8Oabb1KtWjU2bdqU2zmKiMhNkKGNCcDZBA5+9Bj1E9deXr67G597P00b9yIZRubUykQk77EYY3LlN3LDhg3MmDGDyZMn58bhCqWkpCT8/f1JTEzEz8/P2emISCGQfks1JrIS0aH7uTD3SXxST5Pi4o3Hg2/D3V2cnaJIoZeT+iDXCju5cSrsRMQZpsRux3XFGJ5x+x6A40XCKNXrCwjQs3MieUFO6oNsPWPXqlUrVq9efc24s2fPMmHCBN59993sZSoiIs515gD99kfbi7rZ1haUGrhKRZ1IPpWtZ+weffRROnXqhK+vL+3bt6dOnTqEhITg5eXF6dOn2bZtG6tWrWLhwoW0bduWN95442bnLSIiN2r7D/Dts3ApkSTjw4vWp/g+7V5O/npQM11F8qlsFXa9e/cmKiqKr7/+mrlz5/Lhhx9y5swZACwWC+Hh4bRs2ZJ169ZRuXLlm5mviIjcqNRLEDsC/pgGwAZbRTbUnch/H2hCmNqYiORr2Z4V6+HhQbdu3ejWrRsAiYmJXLx4kZIlS+Lu7n7TEhQRkVx0ci981RMSLncw+CDtAVIbv0S/yMtvE/r3e2H/vSwi+cN1tTsB8Pf3x9/fPzdzERGRXDIpdheuLhbHwmzTV/DDQEg5xwW3YvxcaRTJJe9XGxORAuS6CzsREcm7HHrUNQyBn56H9Z8B8LutKltrvUXvNvdlub9G6kTyJxV2IiIFUHph9sPiJXRf9wElL/yNwcI7aR1xbTKU/s2rODlDEbkZsv1KMRERyUeMIbrYbyz0fpmSF/7mqClOt5ThuDQdrqJOpADTiJ2ISEFzKenys3Rb5uEG/Gq7m0EpfTnrWpwvdItVpEDL8YhdhQoVOHnyZIb1Z86coUKFCrmSlIiIXKdDf8EHjWDLPHBxY1X5aHqk/IezrsVJsdqYvGS3szMUkZsoxyN2+/btw2q1ZlifnJzMoUOHciUpERHJXKazXQGM4ddPX+G++Mm4mjTwL8uX5Ubz/BrPy++BjQizvxcWNDlCpKDKdmH33Xff2f/8888/O7Q6sVqtLFmyhHLlyuVqciIi4shhtmt6cXbhFH9/1INGp369vFylLR8Ui2Hc8gR7UffveBV3IgVXtgu7Dh06AJffNNGjRw+Hbe7u7pQrV4633norV5MTERFHGYqzO05wdnYPKqQcxWpxx7X1OLjnSS4s3k1MpJ961IkUMhZjTI5+u8uXL8/atWsJCAi4WTkVWklJSfj7+5OYmIifn5+z0xGRPGzy4p1cWvYmMW5f42axcdq7LMUf/wyC73Z2aiKSy3JSH+T4Gbv4+PjrTkxERHLB2aNEH34e3JcD8K3tfh4cOAc8fZ2bl4g43XW1O1myZAlLlizh2LFj2Gw2h20zZszIlcRERCQTe5fBN0/B+WNcMJ68Yu3FnLSG7F+VQHSECjuRwi7Hhd3o0aN55ZVXqFOnDsHBwVgslpuRl4iI/Js1DZaPg5VvAYYdtjKsvWci49u1IESzXUXk/+W4j93UqVOZNWsWa9asYcGCBcyfP9/h62Z77733KF++PF5eXtSuXZuVK1deNX7FihXUrl0bLy8vKlSowNSpUzPEzJs3j/DwcDw9PQkPD8/0Oq51XmMMo0aNIiQkBG9vb5o0acLWrVtv7GJFpFCZFLsr8z5zif9w6J0IWPkmYJidFsGS+78gql0L4HIxFxNZiYlZ7S8ihUaOC7uUlBQaNGhwM3K5prlz5zJw4EBefPFF1q9fT8OGDWndujUHDhzIND4+Pp42bdrQsGFD1q9fz/Dhw4mOjmbevHn2mLi4ODp37kxUVBQbN24kKiqKTp06sWbNmhyd9/XXX2fixIlMmTKFtWvXEhQURGRkJGfPnr153xARKVDSW5k4FGc7f+LifxtwW9IGkl2L8GPlsZxsOoHnWlR32De9uNNsV5HCLcezYocOHUrRokUZMWLEzcopS3Xr1qVWrVq8//779nVVq1alQ4cOjBs3LkP80KFD+e6779i+fbt9Xd++fdm4cSNxcXEAdO7cmaSkJH766Sd7TKtWrShevDhffPFFts5rjCEkJISBAwcydOhQ4HLD5sDAQCZMmMDTTz+drevTrFgRSW8iPCSiPP2sn8Lv7wFwtGg4gU98DiXKOzlDEbnVbuqs2EuXLjFt2jQWL17MXXfdhbu7u8P2iRMn5vSQ2ZKSksK6desYNmyYw/oWLVqwevXqTPeJi4ujRYsWDutatmzJ9OnTSU1Nxd3dnbi4OAYNGpQh5u233872eePj40lISHA4l6enJ40bN2b16tVZFnbJyckkJyfbl5OSkq7yHRCRwiA6Igz/i/9QY2V3cPkbgPW3daNmr3fAzcPJ2YlIXpfjwm7Tpk3UqFEDgC1btjhsu5kTKU6cOIHVaiUwMNBhfWBgIAkJCZnuk5CQkGl8WloaJ06cIDg4OMuY9GNm57zp/80sZv/+/Vle07hx4xg9enSW20WkENoyjx6bBoDLWU6bogy1PsO0Pi85OysRySdyXNgtW7bsZuSRbVcWj8aYqxaUmcVfuT47x8ytmH974YUXiImJsS8nJSVRpkyZLONFpABLuQCLhsFfHwPwh60yQ6z9OWAtweQluzXbVUSy5br62AHs2bOHvXv30qhRI7y9va9ZxNyogIAAXF1dM4zOHTt2LMNIWbqgoKBM493c3ChZsuRVY9KPmZ3zBgUFAZdH7oKDg7OVG1y+Xevp6ZnldhEpJI7tgK96wvHtGCz8N+1BLE2G8WvzqvZn7kCtTETk2nI8K/bkyZNERERQqVIl2rRpw5EjRwB48sknGTx4cK4nmM7Dw4PatWsTGxvrsD42NjbLWbr169fPEP/LL79Qp04d+7OBWcWkHzM75y1fvjxBQUEOMSkpKaxYscJpM4hFJO/Iso2JMSye/SapUxvB8e2cdy9J95QXoOlL9G9eFVArExHJIZNDUVFRpmXLlubgwYOmaNGiZu/evcYYY37++WcTHh6e08PlyJw5c4y7u7uZPn262bZtmxk4cKApUqSI2bdvnzHGmGHDhpmoqCh7/N9//218fHzMoEGDzLZt28z06dONu7u7+frrr+0xv/32m3F1dTXjx48327dvN+PHjzdubm7m999/z/Z5jTFm/Pjxxt/f33zzzTdm8+bNpmvXriY4ONgkJSVl+/oSExMNYBITE2/k2yQiecw7i3eZ0KE/mHcW7/rfyktJZvuUzsaM9Lv89fGDZuoPqx1jrjjGxF923qKMRSQvyUl9kOPCLjAw0GzYsMEYYxwKu7///tsUKVIkp4fLsXfffdeEhoYaDw8PU6tWLbNixQr7th49epjGjRs7xC9fvtzUrFnTeHh4mHLlypn3338/wzG/+uorU7lyZePu7m6qVKli5s2bl6PzGmOMzWYzI0eONEFBQcbT09M0atTIbN68OUfXpsJOpOByKO4ObzCnxlczZqSfsY4sbsyvbxljtTo7RRHJo3JSH+S4j52vry9//fUXYWFh+Pr6snHjRipUqMDatWtp1aoVJ0+evBkDi4WC+tiJFGyTF+/ixLIpvOQ2Gw9LGmc9A/Ht/gmUrefs1EQkD8tJfZDjZ+waNWrEJ598Yl+2WCzYbDbeeOMNmjZtmvNsRUQKg4uniT4xmlfcP8bDksZiW218B/yuok5EclWOZ8W+8cYbNGnShD///JOUlBSef/55tm7dyqlTp/jtt99uRo4iIvnbwT/g6ycg8SApxpXXbd35KLUlMXEniY4o4ezsRKQAyXFhFx4ezqZNm3j//fdxdXXl/PnzPPTQQzz33HMOrT5ERAo9mw1WvwNLXgVjZZ8tkLhab/JSh/b4qY2JiNwEOSrsUlNTadGiBR988IHemCAiwuVWJq4ulozF2bnj7J8eRejpy++l/s5an8P3j6Nvy5rA/4o5FXcikptyVNi5u7uzZcuWm9qIWEQkP3F1sWQszv5ewfkvehGaepJUF0+Wlx/CvuAORDev5LBverzVlqM5bCIiWcrxrNjBgwfj7u7O+PHjb1ZOhZZmxYrkT+lvhxjcvAL9Xb7B/PoGFgwnfSpQssdsCAx3dooiko/lpD7I8TN2KSkpfPTRR8TGxlKnTh2KFCnisH3ixIk5PaSISL4WHRFGkeRjVPu1J7jswAJsDXyQO3tPBQ8fZ6cnIoVIjgu7LVu2UKtWLQB27drlsE23aEWkUNr1M7039wWXU5wzXrxsfZKJz4xxdlYiUgjlqLCzWq2MGjWK6tWrU6KEpuiLSCGXlgJLRkPcFAC22MoxyDqA3dZAyi3ZrQkRInLL5aiwc3V1pWXLlmzfvl2FnYgUbqf3Xe5Nd2gdADPTWnKh8UhiI++0P3MHmu0qIrdWjm/FVq9enb///pvy5cvfjHxERPKELNuYAAvnvE+z3a/hZT3HJTdfoi88SbWI7vZYtTIREWfJcWE3ZswYhgwZwquvvkrt2rUzTJ7QbE4RKQgybWOSepFN05+jTcK8y8tl6vJZ6Rep5hOSoXhTKxMRcYYctztxcfnf62X/PVnCGIPFYsFqteZedoWM2p2I5C3pt1RjIisRfZfh+KxulDq/+/LG+2Og6XBwdXdukiJS4N3UdifLli277sRERPKT9FG3fUs+4sKvMyllSeaCewl8On8Ed0Q4OTsRkYxyXNg1btz4ZuQhIpL3JJ8jOukt8PgCgDjbndSP/gZ8g5ycmIhI5nJc2P36669X3d6oUaPrTkZEJM9I2Axf9YKTu7EaC5Ntj/Lf1PYM/OMs0REq7EQkb8pxYdekSZMM6/79rJ2esRORfM0Y+HM6LBoO1mSOmBL8dvd4Bj3UGVe1MRGRPC7Hhd3p06cdllNTU1m/fj0jRoxgzBh1WheRvC/LViYXz7Dno17ccXIpAEusNdl73xs81eoeQG1MRCTvy3Fh5+/vn2FdZGQknp6eDBo0iHXr1uVKYiIiN0umrUz+WUfip49xR/JhrBY3VoX2Y+vt3YluXslhX7UxEZG8LMeFXVZKlSrFzp07c+twIiI3jcPIm7ER7fML1tiR+BsriV634f/YpzS+vTZZTRXTSJ2I5FU5Luw2bdrksGyM4ciRI4wfP56777471xITEbmZoiPC8Eo9wx0rngLXDbgCuwOaE/bkDPDKeGdCRCQ/yHFhV6NGDSwWC1f2Na5Xrx4zZszItcRERG6qfat4auuT4HqEZOPOa7YevPrcm/CvyWAiIvlNjgu7+Ph4h2UXFxdKlSqFl5dXriUlInLT2Kzw65uwYjwYG3tsIQyyDmCztQyllu7RbVYRyddyXNiFhobejDxERG6+pCPwTR/YtxKAr9IacbzRGL5vcZf99WGgZ+hEJP9yuXbIZUuXLiU8PJykpKQM2xITE7nzzjtZuXJlriYnIpJTk2J3MXnJ7owbdi/mwn/rw76VpLh4MyjlGY40ncizLe4CLhdzMZGVmJjV/iIi+UC2R+zefvtt+vTpk+nLZ/39/Xn66aeZOHEiDRs2zNUERURyIkMrE2sqLH0VfnsHH+B4kUr8WGkM5YuWzzAyp1YmIpLfWcyVsyCyEBoayqJFi6hatWqm23fs2EGLFi04cOBAriZYmCQlJeHv709iYmKmBbSIZE/6bdWRDYvS68ir8M9aADYFP8pdT0wBdz0TLCL5R07qg2yP2B09ehR3d/esD+TmxvHjx7OfpYjITRIdEUaFE8to+MdIsFwgyfiwMnwkD3Tu6+zURERuqmw/Y3fbbbexefPmLLdv2rSJ4ODgXElKROS6pV6Chf+h7fb/4G+5wHrbHTyYNl5FnYgUCtku7Nq0acPLL7/MpUuXMmy7ePEiI0eOpG3btrmanIhIjpzYA9Obwx/TAJia1pbu1lHEWwM0IUJECoVsP2N39OhRatWqhaurK/369aNy5cpYLBa2b9/Ou+++i9Vq5a+//iIwMPBm51xg6Rk7kRuw6Uv4YRCknOOiWzGeufAUtSI6ER0RZn/mLiayklqZiEi+c1OesQsMDGT16tU888wzvPDCC/Y3T1gsFlq2bMl7772nok5EbppJsbtwdbFkLMxSzrN1+tPcefR7AP7xr81DR5/gsci69liHd8OiPnUiUnDlqEFxaGgoCxcu5PTp0+zZswdjDGFhYRQvXvxm5SciAmTSxgTg6DZOzerGnRfjseGCS5OhfJ3yII/d5aZWJiJSKOX4zRMAxYsX55577sntXEREsuQw6mYM0cVXk/bj85SwJXPOI4CiXWdB+YYMzMYxREQKqusq7EREnCE6IgyPtHPctqI/uP6OG7CveAPKPfkpFAlwdnoiIk6nwk5E8o9Df9F3Ry9w3UeqcWWirQtD+08Bl2xP8BcRKdD0f0MRyfuMgbh3YXoLOL2Pf0wA3ayjeD/1ASYv2+vs7ERE8gyN2IlI3nbhFCx4Fnb9BMBP1ns4cN8EvmpV297GBPT8nIgIqLATkTwgy1Ym++M4O7sHvilHsVrcGZnyGKWbPkt080qA2piIiFxJhZ2IOF2GViY2K6yaiG3pOHyxctq7LD9WGktpvypqYyIichX55hm706dPExUVhb+/P/7+/kRFRXHmzJmr7mOMYdSoUYSEhODt7U2TJk3YunWrQ0xycjL9+/cnICCAIkWK0L59e/75558cn9tisWT4mjp1am5cukiBFx0RRkxkJSbG7uKjn+Lg046w9DVcsLKjdGuKD1zNYx3bZTkiFx0RxqDISrc4axGRvCffFHbdunVjw4YNLFq0iEWLFrFhwwaioqKuus/rr7/OxIkTmTJlCmvXriUoKIjIyEjOnj1rjxk4cCDz589nzpw5rFq1inPnztG2bVusVmuOzz1z5kyOHDli/+rRo0fufQNECrjoiDAm1TnFg793gvgVXDCexFYaSZVnvgBPX2enJyKSP5h8YNu2bQYwv//+u31dXFycAcyOHTsy3cdms5mgoCAzfvx4+7pLly4Zf39/M3XqVGOMMWfOnDHu7u5mzpw59phDhw4ZFxcXs2jRohydGzDz58+/oetMTEw0gElMTLyh44jkO2mpxsSOMmakvzEj/cy2EXeaVsM/dHZWIiJ5Qk7qg3wxYhcXF4e/vz9169a1r6tXrx7+/v6sXr06033i4+NJSEigRYsW9nWenp40btzYvs+6detITU11iAkJCaFatWr2mJycu1+/fgQEBHDPPfcwdepUbDbbjV+8SEF35iDMagOrJgKG2WkRdLKOYbs1mMlLdjs7OxGRfCVfTJ5ISEigdOnSGdaXLl2ahISELPcBCAwMdFgfGBjI/v377TEeHh4Z3nUbGBho3z+753711VeJiIjA29ubJUuWMHjwYE6cOMFLL72U5XUlJyeTnJxsX05KSsoyVqRA2rEQFjwDl86Q7FqEQRd7UyXicTZHhKmViYjIdXBqYTdq1ChGjx591Zi1a9cClycnXMkYk+n6f7tye3b2uTImO+f+dwFXo0YNAF555ZWrFnbjxo275vWL5GdZtjFJS2b9jAHUPPwFAEeLhvPIyT482ryhPVatTEREcs6phV2/fv3o0qXLVWPKlSvHpk2bOHr0aIZtx48fzzAily4oKAi4POIWHBxsX3/s2DH7PkFBQaSkpHD69GmHUbtjx47RoEEDe0xOzw2Xb9cmJSVx9OjRLONeeOEFYmJi7MtJSUmUKVMmy2OK5DcZ2pgAnNzLsZndqHlux+Xl+v2YQzcedfNQKxMRkRvk1MIuICCAgIBrv7i7fv36JCYm8scff3DvvfcCsGbNGhITE+0F2JXKly9PUFAQsbGx1KxZE4CUlBRWrFjBhAkTAKhduzbu7u7ExsbSqVMnAI4cOcKWLVt4/fXXr/vcAOvXr8fLy4tixYplGePp6Ymnp+c1r18kv8ow6lZ6IykLoiltPc9FN3+8H50GlVsxIBvHEBGRa8sXz9hVrVqVVq1a0adPHz744AMAnnrqKdq2bUvlypXtcVWqVGHcuHF07NgRi8XCwIEDGTt2LGFhYYSFhTF27Fh8fHzo1q0bAP7+/vTu3ZvBgwdTsmRJSpQowZAhQ6hevTrNmzfP9rm///57EhISqF+/Pt7e3ixbtowXX3yRp556SoWbFHrREWG4WS9RbNl/wG0ZHsAhv5rc1ns2+N/m7PRERAqUfFHYAcyePZvo6Gj7DNb27dszZcoUh5idO3eSmJhoX37++ee5ePEizz77LKdPn6Zu3br88ssv+Pr+ryfWpEmTcHNzo1OnTly8eJGIiAhmzZqFq6trts/t7u7Oe++9R0xMDDabjQoVKvDKK6/w3HPP3ZTvhUi+cmwHz+7uA27bsRkL79s68NyAj8A13/zvR0Qk37AYY/TwSh6RlJSEv78/iYmJ+Pn5OTsdkRtjDKz/DBb+B9Iuctz4M8TajxVpdxITWUm3WEVEsikn9YH+ySwiuS/5LPwwCDZ/BcCv1ursbPAmH7eupzYmIiI3kQo7EbkuWbYyObyBM588RrFLB7Hhyhupj+LdNIbo5pefSVUbExGRm0eFnYhclwytTIyBNR9g/fkliplUznoGsrDSGLyL1VAbExGRW0SFnYhcl3+PvHmmJvL0mUmw4wdcgb0lG1Ox9yw6+5S45v4iIpJ7VNiJyHWLjggjKGkjDeI6g+UEycaNuIoDaRL1ElzjDS8iIpL7VNiJyPWx2eC3t+m06TWwWIm3BTLIOpAFjz/r7MxERAotF2cnICL50LljMPthWDIajJVvrQ14yDqODdZQJi/Z7ezsREQKLY3YiUjO/L0CvukD546S6uLJi8mPc3vTp1jfvJJamYiIOJkKOxFxkGUbE2saa2YN5d6D07FgOOlTgS6n+9KueTN7rFqZiIg4lwo7EXGQoY0JQOIhDk3vTt2k9ZeXaz3O5x59aOfurVYmIiJ5iAo7EXGQYdStzN9c/LIPt6UlkuLqg0eH/0L1R+ifjWOIiMitpcJORDKIjgjDxZaK57KXwW0h3sCxolUo3etzKFnR2emJiEgWVNiJSEan4ukX/xy4/QXAx9ZW9Bj4Cbh5OjkxERG5GhV2IuJo63z4LhqSkzhjijDc2peFabVJXHFAt1hFRPI4FXYiclnqRVj0AqybCcCftkpsqjuR9x5oqDYmIiL5hAo7kUIky1Ymx3dxYmZXAi7swWDhvbR2mMbD6RdZFVAbExGR/EKFnUghkmkrkw2fk/rdIAJsl7jgXoKfwkZjLVlfbUxERPIhFXYihci/R97crRd45tx7sGkO7sCBYvdStvdnPOwbeM39RUQkb1JhJ1LIREeEEXBuJ3VXdQaXI1iNhTXl+tKgxxhwcXV2eiIicgNU2IkUJsbA2o/otulFcEnmsCnBEGt/Pu8V4+zMREQkF7g4OwERuUUunoEvH4eFQ8CazGJrTTqkTWB1WmUmL9nt7OxERCQXaMROpDD450/4uhecOYDV4saYlK4UaxrNH80rqZWJiEgBosJOpADIso2JzcbKT0fRYN8UXI2VRK/beCzxGSKbt7LHqpWJiEjBocJOpADItI3J+RPs+yiKhqdXX16+syOfFY0m0sNXrUxERAooFXYiBUCGUbeKRzn3eU/KpRwnzcUTtzYToHZPnrNYrnkMERHJv1TYiRQQ0RFhWIyVtKXjsP76DUUthlPe5SjR83MIvNPZ6YmIyC2gwk6koEg6Qv9/hoD7SgC+tjXhkUGfg0cRJycmIiK3itqdiBQEuxfD1Ptg30rOG0+GpD3HkJSnmLzysLMzExGRW0gjdiL5mTUVlrwCqycDsNUWyrp7JvJmu+aUVRsTEZFCR4WdSB6XZSuT0/s5MqM7wWc3A/BxWiTnGo/iuchqgNqYiIgURirsRPK4TFuZbPuOS988S3DaWZJdi/JL2AgSS0WojYmISCGnwk4kj/v3yJurLZnnUmbB2g/xAo74Vif4idm0Kx56zf1FRKTgU2Enkg9ER4RR7MJ+aq/sBi77Afjz9sep02siuLo7OTsREckrNCtWJD/YOIfHNz3OnS77OWH86J02jDpP/ldFnYiIOFBhJ5KXpZyHBc/C/Kch9Txx1nA6pI1nSdpdTF6y29nZiYhIHqNbsSJ5VcIW+LoXnNiFDRfeTn0I96b/YVXzKkxWKxMREcmECjsRJ8myjYkxLJ09gUZ/T8TNlsw5j1L0PtuX+5o/SP//j1UrExERyYwKOxEnybSNyaVEdn3Um2YnYi8v3xHJZyX+w32eJdTKRERErkmFnYiTZBh1q5JE4qdRVLp0CKvFFdfmo6B+P/q6ZP0orEbqRETk31TYiThRdEQYGEPi0rdJ+fUL/C1WEj1D8I/6FG6v4+z0REQkn1FhJ+JMF04RfWwEuC8CYJHtXloN/Aq8izk3LxERyZfyTbuT06dPExUVhb+/P/7+/kRFRXHmzJmr7mOMYdSoUYSEhODt7U2TJk3YunWrQ0xycjL9+/cnICCAIkWK0L59e/755x+HmDFjxtCgQQN8fHwoVqxYpuc6cOAA7dq1o0iRIgQEBBAdHU1KSsqNXLIUdPtXw9T7Ydciko07o6xP0DdlAJNXH3d2ZiIikk/lm8KuW7dubNiwgUWLFrFo0SI2bNhAVFTUVfd5/fXXmThxIlOmTGHt2rUEBQURGRnJ2bNn7TEDBw5k/vz5zJkzh1WrVnHu3Dnatm2L1Wq1x6SkpPDoo4/yzDPPZHoeq9XKAw88wPnz51m1ahVz5sxh3rx5DB48OHcuXgoWmxVWvAGzHoCkQ+y1BTOv1ixGvTqJmMjKTIzdpR51IiJyfUw+sG3bNgOY33//3b4uLi7OAGbHjh2Z7mOz2UxQUJAZP368fd2lS5eMv7+/mTp1qjHGmDNnzhh3d3czZ84ce8yhQ4eMi4uLWbRoUYZjzpw50/j7+2dYv3DhQuPi4mIOHTpkX/fFF18YT09Pk5iYmO3rTExMNECO9pG8aeIvO807i3dl3JB0xOyf1NyYkX7GjPQzX7/Uxrz38waHkHcW7zKhQ3/IfH8RESl0clIf5IsRu7i4OPz9/albt659Xb169fD392f16tWZ7hMfH09CQgItWrSwr/P09KRx48b2fdatW0dqaqpDTEhICNWqVcvyuFnlV61aNUJCQuzrWrZsSXJyMuvWrctyv+TkZJKSkhy+pGBIb2XiMPK2ZwkXJten7Jk/SHXx4uewkRxq8jbPtLjbYd/oiDBiIiupjYmIiORYvpg8kZCQQOnSpTOsL126NAkJCVnuAxAYGOiwPjAwkP3799tjPDw8KF68eIaYrI6b1bmuPE/x4sXx8PC46nHGjRvH6NGjs30eyT/+3crExZZGP8tcWDUJH+B4kTBK9fyclqUq0fIa+4uIiOSEU0fsRo0ahcViuerXn3/+CYDFYsmwvzEm0/X/duX27OyTnZhrnSc7x3nhhRdITEy0fx08eDBH55S8LToijJcb+lFv5eOwahIAm4IeptTAlVCqkpOzExGRgsipI3b9+vWjS5cuV40pV64cmzZt4ujRoxm2HT9+PMNIWbqgoCDg8mhacHCwff2xY8fs+wQFBZGSksLp06cdRu2OHTtGgwYNsn0dQUFBrFmzxmHd6dOnSU1NzTI/uHxr2NPTM9vnkXxmx488sflZcDlDkvHmJevTTO6rEVoREbl5nDpiFxAQQJUqVa765eXlRf369UlMTOSPP/6w77tmzRoSExOzLMDKly9PUFAQsbGx9nUpKSmsWLHCvk/t2rVxd3d3iDly5AhbtmzJUWFXv359tmzZwpEjR+zrfvnlFzw9Paldu3a2jyMFRFoy/DQU5nSDS2fYYKtAx7TxfJd2r2a7iojITZUvnrGrWrUqrVq1ok+fPnzwwQcAPPXUU7Rt25bKlSvb46pUqcK4cePo2LEjFouFgQMHMnbsWMLCwggLC2Ps2LH4+PjQrVs3APz9/enduzeDBw+mZMmSlChRgiFDhlC9enWaN29uP+6BAwc4deoUBw4cwGq1smHDBgDuuOMOihYtSosWLQgPDycqKoo33niDU6dOMWTIEPr06YOfn9+t+0bJTTcpdheuLpZMn4GbvGQ3Rc8f4Ikjo+HIRgCmpT1ASuOXWBIZzuQluzO+G1ZERCQX5YvCDmD27NlER0fbZ7C2b9+eKVOmOMTs3LmTxMRE+/Lzzz/PxYsXefbZZzl9+jR169bll19+wdfX1x4zadIk3Nzc6NSpExcvXiQiIoJZs2bh6upqj3n55Zf5+OOP7cs1a9YEYNmyZTRp0gRXV1d+/PFHnn32We677z68vb3p1q0bb7755k35XojzpM92BcfibPKS3exZMos3vWeC9TwX3fx59sJT1IzobI/L8G5YFXciIpLLLMYY9VTII5KSkvD39ycxMVEjfXlY+shbTGQloiPCeO+XzRT/dQRd3ZZdDijbgA9LD+eid1CWI3tWm2FQpCZQiIjIteWkPlBhl4eosMs/0ou7cNdDTHJ9h8ou/wAWaPQfaDwUXPPNYLiIiORxOakP9LePyHWIbnYHCcumMcJ1Ft6WFCgaCA9NgwpNnJ2aiIgUYvnizRMiecqlJHa+15mxbtPwtqTwq7U6H975iYo6ERFxOo3YieTE4fWc+SSKypcOYsMVl4iX2JD8ABMX7+Gix25NiBAREadSYSeSHcbAmg+w/vwixUwaSZ5B+HX/BMrWJRrA4qLZriIi4nQq7ESu5cIp+LYf7PwRV2BvySZU7D0TfErYQ9KLOatNc5FERMR5VNiJ/L9Mmw8fWANfPwFJ/5Bmccet1Rgq3vsUZPIOYI3UiYiIs2nyhMj/S28+PHnJbrDZYOVbMLM1JP1DvC2Qr2rMhLpPZ1rUiYiI5AUasRP5f+kjbp/EruXBzf0IPfM7AAusDThy/zieaVnDidmJiIhcmwo7kX+JLn+I3r4vUeTMSS4aD15O60mZpn2Ibq63RIiISN6nwk4EwJoGK8bDr29SBMMuczvPpUSz36Usu1TUiYhIPqFn7KRQmJT+7FwmZixcxT/vNIdf3wAMW4I60j75Vfa7lCXFastyPxERkbxGI3ZSKKRPjADH2avffTmDDltHUMJyDjx8+an8MJ7ZWJ6YyEpER4TZ3wl75X4iIiJ5kQo7KRTSizJ7kdY4lL9mDqT9odlgAYLv5uPbRjFy1UV7UZfpfiruREQkD1NhJ4VGelH21eKVNFnRmVouey9vqPsMRI7m1LL9xERaMhRvaj4sIiL5hcUYo7+t8oikpCT8/f1JTEzEz8/P2ekUTFvnc/bLZ/C1XOSMKUKxrh9ClQecnZWIiEiWclIfaPKEFA6pF+H7gfBVT3wtF1lnq0Sb5HFMPqQZryIiUnDoVqwUfMd3wle94NhWbMbCurI9uafnG3RZvk/PzomISIGiwk4KLmNgw+ewcAikXuC48eO36mPp8EgUoIkRIiJS8Kiwk4Ip+Sz8OBg2zQXggP89/FLlVZ5sXd8hTBMjRESkIFFhJwXPkU3wVU84tRcsLtB0OGXvj+FJF9dMwzVSJyIiBYUmT0i+5vBGCWPgjw/howg4tZezHoF8eedUaPQfyKKoExERKUg0Yif5WvobJTxTk3g6cRJs/x6A+BIN6Xj4MZ4oXsvJGYqIiNw6KuwkX4uOCCMoaTMN4jqB5QS4uPNruf48vq02MZGVdZtVREQKFRV2kn/ZbLB6Mp02vwqWNPbbShOTNoB128o7vBZMRESksFBhJ3nWpNhduLpkfMUXwIc//UHznS9T/kzc5RV3dqTjhnacsnrh4eqiok5ERAolTZ6QPCv9+Tn75Ij/N2/eF7T//dHLRZ2bF7R9m8nFh9uLuhSrLcM+IiIihYFG7CTPytBAuGkF1swaSof9H+FqMRBQCR6dxeQtHkxcvMt++3Xykt1qOiwiIoWSCjvJ09ILs9mxv1N3RRR1XbaDBajxGLR5nckrDzMxdpfDM3V6o4SIiBRWKuwkz4suu4/HPF+ghOUs540nRR76L9zdGbj8xojMJkrojRIiIlIYqbCTvMuaCktegdWTKWGBbaYcz6X0p+OJWkT/f8igyEpZ7q6ROhERKWxU2Mktd7XZrpOX7MZqMwyq4wlfPwGH/gRgY/Cj3P3EFDr+elC3WEVERLKgwk5uufTZruBYnKVPeni35kGYOhaSE0kyPqwMH8UDnZ92iFdxJyIikpEKO7nlMivOJi/ZzbuxW/iu/A/ctf1rAA77VufnKmPo9UDjTPfX83MiIiKOLMYY/e2YRyQlJeHv709iYiJ+fn7OTuemSx+h83B14XbbP8wtMY1S5y8Xe9w3AJqNAFd35yYpIiLiZDmpD9SgWJwmOiIMD1cXHjAr+N7jxctFnU8AdJ8Hka+oqBMREckh3YoVp3n/l42MdXmPR1x/BeAf/9rc3ns2+AU7OTMREZH8SSN24hSzv1tI5Koul4s6iwu/l32aRkcHMXntOWenJiIikm9pxE5uLWNYOns8j+x+C0+XVPANhoc/ol65+xmoV4GJiIjcEBV2cutcSoTvomm2Z8Hl14LdEQkdp0KRAECzXUVERG5UvrkVe/r0aaKiovD398ff35+oqCjOnDlz1X2MMYwaNYqQkBC8vb1p0qQJW7dudYhJTk6mf//+BAQEUKRIEdq3b88///zjEDNmzBgaNGiAj48PxYoVy/RcFoslw9fUqVNv5JLzpUmxu5i8ZHfGDYfWkfh2fdi2AFzcIPJV6PalvahLFx0RdtW3SYiIiEjW8k1h161bNzZs2MCiRYtYtGgRGzZsICoq6qr7vP7660ycOJEpU6awdu1agoKCiIyM5OzZs/aYgQMHMn/+fObMmcOqVas4d+4cbdu2xWq12mNSUlJ49NFHeeaZZ656vpkzZ3LkyBH7V48ePW7sovOh9ObD9uLOGFg9BetHLfC/dIhEzxB44me4Lxpc8s2Pn4iISP5g8oFt27YZwPz+++/2dXFxcQYwO3bsyHQfm81mgoKCzPjx4+3rLl26ZPz9/c3UqVONMcacOXPGuLu7mzlz5thjDh06ZFxcXMyiRYsyHHPmzJnG398/0/MBZv78+ddxdf+TmJhoAJOYmHhDx3G2dxbvMqFDfzAf/PSHMZ89asxIP2NG+pndkzsYc+G0s9MTERHJV3JSH+SLIZO4uDj8/f2pW7eufV29evXw9/dn9erVme4THx9PQkICLVq0sK/z9PSkcePG9n3WrVtHamqqQ0xISAjVqlXL8rhX069fPwICArjnnnuYOnUqNpvtqvHJyckkJSU5fBUE0RFhvHnvedrFdYLdP5Ns3FlWcSh39PsGvIs5Oz0REZECK19MnkhISKB06dIZ1pcuXZqEhIQs9wEIDAx0WB8YGMj+/fvtMR4eHhQvXjxDTFbHzcqrr75KREQE3t7eLFmyhMGDB3PixAleeumlLPcZN24co0ePztF58jybFVZO5JHNY8FiY68tmBjrAL6NuvptbBEREblxTh2xGzVqVKaTDv799eeffwKXJydcyRiT6fp/u3J7dvbJTsyVXnrpJerXr0+NGjUYPHgwr7zyCm+88cZV93nhhRdITEy0fx08eDBH58xzzibApx1g2WtgbMyz3s/D1nFstJbNfEKFiIiI5Cqnjtj169ePLl26XDWmXLlybNq0iaNHj2bYdvz48QwjcumCgoKAy6NywcH/e5PBsWPH7PsEBQWRkpLC6dOnHUbtjh07RoMGDXJ8Pf9Wr149kpKSOHr0aJY5enp64unpeUPnyTP2LIZvnoYLJ0h18WLYpZ6ERjzJhogw+zthQf3pREREbianFnYBAQEEBARcM65+/fokJibyxx9/cO+99wKwZs0aEhMTsyzAypcvT1BQELGxsdSsWRO4PLt1xYoVTJgwAYDatWvj7u5ObGwsnTp1AuDIkSNs2bKF119//Yaubf369Xh5eWXZHqXAsKbCsjGwahIAx4uE0eXUUzzYvKm9iEv/r4o7ERGRmytfPGNXtWpVWrVqRZ8+ffjggw8AeOqpp2jbti2VK1e2x1WpUoVx48bRsWNHLBYLAwcOZOzYsYSFhREWFsbYsWPx8fGhW7duAPj7+9O7d28GDx5MyZIlKVGiBEOGDKF69eo0b97cftwDBw5w6tQpDhw4gNVqZcOGDQDccccdFC1alO+//56EhATq16+Pt7c3y5Yt48UXX+Spp54qMCNyk2J34epicSzKzhyAr3vDP39cXq7zBF+49eJBd+8MxZuaD4uIiNwCN3uKbm45efKk6d69u/H19TW+vr6me/fu5vTp0w4xgJk5c6Z92WazmZEjR5qgoCDj6elpGjVqZDZv3uywz8WLF02/fv1MiRIljLe3t2nbtq05cOCAQ0yPHj0MkOFr2bJlxhhjfvrpJ1OjRg1TtGhR4+PjY6pVq2befvttk5qamqNrzMvtTtJbmLyzeNflFdu+N2ZcGWNG+pnElwPNj1+859wERURECqic1AcWY4yGUPKIpKQk/P39SUxMxM/Pz9npZDB5yW6mxG7lywoLqXF4LgAbbBXYcO9EerZt6uTsRERECqac1Af54las5A3RNVzo8tdYSh/eCcC0tAdIafwS/SLDnZyZiIiIgAo7ya5NX8EPAymdco7TxpeY1L78ZqnNLhV1IiIieUa+ePOE3FyT/v1u1yu898smtrwXBd88CSnnOORXi9bJY/nNUpsUq0396URERPIQjdgJri6WTFuRfPbdTzT/czCVXA4BFtaUfZJuuxoxILIq0epPJyIikueosJOMfeaa3cGS2W/wyO7X8XJJhaKBfFN+FDFr/YmJrKT+dCIiInmUCjsB/leUTYvdQIUV0bR1WQ0WoGIz6DiN/XGniYm0qD+diIhIHqZ2J3mI09udHF7P/g86E2o5SppxwS3yZWgwAFz0KKaIiIiz5KQ+0N/YAsbA7+9j/bA5oZajHDIBdEp5mcnJbVXUiYiI5CO6FVvYXTgF3z4HOxfiCuwt2YSKvWfSJO6knp0TERHJZ1TYFWYHfr/8rtekf0g2bsRVHEiTqJfAYiE6ogSgiREiIiL5iQq7wshmg98mwdIxYKyc9irDoqrj6PpgO4cwTYwQERHJX1TYFWCTYnfh6nLFTNZzx+Cbp+DvZZeXqz9K8baT6Orpm+kxNFInIiKSf+jJ+AIsvfGw/e0Qe5fB+/fB38u4aDxYHDYCHvoQsijqREREJH/RiF0Blj7a9k7sdurte597D84ADDtsZfizzls81r6lcxMUERGRXKXCroCLjgjDYtIwv44BF8PnaU050+hVnm1R3dmpiYiISC5TYVcI9G9elYbLorkrbQexlgbsUlEnIiJSIOkZu0Jg8pLdHLQWI9bSgBSr7X/P3ImIiEiBohG7Am7ykt1MjN1FTGQloiPC7MugGa8iIiIFjQq7AuzKog7+V8ypuBMRESl4VNgVYFabcSjq0qnxsIiISMFkMcbob/c8IikpCX9/fxITE/Hz83N2OiIiIpIH5KQ+0OQJERERkQJChZ2IiIhIAaHCTkRERKSAUGEnIiIiUkCosBMREREpIFTYiYiIiBQQKuxERERECggVdiIiIiIFhAo7ERERkQJChZ2IiIhIAaHCTkRERKSAcHN2AvI/6a/tTUpKcnImIiIiklek1wXpdcLVqLDLQ86ePQtAmTJlnJyJiIiI5DVnz57F39//qjEWk53yT24Jm83G4cOH8fX1xWKx5Oqxk5KSKFOmDAcPHsTPzy9Xj53XFKZrhcJ1vbrWgknXWjDpWnOPMYazZ88SEhKCi8vVn6LTiF0e4uLiwu23335Tz+Hn51fgf8HSFaZrhcJ1vbrWgknXWjDpWnPHtUbq0mnyhIiIiEgBocJOREREpIBQYVdIeHp6MnLkSDw9PZ2dyk1XmK4VCtf16loLJl1rwaRrdQ5NnhAREREpIDRiJyIiIlJAqLATERERKSBU2ImIiIgUECrsCogxY8bQoEEDfHx8KFasWLb2McYwatQoQkJC8Pb2pkmTJmzdutUhJjk5mf79+xMQEECRIkVo3749//zzz024gpw5ffo0UVFR+Pv74+/vT1RUFGfOnLnqPhaLJdOvN954wx7TpEmTDNu7dOlyk6/m6q7nWnv27JnhOurVq+cQkxc/25xea2pqKkOHDqV69eoUKVKEkJAQHn/8cQ4fPuwQlxc+1/fee4/y5cvj5eVF7dq1Wbly5VXjV6xYQe3atfHy8qJChQpMnTo1Q8y8efMIDw/H09OT8PBw5s+ff7PSz5GcXOs333xDZGQkpUqVws/Pj/r16/Pzzz87xMyaNSvT391Lly7d7Eu5ppxc6/LlyzO9jh07djjE5dXPFXJ2vZn9f8hisXDnnXfaY/LiZ/vrr7/Srl07QkJCsFgsLFiw4Jr75KnfVyMFwssvv2wmTpxoYmJijL+/f7b2GT9+vPH19TXz5s0zmzdvNp07dzbBwcEmKSnJHtO3b19z2223mdjYWPPXX3+Zpk2bmrvvvtukpaXdpCvJnlatWplq1aqZ1atXm9WrV5tq1aqZtm3bXnWfI0eOOHzNmDHDWCwWs3fvXntM48aNTZ8+fRzizpw5c7Mv56qu51p79OhhWrVq5XAdJ0+edIjJi59tTq/1zJkzpnnz5mbu3Llmx44dJi4uztStW9fUrl3bIc7Zn+ucOXOMu7u7+fDDD822bdvMgAEDTJEiRcz+/fszjf/777+Nj4+PGTBggNm2bZv58MMPjbu7u/n666/tMatXrzaurq5m7NixZvv27Wbs2LHGzc3N/P7777fqsjKV02sdMGCAmTBhgvnjjz/Mrl27zAsvvGDc3d3NX3/9ZY+ZOXOm8fPzy/A77Gw5vdZly5YZwOzcudPhOv79O5dXP1djcn69Z86ccbjOgwcPmhIlSpiRI0faY/LiZ7tw4ULz4osvmnnz5hnAzJ8//6rxee33VYVdATNz5sxsFXY2m80EBQWZ8ePH29ddunTJ+Pv7m6lTpxpjLv9Suru7mzlz5thjDh06ZFxcXMyiRYtyPffs2rZtmwEcfiHi4uIMYHbs2JHt4zz44IOmWbNmDusaN25sBgwYkFup3rDrvdYePXqYBx98MMvtefGzza3P9Y8//jCAw182zv5c7733XtO3b1+HdVWqVDHDhg3LNP755583VapUcVj39NNPm3r16tmXO3XqZFq1auUQ07JlS9OlS5dcyvr65PRaMxMeHm5Gjx5tX87u/9dutZxea3phd/r06SyPmVc/V2Nu/LOdP3++sVgsZt++ffZ1efWzTZedwi6v/b7qVmwhFR8fT0JCAi1atLCv8/T0pHHjxqxevRqAdevWkZqa6hATEhJCtWrV7DHOEBcXh7+/P3Xr1rWvq1evHv7+/tnO6+jRo/z444/07t07w7bZs2cTEBDAnXfeyZAhQzh79myu5Z5TN3Kty5cvp3Tp0lSqVIk+ffpw7Ngx+7a8+NnmxucKkJiYiMViyfBIgrM+15SUFNatW+fwvQZo0aJFltcVFxeXIb5ly5b8+eefpKamXjXGmb+b13OtV7LZbJw9e5YSJUo4rD937hyhoaHcfvvttG3blvXr1+da3tfjRq61Zs2aBAcHExERwbJlyxy25cXPFXLns50+fTrNmzcnNDTUYX1e+2xzKq/9vupdsYVUQkICAIGBgQ7rAwMD2b9/vz3Gw8OD4sWLZ4hJ398ZEhISKF26dIb1pUuXznZeH3/8Mb6+vjz00EMO67t370758uUJCgpiy5YtvPDCC2zcuJHY2NhcyT2nrvdaW7duzaOPPkpoaCjx8fGMGDGCZs2asW7dOjw9PfPkZ5sbn+ulS5cYNmwY3bp1c3hfozM/1xMnTmC1WjP9XcvquhISEjKNT0tL48SJEwQHB2cZ48zfzeu51iu99dZbnD9/nk6dOtnXValShVmzZlG9enWSkpJ45513uO+++9i4cSNhYWG5eg3ZdT3XGhwczLRp06hduzbJycl8+umnREREsHz5cho1agRk/dk783OFG/9sjxw5wk8//cTnn3/usD4vfrY5ldd+X1XY5WGjRo1i9OjRV41Zu3YtderUue5zWCwWh2VjTIZ1V8pOzPXI7vVCxrxzmteMGTPo3r07Xl5eDuv79Olj/3O1atUICwujTp06/PXXX9SqVStbx86Om32tnTt3tv+5WrVq1KlTh9DQUH788ccMxWxOjns9btXnmpqaSpcuXbDZbLz33nsO227V53o1Of1dyyz+yvXX8/t7K1xvXl988QWjRo3i22+/dSjy69Wr5zD557777qNWrVr897//ZfLkybmX+HXIybVWrlyZypUr25fr16/PwYMHefPNN+2FXU6Peatdb26zZs2iWLFidOjQwWF9Xv5scyIv/b6qsMvD+vXrd82Ze+XKlbuuYwcFBQGX/6URHBxsX3/s2DH7vyqCgoJISUnh9OnTDiM7x44do0GDBtd13qvJ7vVu2rSJo0ePZth2/PjxDP8iyszKlSvZuXMnc+fOvWZsrVq1cHd3Z/fu3blaANyqa00XHBxMaGgou3fvBm7tZ3srrjU1NZVOnToRHx/P0qVLHUbrMnOzPtfMBAQE4OrqmuFf5v/+XbtSUFBQpvFubm6ULFnyqjE5+bnIbddzrenmzp1L7969+eqrr2jevPlVY11cXLjnnnvsP8/OcCPX+m/16tXjs88+sy/nxc8Vbux6jTHMmDGDqKgoPDw8rhqbFz7bnMpzv6+5/tSeOFVOJ09MmDDBvi45OTnTyRNz5861xxw+fDjPTJ5Ys2aNfd3vv/+e7Yfse/TokWHWZFY2b95sALNixYrrzvdG3Oi1pjtx4oTx9PQ0H3/8sTEmb36213utKSkppkOHDubOO+80x44dy9a5bvXneu+995pnnnnGYV3VqlWvOnmiatWqDuv69u2b4WHs1q1bO8S0atXK6Q/Z5/RajTHm888/N15eXtd8SD2dzWYzderUMb169bqRVG/Y9VzrlR5++GHTtGlT+3Je/VyNuf7rTZ80snnz5mueI698tunI5uSJvPT7qsKugNi/f79Zv369GT16tClatKhZv369Wb9+vTl79qw9pnLlyuabb76xL48fP974+/ubb775xmzevNl07do103Ynt99+u1m8eLH566+/TLNmzZzeEsOYy78Qd911l4mLizNxcXGmevXqGdpiXHm9xhiTmJhofHx8zPvvv5/hmHv27DGjR482a9euNfHx8ebHH380VapUMTVr1nR6C5CcXOvZs2fN4MGDzerVq018fLxZtmyZqV+/vrntttvy/Geb02tNTU017du3N7fffrvZsGGDQ7uE5ORkY0ze+FzT20RMnz7dbNu2zQwcONAUKVLEPjtw2LBhJioqyh6f3j5h0KBBZtu2bWb69OkZ2if89ttvxtXV1YwfP95s377djB8/Pk+0xcjptX7++efGzc3NvPvuu1m2oxk1apRZtGiR2bt3r1m/fr3p1auXcXNzc/hHgDPk9FonTZpk5s+fb3bt2mW2bNlihg0bZgAzb948e0xe/VyNyfn1pnvsscdM3bp1Mz1mXvxsz549a/87FDATJ04069evt8+0z+u/ryrsCogePXoYIMPXsmXL7DGAmTlzpn3ZZrOZkSNHmqCgIOPp6WkaNWqU4V9UFy9eNP369TMlSpQw3t7epm3btubAgQO36KqydvLkSdO9e3fj6+trfH19Tffu3TO0ELjyeo0x5oMPPjDe3t6Z9jA7cOCAadSokSlRooTx8PAwFStWNNHR0Rn6v91qOb3WCxcumBYtWphSpUoZd3d3U7ZsWdOjR48Mn1te/Gxzeq3x8fGZ/tz/+2c/r3yu7777rgkNDTUeHh6mVq1aDqOFPXr0MI0bN3aIX758ualZs6bx8PAw5cqVy/QfI1999ZWpXLmycXd3N1WqVHEoEJwpJ9fauHHjTD+/Hj162GMGDhxoypYtazw8PEypUqVMixYtzOrVq2/hFWUtJ9c6YcIEU7FiRePl5WWKFy9u7r//fvPjjz9mOGZe/VyNyfnP8ZkzZ4y3t7eZNm1apsfLi59t+ghjVj+Tef331WLM/z/hJyIiIiL5mvrYiYiIiBQQKuxERERECggVdiIiIiIFhAo7ERERkQJChZ2IiIhIAaHCTkRERKSAUGEnIiIiUkCosBMREREpIFTYiYjkgnLlyvH22287O42rmj59Oi1atMjRPseOHaNUqVIcOnToJmUlIrlJhZ2IFGrt2rWjefPmmW6Li4vDYrHw119/5fi4FouFBQsW3GB2uSc5OZmXX36ZESNGANC/f3/CwsIyjT106BCurq588803lC5dmqioKEaOHHkr0xWR66TCTkQKtd69e7N06VL279+fYduMGTOoUaMGtWrVckJmuWvevHkULVqUhg0bApeve8+ePaxcuTJD7KxZsyhZsiTt2rUDoFevXsyePZvTp0/f0pxFJOdU2IlIoda2bVtKly7NrFmzHNZfuHCBuXPn0rt3b+ByYXTnnXfi6elJuXLleOutt7I8Zrly5QDo2LEjFovFvrx3714efPBBAgMDKVq0KPfccw+LFy922PfIkSM88MADeHt7U758eT7//PMMt3kTExN56qmnKF26NH5+fjRr1oyNGzde9TrnzJlD+/bt7cvpBeuMGTMyxM6aNYvHH38cd3d3AKpXr05QUBDz58+/6jlExPlU2IlIoebm5sbjjz/OrFmzMMbY13/11VekpKTQvXt31q1bR6dOnejSpQubN29m1KhRjBgxIkMxmG7t2rUAzJw5kyNHjtiXz507R5s2bVi8eDHr16+nZcuWtGvXjgMHDtj3ffzxxzl8+DDLly9n3rx5TJs2jWPHjtm3G2N44IEHSEhIYOHChaxbt45atWoRERHBqVOnsrzOlStXUqdOHYd1vXv35quvvuLcuXP2dStWrGDPnj088cQTDrH33ntvpqN7IpLHGBGRQm779u0GMEuXLrWva9Sokenatasxxphu3bqZyMhIh33+85//mPDwcPtyaGiomTRpkn0ZMPPnz7/mucPDw81///tfhzzWrl1r3757924D2I+9ZMkS4+fnZy5duuRwnIoVK5oPPvgg03OcPn3aAObXX3/NsN7Ly8vMmDHDvu7xxx839evXz3CMQYMGmSZNmlzzekTEuTRiJyKFXpUqVWjQoIH9tuTevXtZuXKlfdRq+/bt3HfffQ773HfffezevRur1Zrt85w/f57nn3+e8PBwihUrRtGiRdmxY4d9xG7nzp24ubk5PNN3xx13ULx4cfvyunXrOHfuHCVLlqRo0aL2r/j4ePbu3ZvpeS9evAiAl5eXw/pixYrx0EMP2a/77NmzzJs3L8NoHYC3tzcXLlzI9rWKiHO4OTsBEZG8oHfv3vTr1493332XmTNnEhoaSkREBHD59qfFYnGIN/+6bZtd//nPf/j555958803ueOOO/D29uaRRx4hJSXlqsf893qbzUZwcDDLly/PEFesWLFM9y9ZsiQWiyXTyQ+9e/cmIiKC3bt3s2LFCgA6d+6cIe7UqVOUKlXqWpcoIk6mETsREaBTp064urry+eef8/HHH9OrVy97MRceHs6qVasc4levXk2lSpVwdXXN9Hju7u4ZRvNWrlxJz5496dixo31Cwr59++zbq1SpQlpaGuvXr7ev27NnD2fOnLEv16pVi4SEBNzc3LjjjjscvgICAjLNxcPDg/DwcLZt25ZhW9OmTalQoQKzZs1ixowZdOrUCV9f3wxxW7ZsoWbNmpkeX0TyDhV2IiJA0aJF6dy5M8OHD+fw4cP07NnTvm3w4MEsWbKEV199lV27dvHxxx8zZcoUhgwZkuXxypUrx5IlS0hISLCPlN1xxx188803bNiwgY0bN9KtWzdsNpt9nypVqtC8eXOeeuop/vjjD9avX89TTz2Ft7e3vchs3rw59evXp0OHDvz888/s27eP1atX89JLL/Hnn39mmU/Lli0zFKdwud9er169eP/994mLi7PPAv63CxcusG7duhw3NxYRJ3DuI34iInnH6tWrDWBatGiRYdvXX39twsPDjbu7uylbtqx54403HLZfOXniu+++M3fccYdxc3MzoaGhxhhj4uPjTdOmTY23t7cpU6aMmTJlimncuLEZMGCAfb/Dhw+b1q1bG09PTxMaGmo+//xzU7p0aTN16lR7TFJSkunfv78JCQkx7u7upkyZMqZ79+7mwIEDWV7b9u3bjbe3tzlz5kyGbQcPHjQuLi6mcuXKme77+eefZ7lNRPIWizHX8aCIiIjcEv/88w9lypRh8eLF9mf+rlenTp2oWbMmL7zwQo72u/feexk4cCDdunW7ofOLyM2nW7EiInnI0qVL+e6774iPj2f16tV06dKFcuXK0ahRoxs+9htvvEHRokVztM+xY8d45JFH6Nq16w2fX0RuPo3YiYjkIT///DODBw/m77//xtfXlwYNGvD2228TGhrq7NREJB9QYSciIiJSQOhWrIiIiEgBocJOREREpIBQYSciIiJSQKiwExERESkgVNiJiIiIFBAq7EREREQKCBV2IiIiIgWECjsRERGRAkKFnYiIiEgB8X/ECn+F9xcYQQAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plt.figure()\n", "plt.plot(SMU_mesV1, SMU_mesI1, 'x',label=\"measured\")\n", "plt.plot(SMU_mesV1, a * SMU_mesV1 + b, label=\"fit\")\n", "plt.xlabel(\"Voltage (V)\")\n", "plt.ylabel(\"Current (A)\")\n", "plt.title(f\"Resistance = {1/a:.1f} Ohm\")\n", "plt.legend()\n", "plt.tight_layout()\n" ] }, { "cell_type": "code", "execution_count": null, "id": "33405deb-b38e-47b6-8382-ffbef50e1415", "metadata": {}, "outputs": [], "source": [] } ], "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.10.10" } }, "nbformat": 4, "nbformat_minor": 5 }