{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "760dd9eb-df92-42e2-801a-851ef5a0162b",
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "import gasrari as GR\n",
    "import scipy.optimize as opt"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "85001e4b-52c7-49ef-9bb3-c3ab0326b6dd",
   "metadata": {},
   "outputs": [],
   "source": [
    "en_Ne=GR.energy(\"Ne\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "2025e60f-a86c-4f05-bfd8-0c30f2ff76b2",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "[<matplotlib.lines.Line2D at 0x7f90294d0950>]"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAlwAAAGdCAYAAADHQK08AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAABUEUlEQVR4nO3de3xT9f0/8FeSNklvSVpKmxYKFLkpFIpgsyKKjM5O+/1+6eZXWMekQxyIdRPrxsVJ6/ZVy0B/OhwbXkGnXHXqBAQr1wG1QKFACyJgsVBIuZQkvadNPr8/2h6IpBegp2na1/PxOI8057zPyTvny3d5ec7nnKMQQggQERERkWyUnm6AiIiIqKtj4CIiIiKSGQMXERERkcwYuIiIiIhkxsBFREREJDMGLiIiIiKZMXARERERyYyBi4iIiEhmPp5uoCtzOp04d+4cgoKCoFAoPN0OERERtYEQAuXl5YiMjIRS2T7Hphi4ZHTu3DlERUV5ug0iIiK6CWfOnEHv3r3bZVsMXDIKCgoC0PB/MJ1O5+FuiIiIqC1sNhuioqKk3/H2wMAlo6bTiDqdjoGLiIjIy7TncCAOmiciIiKSGQMXERERkcwYuIiIiIhkxsBFREREJDMGLiIiIiKZMXARERERyYyBi4iIiEhmDFxEREREMmPgIiIiIpIZAxcRERGRzBi4iIiIiGTGwEVEREQkMz682gt9W1qOfx0oQbC/L2aOu83T7RAREVEreITLC529UoVlO07ho7yznm6FiIiI2oCBywvF9DIAAE5erEBlbb1nmyEiIqJWMXB5oZ5BGkTqtRACKDxn83Q7RERE1AoGLi8V01sPADh81uLZRoiIiKhVDFxeanhvAwDg8FmrZxshIiKiVjFweamYXg1HuI6UMHARERF1dgxcXmp44ynFokuVsFbXebgbIiIiagkDl5cy+KvRJ8QfAFDAo1xERESdGgOXF7s6cJ6Bi4iIqDNj4PJiw6VxXBbPNkJEREQtkj1wLV26FP369YNWq4XJZMLevXtbrF+3bh2GDBkCrVaLmJgYbNy40WW5EAIZGRmIiIiAn58fEhIScOLECZeaF198EWPGjIG/vz8MBsN1n7FixQooFAq304ULFwAA27dvd7vcbDbf2g5pR01XKh46wyNcREREnZmsgWvNmjVIT09HZmYmDhw4gBEjRiAxMVEKNT+0Z88epKSkYPr06Th48CCSk5ORnJyMgoICqWbRokVYsmQJli1bhtzcXAQEBCAxMRE1NTVSjd1ux8MPP4xZs2a5/ZzJkyfj/PnzLlNiYiLGjRuHsLAwl9rjx4+71P1wuScN66UDAJRYqnG5otbD3RAREVGzhIzi4uJEWlqa9N7hcIjIyEiRlZXltn7SpEkiKSnJZZ7JZBIzZ84UQgjhdDqF0WgUixcvlpZbLBah0WjEqlWrrtve8uXLhV6vb7XPCxcuCF9fX/H+++9L87Zt2yYAiCtXrrS6fnOsVqsAIKxW601vozXjX94m+s5dL7Z9UyrbZxAREXUncvx+y3aEy263Iy8vDwkJCdI8pVKJhIQE5OTkuF0nJyfHpR4AEhMTpfqioiKYzWaXGr1eD5PJ1Ow22+L999+Hv78//vd///e6ZbGxsYiIiMBPfvIT7N69u8Xt1NbWwmazuUxyk8ZxceA8ERFRpyVb4Lp06RIcDgfCw8Nd5oeHhzc7DspsNrdY3/R6I9tsi3feeQe//OUv4efnJ82LiIjAsmXL8PHHH+Pjjz9GVFQU7rvvPhw4cKDZ7WRlZUGv10tTVFTUTffUVtI4LgYuIiKiTsvH0w14Wk5ODo4dO4Z//vOfLvMHDx6MwYMHS+/HjBmDU6dO4dVXX72utsn8+fORnp4uvbfZbLKHrqYboPJKRSIios5LtiNcoaGhUKlUKC0tdZlfWloKo9Hodh2j0dhifdPrjWyzNW+//TZiY2MxatSoVmvj4uJw8uTJZpdrNBrodDqXSW53ROqgVACltlqU2mpaX4GIiIg6nGyBS61WY9SoUdiyZYs0z+l0YsuWLYiPj3e7Tnx8vEs9AGRnZ0v10dHRMBqNLjU2mw25ubnNbrMlFRUVWLt2LaZPn96m+vz8fERERNzw58jJX+2DgWFBADiOi4iIqLOS9ZRieno6UlNTMXr0aMTFxeG1115DZWUlpk2bBgCYOnUqevXqhaysLADAU089hXHjxuGVV15BUlISVq9ejf379+PNN98EACgUCsyePRsvvPACBg4ciOjoaCxYsACRkZFITk6WPre4uBhlZWUoLi6Gw+FAfn4+AGDAgAEIDAyU6tasWYP6+nr86le/uq731157DdHR0Rg6dChqamrw9ttvY+vWrfjyyy9l2ls3b3hvPY6XluPwWQsS7ghvfQUiIiLqULIGrsmTJ+PixYvIyMiA2WxGbGwsNm3aJA16Ly4uhlJ59SDbmDFjsHLlSjz33HN49tlnMXDgQHz66acYNmyYVDNnzhxUVlZixowZsFgsGDt2LDZt2gStVivVZGRk4L333pPejxw5EgCwbds23HfffdL8d955Bz//+c/d3hzVbrfjmWeeQUlJCfz9/TF8+HB89dVXGD9+fHvtnnYzvLce6/LO4jCfqUhERNQpKYQQwtNNdFU2mw16vR5Wq1XW8Vz5ZyxIXrobPQLU2P9cAhQKhWyfRURE1NXJ8fvNZyl2AUOMQfBRKnC50o5zVg6cJyIi6mwYuLoAra8KQyIaBs4fPmPxbDNERER0HQauLiKmlwEAOI6LiIioE2Lg6iKkG6Dy1hBERESdDgNXFxHT+EzFw2ct4HUQREREnQsDVxcx2BgEtY8Stpp6fH+5ytPtEBER0TUYuLoIX5USd0Q0XLrKcVxERESdCwNXF3J1HJfFs40QERGRCwauLuTqOC4e4SIiIupMGLi6kBFRBgBAQYkVDicHzhMREXUWDFxdyG09A+Hnq0Kl3YGiSxWeboeIiIgaMXB1ISqlAsN6NQ6c52lFIiKiToOBq4uR7jjPwEVERNRpMHB1MSOirt4AlYiIiDoHBq4upulKxcJzNtQ7nB7uhoiIiAAGri6nX48ABGl8UFvvxIkLHDhPRETUGTBwdTFKpQLDevFB1kRERJ0JA1cXNLxxHNchjuMiIiLqFBi4uqDhjVcqHuEzFYmIiDoFBq4uqOmZisfO21Bb7/BwN0RERMTA1QX1DvZDsL8v6hwC35o5cJ6IiMjTGLi6IIVCgZjeBgAcx0VERNQZMHB1UcN5pSIREVGnwcDVRcU0juM6zIHzREREHsfA1UU1DZz/trQcNXUcOE9ERORJDFxdlFGnRc8gDRxOgcJzNk+3Q0RE1K0xcHVRCoXimnFcFs82Q0RE1M0xcHVhHMdFRETUOTBwdWFN47gO80pFIiIij2Lg6sJiGh/xc+piBSpq6z3bDBERUTfGwNWF9QzSIFKvhRBAIU8rEhEReYzsgWvp0qXo168ftFotTCYT9u7d22L9unXrMGTIEGi1WsTExGDjxo0uy4UQyMjIQEREBPz8/JCQkIATJ0641Lz44osYM2YM/P39YTAY3H6OQqG4blq9erVLzfbt23HnnXdCo9FgwIABWLFixQ1/f09rGsfFB1kTERF5jqyBa82aNUhPT0dmZiYOHDiAESNGIDExERcuXHBbv2fPHqSkpGD69Ok4ePAgkpOTkZycjIKCAqlm0aJFWLJkCZYtW4bc3FwEBAQgMTERNTU1Uo3dbsfDDz+MWbNmtdjf8uXLcf78eWlKTk6WlhUVFSEpKQnjx49Hfn4+Zs+ejcceewybN2++tZ3SwYZLj/hh4CIiIvIYIaO4uDiRlpYmvXc4HCIyMlJkZWW5rZ80aZJISkpymWcymcTMmTOFEEI4nU5hNBrF4sWLpeUWi0VoNBqxatWq67a3fPlyodfr3X4WAPHJJ5802/ucOXPE0KFDXeZNnjxZJCYmNrvOD1mtVgFAWK3WNq/T3nZ+e0H0nbtejFu01WM9EBEReRM5fr9lO8Jlt9uRl5eHhIQEaZ5SqURCQgJycnLcrpOTk+NSDwCJiYlSfVFREcxms0uNXq+HyWRqdpstSUtLQ2hoKOLi4vDuu+9CCNHmXtypra2FzWZzmTwtpvFeXKcvV8FaVefhboiIiLon2QLXpUuX4HA4EB4e7jI/PDwcZrPZ7Tpms7nF+qbXG9lmc/785z9j7dq1yM7OxkMPPYQnnngCr7/+equ92Gw2VFdXu91mVlYW9Hq9NEVFRd1QT3Iw+KvRJ8QfAFBwjqcViYiIPMHH0w14yoIFC6S/R44cicrKSixevBi/+93vbnqb8+fPR3p6uvTeZrN1itAV01uP4rIqHDprwd0DQj3dDhERUbcj2xGu0NBQqFQqlJaWuswvLS2F0Wh0u47RaGyxvun1RrbZViaTCWfPnkVtbW2Lveh0Ovj5+bndhkajgU6nc5k6gxFNVypy4DwREZFHyBa41Go1Ro0ahS1btkjznE4ntmzZgvj4eLfrxMfHu9QDQHZ2tlQfHR0No9HoUmOz2ZCbm9vsNtsqPz8fwcHB0Gg0berFmzTdAJV3nCciIvIMWU8ppqenIzU1FaNHj0ZcXBxee+01VFZWYtq0aQCAqVOnolevXsjKygIAPPXUUxg3bhxeeeUVJCUlYfXq1di/fz/efPNNAA33zpo9ezZeeOEFDBw4ENHR0ViwYAEiIyNdbulQXFyMsrIyFBcXw+FwID8/HwAwYMAABAYG4vPPP0dpaSl+9KMfQavVIjs7Gy+99BJ+//vfS9t4/PHH8be//Q1z5szBo48+iq1bt2Lt2rXYsGGDnLtMFsN6NRxpK7FU43JFLXoEajzcERERUTfTbtc7NuP1118Xffr0EWq1WsTFxYmvv/5aWjZu3DiRmprqUr927VoxaNAgoVarxdChQ8WGDRtcljudTrFgwQIRHh4uNBqNmDBhgjh+/LhLTWpqqgBw3bRt2zYhhBBffPGFiI2NFYGBgSIgIECMGDFCLFu2TDgcDpftbNu2TcTGxgq1Wi369+8vli9ffkPfvTPcFqLJ+Je3ib5z14ut35R6uhUiIqJOTY7fb4UQ19wLgdqVzWaDXq+H1Wr1+Hiup9fk45ODJUj/ySD8bsJAj/ZCRETUmcnx+81nKXYTTffj4jguIiKijsfA1U0Ml56paPFsI0RERN0QA1c3cUekDkoFUGqrRamtpvUViIiIqN0wcHUT/mofDAoPAsDTikRERB2NgasbaRrHdeSsxbONEBERdTMMXN1I0ziuwyU8wkVERNSRGLi6kZjeBgANpxR5NxAiIqKOw8DVjdweEQRflQJllXaUWKo93Q4REVG3wcDVjWh8VBhsbBg4zwdZExERdRwGrm5GepA1x3ERERF1GAaubkYaOM8rFYmIiDoMA1c3czVwceA8ERFRR2Hg6mYGhQdB7aNEeU09vr9c5el2iIiIugUGrm7GV6XEHRENTz7nOC4iIqKOwcDVDUmnFc9YPNsIERFRN8HA1Q0Nb7oBKo9wERERdQgGrm6o6QhXYYkVDicHzhMREcmNgasbuq1nIPx8Vai0O1B0qcLT7RAREXV5DFzdkEqpwLBeDQPnD53haUUiIiK5MXB1U03juI5wHBcREZHsGLi6Kd5xnoiIqOMwcHVTMb0aB86fs6HO4fRwN0RERF0bA1c31a9HAII0Pqitd+JEKQfOExERyYmBq5tSKhWIaTyteKTE4tlmiIiIujgGrm4s5poHWRMREZF8GLi6seG9DAAYuIiIiOTGwNWNNV2p+I3Zhtp6h4e7ISIi6roYuLqx3sF+CPb3RZ1D4Li53NPtEBERdVkMXN2YQqFATNODrHlakYiISDYMXN3c8F68ASoREZHcGLi6OV6pSEREJD/ZA9fSpUvRr18/aLVamEwm7N27t8X6devWYciQIdBqtYiJicHGjRtdlgshkJGRgYiICPj5+SEhIQEnTpxwqXnxxRcxZswY+Pv7w2AwXPcZhw4dQkpKCqKiouDn54fbb78df/3rX11qtm/fDoVCcd1kNptvbkd0UiMaTymeuFCBajsHzhMREclB1sC1Zs0apKenIzMzEwcOHMCIESOQmJiICxcuuK3fs2cPUlJSMH36dBw8eBDJyclITk5GQUGBVLNo0SIsWbIEy5YtQ25uLgICApCYmIiamhqpxm634+GHH8asWbPcfk5eXh7CwsLwwQcfoLCwEH/84x8xf/58/O1vf7uu9vjx4zh//rw0hYWF3eJe6VzCdRr0DNLA4RQ4et7m6XaIiIi6JiGjuLg4kZaWJr13OBwiMjJSZGVlua2fNGmSSEpKcplnMpnEzJkzhRBCOJ1OYTQaxeLFi6XlFotFaDQasWrVquu2t3z5cqHX69vU6xNPPCHGjx8vvd+2bZsAIK5cudKm9d2xWq0CgLBarTe9jY7w6PK9ou/c9eLdXd95uhUiIiKPk+P3W7YjXHa7HXl5eUhISJDmKZVKJCQkICcnx+06OTk5LvUAkJiYKNUXFRXBbDa71Oj1ephMpma32VZWqxUhISHXzY+NjUVERAR+8pOfYPfu3S1uo7a2FjabzWXyBtIjfjiOi4iISBayBa5Lly7B4XAgPDzcZX54eHiz46DMZnOL9U2vN7LNttizZw/WrFmDGTNmSPMiIiKwbNkyfPzxx/j4448RFRWF++67DwcOHGh2O1lZWdDr9dIUFRV10z11pKZxXIdLGLiIiIjk4OPpBjytoKAAEydORGZmJu6//35p/uDBgzF48GDp/ZgxY3Dq1Cm8+uqr+Oc//+l2W/Pnz0d6err03mazeUXoGtZ4a4hTFytQUVuPQE23/2dBRETUrmQ7whUaGgqVSoXS0lKX+aWlpTAajW7XMRqNLdY3vd7INlty9OhRTJgwATNmzMBzzz3Xan1cXBxOnjzZ7HKNRgOdTucyeYOeQRpE6rUQAijgUS4iIqJ2J1vgUqvVGDVqFLZs2SLNczqd2LJlC+Lj492uEx8f71IPANnZ2VJ9dHQ0jEajS43NZkNubm6z22xOYWEhxo8fj9TUVLz44ottWic/Px8RERE39DneguO4iIiI5CPruaP09HSkpqZi9OjRiIuLw2uvvYbKykpMmzYNADB16lT06tULWVlZAICnnnoK48aNwyuvvIKkpCSsXr0a+/fvx5tvvgmg4VE0s2fPxgsvvICBAwciOjoaCxYsQGRkJJKTk6XPLS4uRllZGYqLi+FwOJCfnw8AGDBgAAIDA1FQUIAf//jHSExMRHp6ujT+S6VSoWfPngCA1157DdHR0Rg6dChqamrw9ttvY+vWrfjyyy/l3GUeM7y3AZsLSzmOi4iISAayBq7Jkyfj4sWLyMjIgNlsRmxsLDZt2iQNei8uLoZSefUg25gxY7By5Uo899xzePbZZzFw4EB8+umnGDZsmFQzZ84cVFZWYsaMGbBYLBg7diw2bdoErVYr1WRkZOC9996T3o8cORIAsG3bNtx333346KOPcPHiRXzwwQf44IMPpLq+ffvi9OnTABqusnzmmWdQUlICf39/DB8+HF999RXGjx8vy77ytOHSES6LZxshIiLqghRCCOHpJroqm80GvV4Pq9Xa6cdzWarsiP1zNgDgUMb90Pv7ergjIiIiz5Dj95vPUiQAgMFfjT4h/gCAfB7lIiIialcMXCQxRTfc+HXXiYse7oSIiKhrYeAiyb2DGi4Y2PntJQ93QkRE1LUwcJFk7IBQKBTA8dJymK01ra9AREREbcLARZLgADWGN951/j88rUhERNRuGLjIhXRa8QRPKxIREbUXBi5y0RS4dp24CIeTdwwhIiJqDwxc5CI2yoAgjQ+uVNXxuYpERETthIGLXPiqlIi/rQcAjuMiIiJqLwxcdB3eHoKIiKh9MXDRdcY1Bq4DxVdQXlPn4W6IiIi8HwMXXScqxB/RoQGodwrsOXXZ0+0QERF5PQYucuuegaEAOI6LiIioPTBwkVv3DuQ4LiIiovbCwEVuxd/WA74qBYrLqnD6UqWn2yEiIvJqDFzkVoDGB6P6BgMAdvK0IhER0S1h4KJm3SOdVmTgIiIiuhUMXNSspttD5Jy6DHu908PdEBEReS8GLmrWHRE69AhQo9LuwIHiK55uh4iIyGsxcFGzlEqFdHsInlYkIiK6eQxc1CLpMT8cOE9ERHTTGLioRWMbj3AVlNhwuaLWw90QERF5JwYualFYkBa3R+gAALtO8iaoREREN4OBi1p176CGo1w7OI6LiIjopjBwUavGNd6P6z8nLkEI4eFuiIiIvA8DF7VqVL9g+PmqcLG8Ft+Yyz3dDhERkddh4KJWaXxU+FH/EAC8PQQREdHNYOCiNuHtIYiIiG4eAxe1SVPg2ld0BVX2eg93Q0RE5F0YuKhN+ocGoJfBD3aHE7nflXm6HSIiIq/CwEVtolAopNtD8LQiERHRjZE9cC1duhT9+vWDVquFyWTC3r17W6xft24dhgwZAq1Wi5iYGGzcuNFluRACGRkZiIiIgJ+fHxISEnDixAmXmhdffBFjxoyBv78/DAaD288pLi5GUlIS/P39ERYWhj/84Q+or3c9VbZ9+3bceeed0Gg0GDBgAFasWHHD378rubfx9hAcOE9ERHRjZA1ca9asQXp6OjIzM3HgwAGMGDECiYmJuHDhgtv6PXv2ICUlBdOnT8fBgweRnJyM5ORkFBQUSDWLFi3CkiVLsGzZMuTm5iIgIACJiYmoqamRaux2Ox5++GHMmjXL7ec4HA4kJSXBbrdjz549eO+997BixQpkZGRINUVFRUhKSsL48eORn5+P2bNn47HHHsPmzZvbae94nzEDQqFUAKcuVqLEUu3pdoiIiLyHkFFcXJxIS0uT3jscDhEZGSmysrLc1k+aNEkkJSW5zDOZTGLmzJlCCCGcTqcwGo1i8eLF0nKLxSI0Go1YtWrVddtbvny50Ov1183fuHGjUCqVwmw2S/P+8Y9/CJ1OJ2pra4UQQsyZM0cMHTrUZb3JkyeLxMTEVr71VVarVQAQVqu1zet0dj9bukv0nbterMz93tOtEBERyUKO32/ZjnDZ7Xbk5eUhISFBmqdUKpGQkICcnBy36+Tk5LjUA0BiYqJUX1RUBLPZ7FKj1+thMpma3WZznxMTE4Pw8HCXz7HZbCgsLGxTL+7U1tbCZrO5TF1N09WK/+E4LiIiojaTLXBdunQJDofDJdQAQHh4OMxms9t1zGZzi/VNrzeyzRv5nGs/o7kam82G6mr3p9OysrKg1+ulKSoqqs09eYumwLXrxCXUO5we7oaIiMg78CrFdjR//nxYrVZpOnPmjKdbancjehug0/rAVlOPQ2etnm6HiIjIK8gWuEJDQ6FSqVBaWuoyv7S0FEaj0e06RqOxxfqm1xvZ5o18zrWf0VyNTqeDn5+f2+1qNBrodDqXqatRKRUYO7Dx9hC8WpGIiKhNZAtcarUao0aNwpYtW6R5TqcTW7ZsQXx8vNt14uPjXeoBIDs7W6qPjo6G0Wh0qbHZbMjNzW12m819zpEjR1yulszOzoZOp8Mdd9zRpl66s6bbQ3AcFxERUdvIekoxPT0db731Ft577z0cO3YMs2bNQmVlJaZNmwYAmDp1KubPny/VP/XUU9i0aRNeeeUVfPPNN3j++eexf/9+PPnkkwAabr45e/ZsvPDCC/j3v/+NI0eOYOrUqYiMjERycrK0neLiYuTn56O4uBgOhwP5+fnIz89HRUUFAOD+++/HHXfcgUceeQSHDh3C5s2b8dxzzyEtLQ0ajQYA8Pjjj+O7777DnDlz8M033+Dvf/871q5di6efflrOXeYVmsZx5Z+xwFpV5+FuiIiIvEC7Xe/YjNdff1306dNHqNVqERcXJ77++mtp2bhx40RqaqpL/dq1a8WgQYOEWq0WQ4cOFRs2bHBZ7nQ6xYIFC0R4eLjQaDRiwoQJ4vjx4y41qampAsB107Zt26Sa06dPiwceeED4+fmJ0NBQ8cwzz4i6ujqX7Wzbtk3ExsYKtVot+vfvL5YvX35D370r3haiyYRXtou+c9eLDYfPeboVIiKidiXH77dCCCE8mPe6NJvNBr1eD6vV2uXGc/3586N4d3cRfnFXFBY+NNzT7RAREbUbOX6/eZUi3ZR7Gp+r+J8Tl8DMTkRE1DIGLropP4ruAbWPEiWWapy6WOnpdoiIiDo1Bi66KX5qFeL6hQDg7SGIiIhaw8BFN+3extOKO3l7CCIiohYxcNFNu6fxflxff3cZNXUOD3dDRETUeTFw0U0bYgxCWJAGNXVO5H1/xdPtEBERdVoMXHTTFAqFdJSL47iIiIiax8BFt6RpHNcOBi4iIqJmMXDRLRk7IBQKBfCNuRwXbDWeboeIiKhTYuCiW9IjUINhkXoADTdBJSIiousxcNEt4+0hiIiIWsbARbfs3saB8/85cQlOJx/zQ0RE9EMMXHTLRvYJRoBahbJKOwrP2TzdDhERUafDwEW3TO2jRPxtPK1IRETUHAYuahfjmsZx8fYQRERE12HgonZx76CGcVx5319BRW29h7shIiLqXBi4qF307RGAPiH+qHcK5Jy67Ol2iIiIOhUGLmo39/K0IhERkVsMXNRurt4egoGLiIjoWgxc1G7ib+sBH6UCpy9XofhylafbISIi6jQYuKjdBGl9cWffYADADh7lIiIikjBwUbu6dyDHcREREf0QAxe1q6bbQ+Scuow6h9PD3RAREXUODFzUroZF6hESoEZFbT0OFls83Q4REVGnwMBF7UqpVGDsAJ5WJCIiuhYDF7W7ewbyuYpERETXYuCidtc0jutIiRVllXYPd0NEROR5DFzU7sJ1WgwxBkEIYNfJS55uh4iIyOMYuEgWTUe5OI6LiIiIgYtk0jSO6z8nLkII4eFuiIiIPIuBi2RxV78QaH2VKLXV4tvSCk+3Q0RE5FGyB66lS5eiX79+0Gq1MJlM2Lt3b4v169atw5AhQ6DVahETE4ONGze6LBdCICMjAxEREfDz80NCQgJOnDjhUlNWVoYpU6ZAp9PBYDBg+vTpqKi4+qP//PPPQ6FQXDcFBARINStWrLhuuVarbYc90j1ofVUwRfcAwNOKREREsgauNWvWID09HZmZmThw4ABGjBiBxMREXLhwwW39nj17kJKSgunTp+PgwYNITk5GcnIyCgoKpJpFixZhyZIlWLZsGXJzcxEQEIDExETU1NRINVOmTEFhYSGys7Oxfv167Ny5EzNmzJCW//73v8f58+ddpjvuuAMPP/ywSz86nc6l5vvvv2/nPdS1SeO4eHsIIiLq7oSM4uLiRFpamvTe4XCIyMhIkZWV5bZ+0qRJIikpyWWeyWQSM2fOFEII4XQ6hdFoFIsXL5aWWywWodFoxKpVq4QQQhw9elQAEPv27ZNqvvjiC6FQKERJSYnbz83PzxcAxM6dO6V5y5cvF3q9/sa+8A9YrVYBQFit1lvajrf61mwTfeeuFwP/uFFU1dZ7uh0iIqI2keP3W7YjXHa7HXl5eUhISJDmKZVKJCQkICcnx+06OTk5LvUAkJiYKNUXFRXBbDa71Oj1ephMJqkmJycHBoMBo0ePlmoSEhKgVCqRm5vr9nPffvttDBo0CPfcc4/L/IqKCvTt2xdRUVGYOHEiCgsLW/zOtbW1sNlsLlN3NiAsEBF6Lez1TuQWXfZ0O0RERB4jW+C6dOkSHA4HwsPDXeaHh4fDbDa7XcdsNrdY3/TaWk1YWJjLch8fH4SEhLj93JqaGnz44YeYPn26y/zBgwfj3XffxWeffYYPPvgATqcTY8aMwdmzZ5v9zllZWdDr9dIUFRXVbG13oFAocO/AhtOK/znB+3EREVH31e2vUvzkk09QXl6O1NRUl/nx8fGYOnUqYmNjMW7cOPzrX/9Cz5498cYbbzS7rfnz58NqtUrTmTNn5G6/0xs3uCFwbSoww+nk7SGIiKh7ki1whYaGQqVSobS01GV+aWkpjEaj23WMRmOL9U2vrdX8cFB+fX09ysrK3H7u22+/jf/6r/+67qjZD/n6+mLkyJE4efJkszUajQY6nc5l6u5+PCQMQVoflFiqsecUTysSEVH3JFvgUqvVGDVqFLZs2SLNczqd2LJlC+Lj492uEx8f71IPANnZ2VJ9dHQ0jEajS43NZkNubq5UEx8fD4vFgry8PKlm69atcDqdMJlMLtsuKirCtm3brjud6I7D4cCRI0cQERHRai1dpfVVYWJsJABg7X4e8SMiom6q3Ybfu7F69Wqh0WjEihUrxNGjR8WMGTOEwWAQZrNZCCHEI488IubNmyfV7969W/j4+IiXX35ZHDt2TGRmZgpfX19x5MgRqWbhwoXCYDCIzz77TBw+fFhMnDhRREdHi+rqaqnmpz/9qRg5cqTIzc0Vu3btEgMHDhQpKSnX9ffcc8+JyMhIUV9//RV0f/rTn8TmzZvFqVOnRF5envjFL34htFqtKCwsbPP37+5XKTY5dOaKdLWipdLu6XaIiIhaJMfvt4+cYW7y5Mm4ePEiMjIyYDabERsbi02bNkmn74qLi6FUXj3INmbMGKxcuRLPPfccnn32WQwcOBCffvophg0bJtXMmTMHlZWVmDFjBiwWC8aOHYtNmza53JT0ww8/xJNPPokJEyZAqVTioYcewpIlS1x6czqdWLFiBX79619DpVJd1/uVK1fwm9/8BmazGcHBwRg1ahT27NmDO+64o713U5cX00uPIcYgfGMux78PleCR+H6ebomIiKhDKYTgg+7kYrPZoNfrYbVau/14rnd3FeHP648ippcen/92rKfbISIiapYcv9/d/ipF6hjJI3vBV6XAkRIrjp7r3vcnIyKi7oeBizpESIAaP7mj4VQyB88TEVF3w8BFHebh0Q03gv00vwS19Q4Pd0NERNRxGLiow9w7sCeMOi0sVXX46qj7B5gTERF1RQxc1GFUSgX+d1RvADytSERE3QsDF3WopsC188RFnLNUe7gbIiKijsHARR2qX2gATNEhEAL4OK/5B4ETERF1JQxc1OEmNQ6eX5d3lg+0JiKiboGBizrcAzFGBGp8UFxWhdyiMk+3Q0REJDsGLupw/mof/PeIhgdar+PgeSIi6gYYuMgjJo1uGDy/seA8bDV1Hu6GiIhIXgxc5BGxUQYMDAtETZ0Tnx865+l2iIiIZMXARR6hUCikwfNr9/NqRSIi6toYuMhjfnZnL/goFTh0xoLj5nJPt0NERCQbBi7ymNBADSbcHgaAg+eJiKhrY+Aij2o6rfjJwRLY650e7oaIiEgeDFzkUeMG9UTPIA0uV9qx9Rs+0JqIiLomBi7yKB+VEg/dyQdaExFR18bARR73cOM9ubYfv4BSW42HuyEiImp/DFzkcbf1DMTovsFwCuDjA7xFBBERdT0MXNQpTLqr8YHW+89CCD7QmoiIuhYGLuoUkmIi4K9WoehSJfZ/f8XT7RAREbUrBi7qFAI0Pviv4REAgLX7OHieiIi6FgYu6jSa7sm14ch5VNTWe7gbIiKi9sPARZ3GqL7B6N8zAFV2BzYc5gOtiYio62Dgok6DD7QmIqKuioGLOpWfj+wFlVKBvO+v4OSFCk+3Q0RE1C4YuKhTCdNpMX5wTwDAujwOnicioq6BgYs6nYcbTyt+nFeCOgcfaE1ERN6PgYs6nR8PCUNooBqXKmqx/fhFT7dDRER0yxi4qNPxVSnxcz7QmoiIuhDZA9fSpUvRr18/aLVamEwm7N27t8X6devWYciQIdBqtYiJicHGjRtdlgshkJGRgYiICPj5+SEhIQEnTpxwqSkrK8OUKVOg0+lgMBgwffp0VFRcHYB9+vRpKBSK66avv/76hnoh+Tw8qiFwbf3mAi6U84HWRETk3WQNXGvWrEF6ejoyMzNx4MABjBgxAomJibhw4YLb+j179iAlJQXTp0/HwYMHkZycjOTkZBQUFEg1ixYtwpIlS7Bs2TLk5uYiICAAiYmJqKm5+qM8ZcoUFBYWIjs7G+vXr8fOnTsxY8aM6z7vq6++wvnz56Vp1KhRN9QLyWdgeBBG9jHA4RT49GCJp9shIiK6NUJGcXFxIi0tTXrvcDhEZGSkyMrKcls/adIkkZSU5DLPZDKJmTNnCiGEcDqdwmg0isWLF0vLLRaL0Gg0YtWqVUIIIY4ePSoAiH379kk1X3zxhVAoFKKkpEQIIURRUZEAIA4ePNhs76310hZWq1UAEFartc3r0FUrc78XfeeuFxNe2S6cTqen2yEiom5Cjt9v2Y5w2e125OXlISEhQZqnVCqRkJCAnJwct+vk5OS41ANAYmKiVF9UVASz2exSo9frYTKZpJqcnBwYDAaMHj1aqklISIBSqURubq7Ltv/nf/4HYWFhGDt2LP7973/fUC/u1NbWwmazuUx08/5reAS0vkqcvFCBA8UWT7dDRER002QLXJcuXYLD4UB4eLjL/PDwcJjNZrfrmM3mFuubXlurCQsLc1nu4+ODkJAQqSYwMBCvvPIK1q1bhw0bNmDs2LFITk52CV2t9eJOVlYW9Hq9NEVFRTVbS60L0vriwZiGB1qv4+B5IiLyYt3yKsXQ0FCkp6fDZDLhrrvuwsKFC/GrX/0KixcvvqXtzp8/H1arVZrOnGFIuFWTG+/J9fmhc6iy84HWRETknWQLXKGhoVCpVCgtLXWZX1paCqPR6HYdo9HYYn3Ta2s1PxyUX19fj7KysmY/FwBMJhNOnjzZ5l7c0Wg00Ol0LhPdmrjoEPTr4Y9KuwMbjzR/dJGIiKgzky1wqdVqjBo1Clu2bJHmOZ1ObNmyBfHx8W7XiY+Pd6kHgOzsbKk+OjoaRqPRpcZmsyE3N1eqiY+Ph8ViQV5enlSzdetWOJ1OmEymZvvNz89HREREm3uhjqFQKKQ7z/OeXERE5LXabfi9G6tXrxYajUasWLFCHD16VMyYMUMYDAZhNpuFEEI88sgjYt68eVL97t27hY+Pj3j55ZfFsWPHRGZmpvD19RVHjhyRahYuXCgMBoP47LPPxOHDh8XEiRNFdHS0qK6ulmp++tOfipEjR4rc3Fyxa9cuMXDgQJGSkiItX7FihVi5cqU4duyYOHbsmHjxxReFUqkU77777g310hpepdg+zlmqRPS89aLv3PXiu4sVnm6HiIi6ODl+v2UNXEII8frrr4s+ffoItVot4uLixNdffy0tGzdunEhNTXWpX7t2rRg0aJBQq9Vi6NChYsOGDS7LnU6nWLBggQgPDxcajUZMmDBBHD9+3KXm8uXLIiUlRQQGBgqdTiemTZsmysvLpeUrVqwQt99+u/D39xc6nU7ExcWJdevWXdd7a720hoGr/aS+myv6zl0v/vLFMU+3QkREXZwcv98KIYTw7DG2rstms0Gv18NqtXI81y364sh5zPrwAMJ1Guye+2P4qLrl9R5ERNQB5Pj95q8WeYUJt4cjJECNUlst/nPikqfbISIiuiEMXOQV1D5KJMf2AsDB80RE5H0YuMhrTLqr4YHWXx0rxeWKWg93Q0RE1HYMXOQ1hhh1GN5bjzqHwKf55zzdDhERUZsxcJFXmdR0T659Z8DrPYiIyFswcJFX+e8RkdD4KHG8tByHz1o93Q4REVGbMHCRV9H7+eKBYQ2PV+LgeSIi8hYMXOR1mk4r/jv/HKrtDg93Q0RE1DoGLvI6P+rfA31C/FFeW4/3c057uh0iIqJWMXCR11EqFfjdhIEAgL9tO4krlXYPd0RERNQyBi7ySj8b2Qu3R+hQXlOPv2076el2iIiIWsTARV5JpVRg/gNDAADv55zGmbIqD3dERETUPAYu8lr3DuqJewaGos4hsGjzcU+3Q0RE1CwGLvJq8x+4HQoF8Pmhczh0xuLpdoiIiNxi4CKvdkekDj8f2fCMxZc2HuPd54mIqFNi4CKv98z9g6DxUSK3qAxbjl3wdDtERETXYeAirxdp8MOjY6MBAAs3fYN6h9PDHREREbli4KIuYdZ9tyHY3xcnL1Rg7f6znm6HiIjIBQMXdQk6ra90M9RXv/oWlbX1Hu6IiIjoKgYu6jKmmPqibw9/XCyvxVv/+c7T7RAREUkYuKjLUPsoMSex4Waob+78DhfKazzcERERUQMGLupSHowxIjbKgCq7A699dcLT7RAREQFg4KIuRqFQ4I9JtwMA1uw7g5MXyj3cEREREQMXdUF39QvB/XeEw+EUWPgFH/lDRESex8BFXdLcB4ZApVTgq2OlyP3usqfbISKibo6Bi7qk23oGIiUuCgAf+UNERJ7HwEVd1lMTBiFArcKhs1asP3ze0+0QEVE3xsBFXVbPIA1mjrsNALBo8zeorXd4uCMiIuquGLioS3vsnmiEBWlwpqwaH3xd7Ol2iIiom2Lgoi7NX+2D9J8MAgC8vvUErNV1Hu6IiIi6IwYu6vL+d1RvDAwLhKWqDn/fftLT7RARUTcke+BaunQp+vXrB61WC5PJhL1797ZYv27dOgwZMgRarRYxMTHYuHGjy3IhBDIyMhAREQE/Pz8kJCTgxAnXO4qXlZVhypQp0Ol0MBgMmD59OioqKqTl27dvx8SJExEREYGAgADExsbiww8/dNnGihUroFAoXCatVnuLe4M8wUelxPwHGx75s3z3aZy9UuXhjoiIqLuRNXCtWbMG6enpyMzMxIEDBzBixAgkJibiwoULbuv37NmDlJQUTJ8+HQcPHkRycjKSk5NRUFAg1SxatAhLlizBsmXLkJubi4CAACQmJqKm5upz86ZMmYLCwkJkZ2dj/fr12LlzJ2bMmOHyOcOHD8fHH3+Mw4cPY9q0aZg6dSrWr1/v0o9Op8P58+el6fvvv2/nPUQdZfzgMPyofwjs9U78vy+/9XQ7RETU3QgZxcXFibS0NOm9w+EQkZGRIisry239pEmTRFJSkss8k8kkZs6cKYQQwul0CqPRKBYvXiwtt1gsQqPRiFWrVgkhhDh69KgAIPbt2yfVfPHFF0KhUIiSkpJme33wwQfFtGnTpPfLly8Xer2+7V/WDavVKgAIq9V6S9uh9nH4jEX0nbte9Ju3Xhw5a/F0O0RE1EnJ8fst2xEuu92OvLw8JCQkSPOUSiUSEhKQk5Pjdp2cnByXegBITEyU6ouKimA2m11q9Ho9TCaTVJOTkwODwYDRo0dLNQkJCVAqlcjNzW22X6vVipCQEJd5FRUV6Nu3L6KiojBx4kQUFha2+J1ra2ths9lcJuo8YnrrMTE2EkIAWV/wZqhERNRxZAtcly5dgsPhQHh4uMv88PBwmM1mt+uYzeYW65teW6sJCwtzWe7j44OQkJBmP3ft2rXYt28fpk2bJs0bPHgw3n33XXz22Wf44IMP4HQ6MWbMGJw9e7bZ75yVlQW9Xi9NUVFRzdaSZ/z+/sFQq5TYffIydnx70dPtEBFRN9Htr1Lctm0bpk2bhrfeegtDhw6V5sfHx2Pq1KmIjY3FuHHj8K9//Qs9e/bEG2+80ey25s+fD6vVKk1nzpzpiK9ANyAqxB+pY/oCALI2fgOHk0e5iIhIfrIFrtDQUKhUKpSWlrrMLy0thdFodLuO0Whssb7ptbWaHw7Kr6+vR1lZ2XWfu2PHDvz3f/83Xn31VUydOrXF7+Pr64uRI0fi5Mnmbyug0Wig0+lcJup80sYPgE7rg+Ol5fj4QPNHLImIiNqLbIFLrVZj1KhR2LJlizTP6XRiy5YtiI+Pd7tOfHy8Sz0AZGdnS/XR0dEwGo0uNTabDbm5uVJNfHw8LBYL8vLypJqtW7fC6XTCZDJJ87Zv346kpCT85S9/cbmCsTkOhwNHjhxBREREG749dWYGfzV+++OBAIBXvjyOajsf+UNERDJrt+H3bqxevVpoNBqxYsUKcfToUTFjxgxhMBiE2WwWQgjxyCOPiHnz5kn1u3fvFj4+PuLll18Wx44dE5mZmcLX11ccOXJEqlm4cKEwGAzis88+E4cPHxYTJ04U0dHRorq6Wqr56U9/KkaOHClyc3PFrl27xMCBA0VKSoq0fOvWrcLf31/Mnz9fnD9/XpouX74s1fzpT38SmzdvFqdOnRJ5eXniF7/4hdBqtaKwsLDN359XKXZe1fZ6MSZri+g7d73429YTnm6HiIg6ETl+v2UNXEII8frrr4s+ffoItVot4uLixNdffy0tGzdunEhNTXWpX7t2rRg0aJBQq9Vi6NChYsOGDS7LnU6nWLBggQgPDxcajUZMmDBBHD9+3KXm8uXLIiUlRQQGBgqdTiemTZsmysvLpeWpqakCwHXTuHHjpJrZs2dLfYeHh4sHH3xQHDhw4Ia+OwNX5/bpwbOi79z1YmjGJnGxvMbT7RARUSchx++3QgheGy8Xm80GvV4Pq9XK8VydkNMp8D9Ld6GgxIbU+L7408Rhnm6JiIg6ATl+v7v9VYrUfSmVCjz74O0AgA9zi/HdxYpW1iAiIro5DFzUrY25LRTjB/dEvVNg8ebjnm6HiIi6KAYu6vbmPXA7lArgiwIz8r4v83Q7RETUBTFwUbc32BiEh0c1PBXgxQ185A8REbU/Bi4iAOn3D4KfrwoHii3YVOD+EVBEREQ3i4GLCEC4Tovf3BMNAPjjpwU4U1bl4Y6IiKgrYeAiajTrvgEY1kuHsko7fvP+flTW1nu6JSIi6iIYuIga+alVePOR0QgNVOMbczmeWXsITj7cmoiI2gEDF9E1Ig1+eOORUfBVKbCp0IwlW094uiUiIuoCGLiIfmBU3xC8mBwDAHjtqxP44sh5D3dERETejoGLyI1Jd0Vh2t39AADpaw/h6DmbZxsiIiKvxsBF1Iw/Png7xg4IRXWdA795fz8uV9R6uiUiIvJSDFxEzfBRKfG3X45Evx7+KLFUY9aHB2Cvd3q6LSIi8kIMXEQtMPir8XbqaARqfLC3qAzPf17o6ZaIiMgLMXARtWJAWBD++otYKBTAytxi/PPr7z3dEhEReRkGLqI2mHB7OOYkDgEA/Onfhdhz6pKHOyIiIm/CwEXURo+P64+JsZGodwqkfXiAj/8hIqI2Y+AiaiOFQoG/PDQcMb30uFJVh8fe248KPv6HiIjagIGL6AZofVV4c+oo9AzS4HhpOdLX5PPxP0RE1CoGLqIbFKFvePyPWqXEl0dL8dpX33q6JSIi6uQYuIhuwp19gvHSzxse/7Nk60lsOMzH/xARUfMYuIhu0v+O6o3HxkYDAJ5Zl4+CEquHOyIios6KgYvoFsx7YAjuGRiKmjonZry/H5f4+B8iInKDgYvoFviolPhbyp2IDg3AOWsNZn2Qx8f/EBHRdRi4iG6R3t8Xb00djSCND/advoKMzwogBK9cJCKiqxi4iNrBgLBALPnlSCgUwOp9Z/B+Dh//Q0REVzFwEbWT8YPDMO+nDY//+fP6o9hzko//ISKiBgxcRO1oxr398bORveBwCjyx8gC+v1zp6ZaIiKgTYOAiakcKhQJZP4/BiCgDLFV1+M37fPwPERExcBG1O62vCm8+MgphQRp8W1qB2av5+B8iou6OgYtIBuE6bcPjf3yU+OpYKV7JPu7ploiIyINkD1xLly5Fv379oNVqYTKZsHfv3hbr161bhyFDhkCr1SImJgYbN250WS6EQEZGBiIiIuDn54eEhAScOHHCpaasrAxTpkyBTqeDwWDA9OnTUVFR4VJz+PBh3HPPPdBqtYiKisKiRYtuuBeilozsE4yFjY//WbrtFP74yRHU1js83BUREXmCrIFrzZo1SE9PR2ZmJg4cOIARI0YgMTERFy5ccFu/Z88epKSkYPr06Th48CCSk5ORnJyMgoICqWbRokVYsmQJli1bhtzcXAQEBCAxMRE1NTVSzZQpU1BYWIjs7GysX78eO3fuxIwZM6TlNpsN999/P/r27Yu8vDwsXrwYzz//PN58880b6oWoNT+/szfm/nQIFArgw9xiPLwsB2fKqjzdFhERdTQho7i4OJGWlia9dzgcIjIyUmRlZbmtnzRpkkhKSnKZZzKZxMyZM4UQQjidTmE0GsXixYul5RaLRWg0GrFq1SohhBBHjx4VAMS+ffukmi+++EIoFApRUlIihBDi73//uwgODha1tbVSzdy5c8XgwYPb3EtbWK1WAUBYrdY2r0Nd07ZvSsWIP20WfeeuFyP+tFls/abU0y0REVEz5Pj9lu0Il91uR15eHhISEqR5SqUSCQkJyMnJcbtOTk6OSz0AJCYmSvVFRUUwm80uNXq9HiaTSarJycmBwWDA6NGjpZqEhAQolUrk5uZKNffeey/UarXL5xw/fhxXrlxpUy/u1NbWwmazuUxEAHDf4DCs/+1YjOith6WqDo+u2If/9+VxODiYnoioW5AtcF26dAkOhwPh4eEu88PDw2E2m92uYzabW6xvem2tJiwszGW5j48PQkJCXGrcbePaz2itF3eysrKg1+ulKSoqqtla6n56B/tj7ePxeORHfSEEsGTrSfx6+V6UVdo93RoREcmMVym2o/nz58NqtUrTmTNnPN0SdTIaHxX+L3kYXpscCz9fFf5z4hKSlvwHB4qveLo1IiKSkWyBKzQ0FCqVCqWlpS7zS0tLYTQa3a5jNBpbrG96ba3mh4Py6+vrUVZW5lLjbhvXfkZrvbij0Wig0+lcJiJ3kkf2wqdpd6N/aADOW2sw+Y0crNhdxIdeExF1UbIFLrVajVGjRmHLli3SPKfTiS1btiA+Pt7tOvHx8S71AJCdnS3VR0dHw2g0utTYbDbk5uZKNfHx8bBYLMjLy5Nqtm7dCqfTCZPJJNXs3LkTdXV1Lp8zePBgBAcHt6kXols12BiEz568Gw/GGFHnEHj+86P43ep8VPLO9EREbVLvcKKmzjtut6MQMv4n9Zo1a5Camoo33ngDcXFxeO2117B27Vp88803CA8Px9SpU9GrVy9kZWUBaLgVw7hx47Bw4UIkJSVh9erVeOmll3DgwAEMGzYMAPCXv/wFCxcuxHvvvYfo6GgsWLAAhw8fxtGjR6HVagEADzzwAEpLS7Fs2TLU1dVh2rRpGD16NFauXAkAsFqtGDx4MO6//37MnTsXBQUFePTRR/Hqq69Kt49oSy+tsdls0Ov1sFqtPNpFzRJC4N3dp5G18RjqnQIDwgKx7Fd3YkBYkKdbIyK6afUOJ6rqHKisrUdlrQNV9npU2x2ornOgpq7p1ek6r/HvH76vqXNK61xbX+cQSLg9DG+n3tWuvcvx++3TLltpxuTJk3Hx4kVkZGTAbDYjNjYWmzZtkgajFxcXQ6m8epBtzJgxWLlyJZ577jk8++yzGDhwID799FOXgDNnzhxUVlZixowZsFgsGDt2LDZt2iSFLQD48MMP8eSTT2LChAlQKpV46KGHsGTJEmm5Xq/Hl19+ibS0NIwaNQqhoaHIyMhwuVdXW3ohag8KhQLTx0ZjRG890lYewMkLFfifv+3GwoeG439GRHq6PSLqJhxOgYqaethq6lBRW98QlOxNgcn1fZXdgYraelTZ61FR60BVbX3j+8Z6ez1q6pwd0nc1j3ARj3DRjbpYXovfrTqInO8uAwB+PaYfnn3wdqh9eH0LEbWsps6B8sbAVF5TD1t1XbN/N9XZqutRXlMHW01DYJKDj1KBAI0P/NUq+KlV8PNtnNQqaHya5inh56uC9gfLtb5X32t9VfBTK6/OU1+dr/VVtWvPcvx+M3DJiIGLbka9w4n/l/0t/r79FABgZB8D/j7lTkTo/TzcGRF1hNp6ByxVdbhSZceVyjpYquwoq7I3zKu040rV1XnW6obQZKupg72+fY4oaX2VCNT4IlCjgr/aB4EaH/hrVAjQ+CBA3fTq0/CqUTX+3VDrOq/hb7VKCYVC0S69dRQGLi/DwEW34qujpXh6bT7Ka+oREqDGkl+MxNiBoZ5ui4huQJ3DicsVdlyqqEVZpb0xRF0NTVeagtU14arSfvOnyBQKIFDjA53WFzo/XwRpm/5ufNX6/GD+9X/7qnhEnYHLyzBw0a0qvlyFWR/mofCcDQoFkJ4wCGnjB0Cp9K7/WiTqSq4NURcranGxvBaXKmpxqdyOixW1uNT0vqIWV6rqWt+gG0oFYPBXw+DvixB/NQz+agT7+yI4QI3gxr+bluubQpOfLwLVPvzfh3bAwOVlGLioPdTUOZD5WSHW7G+4ke74wT3x6uRYGPzVraxJRG0lhMCVqjqYrTUoLa9p1xClUioQEqBGj4DGABVwTYDybwxQAb6N89QI8VcjSMvg5EkMXF6GgYva09r9Z7Dg0wLU1jvRy+CHf/zqTgzvbfB0W0SdXp3DiYvltThvrUGprQZmaw3MTa9Nf9tqbmgMlEqpQI8ANUIDNQgN0iA0UI2eQRr0DNQ0zAvUoGfj/GB/NcOTl2Hg8jIMXNTeCs9ZMeuDAyguq4JapcRj90QjdUw/hOu0ra9M1AVV2etdglNTqLo2XF2sqEVbf+l6BKgRptNKYenaENUwjyGqO2Dg8jIMXCQHa3Udnll7CF8da3j0lK9Kgf8eHolHx0ZjWC+9h7sjaj9Np/lKrlSjxFKFs1eqcfZKNUosja9XqmCradutDHyUCoTrtDDqtTBe8xqu1yKi8e8wnQYan/a9vQB5JwYuL8PARXIRQmBzYSne2fUd9p2++uDrH/UPwfSx/TFhSBj/65s6PadT4GJF7TUhqqoxXFVLr1VtuGIvQK1qCFB6LYw6Pxj1msZQ5SeFqx4BPCJFbcfA5WUYuKgjHDpjwTu7irDhyHk4nA3/7xwdGoBpd/fD/47qDX+1rA+UIGpWU6AqLqtC8eWqq6GqMVCds9TA7mh93FTPIA16GfzQO9gPvYL90NvQ8NrL4I9IgxZBWt8O+DbUnTBweRkGLupI5yzVeC/nNFblFkunWXRaH/zS1BepY/ryxqkki5o6B86UVTWEqrIqfH+5yuV9bSsD0ZUKIELvh16NIap38LV/+yNCr233u4gTtYaBy8swcJEnVNbW46O8s1i+uwinL1cBaBi/kjQ8AtPHRvPKRrohQghcqrA3BqhKFF+uvvp3WRVKbbUtrq9SKhBp0CIq2L8xTPlLR6p6Gfxg1Gt5o03qdBi4vAwDF3mSwymw5Vgp3tlVhNyiMml+XL8QPDo2Gj+5IxwqjmkhAOU1dQ2n+8oaxk01HZ1qOlLV2jiqII0P+vTwR5+QxumavyMNfgxU5HUYuLwMAxd1FgUlVryzqwifHzqH+sZxXn1C/DHt7n54eHQUAjUc59VVCSFgq67H2car/Eoar/RrGkt19ko1rNUt38hToQAi9X6ICvFD35AA9Onhj6gQf/RtDFUGf1+ve1YeUUsYuLwMAxd1NmZrDd7POY0Pc4ulH9kgrQ9S4vogdUw/9DJwnJe3ufbWCWevVLlc8dcUsMprW791gsHfVxo/FRXsj76NoapPiD96BfvxdgnUrTBweRkGLuqsquz1+PhACZbvKsJ3lyoBNIy1eWCYEb+M64PYPgZe3dhJNI2huvaI1A/DVVtundAjQI3ejQPRrx2c3vSeRzmJrmLg8jIMXNTZOZ0C245fwDu7irDn1GVpvkqpwKDwIMRG6TGitwEjogwYFB7EMV8yuHovqqs39vzhfalau9IPAMKCNNKVfU23UGiaIg1+DNBEN4CBy8swcJE3KTxnxYrdp7HzxEW3V575q1UY1kuP2ChDYwjTo5fBj2N3WuFwCphtNdJRKWkMlaWqzfeiUigAo057TZByPUoVafDjrROI2hEDl5dh4CJvZbbWIP+MBYfOWpBfbMGREisq3IwDCg3UuBwFG9HbAL1/97oJZZ3DCbO1RhqI7nJ0ylKN85Ya6UKF5qiUChh12qs39gz2R+9rwpVRr4Xah1f6EXUUBi4vw8BFXYXDKfDdxYqrIeyMBd+cL3cbJPqHBjSGLz1GRBlwR6TOqwdc19Y7cN5S0xikrr/Sz2yrQSt5Cr4qBSL0ruOmroarhsfP+PDWCUSdBgOXl2Hgoq6sps6BwnM2HLomhH3feKPVa/mqFBhi1CE0UI1ArS8CNT7QaX0QqPFBYONrkNYHQY3LArU+CGp89fNV3fIpSyEEauudqLY7UGmvR7XdgSppqkeV3dE4rx6VdgfKa+pxznL1KNWF8lq09r+Sah+ly7ipawej9w72Q1iQluPfiLyIHL/fHEVJRDdF66vCqL7BGNU3WJp3pdKOQ2ctOHTGKoWwsko7jpRYb+ozVEpFQwhrDGXSa2M4UypwXXCqtDtQ7fK+vtUjUK1/V6XLYHTptF/jc/1CAzV8MDIRtYiBi4jaTXCAGvcNDsN9g8MANBxdOnulGkfP22CtrkNFTT3Ka+pRUVuHitqGvxve16Oi8dVW07BMiIZTmdbqulZvzNlWah8l/NUq+Puq4K/xgb9aBT9fVcM8jQ/8fVUI0Pgg0qB1eQRNjwA1Lw4golvCwEVEslEoFIgKabiB5o0QQqDK7mgMZXUuoay8ph7ljfOFgEtYkv5uDFIBTX83hiyOkyIiT2HgIqJOR6FQIEDjgwCND8J1Wk+3Q0R0y/ife0REREQyY+AiIiIikhkDFxEREZHMGLiIiIiIZMbARURERCQzBi4iIiIimTFwEREREclMtsBVVlaGKVOmQKfTwWAwYPr06aioqGhxnZqaGqSlpaFHjx4IDAzEQw89hNLSUpea4uJiJCUlwd/fH2FhYfjDH/6A+vp6l5rt27fjzjvvhEajwYABA7BixQqX5VlZWbjrrrsQFBSEsLAwJCcn4/jx4y419913HxQKhcv0+OOP3/wOISIiom5LtsA1ZcoUFBYWIjs7G+vXr8fOnTsxY8aMFtd5+umn8fnnn2PdunXYsWMHzp07h5///OfScofDgaSkJNjtduzZswfvvfceVqxYgYyMDKmmqKgISUlJGD9+PPLz8zF79mw89thj2Lx5s1SzY8cOpKWl4euvv0Z2djbq6upw//33o7Ky0qWf3/zmNzh//rw0LVq0qJ32DhEREXUrQgZHjx4VAMS+ffukeV988YVQKBSipKTE7ToWi0X4+vqKdevWSfOOHTsmAIicnBwhhBAbN24USqVSmM1mqeYf//iH0Ol0ora2VgghxJw5c8TQoUNdtj158mSRmJjYbL8XLlwQAMSOHTukeePGjRNPPfVU27+0G1arVQAQVqv1lrZDREREHUeO329ZjnDl5OTAYDBg9OjR0ryEhAQolUrk5ua6XScvLw91dXVISEiQ5g0ZMgR9+vRBTk6OtN2YmBiEh4dLNYmJibDZbCgsLJRqrt1GU03TNtyxWq0AgJCQEJf5H374IUJDQzFs2DDMnz8fVVVVLX7v2tpa2Gw2l4mIiIhIlmcpms1mhIWFuX6Qjw9CQkJgNpubXUetVsNgMLjMDw8Pl9Yxm80uYatpedOylmpsNhuqq6vh5+fnsszpdGL27Nm4++67MWzYMGn+L3/5S/Tt2xeRkZE4fPgw5s6di+PHj+Nf//pXs987KysLf/rTn5pdTkRERN3TDQWuefPm4S9/+UuLNceOHbulhjpaWloaCgoKsGvXLpf51443i4mJQUREBCZMmIBTp07htttuc7ut+fPnIz09XXpvs9kQFRUlT+NERETkNW4ocD3zzDP49a9/3WJN//79YTQaceHCBZf59fX1KCsrg9FodLue0WiE3W6HxWJxOcpVWloqrWM0GrF3716X9ZquYry25odXNpaWlkKn0113dOvJJ5+UBvT37t27xe9lMpkAACdPnmw2cGk0Gmg0Gum9EAIAeGqRiIjIizT9bjf9jreLdhsNdo2mQfP79++X5m3evLlNg+Y/+ugjad4333zjdtB8aWmpVPPGG28InU4nampqhBANg+aHDRvmsu2UlBSXQfNOp1OkpaWJyMhI8e2337bpO+3atUsAEIcOHWpTvRBCnDlzRgDgxIkTJ06cOHnhdObMmTb/5rdGIUR7xrerHnjgAZSWlmLZsmWoq6vDtGnTMHr0aKxcuRIAUFJSggkTJuD9999HXFwcAGDWrFnYuHEjVqxYAZ1Oh9/+9rcAgD179gBouC1EbGwsIiMjsWjRIpjNZjzyyCN47LHH8NJLLwFouC3EsGHDkJaWhkcffRRbt27F7373O2zYsAGJiYkAgCeeeAIrV67EZ599hsGDB0s96/V6+Pn54dSpU1i5ciUefPBB9OjRA4cPH8bTTz+N3r17Y8eOHW3eB06nE+fOnUNQUBAUCsWt71Qv1XRq9cyZM9DpdJ5up8vgfm1/3Kfy4H6VB/erPJr269GjRzF48GAole10fWG7RbcfuHz5skhJSRGBgYFCp9OJadOmifLycml5UVGRACC2bdsmzauurhZPPPGECA4OFv7+/uJnP/uZOH/+vMt2T58+LR544AHh5+cnQkNDxTPPPCPq6upcarZt2yZiY2OFWq0W/fv3F8uXL3dZjmaSbFNdcXGxuPfee0VISIjQaDRiwIAB4g9/+ANv73CTeHsMeXC/tj/uU3lwv8qD+1Uecu1X2Y5wETWx2WzQ6/WwWq38r7B2xP3a/rhP5cH9Kg/uV3nItV/5LEUiIiIimTFwkew0Gg0yMzNdruCkW8f92v64T+XB/SoP7ld5yLVfeUqRiIiISGY8wkVEREQkMwYuIiIiIpkxcBERERHJjIGLiIiISGYMXHRLsrKycNdddyEoKAhhYWFITk7G8ePH27z+6tWroVAokJycLF+TXuhm96vFYkFaWhoiIiKg0WgwaNAgbNy4sQM67vxudp++9tprGDx4MPz8/BAVFYWnn34aNTU1HdCxd/jHP/6B4cOHQ6fTQafTIT4+Hl988UWL66xbtw5DhgyBVqtFTEwM/426caP79a233sI999yD4OBgBAcHIyEh4bpnD9PN/Xttcqu/VwxcdEt27NiBtLQ0fP3118jOzkZdXR3uv/9+VFZWtrru6dOn8fvf/x733HNPB3TqXW5mv9rtdvzkJz/B6dOn8dFHH+H48eN466230KtXrw7svPO6mX26cuVKzJs3D5mZmTh27BjeeecdrFmzBs8++2wHdt659e7dGwsXLkReXh7279+PH//4x5g4cSIKCwvd1u/ZswcpKSmYPn06Dh48iOTkZCQnJ6OgoKCDO+/cbnS/bt++HSkpKdi2bRtycnIQFRWF+++/HyUlJR3ceed2o/u1Sbv8XrXrfeup27tw4YIAIHbs2NFiXX19vRgzZox4++23RWpqqpg4cWLHNOil2rJf//GPf4j+/fsLu93egZ15r7bs07S0NPHjH//YZV56erq4++675W7PqwUHB4u3337b7bJJkyaJpKQkl3kmk0nMnDmzI1rzai3t1x+qr68XQUFB4r333pO5K+/X2n5tr98rHuGidmW1WgEAISEhLdb9+c9/RlhYGKZPn94RbXm9tuzXf//734iPj0daWhrCw8MxbNgwvPTSS3A4HB3Vpldpyz4dM2YM8vLypFMz3333HTZu3IgHH3ywQ3r0Ng6HA6tXr0ZlZSXi4+Pd1uTk5CAhIcFlXmJiInJycjqiRa/Ulv36Q1VVVairq2v1f4u7s7bu1/b6vfK5pbWJruF0OjF79mzcfffdGDZsWLN1u3btwjvvvIP8/PyOa86LtXW/fvfdd9i6dSumTJmCjRs34uTJk3jiiSdQV1eHzMzMDuy482vrPv3lL3+JS5cuYezYsRBCoL6+Ho8//jhPKf7AkSNHEB8fj5qaGgQGBuKTTz7BHXfc4bbWbDYjPDzcZV54eDjMZnNHtOpVbmS//tDcuXMRGRl5XbilG9uv7fl7xcBF7SYtLQ0FBQXYtWtXszXl5eV45JFH8NZbbyE0NLQDu/NebdmvQEOICAsLw5tvvgmVSoVRo0ahpKQEixcvZuD6gbbu0+3bt+Oll17C3//+d5hMJpw8eRJPPfUU/u///g8LFizooG47v8GDByM/Px9WqxUfffQRUlNTsWPHjjaHA3LvZvfrwoULsXr1amzfvh1arbaDuvUebd2v7f57dXNnPIlcpaWlid69e4vvvvuuxbqDBw8KAEKlUkmTQqEQCoVCqFQqcfLkyQ7q2Du0db8KIcS9994rJkyY4DJv48aNAoCora2Vq0WvcyP7dOzYseL3v/+9y7x//vOfws/PTzgcDrla9HoTJkwQM2bMcLssKipKvPrqqy7zMjIyxPDhwzugM+/W0n5tsnjxYqHX68W+ffs6qCvv19x+be/fKx7holsihMBvf/tbfPLJJ9i+fTuio6NbrB8yZAiOHDniMu+5555DeXk5/vrXvyIqKkrOdr3Gje5XALj77ruxcuVKOJ1OKJUNwzO//fZbREREQK1Wy91yp3cz+7Sqqkral01UKpW0PXLP6XSitrbW7bL4+Hhs2bIFs2fPluZlZ2e3eWxSd9bSfgWARYsW4cUXX8TmzZsxevToDuzMuzW3X9v996o90iF1X7NmzRJ6vV5s375dnD9/XpqqqqqkmkceeUTMmzev2W3wKsXr3cx+LS4uFkFBQeLJJ58Ux48fF+vXrxdhYWHihRde8MRX6HRuZp9mZmaKoKAgsWrVKvHdd9+JL7/8Utx2221i0qRJnvgKndK8efPEjh07RFFRkTh8+LCYN2+eUCgU4ssvvxRCXL9Pd+/eLXx8fMTLL78sjh07JjIzM4Wvr684cuSIp75Cp3Sj+3XhwoVCrVaLjz76yOXfd3l5uae+Qqd0o/v1h27l94qBi24JALfT8uXLpZpx48aJ1NTUZrfBwHW9m92ve/bsESaTSWg0GtG/f3/x4osvivr6+o5tvpO6mX1aV1cnnn/+eXHbbbcJrVYroqKixBNPPCGuXLnS4f13Vo8++qjo27evUKvVomfPnmLChAnSj5cQ7v+drl27VgwaNEio1WoxdOhQsWHDhg7uuvO70f3at29ft/++MzMzO775Tuxm/r1e61Z+rxRC8Lg4ERERkZx4Hy4iIiIimTFwEREREcmMgYuIiIhIZgxcRERERDJj4CIiIiKSGQMXERERkcwYuIiIiIhkxsBFREREJDMGLiIiIiKZMXARERERyYyBi4iIiEhmDFxEREREMvv/Z9jHtFpCnAUAAAAASUVORK5CYII=",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "en_Ne2=GR.energy_2(\"Ne\")\n",
    "dist_Ne2=GR.distanze(\"Ne\")\n",
    "en_Ne2_corr=[x-2*en_Ne for x in en_Ne2]\n",
    "\n",
    "plt.plot(dist_Ne2, en_Ne2_corr)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "ab84004c-12a9-42ad-ad91-8e9ee9873124",
   "metadata": {},
   "outputs": [],
   "source": [
    "def V_morse(r, D, a, r0):\n",
    "    er=np.exp(-a*(r-r0))\n",
    "    return D*(er*er-2*er)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "id": "f58b1bf1-ec0d-43a2-83e2-758a73ef1de1",
   "metadata": {},
   "outputs": [],
   "source": [
    "p_morse, covar = opt.curve_fit(V_morse, dist_Ne2, en_Ne2_corr, p0=[0.0001, 1.0, 2.7])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "id": "6a31f51b-7f34-4e35-98a6-c3721977bc76",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([1.48435292e-04, 3.17996133e+00, 2.78537637e+00])"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "p_morse"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "id": "bce601b8-487d-4b01-8cc7-7e006f985a8e",
   "metadata": {},
   "outputs": [],
   "source": [
    "en_sim=V_morse(dist_Ne2, *p_morse)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "id": "4db5ec80-f6fe-40db-84cf-8eb38b6da95e",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAlwAAAGdCAYAAADHQK08AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAABov0lEQVR4nO3de1hUdf4H8PfMwMxwmwFEGFBUzFteMQ1CLS1JKreibbNYV81sbQ3bjMq0TanWorR+tW4Xu6q7ZV66b5rm4m1TQkURFe+ieBsQkRnuc/v+/hg4OjrclGEYeL+e5zzDnPM5Zz5zVObtOWe+RyaEECAiIiIil5G7uwEiIiKito6Bi4iIiMjFGLiIiIiIXIyBi4iIiMjFGLiIiIiIXIyBi4iIiMjFGLiIiIiIXIyBi4iIiMjFvNzdQFtms9lw9uxZBAQEQCaTubsdIiIiagQhBEpLSxEREQG5vHmOTTFwudDZs2cRGRnp7jaIiIjoGpw6dQqdO3dulm0xcLlQQEAAAPsfmEajcXM3RERE1BhGoxGRkZHS53hzYOByodrTiBqNhoGLiIjIwzTn5UC8aJ6IiIjIxRi4iIiIiFyMgYuIiIjIxXgNFxERkYtYrVaYzWZ3t0FXUCgU8PLyatEhmxi4iIiIXKCsrAynT5+GEMLdrZATvr6+CA8Ph1KpbJHXY+AiIiJqZlarFadPn4avry86duzIwa9bESEETCYTzp8/j7y8PPTs2bPZBjetDwMXERFRMzObzRBCoGPHjvDx8XF3O3QFHx8feHt74+TJkzCZTFCr1S5/TV40T0RE5CI8stV6tcRRLYfXa9FXIyIiImqHGLiIiIiIXIyBi4iIiFzmxIkTkMlkyM7OBgBs2rQJMpkMJSUlbu2rpTFwERERUYsZNmwYzp07B61W6+5WWhS/peiJCg8Ae74CfIKBETPc3Q0REbVRJpOp2cepUiqV0Ol0zbpNT8AjXJ6o5BSw9R/20EVERK2eEAIVJotbpqYMvDpq1ChMnz4dM2bMQEhICBISErBv3z7cfffd8Pf3R1hYGCZMmICioiJpnbVr12LEiBEIDAxEhw4d8Lvf/Q7Hjh2r8zXa6ylFHuHyRLoB9seiw4CpAlD6urcfIiKqV6XZir5z17nltXNfTYCvsvEf90uXLsW0adOwdetWlJSU4I477sDjjz+Od955B5WVlXjhhRcwbtw4bNiwAQBQXl6OlJQUDBw4EGVlZZg7dy4eeOABZGdnt/jQC60ZA5cnCtABfh2B8vNAYS7Qeai7OyIiojaiZ8+emD9/PgBg3rx5GDx4MF5//XVp+eeff47IyEgcPnwYvXr1woMPPuiw/ueff46OHTsiNzcX/fv3b9HeWzMGLk8kkwG6gcCxdECfw8BFRNTK+XgrkPtqgtteuymGDBki/bxnzx5s3LgR/v7+V9UdO3YMvXr1wpEjRzB37lxkZmaiqKgINpsNAJCfn8/AdRkGLg9UWmVGieoGRCIdOJfj7naIiKgBMpmsSaf13MnPz0/6uaysDPfeey/efPPNq+rCw8MBAPfeey+6du2KTz75BBEREbDZbOjfvz9MJlOL9ewJPONPnxxsOFiI9dkqvKeE/QgXERGRC9x000345ptv0K1bN3h5XR0ZLly4gEOHDuGTTz7BrbfeCgD49ddfW7pNj8Cr2TzQwM6B2C+6AQBEwX7AanFvQ0RE1CYlJyejuLgYSUlJ2LFjB44dO4Z169Zh8uTJsFqtCAoKQocOHfDxxx/j6NGj2LBhA1JSUtzddqvEwOWBunXwRbGqE8qFCjJLFXDhqLtbIiKiNigiIgJbt26F1WrFmDFjMGDAAMyYMQOBgYGQy+WQy+VYvnw5srKy0L9/fzzzzDNYsGCBu9tulXhK0QPJZDIM6ByMA/ldMVR22H5aMbSPu9siIiIPt2nTpqvm9ezZE99++22d68THxyM3N9dh3uVjf3Xr1s3h+ahRo5o0NlhbwSNcHmpgZy3227ran/A6LiIiolaNgctDDYq8dB0Xv6lIRETUujFweahBnQORW3OES+hzgHZ4eJaIiMhTMHB5KJ1WjYt+PWAWCsgqLwKG0+5uiYiIiOrAwOXBbozsiKOik/2Jfq97myEiIqI6MXB5sEGdtcgVvHCeiIiotXN54Hr//ffRrVs3qNVqxMbGYvv27fXWr1q1Cn369IFarcaAAQOwZs0ah+VCCMydOxfh4eHw8fFBfHw8jhw54lDz2muvYdiwYfD19UVgYOBVr7FkyRLIZDKnU2FhIQD7V2OdLdfr9de3Q5rRwMhA7Ld1sz/hhfNEREStlksD14oVK5CSkoLU1FTs2rULgwYNQkJCghRqrrRt2zYkJSVhypQp2L17NxITE5GYmIh9+/ZJNfPnz8fChQuxaNEiZGZmws/PDwkJCaiqqpJqTCYTHnroIUybNs3p6zz88MM4d+6cw5SQkICRI0ciNDTUofbQoUMOdVcud6eBnbRS4LIxcBEREbVewoViYmJEcnKy9NxqtYqIiAiRlpbmtH7cuHFi7NixDvNiY2PFE088IYQQwmazCZ1OJxYsWCAtLykpESqVSnz11VdXbW/x4sVCq9U22GdhYaHw9vYW//rXv6R5GzduFADExYsXG1y/LgaDQQAQBoPhmrfRkLve+I8QqRr7VH7BZa9DRESNV1lZKXJzc0VlZaW7W2mVRo4cKZ5++ulG11/5eZ6amioGDRp0XT3U92fkis9vlx3hMplMyMrKQnx8vDRPLpcjPj4eGRkZTtfJyMhwqAeAhIQEqT4vLw96vd6hRqvVIjY2ts5tNsa//vUv+Pr64g9/+MNVy6KjoxEeHo4777wTW7durXc71dXVMBqNDpOr3RAZgXxbR/sTXjhPRETUKrkscBUVFcFqtSIsLMxhflhYWJ3XQen1+nrrax+bss3G+Oyzz/DHP/4RPj4+0rzw8HAsWrQI33zzDb755htERkZi1KhR2LVrV53bSUtLg1arlabIyMhr7qmxBl12I2sGLiIiotap3X9LMSMjAwcOHMCUKVMc5vfu3RtPPPEEhgwZgmHDhuHzzz/HsGHD8M4779S5rdmzZ8NgMEjTqVOnXN1+zS1+utmf8JuKRER0HUaNGoWnnnoKM2bMQFBQEMLCwvDJJ5+gvLwckydPRkBAAHr06IGff/5ZWmfz5s2IiYmBSqVCeHg4Zs2aBYvFIi0vLy/HxIkT4e/vj/DwcLz99ttXvW51dTWee+45dOrUCX5+foiNjXV6X8e62Gw2vPrqq+jcuTNUKhWio6Oxdu3a69oXzc1lgSskJAQKhQIFBQUO8wsKCqDT6Zyuo9Pp6q2vfWzKNhvy6aefIjo6GkOGDGmwNiYmBkePHq1zuUqlgkajcZhcrX8nLQ7CPjSE5Uy2y1+PiIiugRCAqdw9UxPvRLJ06VKEhIRg+/bteOqppzBt2jQ89NBDGDZsGHbt2oUxY8ZgwoQJqKiowJkzZ3DPPffg5ptvxp49e/Dhhx/is88+w7x586TtPf/889i8eTN++OEH/PLLL9i0adNVZ4umT5+OjIwMLF++HDk5OXjooYdw1113XTUKQV3+8Y9/4O2338Zbb72FnJwcJCQk4L777mv0+i3By1UbViqVGDJkCNLT05GYmAjAnkDT09Mxffp0p+vExcUhPT0dM2bMkOatX78ecXFxAICoqCjodDqkp6cjOjoaAGA0GpGZmVnnNxLrU1ZWhpUrVyItLa1R9dnZ2QgPD2/y67iSn8oLFcH9gVJAXnwUMFcC3j4Nr0hERC3HXAG8HuGe137xLKD0a3T5oEGD8NJLLwGwn7l54403EBISgj//+c8AgLlz5+LDDz9ETk4O/vOf/yAyMhLvvfceZDIZ+vTpg7Nnz+KFF17A3LlzUVFRgc8++wxffPEFRo8eDcAe6Dp37iy9Xn5+PhYvXoz8/HxERNj30XPPPYe1a9di8eLFeP311xvs+a233sILL7yARx55BADw5ptvYuPGjXj33Xfx/vvvN/q9u5LLAhcApKSkYNKkSRg6dChiYmLw7rvvSoclAWDixIno1KmTFHiefvppjBw5Em+//TbGjh2L5cuXY+fOnfj4448BADKZDDNmzMC8efPQs2dPREVFYc6cOYiIiJBCHWD/wysuLkZ+fj6sViuys7MBAD169IC/v79Ut2LFClgsFvzpT3+6qvd3330XUVFR6NevH6qqqvDpp59iw4YN+OWXX1y0t65dp8goFO3XIARGoDAX6NTw0ToiIiJnBg4cKP2sUCjQoUMHDBgwQJpXex11YWEhDhw4gLi4OMhkMmn58OHDUVZWhtOnT+PixYswmUyIjY2VlgcHB6N3797S871798JqtaJXr14OfVRXV6NDhw4N9ms0GnH27FkMHz7cYf7w4cOxZ8+eRr5r13Np4Hr44Ydx/vx5zJ07F3q9XjqnWvuHlZ+fD7n80lnNYcOGYdmyZXjppZfw4osvomfPnvj+++/Rv39/qWbmzJkoLy/H1KlTUVJSghEjRmDt2rVQq9VSzdy5c7F06VLp+eDBgwEAGzduxKhRo6T5n332GX7/+987HRzVZDLh2WefxZkzZ+Dr64uBAwfiv//9L26//fbm2j3NZmCXIOTu7YrbFHvtA6AycBERtS7evvYjTe567aaUe3s7PJfJZA7zasOVzWa7/t5gP9ukUCiQlZUFhULhsOzygySezqWBC7Cfl63rFKKzC+IeeughPPTQQ3VuTyaT4dVXX8Wrr75aZ82SJUuwZMmSBnvbtm1bnctmzpyJmTNnNriN1mBQZy22iW64DXshzuVA1vAqRETUkmSyJp3W8xQ33ngjvvnmGwghpCC2detWBAQEoHPnzggODoa3tzcyMzPRpUsXAMDFixdx+PBhjBw5EoD9oIjVakVhYSFuvfXWJveg0WgQERGBrVu3Stus7SMmJqYZ3mXzaPffUmwL+ug0OIQoAIDpTOs5fEpERG3bk08+iVOnTuGpp57CwYMH8cMPPyA1NRUpKSmQy+Xw9/fHlClT8Pzzz2PDhg3Yt28fHn30UYezW7169cL48eMxceJEfPvtt8jLy8P27duRlpaG1atXN6qP559/Hm+++SZWrFiBQ4cOYdasWcjOzsbTTz/tqrfeZC4/wkWup/SSozqkH1ACKM7nAjYrIFc0uB4REdH16NSpE9asWYPnn38egwYNQnBwMKZMmSJddA8ACxYsQFlZGe69914EBATg2WefhcFgcNjO4sWLMW/ePOlSnpCQENxyyy343e9+16g+/vrXv8JgMODZZ59FYWEh+vbtix9//BE9e/Zs1vd7PWRCNPH7otRoRqMRWq0WBoPB5UNEzP1uD2Zl3wlfWTWQvB3o2LvhlYiIyCWqqqqQl5eHqKgoh2uMqfWo78/IFZ/fPKXYRgyIDMYBYT8/Dt7ImoiIqFVh4GojBkUGSiPO2xi4iIiIWhUGrjbiho7+OCq3Xzhfmb/bzd0QERHR5Ri42giFXAZTR/t4ZV6Fe5t8KwciIiJyHQauNiSw20BYhBwqcwlgdNMAe0RERHQVBq42pH/XMBwVnexP9LyOi4jI3TgQQOvV0n82DFxtyKDOgdgvugIALGey3dsMEVE7VnuLGpPJ5OZOqC4VFRUArr6Vkatw4NM2pHOQD1Z43QCIX1F2cjcC3d0QEVE75eXlBV9fX5w/fx7e3t4OI6uTewkhUFFRgcLCQgQGBl51/0ZXYeBqQ2QyGcwdBwCFNRfOExGRW8hkMoSHhyMvLw8nT550dzvkRGBgIHQ6XYu9HgNXG6PpFg0UAv6VZ4HKi4BPkLtbIiJql5RKJXr27MnTiq2Qt7d3ix3ZqsXA1cb07tYFp37riEj5eUC/F4i6zd0tERG1W3K5nLf2IQC8aL7NGRipxX7RDQBQfXqPe5shIiIiAAxcbU5ogBqnlDcAAIx5WW7uhoiIiAAGrjbJ1LEfAEBRwAvniYiIWgMGrjbIr9sQAIC2/DhgrnJzN0RERMTA1Qb16N4LxcIfCtiAwlx3t0NERNTuMXC1QQMiA7Hf1g0AUH5yt3ubISIiIgautkjr440z6p4AgBJeOE9EROR2DFxtlKljfwCAXM8L54mIiNyNgauN8us6GAAQXHYYsFnd3A0REVH7xsDVRnXrNRAVQgWVqIK4cNTd7RAREbVrDFxtVN9OwTgougAASvJ2ubkbIiKi9o2Bq43yUSpwTt0DAHDx6E43d0NERNS+MXC1YdU1F87LCnLc3AkREVH7xsDVhvl2uQkAEFx6CBDCzd0QERG1XwxcbVhknyGwCDm0NgNshrPuboeIiKjdYuBqw3p17ojjiAAAFB7Z7uZuiIiI2i+XB673338f3bp1g1qtRmxsLLZvr/+Df9WqVejTpw/UajUGDBiANWvWOCwXQmDu3LkIDw+Hj48P4uPjceTIEYea1157DcOGDYOvry8CAwOdvo5MJrtqWr58uUPNpk2bcNNNN0GlUqFHjx5YsmRJk9+/O3kr5DhbM+J88TF+U5GIiMhdXBq4VqxYgZSUFKSmpmLXrl0YNGgQEhISUFhY6LR+27ZtSEpKwpQpU7B7924kJiYiMTER+/btk2rmz5+PhQsXYtGiRcjMzISfnx8SEhJQVVUl1ZhMJjz00EOYNm1avf0tXrwY586dk6bExERpWV5eHsaOHYvbb78d2dnZmDFjBh5//HGsW7fu+nZKC6sOsV84Dz0vnCciInIb4UIxMTEiOTlZem61WkVERIRIS0tzWj9u3DgxduxYh3mxsbHiiSeeEEIIYbPZhE6nEwsWLJCWl5SUCJVKJb766qurtrd48WKh1WqdvhYA8d1339XZ+8yZM0W/fv0c5j388MMiISGhznWuZDAYBABhMBgavU5z27JulRCpGqF/tZfbeiAiIvIkrvj8dtkRLpPJhKysLMTHx0vz5HI54uPjkZGR4XSdjIwMh3oASEhIkOrz8vKg1+sdarRaLWJjY+vcZn2Sk5MREhKCmJgYfP755xCXfZOvoV48RecbYwEAYVY9LOUX3dwNERFR++Tlqg0XFRXBarUiLCzMYX5YWBgOHjzodB29Xu+0Xq/XS8tr59VV01ivvvoq7rjjDvj6+uKXX37Bk08+ibKyMvz1r3+ttxej0YjKykr4+Phctc3q6mpUV1dLz41GY5N6coWunTrjjOiITrLzOH1gO7oNTXB3S0RERO2OywJXazdnzhzp58GDB6O8vBwLFiyQAte1SEtLwyuvvNIc7TUbuVyGcz490KnqPIqP7WTgIiIicgOXnVIMCQmBQqFAQUGBw/yCggLodDqn6+h0unrrax+bss3Gio2NxenTp6UjVHX1otFonB7dAoDZs2fDYDBI06lTp66rp+ZS2YEXzhMREbmTywKXUqnEkCFDkJ6eLs2z2WxIT09HXFyc03Xi4uIc6gFg/fr1Un1UVBR0Op1DjdFoRGZmZp3bbKzs7GwEBQVBpVI1qhdnVCoVNBqNw9Qa+HYZDAAINDo/lUtERESu5dJTiikpKZg0aRKGDh2KmJgYvPvuuygvL8fkyZMBABMnTkSnTp2QlpYGAHj66acxcuRIvP322xg7diyWL1+OnTt34uOPPwZgHztrxowZmDdvHnr27ImoqCjMmTMHERERDkM65Ofno7i4GPn5+bBarcjOzgYA9OjRA/7+/vjPf/6DgoIC3HLLLVCr1Vi/fj1ef/11PPfcc9I2/vKXv+C9997DzJkz8dhjj2HDhg1YuXIlVq9e7cpd5hKdbowBMoBIyylUVVZA7ePr7paIiIjal2b7vmMd/vnPf4ouXboIpVIpYmJixG+//SYtGzlypJg0aZJD/cqVK0WvXr2EUqkU/fr1E6tXr3ZYbrPZxJw5c0RYWJhQqVRi9OjR4tChQw41kyZNEgCumjZu3CiEEOLnn38W0dHRwt/fX/j5+YlBgwaJRYsWCavV6rCdjRs3iujoaKFUKkX37t3F4sWLm/TeW8OwEEIIYbNaxcXUCCFSNSJ352a39kJERNTaueLzWyYE72rsKkajEVqtFgaDwe2nF3PfGIW+Vbvxa99UjBiX4tZeiIiIWjNXfH7zXortREWHfgAAcY4XzhMREbU0Bq52wjcyGgAQZDjg3kaIiIjaIQaudiKij33E+ShrHowVVQ1UExERUXNi4GonArv0QxWU8JNV4+jBve5uh4iIqF1h4Gov5AroVd0BAEVHd7i5GSIiovaFgasdkS6cP7vHzZ0QERG1Lwxc7Ygq0j7ivNbAEeeJiIhaEgNXOxLe+2YAwA224zhv5IXzRERELYWBqx3x7TwQVsjRUWbEwaOH3d0OERFRu8HA1Z4ofXFe1QUAUHQky83NEBERtR8MXO1MeXBfAIDtXLZ7GyEiImpHGLjaGXXnaACAtuQgeBtNIiKilsHA1c507Gm/cL6n7ThOX6x0czdERETtAwNXO6OsOcLVVV6I3LzT7m2GiIionWDgam98g1HiHQoAOH9kp5ubISIiah8YuNqh8mD7iPNWjjhPRETUIhi42iFVzWlFjeEgrDZeOE9ERORqDFztUFD3IQCA3iIPx8+XubkbIiKito+Bqx1SdBoEAOgpO429J8+7uRsiIqK2j4GrPdJGolIRAG+ZFQXHdru7GyIiojaPgas9kslQFmQfcd58JsfNzRAREbV9DFztVO14XEGGAzBZbO5thoiIqI1j4GqnNFE3AQD6yPJwSF/q5m6IiIjaNgaudkoWbr9w/kZZPvacKnZzN0RERG0bA1d71aEnzDIV/GVVOHt8v7u7ISIiatMYuNorhRcqAnsBAMxnOOI8ERGRKzFwtWPetRfOGw+iwmRxbzNERERtGANXO+bbJRoA0Fd2EvvPGt3bDBERURvGwNWe6ewXzveTn8CeUyXu7YWIiKgNY+Bqz8L6wQY5OsoMyDtx3N3dEBERtVkMXO2Z0heVmigAgPl0tnt7ISIiasNcHrjef/99dOvWDWq1GrGxsdi+fXu99atWrUKfPn2gVqsxYMAArFmzxmG5EAJz585FeHg4fHx8EB8fjyNHjjjUvPbaaxg2bBh8fX0RGBh41Wvs2bMHSUlJiIyMhI+PD2688Ub84x//cKjZtGkTZDLZVZNer7+2HdFKedfcyDqk7BBKKkxu7oaIiKhtcmngWrFiBVJSUpCamopdu3Zh0KBBSEhIQGFhodP6bdu2ISkpCVOmTMHu3buRmJiIxMRE7Nu3T6qZP38+Fi5ciEWLFiEzMxN+fn5ISEhAVVWVVGMymfDQQw9h2rRpTl8nKysLoaGh+OKLL7B//3787W9/w+zZs/Hee+9dVXvo0CGcO3dOmkJDQ69zr7Qutbf46Sc/gZzTBvc2Q0RE1FYJF4qJiRHJycnSc6vVKiIiIkRaWprT+nHjxomxY8c6zIuNjRVPPPGEEEIIm80mdDqdWLBggbS8pKREqFQq8dVXX121vcWLFwutVtuoXp988klx++23S883btwoAIiLFy82an1nDAaDACAMBsM1b8Pljm4QIlUjjs/pKf6Zftjd3RAREbmdKz6/XXaEy2QyISsrC/Hx8dI8uVyO+Ph4ZGRkOF0nIyPDoR4AEhISpPq8vDzo9XqHGq1Wi9jY2Dq32VgGgwHBwcFXzY+OjkZ4eDjuvPNObN26td5tVFdXw2g0Okytnm4gACBKXoBD+Wfd3AwREVHb5LLAVVRUBKvVirCwMIf5YWFhdV4Hpdfr662vfWzKNhtj27ZtWLFiBaZOnSrNCw8Px6JFi/DNN9/gm2++QWRkJEaNGoVdu3bVuZ20tDRotVppioyMvOaeWoxfB1T7hgMATKc54jwREZEreLm7AXfbt28f7r//fqSmpmLMmDHS/N69e6N3797S82HDhuHYsWN455138O9//9vptmbPno2UlBTpudFo9IjQpeg0CDhyDrrKoygwViFMo3Z3S0RERG2Ky45whYSEQKFQoKCgwGF+QUEBdDqd03V0Ol299bWPTdlmfXJzczF69GhMnToVL730UoP1MTExOHr0aJ3LVSoVNBqNw+QJvMLtpxX7yTgAKhERkSu4LHAplUoMGTIE6enp0jybzYb09HTExcU5XScuLs6hHgDWr18v1UdFRUGn0znUGI1GZGZm1rnNuuzfvx+33347Jk2ahNdee61R62RnZyM8PLxJr+MRagMXv6lIRETkEi49pZiSkoJJkyZh6NChiImJwbvvvovy8nJMnjwZADBx4kR06tQJaWlpAICnn34aI0eOxNtvv42xY8di+fLl2LlzJz7++GMAgEwmw4wZMzBv3jz07NkTUVFRmDNnDiIiIpCYmCi9bn5+PoqLi5Gfnw+r1Yrs7GwAQI8ePeDv7499+/bhjjvuQEJCAlJSUqTrvxQKBTp27AgAePfddxEVFYV+/fqhqqoKn376KTZs2IBffvnFlbvMPWounO8pO40Fp84D6F1/PRERETWJSwPXww8/jPPnz2Pu3LnQ6/WIjo7G2rVrpYve8/PzIZdfOsg2bNgwLFu2DC+99BJefPFF9OzZE99//z369+8v1cycORPl5eWYOnUqSkpKMGLECKxduxZq9aXrjubOnYulS5dKzwcPHgwA2LhxI0aNGoWvv/4a58+fxxdffIEvvvhCquvatStOnDgBwP4ty2effRZnzpyBr68vBg4ciP/+97+4/fbbXbKv3CqwC6xKDZQmIypO74cQwyGTydzdFRERUZshE0IIdzfRVhmNRmi1WhgMhlZ/PZdt8VjIT/6K581TkfzMy+gW4ufuloiIiNzCFZ/fvJciAQDk4fZb/PSVncT2vGI3d0NERNS2MHCRXc2F833lJ7H5yHk3N0NERNS2MHCRnW4AAPsRrq2HC2G18UwzERFRc2HgIruQXhAKFQJkldBWn8Ge0yXu7oiIiKjNYOAiO4U3ZDr7t0GHyA5j8yGeViQiImouDFx0SdRIAMCtir3YfJiBi4iIqLkwcNElN9wBABgh34uc08W4WG5yc0NERERtAwMXXRIZCyj90VFmxI3Ix/+OFrm7IyIiojaBgYsu8VIC3W4FANwmz+F1XERERM2EgYsc1ZxWvFWeg82Hz8PG4SGIiIiuGwMXOeoxGgAwVH4IFWUGHNAb3dwQERGR52PgIkfB3YHALlDKrIiVH+C3FYmIiJoBAxc5ksmAG+xHuW6T52ALAxcREdF1Y+Ciq/W4FLh2nriIsmqLmxsiIiLybAxcdLWo2wCZAjfIz0EnCrGNw0MQERFdFwYuuppaC3S+GcClbysSERHRtWPgIuek4SHst/kRgsNDEBERXSsGLnKu5jquEfJ9OHexDMeLyt3cEBERkedi4CLnIgYD6kBoZBUYJDvGUeeJiIiuAwMXOSdXAN1HAQBuU/A6LiIiouvBwEV1u2x4iMy8C6gyW93cEBERkWdi4KK61Vw4P0h+HEqzEdvzit3cEBERkWdi4KK6aTsDIb2hgA3D5ft5WpGIiOgaMXBR/WpOK3I8LiIiomvHwEX1qzmtOFKxF0cLS3H6YoWbGyIiIvI8DFxUv67DAYUSnWRFuEF2FlsO8zY/RERETcXARfVT+gJd4gDUjjpf6OaGiIiIPA8DFzXssuEhth69ALPV5uaGiIiIPAsDFzWs5jquOMUBmKorsevkRTc3RERE5FkYuKhhYf0B/zD4oBpD5Iex5Qi/rUhERNQULg9c77//Prp16wa1Wo3Y2Fhs37693vpVq1ahT58+UKvVGDBgANasWeOwXAiBuXPnIjw8HD4+PoiPj8eRI0ccal577TUMGzYMvr6+CAwMdPo6+fn5GDt2LHx9fREaGornn38eFovFoWbTpk246aaboFKp0KNHDyxZsqTJ779NkMkufVuRw0MQERE1mUsD14oVK5CSkoLU1FTs2rULgwYNQkJCAgoLnV94vW3bNiQlJWHKlCnYvXs3EhMTkZiYiH379kk18+fPx8KFC7Fo0SJkZmbCz88PCQkJqKqqkmpMJhMeeughTJs2zenrWK1WjB07FiaTCdu2bcPSpUuxZMkSzJ07V6rJy8vD2LFjcfvttyM7OxszZszA448/jnXr1jXT3vEwNYHrVnkO9p0x4nxptZsbIiIi8iDChWJiYkRycrL03Gq1ioiICJGWlua0fty4cWLs2LEO82JjY8UTTzwhhBDCZrMJnU4nFixYIC0vKSkRKpVKfPXVV1dtb/HixUKr1V41f82aNUIulwu9Xi/N+/DDD4VGoxHV1dVCCCFmzpwp+vXr57Deww8/LBISEhp415cYDAYBQBgMhkav02qVFgqRqhEiVSOGvPCl+CbrlLs7IiIicglXfH677AiXyWRCVlYW4uPjpXlyuRzx8fHIyMhwuk5GRoZDPQAkJCRI9Xl5edDr9Q41Wq0WsbGxdW6zrtcZMGAAwsLCHF7HaDRi//79jerFmerqahiNRoepzfDvCOgGAgBGyPfytCIREVETuCxwFRUVwWq1OoQaAAgLC4Ner3e6jl6vr7e+9rEp22zK61z+GnXVGI1GVFZWOt1uWloatFqtNEVGRja6J49QOzyEIgdbDp+H1Sbc3BAREZFn4LcUm9Hs2bNhMBik6dSpU+5uqXndUDse116UVFRj3xmDmxsiIiLyDC4LXCEhIVAoFCgoKHCYX1BQAJ1O53QdnU5Xb33tY1O22ZTXufw16qrRaDTw8fFxul2VSgWNRuMwtSmRsYC3H0JkBtwoy+dpRSIiokZyWeBSKpUYMmQI0tPTpXk2mw3p6emIi4tzuk5cXJxDPQCsX79eqo+KioJOp3OoMRqNyMzMrHObdb3O3r17Hb4tuX79emg0GvTt27dRvbRLXkog6lYA9lHntzBwERERNYpLTymmpKTgk08+wdKlS3HgwAFMmzYN5eXlmDx5MgBg4sSJmD17tlT/9NNPY+3atXj77bdx8OBBvPzyy9i5cyemT58OAJDJZJgxYwbmzZuHH3/8EXv37sXEiRMRERGBxMREaTv5+fnIzs5Gfn4+rFYrsrOzkZ2djbKyMgDAmDFj0LdvX0yYMAF79uzBunXr8NJLLyE5ORkqlQoA8Je//AXHjx/HzJkzcfDgQXzwwQdYuXIlnnnmGVfustZPGh5iL3blX4ShwuzmhoiIiDxAs33fsQ7//Oc/RZcuXYRSqRQxMTHit99+k5aNHDlSTJo0yaF+5cqVolevXkKpVIp+/fqJ1atXOyy32Wxizpw5IiwsTKhUKjF69Ghx6NAhh5pJkyYJAFdNGzdulGpOnDgh7r77buHj4yNCQkLEs88+K8xms8N2Nm7cKKKjo4VSqRTdu3cXixcvbtJ7b1PDQtQ6f0SIVI2oTg0WfV74WqzOOevujoiIiJqVKz6/ZUIIftXMRYxGI7RaLQwGQ9u5nksI4B8DgZJ8PGp6HmE33Yc3/zDQ3V0RERE1G1d8fvNbitQ0Tm7zw8xORERUPwYuarqa4SFGKvZCb6zC4YIyNzdERETUujFwUdNF3QbIFOguO4tOOI/Nh53fG5OIiIjsGLio6XwCgc5DAQC3KnibHyIiooYwcNG1qTmteKs8BzvyLqLCZHFzQ0RERK0XAxddm5oL529T7IfVasZvxy+4uSEiIqLWi4GLrk2nmwC1FgEoxyDZMWw+xNOKREREdWHgomsjVwDdRwGwjzrP67iIiIjqxsBF104aHiIHJy5U4ERRuZsbIiIiap0YuOja1VzHFS0/Bg3KseUIj3IRERE5w8BF1y4wEgjpBTlsGCbfz+u4iIiI6sDARden5rTibfIcbDt2AdUWq5sbIiIian0YuOj61JxWHOW1F5VmC3aeuOjmhoiIiFofBi66Pt2GAwolInAe3WXnsIXfViQiIroKAxddH6Uf0CUOAIeHICIiqgsDF12/2lHn5Tk4qC+F3lDl5oaIiIhaFwYuun497BfOD1fkQgkzTysSERFdgYGLrl9oP8AvFGpUY4j8ME8rEhERXYGBi66fXO5wWvF/R87DYrW5uSkiIqLWg4GLmkfNacVRXvtgrLJgz+kS9/ZDRETUijBwUfPofjsA4EbkIQQGjjpPRER0GQYuah7+HQHdQADAcA4PQURE5ICBi5pP7XVcihzknDGguNzk5oaIiIhaBwYuaj4113Hd7rUPEDb87wiPchEREQEMXNScImMBb18EixL0kZ3iaUUiIqIaDFzUfLxUQLdbAdiHh9hyuAg2m3BzU0RERO7HwEXNSxoeYi+KyqqRe87o5oaIiIjcj4GLmlfNhfNDZYfggyqeViQiIgIDFzW3Dj0AbRd4w4xY+QEGLiIiIjBwUXOTyYAb7IOg3ibfi10nL6K0yuzmpoiIiNyLgYuaX811XKO998FiE9h27IKbGyIiInIvlweu999/H926dYNarUZsbCy2b99eb/2qVavQp08fqNVqDBgwAGvWrHFYLoTA3LlzER4eDh8fH8THx+PIkSMONcXFxRg/fjw0Gg0CAwMxZcoUlJWVSctffvllyGSyqyY/Pz+pZsmSJVctV6vVzbBH2oGokYBMjq7iNCJQxNOKRETU7rk0cK1YsQIpKSlITU3Frl27MGjQICQkJKCwsNBp/bZt25CUlIQpU6Zg9+7dSExMRGJiIvbt2yfVzJ8/HwsXLsSiRYuQmZkJPz8/JCQkoKqqSqoZP3489u/fj/Xr1+Onn37Cli1bMHXqVGn5c889h3PnzjlMffv2xUMPPeTQj0ajcag5efJkM++hNsonEOg0FABwq2IvNh86DyE4PAQREbVjwoViYmJEcnKy9NxqtYqIiAiRlpbmtH7cuHFi7NixDvNiY2PFE088IYQQwmazCZ1OJxYsWCAtLykpESqVSnz11VdCCCFyc3MFALFjxw6p5ueffxYymUycOXPG6etmZ2cLAGLLli3SvMWLFwutVtu0N3wFg8EgAAiDwXBd2/FIG9OESNWI1XPGiK4v/CSOFJS6uyMiIqJGccXnt8uOcJlMJmRlZSE+Pl6aJ5fLER8fj4yMDKfrZGRkONQDQEJCglSfl5cHvV7vUKPVahEbGyvVZGRkIDAwEEOHDpVq4uPjIZfLkZmZ6fR1P/30U/Tq1Qu33nqrw/yysjJ07doVkZGRuP/++7F///5633N1dTWMRqPD1G7V3lfRaz/ksPG0IhERtWsuC1xFRUWwWq0ICwtzmB8WFga9Xu90Hb1eX2997WNDNaGhoQ7Lvby8EBwc7PR1q6qq8OWXX2LKlCkO83v37o3PP/8cP/zwA7744gvYbDYMGzYMp0+frvM9p6WlQavVSlNkZGSdtW1exE2AWgt/WykGyY4h/UCBuzsiIiJym3b/LcXvvvsOpaWlmDRpksP8uLg4TJw4EdHR0Rg5ciS+/fZbdOzYER999FGd25o9ezYMBoM0nTp1ytXtt14KL6D7KAD267i2HbuAU8UV7u2JiIjITVwWuEJCQqBQKFBQ4Hhko6CgADqdzuk6Op2u3vrax4Zqrrwo32KxoLi42Onrfvrpp/jd73531VGzK3l7e2Pw4ME4evRonTUqlQoajcZhatdqTiv+zjcXALBqZzsOoERE1K65LHAplUoMGTIE6enp0jybzYb09HTExcU5XScuLs6hHgDWr18v1UdFRUGn0znUGI1GZGZmSjVxcXEoKSlBVlaWVLNhwwbYbDbExsY6bDsvLw8bN2686nSiM1arFXv37kV4eHiDtVSjJnD1MB+CBuVYlXUaVt7MmoiI2iGXnlJMSUnBJ598gqVLl+LAgQOYNm0aysvLMXnyZADAxIkTMXv2bKn+6aefxtq1a/H222/j4MGDePnll7Fz505Mnz4dACCTyTBjxgzMmzcPP/74I/bu3YuJEyciIiICiYmJAIAbb7wRd911F/785z9j+/bt2Lp1K6ZPn45HHnkEERERDv19/vnnCA8Px913331V76+++ip++eUXHD9+HLt27cKf/vQnnDx5Eo8//riL9lYbFNgF6NATcmFFvM8hnDNUYcsRXjxPRETtj5crN/7www/j/PnzmDt3LvR6PaKjo7F27Vrp9F1+fj7k8kuZb9iwYVi2bBleeuklvPjii+jZsye+//579O/fX6qZOXMmysvLMXXqVJSUlGDEiBFYu3atw6CkX375JaZPn47Ro0dDLpfjwQcfxMKFCx16s9lsWLJkCR599FEoFIqrer948SL+/Oc/Q6/XIygoCEOGDMG2bdvQt2/f5t5NbVuP0cCFIxjf4Qi+PX0TVmw/hdt7hza8HhERURsiE4IjUrqK0WiEVquFwWBov9dzHf0v8MWDsKoC0dfwLqxyFTJmj0bHAJW7OyMiInLKFZ/f7f5biuRi3W8HNJ2hqC7BEx33w2IT+G533UNrEBERtUUMXORacgUwxD7kxgQv+5cdlu84xVv9EBFRu8LARa43eAIgU6DjxV0YqDyL4+fLsfPkRXd3RURE1GIYuMj1NOFAn3sAAC+E2G/BtHw7x+QiIqL2g4GLWsYQ+1AgsWW/wAdVWLP3HIxVZjc3RURE1DIYuKhldL8dCOoGL1MppgTuRqXZiv/sOevuroiIiFoEAxe1DLlcOso1UbkRALBiB08rEhFR+8DARS0nejwg90aocR8GKU4g57QBuWeN7u6KiIjI5Ri4qOX4dwT63gcAeL7DNgDASt7QmoiI2gEGLmpZQx8DANxSsRF+qMS3u06jymx1c1NERESuxcBFLavrcCCkF7ws5XjUfzuMVRas2693d1dEREQuxcBFLUsmky6ef1S1AYDgmFxERNTmMXBRyxv0CKBQoWP5EQyWH0PG8Qs4eaHc3V0RERG5DAMXtTzfYKD/7wEAKUG/AuDF80RE1LYxcJF71Fw8P6xqCzQow6qdp2Gx2tzcFBERkWswcJF7dL4ZCO0HhbUKE3wyUFhajU2Hzru7KyIiIpdg4CL3kMmAobUXz28EILCCpxWJiKiNYuAi9xn4MODth45VJ3Cz7BA2HCxEobHK3V0RERE1OwYuch+1BhjwIADgKe0WWG0CX+867eamiIiImh8DF7lXzcXzw01bEQQjVuw4BSGEm5siIiJqXgxc5F4Rg4HwaChsZoxX/YqTFyrw2/Fid3dFRETUrBi4yP1qjnI9qtoEGWwck4uIiNocBi5yv/4PAioNQkynESfPxZq952CoMLu7KyIiombDwEXup/IHBo4DAEzz24xqiw0/7Dnj5qaIiIiaDwMXtQ41N7QebslER5TwhtZERNSmMHBR66DrD0TGQi4sSPLejNxzRuw7Y3B3V0RERM2CgYtaj5qjXJNUmyCHDct35Lu5ISIioubBwEWtR79EQB2IDpYC3CbPwQ+7z6LSZHV3V0RERNeNgYtaD28fIHo8AGCKz0aUVluwZu85NzdFRER0/Ri4qHUZ8igAYLg1C+G4wBtaExFRm+DywPX++++jW7duUKvViI2Nxfbt2+utX7VqFfr06QO1Wo0BAwZgzZo1DsuFEJg7dy7Cw8Ph4+OD+Ph4HDlyxKGmuLgY48ePh0ajQWBgIKZMmYKysjJp+YkTJyCTya6afvvttyb1Qi7QsRfQ7VbIYcMjXhuxPa8Yx8+XNbweERFRK+bSwLVixQqkpKQgNTUVu3btwqBBg5CQkIDCwkKn9du2bUNSUhKmTJmC3bt3IzExEYmJidi3b59UM3/+fCxcuBCLFi1CZmYm/Pz8kJCQgKqqKqlm/Pjx2L9/P9avX4+ffvoJW7ZswdSpU696vf/+9784d+6cNA0ZMqRJvZCL1BzlmqjaDAWsPMpFRESeT7hQTEyMSE5Olp5brVYREREh0tLSnNaPGzdOjB071mFebGyseOKJJ4QQQthsNqHT6cSCBQuk5SUlJUKlUomvvvpKCCFEbm6uACB27Ngh1fz8889CJpOJM2fOCCGEyMvLEwDE7t276+y9oV4aw2AwCADCYDA0eh0SQpirhHizuxCpGvHn2aliyN9/ESaL1d1dERFRO+GKz2+XHeEymUzIyspCfHy8NE8ulyM+Ph4ZGRlO18nIyHCoB4CEhASpPi8vD3q93qFGq9UiNjZWqsnIyEBgYCCGDh0q1cTHx0MulyMzM9Nh2/fddx9CQ0MxYsQI/Pjjj03qxZnq6moYjUaHia6BlwoY/CcA9iEiispM2HDQ+VFRIiIiT+CywFVUVASr1YqwsDCH+WFhYdDr9U7X0ev19dbXPjZUExoa6rDcy8sLwcHBUo2/vz/efvttrFq1CqtXr8aIESOQmJjoELoa6sWZtLQ0aLVaaYqMjKyzlhowZBIAYJjIRmdZIVbs4GlFIiLyXF7ubsAdQkJCkJKSIj2/+eabcfbsWSxYsAD33XffNW939uzZDts1Go0MXdcquDtwwx2QHduAJMUGvH0oFOcMlQjX+ri7MyIioiZz2RGukJAQKBQKFBQUOMwvKCiATqdzuo5Op6u3vvaxoZorL8q3WCwoLi6u83UBIDY2FkePHm10L86oVCpoNBqHia7D0McAAOOVW6AQFny987SbGyIiIro2LgtcSqUSQ4YMQXp6ujTPZrMhPT0dcXFxTteJi4tzqAeA9evXS/VRUVHQ6XQONUajEZmZmVJNXFwcSkpKkJWVJdVs2LABNpsNsbGxdfabnZ2N8PDwRvdCLaDXXYC/DoG2EoyR78TKrFOw2YS7uyIiImq6Zrv83only5cLlUollixZInJzc8XUqVNFYGCg0Ov1QgghJkyYIGbNmiXVb926VXh5eYm33npLHDhwQKSmpgpvb2+xd+9eqeaNN94QgYGB4ocffhA5OTni/vvvF1FRUaKyslKqueuuu8TgwYNFZmam+PXXX0XPnj1FUlKStHzJkiVi2bJl4sCBA+LAgQPitddeE3K5XHz++edN6qUh/JZiM0ifJ0SqRmTMHSa6vvCT+PXIeXd3REREbZwrPr9dGriEEOKf//yn6NKli1AqlSImJkb89ttv0rKRI0eKSZMmOdSvXLlS9OrVSyiVStGvXz+xevVqh+U2m03MmTNHhIWFCZVKJUaPHi0OHTrkUHPhwgWRlJQk/P39hUajEZMnTxalpaXS8iVLlogbb7xR+Pr6Co1GI2JiYsSqVauu6r2hXhrCwNUMLuYL8XKgEKkacfusj8X0Zbvc3REREbVxrvj8lgkheI7GRYxGI7RaLQwGA6/nuh7LHgYOr8UnlnuwQExE5oujEeSndHdXRETURrni85v3UqTWr+bi+Ye9/weZtQrfZ59xc0NERERNw8BFrV+PeEAbCY0oxd3y7Vi+/RR4YJaIiDwJAxe1fnIFcJN9INQJ3uk4VFCKPacNbm6KiIio8Ri4yDPcNAGQKTBEdgi9ZKewYke+uzsiIiJqNAYu8gwBOqDPPQCAPyrS8WP2WZRXW9zcFBERUeMwcJHnqLl4/g9ev8JmKsfqvefc3BAREVHjMHCR54gaBQRFwR8V+J3iN97QmoiIPAYDF3kOuRwY8igA4E+KdGSdvIh9Z3jxPBERtX4MXORZBv8JkHtjkPwY+snykPbzAQ4RQURErR4DF3kWvxCg730AgAleG7D16AVsOnTezU0RERHVj4GLPE/NxfO/996GAFTgtTUHYLHa3NwUERFR3Ri4yPN0HQ50vBFKWyVmqr/H0cIyLOcF9ERE1IoxcJHnkcmAMfMAAONlP6O3LB/vrD+M0iqzmxsjIiJyjoGLPFPPeODGeyEXVsz3+RculFdj0eZj7u6KiIjIKQYu8lx3vQF4+2KQLRe/l/8Pn/4vD2dLKt3dFRER0VUYuMhzaTsDI2cCAOaql0NlKcWCdYfc3BQREdHVGLjIs92SDIT0QqCtBCleq/Dd7jPIOV3i7q6IiIgcMHCRZ/NSAve8BQCY6PVf9JPl4bXVHAyViIhaFwYu8nzdRwL9H4QcNrymXIzteUVYn1vg7q6IiIgkDFzUNox5DVAGIFp2FOMUm/HGzwdh5mCoRETUSjBwUdugCQdunw0AmO29HMVFenz520k3N0VERGTHwEVtR8xUILQvAlGKmV4r8I/0IzBUcjBUIiJyPwYuajsU3sDYtwEAj3htRJfKA/hg41E3N0VERMTARW1N12HAoCTIITDP+3Ms3Xocp4or3N0VERG1cwxc1Pbc+SqESoMB8hP4A9bjzbUH3d0RERG1cwxc1Pb4h0I2ei4A4HmvFcjIOYhd+Rfd3BQREbVnDFzUNg19DNANhFZWgVleX2HeT7kcDJWIiNyGgYvaJrkCGPt/AICHvLZAdioTP+/Tu7kpIiJqrxi4qO2KvBm4aSIAYJ73YixYsx/VFqubmyIiovaIgYvattEvQ6iDcKM8H7cbv8e/MzgYKhERtTyXB673338f3bp1g1qtRmxsLLZv315v/apVq9CnTx+o1WoMGDAAa9ascVguhMDcuXMRHh4OHx8fxMfH48iRIw41xcXFGD9+PDQaDQIDAzFlyhSUlZVJyzdt2oT7778f4eHh8PPzQ3R0NL788kuHbSxZsgQymcxhUqvV17k3qMX5dYDszpcBAM94fY1l6dtxsdzk3p6IiKjdcWngWrFiBVJSUpCamopdu3Zh0KBBSEhIQGFhodP6bdu2ISkpCVOmTMHu3buRmJiIxMRE7Nu3T6qZP38+Fi5ciEWLFiEzMxN+fn5ISEhAVVWVVDN+/Hjs378f69evx08//YQtW7Zg6tSpDq8zcOBAfPPNN8jJycHkyZMxceJE/PTTTw79aDQanDt3TppOnuTREY80eCJExBAEyCrxV+tS/HMDB0MlIqIWJlwoJiZGJCcnS8+tVquIiIgQaWlpTuvHjRsnxo4d6zAvNjZWPPHEE0IIIWw2m9DpdGLBggXS8pKSEqFSqcRXX30lhBAiNzdXABA7duyQan7++Wchk8nEmTNn6uz1nnvuEZMnT5aeL168WGi12sa/WScMBoMAIAwGw3Vth5rBmV3ClqoVIlUjxv9tvsg7X+bujoiIqJVyxee3y45wmUwmZGVlIT4+Xponl8sRHx+PjIwMp+tkZGQ41ANAQkKCVJ+Xlwe9Xu9Qo9VqERsbK9VkZGQgMDAQQ4cOlWri4+Mhl8uRmZlZZ78GgwHBwcEO88rKytC1a1dERkbi/vvvx/79++t9z9XV1TAajQ4TtRIRgyG7eQoAIFW+GAvW7HVzQ0RE1J64LHAVFRXBarUiLCzMYX5YWBj0eudfz9fr9fXW1z42VBMaGuqw3MvLC8HBwXW+7sqVK7Fjxw5MnjxZmte7d298/vnn+OGHH/DFF1/AZrNh2LBhOH36dJ3vOS0tDVqtVpoiIyPrrCU3uOMlWNQd0FN+Bp0PLcH2vGJ3d0RERO1Eu/+W4saNGzF58mR88skn6NevnzQ/Li4OEydORHR0NEaOHIlvv/0WHTt2xEcffVTntmbPng2DwSBNp06daom3QI3lEwSvu+YBAP7q9S0++nEzbDYOhkpERK7nssAVEhIChUKBgoICh/kFBQXQ6XRO19HpdPXW1z42VHPlRfkWiwXFxcVXve7mzZtx77334p133sHEiRPrfT/e3t4YPHgwjh6t+4JrlUoFjUbjMFErMygJpk6x8JNV4/dFH+A/OWfd3REREbUDLgtcSqUSQ4YMQXp6ujTPZrMhPT0dcXFxTteJi4tzqAeA9evXS/VRUVHQ6XQONUajEZmZmVJNXFwcSkpKkJWVJdVs2LABNpsNsbGx0rxNmzZh7NixePPNNx2+wVgXq9WKvXv3Ijw8vBHvnlotmQzKe/8PNigwVrEdm1avQJWZg6ESEZGLNdvl904sX75cqFQqsWTJEpGbmyumTp0qAgMDhV6vF0IIMWHCBDFr1iypfuvWrcLLy0u89dZb4sCBAyI1NVV4e3uLvXv3SjVvvPGGCAwMFD/88IPIyckR999/v4iKihKVlZVSzV133SUGDx4sMjMzxa+//ip69uwpkpKSpOUbNmwQvr6+Yvbs2eLcuXPSdOHCBanmlVdeEevWrRPHjh0TWVlZ4pFHHhFqtVrs37+/0e+f31Jsvcw/zRQiVSOOzeklFqXnursdIiJqRVzx+e3SwCWEEP/85z9Fly5dhFKpFDExMeK3336Tlo0cOVJMmjTJoX7lypWiV69eQqlUin79+onVq1c7LLfZbGLOnDkiLCxMqFQqMXr0aHHo0CGHmgsXLoikpCTh7+8vNBqNmDx5sigtLZWWT5o0SQC4aho5cqRUM2PGDKnvsLAwcc8994hdu3Y16b0zcLVilQZR8foNQqRqxMK5fxZFpVXu7oiIiFoJV3x+y4QQvGrYRYxGI7RaLQwGA6/naoVsOasg//ZxVAlvfNBvGVLGjXF3S0RE1Aq44vO73X9Lkdov+YA/wBB2C9QyMwbsewNHC8saXomIiOgaMHBR+yWTQfvgP2CBAnfKs7Dmm8Xu7oiIiNooBi5q30L7oHTwEwCAxHP/wG+H6h7YloiI6FoxcFG7F3TX32DwDkUX+Xnkffd3DoZKRETNjoGLSOUP+d1pAIDfV36NtVu2urkhIiJqaxi4iAAEDH4Q+UFxUMksiNz4FA6d5Aj0RETUfBi4iABAJkP4+A9QKtdggOw4jEsfRlGJ0d1dERFRG8HARVTDO6Q7MP5rVECNm205OPrhI6g2Vbu7LSIiagMYuIguE3BDLC7etwQmeOGW6q3I/mAyhM3m7raIiMjDMXARXaHTTXfjyK0LYRUyxJasxp4lM9zdEhEReTgGLiIn+o0ej4x+qQCA6PylOPLtPDd3REREnoyBi6gOwx+agTXhTwIAeuYswLmNH7m5IyIi8lQMXER1kMlkuPPx1/C9/zgAQOjmWSjd9bWbuyIiIk/EwEVUD2+FHCP/8h5+9BoDBWxQ//gETIfT3d0WERF5GAYuogYE+avQ9/FPsE7cAm9YIJaPhzi1w91tERGRB2HgImqEHrpAqB/+FP+zDYDKVonqpQ8ChQfc3RYREXkIBi6iRhrZNxLH7/gIu209oLYYUL34PuDiSXe3RUREHoCBi6gJJo7si5/6v4vDtk5QVRbCtOQ+oKzQ3W0REVErx8BF1AQymQwv/H44/k/3Jk7ZOkJpOAHLvx4Aqgzubo2IiFoxBi6iJlJ6yZE2KQEv+L6C80IDr8J9sH35MGCudHdrRETUSjFwEV2DID8lXpl8H6bhJRiFD+SnMiBWTQKsZne3RkRErRADF9E16hkWgOQ/PoA/m59HlfCG7PA64IdkgDe7JiKiKzBwEV2H23uH4s67H8A08wyYhQLIWQGsmw0I4e7WiIioFWHgIrpOU0ZEIfSm+/Cc+Qn7jMxFwOb57m2KiIhaFQYuouskk8nw98T+ONf1Psw1T7LP3PQ6kPmxexsjIqJWg4GLqBkoveRY9Kch2KhNxDvmB+0zf34eyFnl3saIiKhVYOAiaibBfkp8NulmfOY1DostCQAA8f1fgMO/uLkzIiJyNwYuombUKywAC5MG4+/WCfjOOhwymwVYOQE4meHu1oiIyI0YuIia2R19wjD77n543vwENlgHA5Yq4MuHgOyv+O1FIqJ2ioGLyAUevzUKDwzphifNf0Wm6AeYSoHv/wKs+BNQdt7d7RERUQtj4CJyAZlMhnkP9MeAbjr8sXoWPvEeDyH3Bg7+BHwQC+T+6O4WiYioBbk8cL3//vvo1q0b1Go1YmNjsX379nrrV61ahT59+kCtVmPAgAFYs2aNw3IhBObOnYvw8HD4+PggPj4eR44ccagpLi7G+PHjodFoEBgYiClTpqCsrMyhJicnB7feeivUajUiIyMxf/7V4yY11AtRfVReCiz60xBEBPvjtdKx+INlHoyaXkDFBft1Xd9OBSovurtNIiJqAS4NXCtWrEBKSgpSU1Oxa9cuDBo0CAkJCSgsLHRav23bNiQlJWHKlCnYvXs3EhMTkZiYiH379kk18+fPx8KFC7Fo0SJkZmbCz88PCQkJqKqqkmrGjx+P/fv3Y/369fjpp5+wZcsWTJ06VVpuNBoxZswYdO3aFVlZWViwYAFefvllfPzxx03qhaghHfxV+OYvwxAbFYwsUySGFv4Nm0MnQMjk9lHpP4gDjv7X3W0SEZGLyYRw3VW8sbGxuPnmm/Hee+8BAGw2GyIjI/HUU09h1qxZV9U//PDDKC8vx08//STNu+WWWxAdHY1FixZBCIGIiAg8++yzeO655wAABoMBYWFhWLJkCR555BEcOHAAffv2xY4dOzB06FAAwNq1a3HPPffg9OnTiIiIwIcffoi//e1v0Ov1UCqVAIBZs2bh+++/x8GDBxvVS2MYjUZotVoYDAZoNJpr2IP1qzBX1LlMIVdApVA1qlYuk0Ptpb6m2kpLJer6KySTyeDj5XNNtVWWKthE3fck9PX2vabaams1rDZrs9T6ePlAJpMBAExWEyw2S5213jIV3k0/gvc3HgNkFiR2PINX5R/Dy3DCXjB4AjB6LqD0g9pLDbnM/n8hs9UMs63uG2I3pValUEEhVzS91maGuZ6bcisVSnjJvZpca7FZYLKa6qz1VnjDW+7d5FqrzYpqa3XdtXJveCuaXmsTNlRZqpql1kvuBaXC/rtHCIFKS2Wz1Dbl3z1/RzivddfviKbU8neE3eX/7pubKz6/vZplK06YTCZkZWVh9uzZ0jy5XI74+HhkZDj/inxGRgZSUlIc5iUkJOD7778HAOTl5UGv1yM+Pl5artVqERsbi4yMDDzyyCPIyMhAYGCgFLYAID4+HnK5HJmZmXjggQeQkZGB2267TQpbta/z5ptv4uLFiwgKCmqwF2eqq6tRXX3pl7fRaKx7BzWD2GWxdS67tdOt+CD+A+n5qJWj6vxFPTRsKBbftVh6ftc3d+FitfNTXf069MPy3y2Xnid+n4iz5Wed1t6gvQHfJ34vPU/6KQnHDMec1kb4RWDdH9ZJzx9d+yj2X9jvtDZIFYQtj2yRnk/77zTsLNjptNbHywfbx186jf3MxmfwvzP/c1oLAHsn7ZV+nv2/2Vh/cn2dtZl/zJR++b6S8Qp+PFb3dVmbH96M5xP6YGi3YDy9fg7SA7YhHQCCIu0FFzcBX28CAKx9cC06+XcCACzcvRBL9i+pc7vf3fcdegT1AAB8svcTfLjnwzprvxr7FfqH9AcAfHHgC/xf1v/VWft5wue4WXczAODrw1/j9czX66x9f/T7uK3zbQCA1cdXY87WOXXWvjXyLSR0s49Rlp6fjuc2P1dn7d+H/x2JPRIBANvObkNyenKdtS/GvoikPkkAgF2Fu/DYusfqrE0ZkoLJ/ScDAA4UH0DS6qQ6a6cNmoYno58EABwvOY4HfnygztpH+z2KZ4c+CwA4V34Od31zV521D/d+GC/d8hIA4GL1RYxcMbLO2vtuuA+vjXgNgD2Q1Pfv/s6ud+L/Rl36c+XvCDtP+R0RrA4GAMzfMR8rDq2os5a/I+wu/x3hCVx2SrGoqAhWqxVhYWEO88PCwqDX652uo9fr662vfWyoJjQ01GG5l5cXgoODHWqcbePy12ioF2fS0tKg1WqlKTIyss5aap9u7x2K3w2IqL9o81uAue6jI0RE5Hlcdkrx7Nmz6NSpE7Zt24a4uDhp/syZM7F582ZkZmZetY5SqcTSpUuRlHTpf5wffPABXnnlFRQUFGDbtm0YPnw4zp49i/DwcKlm3LhxkMlkWLFiBV5//XUsXboUhw4dcth2aGgoXnnlFUybNg1jxoxBVFQUPvroI2l5bm4u+vXrh9zcXNx4440N9uKMsyNckZGRPKV4DbVt/XRBpdmE/1t/GEu3nQQAxHXywnsdf4DPwa+hFgLykN7AA4tg1g3g6YIm1vKUIk8pAp7/O4KnFHlKsdFCQkKgUCiuCicFBQXQ6XRO19HpdPXW1z4WFBQ4BK6CggJER0dLNVdelG+xWFBcXOywHWevc/lrNNSLMyqVCiqVqs7lze3yf/juqr38F2Bz1l7+C7s5a1UKFaBo/lqlQil9KDa29pV7b0JcVASe/3oPMk5bEF88Dv++NR79s+YARYeAT+Phfdtz8L7teUDR8C8Vb8WlD/xmrZU3/pdaU2q95F7SL9bmrFXIFfCVN+7vcFNq5TJ5o/9tNKVWJpO5pBbg74hrqW1NvyMaw2X/7j3sd4QncNk7USqVGDJkCNLT05GYmAjAftF8eno6pk+f7nSduLg4pKenY8aMGdK89evXS0fIoqKioNPpkJ6eLgUso9GIzMxMTJs2TdpGSUkJsrKyMGTIEADAhg0bYLPZEBsbK9X87W9/g9lshre3t/Q6vXv3RlBQUKN6IWoOd/XX4cbwADz55S7sP2vEvesD8NyILzCt/EPIc78DNr8JHPoZeOAjIKyvu9slImo2NptAtcWGSrMVVWYrKs1WVJrsP1eZ7fMrzVZUmayostiXOcyrqemtC0Dy7T3c/XYa5NJvKa5YsQKTJk3CRx99hJiYGLz77rtYuXIlDh48iLCwMEycOBGdOnVCWloaAPtQDCNHjsQbb7yBsWPHYvny5Xj99dexa9cu9O9vv5DvzTffxBtvvIGlS5ciKioKc+bMQU5ODnJzc6FW2/8Xc/fdd6OgoACLFi2C2WzG5MmTMXToUCxbtgyA/ZuNvXv3xpgxY/DCCy9g3759eOyxx/DOO+9Iw0c0ppeGuPpbitR2VJmt+PtPufgyMx8AENe9AxYNPgHthln2sboUSuCOl4C46YC8kf+lJiJqBkLYg1F5tQXl1VaUmyyoMNX8XG1BucmKCpMFZdUWVNQur7aizGRBRc3y8mqLFJhqw1WVue7TvE1xa88Q/HtK3V8QuRau+Px2aeACgPfeew8LFiyAXq9HdHQ0Fi5cKB1pGjVqFLp164YlS5ZI9atWrcJLL72EEydOoGfPnpg/fz7uueceabkQAqmpqfj4449RUlKCESNG4IMPPkCvXr2kmuLiYkyfPh3/+c9/IJfL8eCDD2LhwoXw9/eXanJycpCcnIwdO3YgJCQETz31FF544QWH3hvqpSEMXNRUP2Sfwexv96LCZEXHABU+vK8ThuakAkdqvp0VeQuQ+AHQ4Qb3NkpEHsFitaGs2gJjpQXGKjNKqy57rLQ/llaZr1rmEJ5MVlhtrr0PrNJLDh9vBdTetY8K+CgVUHvZH328FVDVLPOpXVZT1yXYF3f2DWv4RZrAIwNXe8bARdfiaGEZnvwyC4cLyiCXAc/e2QvTtNsgX/ciYCoDvH2BMX8Hhk4Bai6yJaK2y2y1oaTCjJIKEy5WmHGxwiT9bKg0S8HJWZAqN9V9Uf+18PFWwE/lBT+VAr5KL/jXPPqpFPBTesFP5QVfZU2NUgFflRf8a+b5Ku2Pam+5PVB5XwpNCrkMEAKwVAHVpVdPpjKg2njZvLKaRyOgGwiMeqHh5puAgcvDMHDRtaowWfDS9/vw7a4zAIBRvTviHwnB0P4yAzhRM0ZQ99uBu94AQvu4r1EiajQhBMqqLSipCU0Xa0NU+WU/1ywzVNYEq3IzSqvr/sZiY/l4KxCg9kKA2gsaH28EqL3tP6u9oZHm2ZcHqLzhr/aqCVC1AcsLPrXByBmrxR5+qgyXHquueF4bkK4KTTXzTWVAPd/OrFP324GJ31/X/rkSA5eHYeCi6yGEwKqdpzHnh32ottgQoVXjvT9G46Zzq4D/ptr/JwgAXYcDQx8DbrwX8Gq5b8kSkf2UXXG5CefLqlFUZkJRaTWKymonk8NjSYUJZuu1feTKZIBG7Y0gX28E+ioR5OuNIF8ltL7e9tDkUxug7CFKClQ1870VDQy7aTUDlSX2a0arSmoCk6H+AHX5c1NZ/dtvKmUAoPIHVAGXJqU/oNLUPL9sWWA3oGd8g5tsCgYuD8PARc0h96wRyct2Ia+oHF5yGWbfcyMe622GLP0V4NAaoHZ8Id8OwOA/AUMmA8FR7m2ayINVW6y4IIWlahSV1gYqx1B1odyEixUmNPVTVOUlR5CvEoE1oSnIzzFE1f4ceFmN1se77qNLtWzWmiBU4hieKi/anzv8bHCc31yBycsHUGsBtcb+qNJceq7SXBaYAi4LTprLAlXNo9ylt3puEAOXh2HgouZSWmXGrG/2YvXecwCAu/rpMP+hgdBUFwK7/gXsWgqUnru0wg2j7Ue9et0FKNrOODZE10MIgeJyE/TGKhQYq6A3VNt/NlRdmmesQklF3QNzOiOXAcF+SoT4q2qmmp8DHJ8H+ykR5KuEj7KBbxoLYQ9AFcVAZTFQcQGouFjz8+XzimtCU+1RKSOA6/xIV2kBH21NSAp0DExXBiiH5zU/ezVu/LDWjoHLwzBwUXMSQuBfGScxb3UuzFaBLsG+eOX+fhjZsyPkwgocXgvs/Bw4ln5ppYAI4KaJ9knbyX3NE7lYldlaE6LsoanQaA9TlweqQmM1TNbGDUXgrZChg58KIQH2sFT7c0cpVF1aFuSrrPvok81mD0MVNSGpNjRd/nNlsT1QXT6vnlHeG27eD/AJtAcmn6DLfg50Mr/m0SfIHpj4HzQADFweh4GLXGHPqRI8+eUunCmx3+Kle0c/TB4ehQdv6gRfpRdQfBzIWgrs/rf9FzgAyORAr7uBmx8Dut/h9sP1RE1htQkUllbhbEklTl+sxJmSSpypeawNWE05KhXir0SYRg2dRo0wbc2jRmWfp1UjLECNQF9v6VY7js1Y7KGovAioKKp5vHD1c2neBUBc4zcFvdSATzDgG2wPRL7B9ksHpHk186Up0B6m2shRJndi4PIwDFzkKoYKMxZuOIIVO06hrOYbTBq1F5JiumDisG7oFOgDWKqBA/+xH/U6ufXSyoFdgaGTgeg/Af4d3fQOiC6ptlhxtqQKZy5W2kOVFKgqpFDVmIvNVV5ye2CqCVOOP9sDVWiAGkqvy/7DURugygqB8vP2kFR+vu4wVVVybW9SpbkUkq4KTkHOlykbfwslal4MXB6GgYtcrbTKjK+zTmPJthM4ecF+Q2GFXIa7+ukweXg3DOkaZP9feuFBIGsxkP0VUG2wryz3BvreZ7/Wq+twjulFLlNWbcGp4grpqJRjqKrE+dK6byBeSyGXQadRo1OQDzoH+qBTkA86BfogPNAHYRoVdBo1tD7ekAGAqfyy8FQbpM4DZecv/Vw7VRSj6dc9yewhyS8E8A0B/DrUPIZc9tjh0qNvB36D2MMwcHkYBi5qKVabwMaDhfh8ax62HbsgzR/YWYvHhkfhngHh9v/VmyqA/d/aj3qdybq0gZDe9uA16GH7BwlRE1isNpwzVOFUcQXya6ZTFyvtj8UVKC43NbgNH2+FFKIiAn3QuebnTkE+6KRVI0xZBUV5oT1A1R6NKiu8OkCVnQcslU18B7KagNTRHpJqH52FKd8O9qNQvMVWm8bA5WEYuMgdDuqNWPzrCXyXfQYmi/0C4dAAFSbc0hV/jO2CDv41/9M+m20PXnu/Bszl9nlePkDvu4EucUDkzUBYf0Dh7Z43Qq2GEAKGSvOlMFV8KUzlF1fgbEklLA3c+iXQ1/tSiAr0RYRWiSg/EyKVZdApShFgKYas/PylQFVWG65qglRTLyL38rGfMverY7p8mU8wLxYnBwxcHoaBi9zpQlk1vtqej39lnERhzSkbpZccidERmDw8CjeG1/ydrDIAOSvt4asw13EjXmogYjDQ+eZLkya8hd8JtYQKkwVnLtpP9Z2+WGkPUxdqj1ZVoLSq/hHAlQo5ugQpcWOgBb39qhDlU47OyjKEyY0IEiVQVV1wPDpVfr7po4qrtYB/GOAXWhOYQq8IVaGXjlCp/BveHlEdGLg8DAMXtQYmiw0/7zuHz37NQ85pgzQ/rnsHPDYiCnf0Cb10H7PTO4BjG4HT2+0/Vxmu3qA2Eug8tCaAxQDhA3l9SitXe4Tqym/4SY8llXWe9lPAimAY0VFmQA/fCvT0q0RXdTkivErRUWZEoO0ifM3F8Kosgqyi6NJAvI3lE1QTnEJrjjzV/nzFPL+O/HtGLYaBy8MwcFFrIoTArvyL+HzrCazdp4e15hRQ1w6+mBTXDQ8N7YwA9WWnD202oPgYcKomfJ3eCRTuv/oDVaG03zw2MuZSENNG8iL8FmSzCRSVVTtciH7lo6m6EhpUQCMrv+KxAhqUQyOrQAdFJcKUVQhTlKEDDNBYL0JlLoGsqReV+3a44ijU5cHpsvl+HTmEAbVKDFwehoGLWqszJZX4V8YJLN9+CoZK+7Ux/iovPDS0M5JiuqBHR3/InQ3kWF0KnN1tD2CndtgfK4qurvPXXQpfkTFAeDS/4n4dLFYbzl0woPB8AYqKzqOkuAhlJUWoKC2GqewibJUl8BPlUnByFqjUsusYSFMmb0SIqpl8Q3g9FHk8Bi4Pw8BFrV2FyYJvd53B4q15OHa+XJofoPbCwM5aDOociIGdAzEoUgudRn31QJBCABdP2I9+1Z6G1O+9+tocmfyyD+UwICDM/ugwhQIBOkDp5/o37g6WascbAl92g2BzeQnKDBdQaSyGudweoGRVBniZS6G2lsJfXGdgksgu3ZKl9tYt0s+X3Z7Ft4NjsPLtwG/lUbvCwOVhGLjIU9hsAv87WoQlNcNKVFuuvg4nNECFgZ0DER2pxcDOgRjYWYtAXyeng8yV9m9Ant5xabr8Po8NUfrXBDNd/QHNL6RlQ4Cl2n6vOieByfl0abmoMkBmqWqWNsplfqj2CoBVaQ9OCt9AqPyDoA4IhsIn0DE8XflcGcC7DBA1AgOXh2HgIk9kttpwuKAUe04ZkHO6BNmnSnCksEy65uty3Tr41hwBC8Sgzlr0i9BefWNeIezfSCs9B5QWAGVXTJfPM1c0vlGZ/FKAkHsBMoX9Ua6omS6f14QayOynTq8MUE0e28k5o/CFEb4wCr+aR18Y4YdKuV9NgAqCyj8IvpoO0AaFIKhDCDqGhCK4Q0fI1RoeaSJqAQxcHoaBi9qKSpMV+88akH2qBDmnDdhzukQa2f5yCrkMvcICMKizVjoV2SssAN6KRh5VqS6tGYOpACjVX/r5ynBWfh5NHx28eVQr/FEh94cRvrho80WRRY0Sm2+dQar2UabWIigoCOFB/ugc5HtpUM+aQT61PnXcu4+IWhwDl4dh4KK2rKTCZA9fp0qwpyaEObtFi8pLjr4RGoT4qxCg8oK/2gv+NY+XnnvDX+WFgMuW+au8oPKS130D4YoioLrMfmNgmwWwXfZ4xTyzxQKTyYRqswkmkxlmsxnVJhPMZjPMFgvMZhMsNT+bzBYUmLxxplKJk+XeOFnuheKaQFUGX9hwdXiUyeynXO3hyVcKUrW3oYkI9IGfiheSE3kKBi4Pw8BF7YkQAnpjFfacsoevnNP2o2ENDZhZHy+57FJAcwhk9oAml9mPvlWYrCg3WaSfK81WVJgs9p9N1gZHQW9MHxGBPg5HpS6/p1+41sfxhshE5NFc8fnN/3IRUbOQyWQI19rDx139dQDsF+PnXSjHwXOlMFSaUVZtRlmVBaXVFpRVWVBWbZ9Ka3++bB4AWGwCJRVmlFQ0xzf07MHJR6mAn9ILvkoFfJSKmkcv+Hrbf/ZTeUGnVTvcyy80QG0fHJaI6BoxcBGRy8jlMtzQ0R83dGzabVZsNoFy06UQ5hDQap6XVpkhBOB7eWiq+dn3ikDl6+0FH6WCR6GIyG0YuIio1ZHLZQhQe9tHvte6uxsiouvH/+4RERERuRgDFxEREZGLMXARERERuRgDFxEREZGLMXARERERuRgDFxEREZGLMXARERERuZjLAldxcTHGjx8PjUaDwMBATJkyBWVlZfWuU1VVheTkZHTo0AH+/v548MEHUVBQ4FCTn5+PsWPHwtfXF6GhoXj++edhsTjeOmTTpk246aaboFKp0KNHDyxZssRheVpaGm6++WYEBAQgNDQUiYmJOHTokEPNqFGjIJPJHKa//OUv175DiIiIqN1yWeAaP3489u/fj/Xr1+Onn37Cli1bMHXq1HrXeeaZZ/Cf//wHq1atwubNm3H27Fn8/ve/l5ZbrVaMHTsWJpMJ27Ztw9KlS7FkyRLMnTtXqsnLy8PYsWNx++23Izs7GzNmzMDjjz+OdevWSTWbN29GcnIyfvvtN6xfvx5msxljxoxBeXm5Qz9//vOfce7cOWmaP39+M+0dIiIialeEC+Tm5goAYseOHdK8n3/+WchkMnHmzBmn65SUlAhvb2+xatUqad6BAwcEAJGRkSGEEGLNmjVCLpcLvV4v1Xz44YdCo9GI6upqIYQQM2fOFP369XPY9sMPPywSEhLq7LewsFAAEJs3b5bmjRw5Ujz99NONf9NOGAwGAUAYDIbr2g4RERG1HFd8frvkCFdGRgYCAwMxdOhQaV58fDzkcjkyMzOdrpOVlQWz2Yz4+HhpXp8+fdClSxdkZGRI2x0wYADCwsKkmoSEBBiNRuzfv1+quXwbtTW123DGYDAAAIKDgx3mf/nllwgJCUH//v0xe/ZsVFRU1Pu+q6urYTQaHSYiIiIil9xLUa/XIzQ01PGFvLwQHBwMvV5f5zpKpRKBgYEO88PCwqR19Hq9Q9iqXV67rL4ao9GIyspK+Pj4OCyz2WyYMWMGhg8fjv79+0vz//jHP6Jr166IiIhATk4OXnjhBRw6dAjffvttne87LS0Nr7zySp3LiYiIqH1qUuCaNWsW3nzzzXprDhw4cF0NtbTk5GTs27cPv/76q8P8y683GzBgAMLDwzF69GgcO3YMN9xwg9NtzZ49GykpKdJzo9GIyMhI1zROREREHqNJgevZZ5/Fo48+Wm9N9+7dodPpUFhY6DDfYrGguLgYOp3O6Xo6nQ4mkwklJSUOR7kKCgqkdXQ6HbZv3+6wXu23GC+vufKbjQUFBdBoNFcd3Zo+fbp0QX/nzp3rfV+xsbEAgKNHj9YZuFQqFVQqlfRcCAEAPLVIRETkQWo/t2s/x5tFs10Ndpnai+Z37twpzVu3bl2jLpr/+uuvpXkHDx50etF8QUGBVPPRRx8JjUYjqqqqhBD2i+b79+/vsO2kpCSHi+ZtNptITk4WERER4vDhw416T7/++qsAIPbs2dOoeiGEOHXqlADAiRMnTpw4cfLA6dSpU43+zG+ITIjmjG+X3H333SgoKMCiRYtgNpsxefJkDB06FMuWLQMAnDlzBqNHj8a//vUvxMTEAACmTZuGNWvWYMmSJdBoNHjqqacAANu2bQNgHxYiOjoaERERmD9/PvR6PSZMmIDHH38cr7/+OgD7sBD9+/dHcnIyHnvsMWzYsAF//etfsXr1aiQkJAAAnnzySSxbtgw//PADevfuLfWs1Wrh4+ODY8eOYdmyZbjnnnvQoUMH5OTk4JlnnkHnzp2xefPmRu8Dm82Gs2fPIiAgADKZ7Pp3qoeqPbV66tQpaDQad7fTZnC/Nj/uU9fgfnUN7lfXqN2vubm56N27N+TyZvp+YbNFtytcuHBBJCUlCX9/f6HRaMTkyZNFaWmptDwvL08AEBs3bpTmVVZWiieffFIEBQUJX19f8cADD4hz5845bPfEiRPi7rvvFj4+PiIkJEQ8++yzwmw2O9Rs3LhRREdHC6VSKbp37y4WL17ssBx1JNnauvz8fHHbbbeJ4OBgoVKpRI8ePcTzzz/P4R2uEYfHcA3u1+bHfeoa3K+uwf3qGq7ary47wkVUy2g0QqvVwmAw8H9hzYj7tflxn7oG96trcL+6hqv2K++lSERERORiDFzkciqVCqmpqQ7f4KTrx/3a/LhPXYP71TW4X13DVfuVpxSJiIiIXIxHuIiIiIhcjIGLiIiIyMUYuIiIiIhcjIGLiIiIyMUYuOi6pKWl4eabb0ZAQABCQ0ORmJiIQ4cONXr95cuXQyaTITEx0XVNeqBr3a8lJSVITk5GeHg4VCoVevXqhTVr1rRAx63fte7Td999F71794aPjw8iIyPxzDPPoKqqqgU69gwffvghBg4cCI1GA41Gg7i4OPz888/1rrNq1Sr06dMHarUaAwYM4N9RJ5q6Xz/55BPceuutCAoKQlBQEOLj46+69zBd29/XWtf7ecXARddl8+bNSE5Oxm+//Yb169fDbDZjzJgxKC8vb3DdEydO4LnnnsOtt97aAp16lmvZryaTCXfeeSdOnDiBr7/+GocOHcInn3yCTp06tWDnrde17NNly5Zh1qxZSE1NxYEDB/DZZ59hxYoVePHFF1uw89atc+fOeOONN5CVlYWdO3fijjvuwP3334/9+/c7rd+2bRuSkpIwZcoU7N69G4mJiUhMTMS+fftauPPWran7ddOmTUhKSsLGjRuRkZGByMhIjBkzBmfOnGnhzlu3pu7XWs3yedWs49ZTu1dYWCgAiM2bN9dbZ7FYxLBhw8Snn34qJk2aJO6///6WadBDNWa/fvjhh6J79+7CZDK1YGeeqzH7NDk5Wdxxxx0O81JSUsTw4cNd3Z5HCwoKEp9++qnTZePGjRNjx451mBcbGyueeOKJlmjNo9W3X69ksVhEQECAWLp0qYu78nwN7dfm+rziES5qVgaDAQAQHBxcb92rr76K0NBQTJkypSXa8niN2a8//vgj4uLikJycjLCwMPTv3x+vv/46rFZrS7XpURqzT4cNG4asrCzp1Mzx48exZs0a3HPPPS3So6exWq1Yvnw5ysvLERcX57QmIyMD8fHxDvMSEhKQkZHREi16pMbs1ytVVFTAbDY3+Lu4PWvsfm2uzyuv61qb6DI2mw0zZszA8OHD0b9//zrrfv31V3z22WfIzs5uueY8WGP36/Hjx7FhwwaMHz8ea9aswdGjR/Hkk0/CbDYjNTW1BTtu/Rq7T//4xz+iqKgII0aMgBACFosFf/nLX3hK8Qp79+5FXFwcqqqq4O/vj++++w59+/Z1WqvX6xEWFuYwLywsDHq9viVa9ShN2a9XeuGFFxAREXFVuKWm7dfm/Lxi4KJmk5ycjH379uHXX3+ts6a0tBQTJkzAJ598gpCQkBbsznM1Zr8C9hARGhqKjz/+GAqFAkOGDMGZM2ewYMECBq4rNHafbtq0Ca+//jo++OADxMbG4ujRo3j66afx97//HXPmzGmhblu/3r17Izs7GwaDAV9//TUmTZqEzZs3NzockHPXul/feOMNLF++HJs2bYJarW6hbj1HY/drs39eXdsZTyJHycnJonPnzuL48eP11u3evVsAEAqFQppkMpmQyWRCoVCIo0ePtlDHnqGx+1UIIW677TYxevRoh3lr1qwRAER1dbWrWvQ4TdmnI0aMEM8995zDvH//+9/Cx8dHWK1WV7Xo8UaPHi2mTp3qdFlkZKR45513HObNnTtXDBw4sAU682z17ddaCxYsEFqtVuzYsaOFuvJ8de3X5v684hEuui5CCDz11FP47rvvsGnTJkRFRdVb36dPH+zdu9dh3ksvvYTS0lL84x//QGRkpCvb9RhN3a8AMHz4cCxbtgw2mw1yuf3yzMOHDyM8PBxKpdLVLbd617JPKyoqpH1ZS6FQSNsj52w2G6qrq50ui4uLQ3p6OmbMmCHNW79+faOvTWrP6tuvADB//ny89tprWLduHYYOHdqCnXm2uvZrs39eNUc6pPZr2rRpQqvVik2bNolz585JU0VFhVQzYcIEMWvWrDq3wW8pXu1a9mt+fr4ICAgQ06dPF4cOHRI//fSTCA0NFfPmzXPHW2h1rmWfpqamioCAAPHVV1+J48ePi19++UXccMMNYty4ce54C63SrFmzxObNm0VeXp7IyckRs2bNEjKZTPzyyy9CiKv36datW4WXl5d46623xIEDB0Rqaqrw9vYWe/fudddbaJWaul/feOMNoVQqxddff+3w97u0tNRdb6FVaup+vdL1fF4xcNF1AeB0Wrx4sVQzcuRIMWnSpDq3wcB1tWvdr9u2bROxsbFCpVKJ7t27i9dee01YLJaWbb6VupZ9ajabxcsvvyxuuOEGoVarRWRkpHjyySfFxYsXW7z/1uqxxx4TXbt2FUqlUnTs2FGMHj1a+vASwvnf05UrV4pevXoJpVIp+vXrJ1avXt3CXbd+Td2vXbt2dfr3OzU1teWbb8Wu5e/r5a7n80omBI+LExEREbkSx+EiIiIicjEGLiIiIiIXY+AiIiIicjEGLiIiIiIXY+AiIiIicjEGLiIiIiIXY+AiIiIicjEGLiIiIiIXY+AiIiIicjEGLiIiIiIXY+AiIiIicjEGLiIiIiIX+3+uQ89qC8WuvwAAAABJRU5ErkJggg==",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.plot(dist_Ne2, en_Ne2_corr, label=\"reali\")\n",
    "plt.plot(dist_Ne2, en_sim, label=\"modello\")\n",
    "plt.plot(dist_Ne2, en_sim*0, '--')\n",
    "plt.legend();"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "5c45f1b6-a2c8-4500-a42d-2321d0b81b2b",
   "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.11.2"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}
