{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "e7c9ad38",
   "metadata": {},
   "source": [
    "# Channel Ablation: What the Potassium Current Does\n",
    "\n",
    "To build intuition for what each conductance contributes, we can **ablate** a channel — set its maximal conductance to zero — and watch how the dynamics change. Here we compare an intact HH neuron against one whose delayed-rectifier potassium current has been removed. Potassium repolarizes the membrane after a spike, so removing it should abolish normal, repetitive spiking."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "6dc7c218",
   "metadata": {
    "execution": {
     "iopub.execute_input": "2026-05-25T10:06:14.613330Z",
     "iopub.status.busy": "2026-05-25T10:06:14.613186Z",
     "iopub.status.idle": "2026-05-25T10:06:16.141691Z",
     "shell.execute_reply": "2026-05-25T10:06:16.140594Z"
    }
   },
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "An NVIDIA GPU may be present on this machine, but a CUDA-enabled jaxlib is not installed. Falling back to cpu.\n"
     ]
    }
   ],
   "source": [
    "import brainstate\n",
    "import brainunit as u\n",
    "import matplotlib.pyplot as plt\n",
    "import braincell"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "5ba8dc8c",
   "metadata": {},
   "source": [
    "## A neuron with a tunable potassium conductance\n",
    "\n",
    "We expose the potassium maximal conductance `gK` as a constructor argument so we can instantiate an intact cell and an ablated cell (`gK = 0`) from the same class."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "d0688c90",
   "metadata": {
    "execution": {
     "iopub.execute_input": "2026-05-25T10:06:16.144483Z",
     "iopub.status.busy": "2026-05-25T10:06:16.144098Z",
     "iopub.status.idle": "2026-05-25T10:06:16.149436Z",
     "shell.execute_reply": "2026-05-25T10:06:16.148574Z"
    }
   },
   "outputs": [],
   "source": [
    "class HH(braincell.SingleCompartment):\n",
    "    def __init__(self, size, gK=36. * (u.mS / u.cm ** 2), solver='exp_euler'):\n",
    "        super().__init__(size, V_th=20. * u.mV, solver=solver)\n",
    "        self.na = braincell.ion.SodiumFixed(size, E=50. * u.mV)\n",
    "        self.na.add(INa=braincell.channel.Na_HH1952(size))\n",
    "        self.k = braincell.ion.PotassiumFixed(size, E=-77. * u.mV)\n",
    "        self.k.add(IK=braincell.channel.K_HH1952(size, g_max=gK))\n",
    "        self.IL = braincell.channel.IL(size, E=-54.387 * u.mV,\n",
    "                                       g_max=0.03 * (u.mS / u.cm ** 2))"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "d832161f",
   "metadata": {},
   "source": [
    "## Simulate intact vs. ablated\n",
    "\n",
    "Both cells receive the same 5 uA/cm^2 current. The intact cell uses the default potassium conductance; the ablated cell sets `gK = 0`."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "0eef8f6e",
   "metadata": {
    "execution": {
     "iopub.execute_input": "2026-05-25T10:06:16.151615Z",
     "iopub.status.busy": "2026-05-25T10:06:16.151392Z",
     "iopub.status.idle": "2026-05-25T10:06:19.639136Z",
     "shell.execute_reply": "2026-05-25T10:06:19.638138Z"
    }
   },
   "outputs": [],
   "source": [
    "intact = HH(1)\n",
    "ablated = HH(1, gK=0. * (u.mS / u.cm ** 2))\n",
    "intact.init_state()\n",
    "ablated.init_state()\n",
    "\n",
    "I = 5. * u.uA / u.cm ** 2\n",
    "\n",
    "def step(t):\n",
    "    with brainstate.environ.context(t=t):\n",
    "        intact.update(I)\n",
    "        ablated.update(I)\n",
    "    return intact.V.value, ablated.V.value\n",
    "\n",
    "with brainstate.environ.context(dt=0.01 * u.ms):\n",
    "    times = u.math.arange(0. * u.ms, 80. * u.ms, brainstate.environ.get_dt())\n",
    "    v_intact, v_ablated = brainstate.transform.for_loop(step, times)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "47b2d97e",
   "metadata": {},
   "source": [
    "## Compare the traces\n",
    "\n",
    "The intact neuron fires a clean spike train; without potassium the membrane cannot repolarize and gets stuck in a depolarized state (depolarization block)."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "4a87054d",
   "metadata": {
    "execution": {
     "iopub.execute_input": "2026-05-25T10:06:19.641282Z",
     "iopub.status.busy": "2026-05-25T10:06:19.640962Z",
     "iopub.status.idle": "2026-05-25T10:06:20.458264Z",
     "shell.execute_reply": "2026-05-25T10:06:20.456976Z"
    }
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAxYAAAEiCAYAAABkykQ1AAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjgsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvwVt1zgAAAAlwSFlzAAAPYQAAD2EBqD+naQAAi79JREFUeJzt3Xd4VFX6wPHvpMyk90Z6CL03QUBABcWGolh/CgIqFkABu64g6mJbsSt2XHVlxbLrrmtBBRRERHqvoaXXmfQyc39/3KnJJJlJIe39PM88mblzZ+bk5uTOfU95j0ZRFAUhhBBCCCGEaAaPti6AEEIIIYQQouOTwEIIIYQQQgjRbBJYCCGEEEIIIZpNAgshhBBCCCFEs0lgIYQQQgghhGg2CSyEEEIIIYQQzSaBhRBCCCGEEKLZJLAQQgghhBBCNJsEFkIIIYQQQohmk8BCCNEkJSUl3HrrrcTExKDRaFiwYAEA2dnZXH311YSHh6PRaHjppZfatJzuqO93ao7k5GQuu+yyRvdbt24dGo2GdevWNfsz7Wk0Gh5//PEWfc/mOvfccxkwYEBbF6NDS05OZubMmW1dDCGEcCCBhRDCauXKlWg0mnpvv//+u3XfZcuWsXLlSu68804++ugjpk+fDsDChQv5/vvvefjhh/noo4+46KKLWrycy5Yt41//+lervK+z36m9+9///tfugoeMjAwef/xxduzY0dZFaTFlZWU8/vjjLR78ifbrjTfeYOXKlW1dDCE6DK+2LoAQov154oknSElJqbO9R48e1vs///wzZ599NkuWLHHY5+eff+aKK67gvvvua7XyLVu2jKuvvpqpU6e26PvW9zu1d//73/94/fXXnQYX5eXleHmd+VN9RkYGS5cuJTk5mSFDhpzxz28NZWVlLF26FFB7XdrSwYMH8fCQtsHW9sYbbxARESG9Q0K4SAILIUQdF198MSNGjGhwn5ycHPr16+d0e0hISCuVrHXV9zt1ZD4+Pm1dBNEKdDpdWxfhjCsrK8PPz6/O9pqaGkwmE1qttg1KJYSwJ80dQgi3WOYCpKWl8c0331iHSVmGUSmKwuuvv27dblFUVMSCBQtISEhAp9PRo0cPnn32WUwmk8P7m0wmXn75ZQYOHIiPjw+RkZFcdNFF/Pnnn4A6Z6C0tJQPP/zQ+hmNtSbm5ORwyy23EB0djY+PD4MHD+bDDz9s9Hc6fvx4ve/5wQcfcP755xMVFYVOp6Nfv368+eab9e7/ww8/MGTIEHx8fOjXrx9ffvllg2UG+PXXX7nmmmtITExEp9ORkJDAwoULKS8vt+4zc+ZMXn/9deuxqX3ca8+xePzxx9FoNBw5coSZM2cSEhJCcHAws2bNoqyszOHzy8vLufvuu4mIiCAwMJDLL7+c9PT0RudtrFu3jrPOOguAWbNmOdQRe/v27eO8887Dz8+PuLg4nnvuuTrvVVlZyZIlS+jRo4f1GDzwwANUVlY2evwsczm2bt3KmDFj8PX1JSUlhRUrVtTZt7E6cvz4cSIjIwFYunSp9XeyHIddu3Yxc+ZMunfvjo+PDzExMcyePZv8/HyHzykuLmbBggUkJyej0+mIioriggsuYNu2bdZ9Dh8+zLRp04iJicHHx4f4+Hiuv/569Hq9dZ/acywsf9faLP+X9nXZMu9n3bp1jBgxAl9fXwYOHGgd4vXll19a//+GDx/O9u3bGz3WoP6PL1y40Pq7xcfHM2PGDPLy8uotCzifX2T/txs/fjx+fn488sgjHD9+HI1Gw9/+9jdeeuklUlNT0el07Nu3D4ADBw5w9dVXExYWho+PDyNGjODrr792ekw2btzIokWLiIyMxN/fnyuvvJLc3FyH47R3717Wr19v/Xu3dU+VEO2d9FgIIerQ6/XWiwELjUZDeHg4ffv25aOPPmLhwoXEx8dz7733AjB06FDrvIQLLriAGTNmWF9bVlbGhAkTSE9P5/bbbycxMZHffvuNhx9+mMzMTIcJ3rfccgsrV67k4osv5tZbb6WmpoZff/2V33//nREjRvDRRx9x6623MnLkSObMmQNAampqvb9LeXk55557LkeOHGHevHmkpKSwevVqZs6cSVFREffcc0+9v5PlQtKZN998k/79+3P55Zfj5eXFf/7zH+666y5MJhNz58512Pfw4cNcd9113HHHHdx888188MEHXHPNNXz33XdccMEF9X7G6tWrKSsr48477yQ8PJw//viDV199ldOnT7N69WoAbr/9djIyMlizZg0fffRRve9V27XXXktKSgpPP/0027Zt49133yUqKopnn33Wus/MmTP57LPPmD59OmeffTbr16/n0ksvbfS9+/btyxNPPMHixYuZM2cO48aNA2DMmDHWfQoLC7nooou46qqruPbaa/n888958MEHGThwIBdffDGgBpmXX345GzZsYM6cOfTt25fdu3fz4osvcujQIZfm2RQWFnLJJZdw7bXXcsMNN/DZZ59x5513otVqmT17NuBaHYmMjOTNN9/kzjvv5Morr+Sqq64CYNCgQQCsWbOGY8eOMWvWLGJiYti7dy9vv/02e/fu5ffff7de9N9xxx18/vnnzJs3j379+pGfn8+GDRvYv38/w4YNo6qqismTJ1NZWcn8+fOJiYkhPT2d//73vxQVFREcHOzCX7dxR44c4f/+7/+4/fbbuemmm/jb3/7GlClTWLFiBY888gh33XUXAE8//TTXXntto0OvSkpKGDduHPv372f27NkMGzaMvLw8vv76a06fPk1ERITbZczPz+fiiy/m+uuv56abbiI6Otr63AcffEBFRQVz5sxBp9MRFhbG3r17GTt2LHFxcTz00EP4+/vz2WefMXXqVL744guuvPJKh/efP38+oaGhLFmyhOPHj/PSSy8xb948/vnPfwLw0ksvMX/+fAICAnj00UcBHMoghHBCEUIIsw8++EABnN50Op3DvklJScqll15a5z0AZe7cuQ7bnnzyScXf3185dOiQw/aHHnpI8fT0VE6ePKkoiqL8/PPPCqDcfffddd7XZDJZ7/v7+ys333yzS7/TSy+9pADKxx9/bN1WVVWljB49WgkICFAMBkOjv5MzZWVldbZNnjxZ6d69u8O2pKQkBVC++OIL6za9Xq9069ZNGTp0qHXb2rVrFUBZu3Ztg5/x9NNPKxqNRjlx4oR129y5c5X6TueAsmTJEuvjJUuWKIAye/Zsh/2uvPJKJTw83Pp469atCqAsWLDAYb+ZM2fWeU9ntmzZogDKBx98UOe5CRMmKIDy97//3bqtsrJSiYmJUaZNm2bd9tFHHykeHh7Kr7/+6vD6FStWKICycePGBstg+ZwXXnjB4XOGDBmiREVFKVVVVYqiuF5HcnNz6/3dnf2tPv30UwVQfvnlF+u24ODgOv8f9rZv364AyurVqxv83ZKSkhz+Byx/19os/9NpaWkOrwWU3377zbrt+++/VwDF19fXoW699dZbdeqlM4sXL1YA5csvv6zznOV/11lZFMV53bf87VasWOGwb1pamgIoQUFBSk5OjsNzEydOVAYOHKhUVFQ4fPaYMWOUnj171jkmkyZNcjivLFy4UPH09FSKioqs2/r3769MmDChwd9dCGEjQ6GEEHW8/vrrrFmzxuH27bffNvn9Vq9ezbhx4wgNDSUvL896mzRpEkajkV9++QWAL774Ao1G43TytLNhHq743//+R0xMDDfccIN1m7e3N3fffTclJSWsX7++Se/r6+trvW/p4ZkwYQLHjh1zGLICEBsb69BaGhQUxIwZM9i+fTtZWVkufUZpaSl5eXmMGTMGRVFcHp5SnzvuuMPh8bhx48jPz8dgMADw3XffAVhbri3mz5/frM+1CAgI4KabbrI+1mq1jBw5kmPHjlm3rV69mr59+9KnTx+HenP++ecDsHbt2kY/x8vLi9tvv93hc26//XZycnLYunUr0DJ1xP5vVVFRQV5eHmeffTaAwzCnkJAQNm/eTEZGhtP3sfRIfP/993WGprWkfv36MXr0aOvjUaNGAXD++eeTmJhYZ7v938WZL774gsGDB9fpFYCm/+/qdDpmzZrl9Llp06Y59CgWFBTw888/c+2111JcXGytK/n5+UyePJnDhw+Tnp7u8B5z5sxxKNu4ceMwGo2cOHGiSeUVQshQKCGEEyNHjmx08rY7Dh8+zK5du+odWpSTkwPA0aNHiY2NJSwsrMU++8SJE/Ts2bPOMI6+fftan2+KjRs3smTJEjZt2lTnAlCv1zsMWenRo0edi6tevXoB6tj9mJgYp59x8uRJFi9ezNdff01hYWGdz2gO+4tHgNDQUEAdOhQUFMSJEyfw8PCokx3MPjNYc8THx9c5JqGhoezatcv6+PDhw+zfv7/RetOQ2NhY/P39HbbZH/uzzz67RepIQUEBS5cuZdWqVXXKZf+3eu6557j55ptJSEhg+PDhXHLJJcyYMYPu3bsDkJKSwqJFi1i+fDmffPIJ48aN4/LLL+emm25qsWFQUPfvb3nvhIQEp9tr17/ajh49yrRp01qsfABxcXH1TsiuXS+PHDmCoig89thjPPbYY05fk5OTQ1xcnPVxQ/8DQoimkcBCCNHqTCYTF1xwAQ888IDT5y0Xeh3F0aNHmThxIn369GH58uUkJCSg1Wr53//+x4svvlhnQnpTGI1GLrjgAgoKCnjwwQfp06cP/v7+pKenM3PmzGZ/hqenp9PtiqI0631b8vNNJhMDBw5k+fLlTvetfRHclq699lp+++037r//foYMGUJAQAAmk4mLLrrI4W917bXXMm7cOL766it++OEHnn/+eZ599lm+/PJL69ySF154gZkzZ/Lvf/+bH374gbvvvpunn36a33//nfj4eKefX1+vgNFodLq9vuPfmvXC3TLa9wI19pzlGN93331MnjzZ6WtqB8Vt/T8gRGckgYUQotWlpqZSUlLCpEmTGt3v+++/p6CgoMFeC3eGViQlJbFr1y5MJpNDi/SBAwesz7vrP//5D5WVlXz99dcOrZ71Dc2xtKbal/vQoUOAmnnGmd27d3Po0CE+/PBDh4nwa9asqbNvU4eaNCQpKQmTyURaWho9e/a0bj9y5IhLr2+JMqWmprJz504mTpzY5PfLyMigtLTUodei9rF3tY7UV4bCwkJ++uknli5dyuLFi63bDx8+7HT/bt26cdddd3HXXXeRk5PDsGHD+Otf/2oNLAAGDhzIwIED+ctf/sJvv/3G2LFjWbFiBU899ZTT97S0thcVFTmkez5Tw3pSU1PZs2dPg/vYl9FeS5TR0uPj7e3d6HnGHa3xvyVEZyZzLIQQre7aa69l06ZNfP/993WeKyoqoqamBlDHTSuKYl2EzJ59K6K/v3+di5P6XHLJJWRlZVkzvYCa9/7VV18lICCACRMmuPnb2Fo67cuk1+v54IMPnO6fkZHBV199ZX1sMBj4+9//zpAhQ+odBuXsMxRF4eWXX66zr+Wi2dVj4gpLq+8bb7zhsP3VV1916fUtUaZrr72W9PR03nnnnTrPlZeXU1pa2uh71NTU8NZbb1kfV1VV8dZbbxEZGcnw4cMB1+uIZQ2F2r+Ts78V4JDtDNSW+dpD2KKiooiNjbWmzzUYDNb/B4uBAwfi4eHRYIpdS2Y0y3wlwJqW+UyYNm0aO3fudKjnFpbj4qyMRqORt99+u9mfHxUVxbnnnstbb71FZmZmneft08i6w51zjRBCeiyEEE58++231tZae2PGjLG2DLrj/vvv5+uvv+ayyy5j5syZDB8+nNLSUnbv3s3nn3/O8ePHiYiI4LzzzmP69Om88sorHD582DqM5Ndff+W8885j3rx5AAwfPpwff/yR5cuXExsbS0pKinWSaW1z5szhrbfeYubMmWzdupXk5GQ+//xzNm7cyEsvvURgYKDbv8+FF16IVqtlypQp3H777ZSUlPDOO+8QFRXl9KKmV69e3HLLLWzZsoXo6Gjef/99srOz6w1EAPr06UNqair33Xcf6enpBAUF8cUXXzgd/225QL777ruZPHkynp6eXH/99W7/XrXfc9q0abz00kvk5+db081aWvsba8lNTU0lJCSEFStWEBgYiL+/P6NGjXK6ont9pk+fzmeffcYdd9zB2rVrGTt2LEajkQMHDvDZZ5/x/fffNzoXKDY2lmeffZbjx4/Tq1cv/vnPf7Jjxw7efvttvL29AdfriK+vL/369eOf//wnvXr1IiwsjAEDBjBgwADGjx/Pc889R3V1NXFxcfzwww+kpaU5lKW4uJj4+HiuvvpqBg8eTEBAAD/++CNbtmzhhRdeANTV3+fNm8c111xDr169qKmp4aOPPsLT07PBOQwXXnghiYmJ3HLLLdx///14enry/vvvExkZycmTJ10+5k11//338/nnn3PNNdcwe/Zshg8fTkFBAV9//TUrVqxg8ODB9O/fn7PPPpuHH37Y2iu5atWqOoFUU73++uucc845DBw4kNtuu43u3buTnZ3Npk2bOH36NDt37nT7PYcPH86bb77JU089RY8ePYiKirImDxBCONEWqaiEEO1TQ+lmqZU61J10s4qiKMXFxcrDDz+s9OjRQ9FqtUpERIQyZswY5W9/+5s17aeiKEpNTY3y/PPPK3369FG0Wq0SGRmpXHzxxcrWrVut+xw4cEAZP3684uvrqwCNpp7Nzs5WZs2apURERCharVYZOHCg0zSo7qSb/frrr5VBgwYpPj4+SnJysvLss88q77//vtPUnpdeeqny/fffK4MGDVJ0Op3Sp0+fOulEnaXc3LdvnzJp0iQlICBAiYiIUG677TZl586ddf4WNTU1yvz585XIyEhFo9E4pB2lnnSzubm5Dp/vLBVoaWmpMnfuXCUsLEwJCAhQpk6dqhw8eFABlGeeeabRY/Tvf/9b6devn+Ll5eVQ5gkTJij9+/evs//NN9+sJCUlOWyrqqpSnn32WaV///6KTqdTQkNDleHDhytLly5V9Hp9g59v+Zw///xTGT16tOLj46MkJSUpr732Wp19Xa0jv/32mzJ8+HBFq9U6HNvTp08rV155pRISEqIEBwcr11xzjZKRkeGwT2VlpXL//fcrgwcPVgIDAxV/f39l8ODByhtvvGF9/2PHjimzZ89WUlNTFR8fHyUsLEw577zzlB9//NGhHLXTzSqKmiJ41KhRilarVRITE5Xly5fXm27W1f9dS3rX559/voEjrcrPz1fmzZunxMXFKVqtVomPj1duvvlmJS8vz7rP0aNHlUmTJik6nU6Jjo5WHnnkEWXNmjVO0806qyONlefo0aPKjBkzlJiYGMXb21uJi4tTLrvsMuXzzz+37mM5Jlu2bHF4rbP/waysLOXSSy9VAgMDFUBSzwrRCI2iyCwlIYQQrtmxYwdDhw7l448/5sYbb2zr4jTo3HPPJS8vr9Gx/0IIIVqGzLEQQgjhVHl5eZ1tL730Eh4eHowfP74NSiSEEKI9kzkWQgghnHruuefYunUr5513Hl5eXnz77bd8++23zJkzp12lehVCCNE+SGAhhBDCqTFjxrBmzRqefPJJSkpKSExM5PHHH+fRRx9t66IJIYRoh2SOhRBCCCGEEKLZZI6FEEIIIYQQotkksBBCCCGEEEI0m8yxsGMymcjIyCAwMLDRxZ+EEEIIIYTo7BRFobi4mNjYWDw8Gu6TkMDCTkZGhmQ6EUIIIYQQopZTp04RHx/f4D4SWNgJDAwE1AMXFBTUxqURQgghhBCibRkMBhISEqzXyQ2RwMKOZfhTUFCQBBZCCCGEEEKYuTJNQCZvCyGEEEIIIZpNAgshhBBCCCFEs0lgIYQQQgghhGg2mWPRBEajkerq6rYuRqfj7e2Np6dnWxdDCCGEEEI0QbsJLH755Reef/55tm7dSmZmJl999RVTp061Pq8oCkuWLOGdd96hqKiIsWPH8uabb9KzZ0/rPgUFBcyfP5///Oc/eHh4MG3aNF5++WUCAgJapIyKopCVlUVRUVGLvJ+oKyQkhJiYGFlHRAghhBCig2k3gUVpaSmDBw9m9uzZXHXVVXWef+6553jllVf48MMPSUlJ4bHHHmPy5Mns27cPHx8fAG688UYyMzNZs2YN1dXVzJo1izlz5vCPf/yjRcpoCSqioqLw8/OTi98WpCgKZWVl5OTkANCtW7c2LpEQQgghRAurKoOaCqipBGOl+tNyA0g4y7bv4R+hOANCEqH7uW1SXHdpFEVR2roQtWk0GoceC0VRiI2N5d577+W+++4DQK/XEx0dzcqVK7n++uvZv38//fr1Y8uWLYwYMQKA7777jksuuYTTp08TGxvb6OcaDAaCg4PR6/V10s0ajUYOHTpEVFQU4eHhLfsLC6v8/HxycnLo1auXDIsSQgghRMsx1tgu6mvKzT8r1JunFmIG2vbd+xWUF9r2qa6wvdYvHCbcb9v3X3OhMM38fJU5YDDfD4yB29fb9l1xDmTtdl4+/yi4/7Dt8fsXwclNMGAaXP1+yx4LNzR0fVxbu+mxaEhaWhpZWVlMmjTJui04OJhRo0axadMmrr/+ejZt2kRISIg1qACYNGkSHh4ebN68mSuvvLLO+1ZWVlJZWWl9bDAY6i2DZU6Fn59fS/xKoh6W41tdXS2BhRBCCNEZGWvUC/vqCqguUy/qg8wjFRQFDv7PLgAw/6w2BwIhiTD0Rtt7fTkHKgx196+pUAOFaz+07ftCbyjLc16mbkMcA4A1i6HopPN9w3s6BhYZ2yBnn/N9NbXyJHn5mLd7gpdOvXmaf/pHOO6bMAp8gqHbYOfv3Q51iMAiKysLgOjoaIft0dHR1ueysrKIiopyeN7Ly4uwsDDrPrU9/fTTLF261K2yyPCn1iXHVwghhGgjJhOU5Tte9NeYf1ZXqK3vccPUfavLYeMrdfexPE4cDePVUSbUVMJLA23Pm2olwOlzGVz/iXpfo4HPZoCpxnkZUyY4BhaHvoMKvfN9fUMdH1su6kENZrx8zBf3PhDgeI1J9/PUY+GlAy9f235eurr7TloK1aXq855a276eWvCu1SB983/Bwws8XbgEv8C9a9T2oEMEFq3l4YcfZtGiRdbHliXLhRBCCCHaDUUBxQQe5p786nLI2a9epFeVqRe1VWXmx6UQO8Q2Jr8kB76519zibwkAym23wdfbLmBLc+GFXvWXY/D/wZVvmstkgnXL6t/X/oLaU6uWAyej771867bqJ442P2e5mPexXdRH9nbc98Kn1ONjHyR4m/f3CXHc967f1LJ46sCjkRUXLn+l4eft9brQ9X29fRrfpwPrEIFFTEwMANnZ2Q6TerOzsxkyZIh1H8vEX4uamhoKCgqsr69Np9Oh0+lap9DtyLnnnsuQIUN46aWX2rooQgghROdkMqoX9Z7e4O2rbqvQQ/pW80V/uV0AYP7Z/VzoPkHdN/8o/HeBup9DsGB+3bh74fy/qPsWnYR3zqu/LKPutAUWphrY/3X9+5bl2+5bLnotF/Lefuo2bz/1cWiybV8vXxg+S/1dvX3Vx96+tv3t99Vo4I4N5vf0td3qu8Cf+d/6y1vbsBmu7+sT7Pq+okk6RGCRkpJCTEwMP/30kzWQMBgMbN68mTvvvBOA0aNHU1RUxNatWxk+fDgAP//8MyaTiVGjRrVV0duFL7/8Em9vb5f2PX78OCkpKWzfvt16rFvCypUrWbBggaTqFUII0fZq9wCUF6o9AJUlUFWiBghVpVBVrP7sdTEkmq8lMnbAtw+o2yuLbfvWlKvPX/gUjJmv3s89BB/VneNp5am1BRbGKkj7pf59q8ps93WBEBSvXpxr/cDb3/zTD7T+tuFKoA4HuuRvdhf0fo5Bg7/dMHJdECwubLw1H9R9przU+H4WMQNc31d0WG4FFvv372fVqlX8+uuvnDhxgrKyMiIjIxk6dCiTJ09m2rRpTe4BKCkp4ciRI9bHaWlp7Nixg7CwMBITE1mwYAFPPfUUPXv2tKabjY2NtWaO6tu3LxdddBG33XYbK1asoLq6mnnz5nH99de7lBGqMwsLC2vrIgghhBDNV2EA/alaAYDdz56TIbqfuu+pP2D9c3X3sQQCl78Kw6ar+6ZvhY+n1f+5/lG2wMJYBac217+vfQDgGwpR/R0v+r39bMFAvC3hDMHxcNW7tfb1td3X2WXjCYqFRXtdO2bevjDyNtf21WjUmxBN5FJgsW3bNh544AE2bNjA2LFjGTVqFFdeeSW+vr4UFBSwZ88eHn30UebPn88DDzzAggUL3A4w/vzzT847z9atZ5n7cPPNN7Ny5UoeeOABSktLmTNnDkVFRZxzzjl899131jUsAD755BPmzZvHxIkTrQvkvfKKG2PkOin7oVDJycnMmTOHI0eOsHr1akJDQ/nLX/7CnDlzALV3CGDo0KEATJgwgXXr1rFlyxYeeeQRtm/fTnV1NUOGDOHFF19k2DBbq0hRUREPPvgg//rXv9Dr9fTo0YNnnnmGgIAAZs2aBdgmZy9ZsoTHH3/8DB4FIYQQZ5SiqGP6K4vNN4PtftxwdSIwwOmtsPNT5/tVlcBlL0GfS9R9D/8AX9xS/2f6R9oCiwo9HFlT/75VJbb7vqEQ3kO9gNcGmn/6gy4AtAGOaUgjesJ1H5v3Caj1099xgnBED3Vcvyt0gTDoGtf2FaKdcmkdi5SUFO677z5uvPFGQkJC6t1v06ZNvPzyywwaNIhHHnmkJct5RjSUp7eiooK0tDRSUlKswYyiKJRXG9uiqPh6e7qcQal2YFFcXMyTTz7JhRdeyOeff86jjz7Kvn376N27N1u2bGHkyJH8+OOP9O/fH61WS1hYGD///DMZGRmMGDECRVF44YUX+O9//8vhw4cJDAzEZDIxduxYiouLefHFF0lNTWXfvn14enoyceJE3nzzTRYvXszBgwcBCAgIcLoiurPjLIQQ4gyrrlAvzH2CbePucw/Bqd/tAoBat/MfhVi1UYptH6nzBerL7HPdx9B3inp/9+cNBwtXvGHLAnTkR/jydvMFf6D5p78tIBhyAySfo+5ryICjP5sv+gPsAgVzIOATrE72FUI0qMXXsTh06JBLY/RHjx7N6NGjrWs+dHbl1Ub6Lf6+TT573xOT8dM2bYrMJZdcwl133QXAgw8+yIsvvsjatWvp3bs3kZGRAISHhztMej///PMd3uPtt98mJCSE9evXc9lll/Hjjz/yxx9/sH//fnr1UjNKdO/e3bp/cHAwGo2m3on0QgghWojJpLb6V+htPyvsHve7wtZbsP8/sOVdx+cr9OpwH4Cb/wMp49X7x3+FbxY5/0yA4TNtgYWXj11QoVEv5HWBtpt9xqDo/jD+frvngxz3DbbL1thjEjxw1LXjEBQLQ29ybV8hRItw6crU29ub1157jZtuuqnBHgv7/UX7NWjQIOt9y8V+7YxatWVnZ/OXv/yFdevWkZOTg9FopKysjJMn1cVjduzYQXx8vDWoEEII0UQmo3pxrwtUMwwBZO2B01vMF/5OgoYpL0Ok+fy78UX46Yn63z+qry2wKMmGY+vq2VGjzkWwCOuuzmGwv+i3DwK62b5b6H0xLNpvDiL8G54MHNXXlu1ICNGhudzk/eijj/LAAw8wdepUbr311jot2F2Rr7cn+56Y3Gaf3VS1Az+NRoPJZGrwNTfffDP5+fm8/PLLJCUlodPpGD16NFVVaquWr69vk8sjhBCdjskElXooL4KKItvP1IngYx5KsPdfsOcLx+fLi9SAAdT0nJax/Ye/bzhYKMm2BRaWSb5ePup9n2D1M32CzbcQ2+tSzoUr37Y9r7PbTxvgGBCknqfeXKELUG9CiC7F5cAiKyuL1atX88EHH3DBBReQmJjI7NmzmTlzZpddVE6j0TR5OFJ7pdVqATAaHeeObNy4kTfeeINLLlEn0J06dYq8vDzr84MGDeL06dMcOnTIaa+FVqut855CCNEhlOaD4bTjxX95oe3+eY9CgDqMlF9fUFcjrtDjdDGw23+1teznH2l4fQH71YQj+6gpT2sHCZZAILKPbd9hM2DodNcW4orood6EEKIFuHxV7Ovry4wZM5gxYwbHjh1j5cqVvPfeeyxdupRJkyZxyy23MHXqVBkG1cFFRUXh6+vLd999R3x8PD4+PgQHB9OzZ08++ugjRowYgcFg4P7773fopZgwYQLjx49n2rRpLF++nB49enDgwAE0Gg0XXXQRycnJlJSU8NNPPzF48GD8/Pzw8/NroCRCCNGCTEZ1dV9L0ovMXZC1C8oKoLzA7meh+vOmL9Qx+qAOLfrt1frfe8RsW2ABasBh4eWrZhzyDVF7Cjzsept7TLT1IFiet//pafd92udS9eYKmZAshGgjTWpu7969O0888QRLly7lxx9/ZOXKlcycORN/f/9Gx+qL9s3Ly4tXXnmFJ554gsWLFzNu3DjWrVvHe++9x5w5cxg2bBgJCQksW7aM++67z+G1X3zxBffddx833HADpaWl1nSzAGPGjOGOO+7guuuuIz8/X9LNCiGarqrMFgxE9bVdgB/8DtLWOwkWCtTW/wW7IcTcw777s4aDhdJcW2AREAMB0eYL/tC6QYBfuO11Q6dDn8tszzd0kR871DbZWQghOgGX0s02ZO3atbz33nt8+eWX6HQ6CgsLW6psZ5y76WZFy5PjLEQXYzKpLfyleVCWD2V56gRhL3VYJltXqnMRyvLUAKEsX10bwWLBbghJVO//8Bj81sDaRXPWQ+wQ9f6OT2HP5+AbBn5hdj9D1Z9xI2xzIYQQogtr8XSztZ06dYoPPviAlStXcvLkScaPH88777zDtGkNrFophBCi8zPWqL0EpXlqMGAJGOzH/G98Wb2wtwQLSq35V/bBQsExOLa27ud4eKnBgP0qxykT1KFOvmFqL4JDwGD+aTHkBvUmhBCixbgcWFRVVfHll1/y/vvv8/PPP9OtWzduvvlmZs+e7bBegRBCiE5EUdQehZJcKM2Bkhx1mFBJDoy7F7TmuVI/LoU/33ecX2Cv54UQmqTeLyuA3P2Oz+uC1GDAPwKMdmsh9bsCIvuq2/3MAYNvmJrGtPYioT0nqTchhBBtwuXAIiYmhrKyMi677DL+85//MHnyZDwayksthBCifTKZ1F6FkhxzsGAXNJz7EHibEzN8/yj88bZtsbTahk0HbbJ6XzHZBRUadUiRfwT4RYB/uGMQMOT/oPsE83MRarBQ31yEuOHqTQghRLvncmDxl7/8henTp1tXZhZCCNGOmIzqsKOS7LrBwnmP2noWvnsENq+oO/zIYvhMCEtR73tqbUGFLljNfORvvgVEgaddMDByDgy+Xg0W/MIcsx/VFtlbvQkhhOhUXA4sFi1aVGdbSUlJnYXVGpvUIYQQwg01lWqwUJwNJVlQnKU+HnevrWfhu0dg85tqr4EzZ92irpoM6jwHS1DhG6YGCJZAwT/Ssefg7LvUVKr+kY2viRAcB8Q161cVQgjRsbk9eTstLY158+axbt06KipsmTkURUGj0cgiaEII4YrKEnPAkGUOGLJhxCxbsPDzX2HLO+pCbM4MvgHCU9X7Wj81qNB4mOcpRJl7F6LUgMHLtuYMZ8+Fs25ThyB5NrLuUID0UAshhHCd24HFTTfdhKIovP/++0RHR6OpPXlOCCG6suoKKM4EQ4b6s89ljtmQtn6oBhRVJXVf22uyLVhAsQUVnlp1HYWAaAiMUW/2PQuj7oSzblV7FhoaggTqfAchhBCiFbgdWOzcuZOtW7fSu7eMjxVCdCGK+UK/OBMietla+3eugj1fgCETDOnqpGh787fZgoXKEig4anvO2x8Co9UF2AKj1R4Hi+GzoP9VahDhG1o3A5I9CRaEEEK0A24HFmeddRanTp2SwEII0XkYq0HjCZZMd8fWwZGfzD0P5oChONO2MNu8rRDRQ71fcAwO/+D4fl4+ENhNXbnZVGPbPvh6NRuSJZDQBdZfpuA487wFIYQQomNwO7B49913ueOOO0hPT2fAgAF4ezuO0R00aFCLFU4IIVpMwTHI3KUGCfrTtpshXc2cNO9PW7BwYlP9Kzj7hUOF3va498UQFGe+dVMDivp6GMJT7YY6CSGEEJ2L24FFbm4uR48eZdasWdZtGo1GJm93URMmTCA1NZX333+/rYsiuqoKPejNwYLBEjCYH099HUKT1f12/hPWP1P/+xjSbYFF0mg1I5Kl1yEoVr0f2K1udqTYoepNCCGE6OLcDixmz57N0KFD+fTTT2XydhenKArbt2/nmmuuaeuiiM7KWF2rh+EUDJ2hDiMC+OV5+Pmp+l9feNwWWET2hviREBxvHmaUoPYyBJt7G/wibK/rfq56E0IIIYTL3A4sTpw4wddff02PHj1aozyilZx77rkMGjQIHx8f3n33XbRaLXfccQePP/44AJWVldx///2sWrUKg8HAiBEjePHFFznrrLPqfc/Dhw9TXFzM8OGyKq5oImO1GjAERNsWcNv3Nfz+JhSdhOKMumszxJ9lCywCzD99QyEo3i5oiFcfR/a1vW7AVepNCCGEEK3C7cDi/PPPZ+fOnRJYdEAffvghixYtYvPmzWzatImZM2cyduxYLrjgAh544AG++OILPvzwQ5KSknjuueeYPHkyR44cISwszOn7bd26FU9PTwYPHnyGfxPR4RQcg+Mb1WDB/mYJHGZ8rU5qBnVo08nfbK/11JkDBvPNN9T23ICrYcA00Pqf2d9HCCGEEHW4HVhMmTKFhQsXsnv3bgYOHFhn8vbll1/eYoXrMKpK639O4+k4JrvBfT1si2M1tG8TL6IGDRrEkiVLAOjZsyevvfYaP/30E2PGjOHNN99k5cqVXHzxxQC88847rFmzhvfee4/777/f6ftt27aN3r174+entjQvX76cvLw8li1bRnFxMVdddRXjxo1j8eLFTSqv6ABMRnW9hsI0ddiRfdAwaSkkjlL3O74Rvp7n/D28fBxTtKaMh2nvQUgShCSa12bwcP5aSy+HEEIIIdqc24HFHXfcAcATTzxR57kuO3l7WWz9z/W8EG5cbXv8fA+oLnO+b9I5MOsb2+OXBkJZft39HtfX3eaC2hm7unXrRk5ODkePHqW6upqxY8dan/P29mbkyJHs37+/3vfbtm2bwzCoPXv2MGnSJLKzs5kyZQq33norc+bMaVJZRTtSVQqFJ9TgodtgtdcA1LUbvroDjFXOX5d3yBZYRPWFHpPUQMF6swsc7OdqhSapNyGEEEJ0KG4HFiaTqfGdRLtUu3dJo9E06++5bds2pkyZYn28Z88epkyZwsSJE3nqqaeYOnVqk99btJG8I7B7ta0HovC4ukq0xRWvw9Cb1Pt+EWpQ4eGtBgihyWpAEJKoToxOGGV7XfwIuOmLM/iLCCGEEOJMczuwEE48klH/cxpPx8f3H2lg31rDPRbsbnqZ3JCamopWq2Xjxo0kJaktxdXV1WzZsoUFCxY4fc2xY8coKiqy9lgoisK+ffuYN28eq1atYty4cWek7MJFJpOaijX/qLryc/4xc+CQBuPuhYFXq/vpTzlPyeoTogYOXnbD+uLPUutoUBx4eNZ9jRBCCCG6FJcCi1WrVnH99de79IanTp3i5MmTDsNqOj135jy01r7N4O/vz5133sn9999PWFgYiYmJPPfcc5SVlXHLLbc4fc3WrVvRaDQMGTIEgLS0NGJjYwkNDaWgoMDpa0QrUxQozlIDh8ButoXYTvwGf58Kxkrnr8s9aLsf2QeGzYDQFDWQCDP/tJ8wbaH1A21iC/8SQgghhOioXAos3nzzTZYuXcqsWbOYMmUKffv2dXher9ezceNGPv74Y+uEX9GxPPPMM5hMJqZPn05xcTEjRozg+++/JzTUyQUl6jConj17EhgYCKjDoIYMGcLTTz/NxIkTGTp0KImJctHZasoK4OC35t4H863gGFSbJ/yf+zCc+5B6P7CbGlR4eKuBQlgqhHU3Bw0pEN3f9r5B3eDyV8/87yOEEEKIDk+jKIriyo5ff/01r776Kj///DP+/v5ER0fj4+NDYWEhWVlZREREMHPmTBYuXEh0dHRrl7tVGAwGgoOD0ev1BAUFOTxXUVFBWloaKSkp+Pj41PMOXdeyZcuoqalh8eLFfPzxx7zxxhv88ssveHm5N9pOjrNZdYUaNOQehLzDkH9YXbDNMr8h/yi8Oqzu6zSe6hyH4TPhnAXqNpNRzdIUnACeMvpRCCGEEK5r6Pq4NpevMi6//HIuv/xy8vLy2LBhAydOnKC8vJyIiAiGDh3K0KFD8agvJaTo9Pbs2cPVV6vj9G+66SbWrFnDo48+yrPPPtvGJWvnaqrAS6veL8mFf9+lZlMqPAHUivk1nrbAIiQRup+n9jyEp6q9EOE91O2W97Pw8FR7J4QQQgghWpHLPRZdgfRYtL1OeZxNRig6ofY85B0y3w6rvRF9LlEzLQHUVMJfY2wrTfsEQ0RviOgFET3UydLJ57Td7yGEEEKILqdVeiyEEI0w1qjzHGoqoJt5zZCaSngmCWrKnb8m95DtvpcOrnwLgmLVYKL2+g5CCCGEEO2YBBZCNEX+Ucg9ADn7IOeAej/vkLqug/1Ch146NUAoyYaInuZbL9stvIfj+w669sz/LkIIIYQQLUACCyHqoyjqug45+9XVpwdcZXvug4sdF46z8PZTgwl7t/4I/hGy1oMQQgghOjUJLISwOLUFMrZB1m61JyL3IFSVqM8FxjoGFrHDwJAOUX3VtR+i+kFUHwhOhNpJDAI7ZpY0IYQQQgh3SGAhuhZFUSdSZ+0G/Wk4+07bcz/8BU797ri/h7c6fCmqrzqHwpKu9f9WnbkyCyGEEEJ0AC4FFosWLXL5DZcvX97kwnQEkkSrdbX48c3ZD6f+gOw9kLVH/VlpMD+pUVeZtqxwnnqeusJ0zAB10biofmo6V0/vli2TEEIIIUQn5FJgsX37dpfeTNOJM9h4e6sXl2VlZfj6+rZxaTqvsrIywHa8XVaSA1m71OBh9FxbMLDxFdj5D8d9PbzVYUvRA9W5E5bAwrJStRBCCCGEcJtLgcXatWtbuxztnqenJyEhIeTk5ADg5+fXqQOpM01RFMrKysjJySEkJARPzwYmOhedhJObIXu3OqQpaw+U5tie73khRPdT7yeerc6FiBmo3qIHqNmYai8iJ4QQQgghmkXmWLghJiYGwBpciJYXEhJiPc6UF0H2XnX4Ur+ptknQO/8Ja5+q9UqNmro1ZqDj2g/Db1ZvQgghhBCiVTUpsPjzzz/57LPPOHnyJFVVVQ7Pffnlly1SsPZIo9HQrVs3oqKiqK6ubuvidC5lBXhn/onn/t2wbq/aG1F00vZ8UBz0vUy9HzcMEkbZeiBiBqmTq7V+bVN2IYQQQgjhfmCxatUqZsyYweTJk/nhhx+48MILOXToENnZ2Vx55ZWtUcZ2x9PTs+GhOqJ+lcWQvU8NHJLGqgEBwOHf4Itb6u4fFK8GED52S8j3mKjehBBCCCFEu+F2YLFs2TJefPFF5s6dS2BgIC+//DIpKSncfvvtdOvWrTXKCMDjjz/O0qVLHbb17t2bAwcOAFBRUcG9997LqlWrqKysZPLkybzxxhtER8saAm2mNB+OrFHXg8g7pA5pKjxue/6CJ22BRbfB0G2IOSOT5dYf/MLaouRCCCGEEMJNbgcWR48e5dJLLwVAq9VSWlqKRqNh4cKFnH/++XUu/ltS//79+fHHH62PvbxsxV+4cCHffPMNq1evJjg4mHnz5nHVVVexcePGVitPl6coajamPHPgkHsIUsbbhizpT8FXt9d9XWA3NXAIjrdti+gJt68/M+UWQgghhBAtzu3AIjQ0lOLiYgDi4uLYs2cPAwcOpKioyJoqtLV4eXnZJvba0ev1vPfee/zjH//g/PPPB+CDDz6gb9++/P7775x99tmtWq4upTgbfn7C1gtRoXd83lRtCywiekLSORDZS83EFNVPDSj8w898uYUQQgghRKtyO7AYP348a9asYeDAgVxzzTXcc889/Pzzz6xZs4aJE1t33Pvhw4eJjY3Fx8eH0aNH8/TTT5OYmMjWrVuprq5m0qRJ1n379OlDYmIimzZtksDCXZXFcPRnOPk7nN4CSWPggifU57x9YPvHtn01HhCSBJG91UAiZYLtOa0/zPrmzJZdCCGEEEK0CbcDi9dee42KigoAHn30Uby9vfntt9+YNm0af/nLX1q8gBajRo1i5cqV9O7dm8zMTJYuXcq4cePYs2cPWVlZaLVaQkJCHF4THR1NVlZWve9ZWVlJZWWl9bHBYKh33y4hYztseBEOfgdG23GxLiAH4BOsBhnBCWowEZaqBhtCCCGEEKJLczuwCAuzTab18PDgoYfOzGrFF198sfX+oEGDGDVqFElJSXz22WdNXgn76aefbtU5IR1GTRX88Bf44y3btrBUSD0fEkZC/AjH/cfec2bLJ4QQQggh2j2XAguDwUBQUJD1fkMs+7W2kJAQevXqxZEjR7jggguoqqqiqKjIodciOzvb6ZwMi4cffphFixZZHxsMBhISElqz2O2TRgPFmer9gdfC2LvVuRCysrgQQgghhHCRS4FFaGgomZmZREVFERISgsbJBaeiKGg0GoxGY4sX0pmSkhKOHj3K9OnTGT58ON7e3vz0009MmzYNgIMHD3Ly5ElGjx5d73vodDp0Ot0ZKW+75ukNV78P6dsgcVRbl0YIIYQQQnRALgUWP//8s3UI1Nq1a1u1QPW57777mDJlCklJSWRkZLBkyRI8PT254YYbCA4O5pZbbmHRokWEhYURFBTE/PnzGT16tEzcbojJBB4e6n1PbwkqhBBCCCFEk7kUWEyYYMv0k5KSQkJCQp1eC0VROHXqVMuWzs7p06e54YYbyM/PJzIyknPOOYfff/+dyMhIAF588UU8PDyYNm2awwJ5ogH/vQeqyuD8RyGse1uXRgghhBBCdGAaRVEUd17g6elpHRZlLz8/n6ioqDM2FKo1GAwGgoOD0ev1Z2yuSJupKoXne0B1Gcz+QXorhBBCCCFEHe5cH3u4++aWuRS1lZSU4OMjaUc7jIPfqkFFaIqa+UkIIYQQQohmcDndrCV7kkaj4bHHHsPPz8/6nNFoZPPmzQwZMqTFCyhayZEf1Z/9rpDsT0IIIYQQotlcDiy2b98OqD0Wu3fvRqvVWp/TarUMHjyY++67r+VLKFrHiY3qz5RxbVsOIYQQQgjRKbgcWFiyQc2aNYuXX365889B6Mz06VB0EjQekCBzK4QQQgghRPO5vfL2Bx980BrlEGdShtr7RHR/0AW2bVmEEEIIIUSn4Pbk7dLSUh577DHGjBlDjx496N69u8NNdAAenhA7FOJGtHVJhBBd0J/HC7jqjY38e0d6WxdFdALf7cninGd/5rM/Wy/lveg6vtuTKfWpGdzusbj11ltZv34906dPp1u3bk4zRIl2rvfF6k0IIdrA898fZNvJIvak72Jy/xh8vD3bukiiA1u+5iCnC8t54PNdXDEkFp2X1CfRdC/8cEjqUzO4HVh8++23fPPNN4wdO7Y1yiOEaKe+25PJ51vTefTSvqRE+Ld1cUQHpSgKezMMAFQZTWxOK2BCr8g2LpXoqBRF4XhemfXx78ekPommUxSFEwVSn5rD7aFQoaGhhIWFtUZZxJmgKOpNCDeYTAr3r97Fj/uzWfL13rYujujADBU1lFTWWB9vOJzbhqURHV1eSRVVRpP18Za0gjYsjejo8kqqqKqx1ac/0vLbsDQdk9uBxZNPPsnixYspKytrfGfR/pTmwlPR8OpwMHXcVdLFmXWqsIxi88XghsO5lNpdGArhjsLSKofHO04VtU1BRKeQbahwePzHcQksRNPVrk9bjhe2UUk6LreHQr3wwgscPXqU6OhokpOT8fb2dnh+27ZtLVY40Qr0p8BYCVWl6iRuIVxwLLfUet+kwPaTRZzTM6INSyQ6qoIyx8BiT7oBo0nB00Pm6wn3VdY4NpDtOl0k9Uk0mdSn5nM7sJg6dWorFEOcMfrT6s/g+LYth+hQCmq1Mu88LYGFaBpLj0X/2CCO55VSWmXkSE4JvWMk9bVwX0W1OmylR1QA6YXllFcbScsrpUdUQBuXTHREUp+az+3AYsmSJa1RDnGm6M3pHSWwEG7Ql1c7PN5+sqhtCiI6vOIKdRhdqJ+WgDgvNqcVsOt0kQQWokkqqtUWZj+tJ71jAtlxqoj9mQa5EBRNYumxkPrUdG7PsQAoKiri3Xff5eGHH6agQB3PuG3bNtLTJSd5u2fpsQiKa9tyiA6lyBxY9DFf/O3L0LdlcUQHZvni9vH2oG+3IAAOZRe3ZZFEB1Zpnmjr4+VprU/7Mw1tWSTRgVl6LHReHvSLlfrUFG73WOzatYtJkyYRHBzM8ePHue222wgLC+PLL7/k5MmT/P3vf2+NcoqWYrAMhUpo23KIDkVvHhd/VnIYB7KKydBXUFxRTaCPdyOvFMKR5UJQ5+VJr2g1UD2UXdKWRRIdmKXHQudtuxDcJxeCooks9cnHWwLVpnK7x2LRokXMnDmTw4cP4+PjY91+ySWX8Msvv7Ro4UQrKM1TfwZIXmbhOstQqKRwP6KDdIBcDIqmqbRrEewVrQ4vOCw9FqKJbC3MnvTrpgaqBzKlPommsW/4sNSn/VKf3OJ2YLFlyxZuv/32Otvj4uLIyspqkUKJVhTVD+JHSo+FcIvBPC4+0MfL2sosF4OiKSxDoXTeHvQ01yVLD5gQ7rIfWtcjUq1PWYYKh7VShHCVfQ9YjyipT03hdmCh0+kwGOp2Cx06dIjISGkFb/cu/RvcugYSRrZ1SUQHYvvy9qS3+WLwoAQWognsW5iDfb2tPWCHc6QHTLjPoT75eRMRoNanY7lSn4T7LPXJx3x+igyU+uQutwOLyy+/nCeeeILqarV1SaPRcPLkSR588EGmTZvW4gUUQrQ968nW235cvAQWwn32PRaA9ICJZrGNiVfrU2qkPwBH5UJQNIF9DxhIfWoKtwOLF154gZKSEqKioigvL2fChAn06NGDwMBA/vrXv7ZGGYUQbcx6MejlQU/ruHg50Qr32Y9hBugZZQkspD4J91mzQnmr9cmSFvRoTmm9rxGiPvY9YACpkVKf3OV2Vqjg4GDWrFnDxo0b2blzJyUlJQwbNoxJkya1RvlES9Knw2sjwD8S7tkJGllJUrjGvseie4R6os0prqS0sgZ/ndunEdGF2U/eBuhubhFMy5MvbuE+65h4L0sLs3p+OiJD60QT1O0Bk/rkLrevCP7+979z3XXXMXbsWMaOHWvdXlVVxapVq5gxY0aLFlC0oPJCqC5TbxJUCDfYf3kH+3kT5q+loLSKtLxSBsQFt3HpREdi3/sF0D1CAgvRdPbzvwBSLT0WMnRFNEHtHjCpT+5zeyjUrFmz0OvrLo5VXFzMrFmzWqRQopWUF6o/fUPbthyiw6l9sk2Ri0HRRNahUJa6ZO6xOFlQRrXR1GblEh1TpbU31XFM/PH8UmqkPgk3VdbpAZP65C63AwtFUdA4ae0+ffo0wcHSctmuSWAhmqj2cANLYHFcAgvhptp1KTrQB19vT2pMCqcLy9uyaKIDqrD2gKmBamywLz7eHlQbFU5JfRJuqqjVAxYb7Iuvt6fUJze4PBRq6NChaDQaNBoNEydOxMvL9lKj0UhaWhoXXXRRqxRStBBrYBHWtuUQHY70WIiWYpu8rQYWHh4aUiL82Zdp4FhuibVuCeGKilo9Fh4eGrpHBLAv08CRHKlPwj21e8A8PDR0j/Rnb4bUJ1e5HFhMnToVgB07djB58mQCAgKsz2m1WpKTkyXdbHtXXqD+9A3FaFL4y792k1tcxYvXDSbQx7ttyybaLZNJoarWxaBlXPwxCSyEm2oHqaAOh9qXaZBAVbit9hwLUBMCWAJViG6jkomOqHYPGED3yAD2Zkh9cpXLgcWSJUsASE5O5rrrrsPHx6fVCiVaid1QqF8O5fLpH6cA+PSPk8wZn9qGBRPtWZXduFLLl3dyhG3cqRDuqD15GyRQFU1XUSvLGKgXgiA9qsJ9tXvAwO78lCv1yRVuZ4W6+eabAdi6dSv79+8HoH///gwdOrRlSyZaXlAcJIyC8FTWH8q1bv7lUJ4EFqJeljHxYPvyTg5XT7RFZdUUllYR6q9tk7KJjqeyVp54sBtaJ1/cwk3WOTt2PRaWCbdyISjcZVvA07EHDOBYnmSGcoXbgUVOTg7XX38969atIyQkBICioiLOO+88Vq1aRWRkZEuXUbSUUberN+Dwu79bN+9O19c7KV8Iy9AVLw8NXp5qYOGr9SQ22IcMfQXH8koZLoGFcJEtK5StRVDm7Iimsg6tsx+6Yl5rRy4Ehbuc9YBZ1rKQQNU1bmeFmj9/PsXFxezdu5eCggIKCgrYs2cPBoOBu+++uzXKKFqBfcugvrya3JLKNiyNaM9qZ/GxSJGFzUQTOKtPlsAiy1BBaWVNm5RLdEy2Hgu7+mQ+N+WVVKEvr26TcomOybZAXt0e1fzSKvRlUp8a43Zg8d133/HGG2/Qt29f67Z+/frx+uuv8+2337Zo4UTrKK8ykqGvACDIR+20OpItLTvCOWeTbcE2HEpSzgp32LJC2epTiJ+WMHOvlwSqwh3WMfF29SlA50V0kA7APOFWCNc46wHzt6tPR6UXrFFuBxYmkwlv77oZhLy9vTGZZPGQdu3V4bC8H4XHdwLq5KThSeqaFicKytqyZKIdq7fHIkLGnQr3OZu8DTIcSjSNLSuUY32yDoeS4SvCDc56wEDqkzvcDizOP/987rnnHjIyMqzb0tPTWbhwIRMnTmzRwokWZsgAQzqF1WovRbi/jrhQXwAyimThF+FcfT0W3a1DoSQoFa5RFKXe+iSBhWgKazKAes5P0vAh3GFbx6Ke+iQ9YI1yO7B47bXXMBgMJCcnk5qaSmpqKikpKRgMBl599dXWKKNoCdUVUK1eAOYZ/QAI9fcmNkQNLNIlsBD1sLTgaOu0MKstOMfzSjGZlDNeLtHxVBsVFHNVqd0iKIGFcJfRpFjTYfvUOj91lwm3wk0mqU8twu2sUAkJCWzbto0ff/yRAwcOANC3b18mTZrU4oUTLaiiSP2p8SS3Sh0rGOqnJS5EeixEw+prwYkP9cXLQ0N5tZHs4gq6Bfu2RfFEB2IZtgJ1h0KlSougcJNl4U5oqIVZLgSFayrt6pP0gDWd24EFgEaj4YILLuCCCy5o6fKI1lJmWXU7hEJzVoMwf621xyKjqKKtSibauYp6xsR7e3qQGObHsbxS0nJLJbAQjbL/4tZ6Ou8BO5ZbKumvhUucrbFjYVnULC2/FKNJwdND6pNomH19qt1jkWrpoc8vk/rUCLeHQgH89NNPXHbZZdahUJdddhk//vhjS5dNtCS7VbcLyqoAtcfCElhk6stlOItwqr4eC7AbviIrcAsX2CcCqB04JIX7odFAcWWNpL8WLrE0etivsWMRH+qH1tODqhqT9MgLlzhbs8kiLtQXrZdan9ILpT41xO3A4o033uCiiy4iMDCQe+65h3vuuYegoCAuueQSXn/99dYoo2gJdoFFYakaWIT5a4kO1OGhUcc+55XKl7moq74eC5AVk4V7bKlm69YlH29P4s3JJGT4inBFQ40enh4aksLV+YRHZXidcEF9GRBBrU/Jlvokw6Ea5HZgsWzZMl588UU+/fRT7r77bu6++27+8Y9/8OKLL7Js2bLWKKNoCVo/SBgF0QMosAssvDw9CA9Q51zkGCSwEHU12GMhi+QJN9SXwcdCUjoKdzTU6AEyz0K4p6Km7uJ49uT85Bq3A4uioiIuuuiiOtsvvPBC9Hp9ixRKtILU8+GWH2DKSxSW2QILgKhANbDILZbAQtTVYI9FuAQWwnX1rTlgISkdhTsqGmj0ALtMPtLCLFzQUCMayPnJVW4HFpdffjlfffVVne3//ve/ueyyy1qkUKJ1WXosQv0cA4tsg0zgFnW50mNxsqCMaqMskCka5mzVbXu2C0EJVEXjKutZzMzCMoFbWpiFKxoaCgWSctZVbmeF6tevH3/9619Zt24do0ePBuD3339n48aN3HvvvbzyyivWfe++++6WK6kbXn/9dZ5//nmysrIYPHgwr776KiNHjmyTsrRH9lmhAKICfQDIkR4L4URDPRbRgT74entSXm3kdGG5dc6FEM40NMcCIFXWshBuqHA1UJULQeECa31qrMdCesAa5HZg8d577xEaGsq+ffvYt2+fdXtISAjvvfee9bFGo2mTwOKf//wnixYtYsWKFYwaNYqXXnqJyZMnc/DgQaKios54edqNr+6AtF8wTXycwjL1nyPU3xuA6CDzHIti6bEQdTXUY+HhoSE5wp/9mQbS8koksBANcrVF8GRBGVU1pjqLMgphz1Kf6htaZ1kbJctQQWllDf66JmXYF11EZWP1yTzHIttQSUllDQFSn5xy+6ikpaW1RjlazPLly7ntttuYNWsWACtWrOCbb77h/fff56GHHmrj0rUhQwYY0imrqrGufGsZChUZZO6xkMnbwonKxiZImgOLY7mlnN/nTJZMdDSNDYWKDtLhp/WkrMrIyYIyekQFnMniiQ7GUp986qlPIX5awvy1FJRWkZZXyoC44DNZPNHBVDTSoxrs5024v5b80irScksZGC/1yZlO1RxUVVXF1q1bHVYB9/DwYNKkSWzatKkNS9YOmNPNGjSBAAT6eOFtztNsmWMhQ6GEM41NkEyR4SvCRY2NiddoNNb6JBMkRWMqGqlPYDfPQs5PohG2HjDn33Ugw6Fc0akCi7y8PIxGI9HR0Q7bo6OjycrKqrN/ZWUlBoPB4dZpmQOLIkX9p7DMrwC7wEImbwsnrD0W9Xx5S2AhXNXYHAuQCdzCddahK/X0WIBk8hGua6wHDCTlrCs6VWDhrqeffprg4GDrLSEhoa2L1HrMgUW+yTy/ws8WWESbh0LlllSiKLL6tnBk7bGo52Qra1kIVzU2FArsM/nIhaBomPVCsKEeC5nALVzUWI8q2PdYSH2qT6cKLCIiIvD09CQ7O9the3Z2NjExMXX2f/jhh9Hr9dbbqVOnzlRRz6yaKqhSv6Rzq9WVI+17LCLMC+RVGxVrxighLBrrsbBcCGbqKyirqjlj5RIdT2PrWIAsaiZcZ0sG4EKgKkNXRCNc6rGwBqpSn+rTqQILrVbL8OHD+emnn6zbTCYTP/30kzU1rj2dTkdQUJDDrVOqKDLf0ZBdrQYU9j0WWi8Pa6AhmaFEbZYei/q+vEP8tIT6qRnGjueVnbFyiY6nspG6BJAqQ6GEi9zpsUjLLZUeedGgxrKMga3hIy1P6lN9mhRY/Prrr9x0002MHj2a9PR0AD766CM2bNjQooVrikWLFvHOO+/w4Ycfsn//fu68805KS0utWaK6JGMVJIyCuGEUlKn/OOEBWoddbIvkyQRu4ciVk22yzLMQLmhoTRQLy5ydgtIqisqqzki5RMfkymTbxDA/PD00lFYZ5ftNNMiWDKDx+lRWZSRL5qU65XZg8cUXXzB58mR8fX3Zvn07lZXqP6per2fZsmUtXkB3XXfddfztb39j8eLFDBkyhB07dvDdd9/VmdDdpQTHwy0/wG0/U1CqDnWy77EAiLKmnJV/FOHIlS9vy8Xg8XwJLET9rD0WDQSp/jovYsznI+m1EA2x9abWX5+0Xh4khPoCMnxFNMw2FKr++uTt6UFimDqkXIZrOud2YPHUU0+xYsUK3nnnHby9va3bx44dy7Zt21q0cE01b948Tpw4QWVlJZs3b2bUqFFtXaR2o9DcAhjm7+2wXVLOivo0lm4W7CfcyolW1M86x6KBoVAg8yyEa2zzvxqrT+pwqKMSqIoGuNJjAZJgojFuBxYHDx5k/PjxdbYHBwdTVFTUEmUSraigVA0s6vRYmAOLXAksRC2uDIVKMafgS5MJkqIBrgSpgKxlIVzian2SC0HhCld6wMDW8HFUGj6ccjuwiImJ4ciRI3W2b9iwge7du7dIoUQL+30FvNAXfnzcrseivjkWMhRKOKpwIVe8rGUhXFFe1XiQCpIiVLjGlhVK6pNoPlvWOtd6wGSopnNuBxa33XYb99xzD5s3b0aj0ZCRkcEnn3zCfffdx5133tkaZRTNVZIFxRlQXWHrsagVWFjWspDAQtRWUdN4q2ByhDrmtLCs2lrHhKitwuWhK5IiVDTOlXMTSH0SrpEesJbh5e4LHnroIUwmExMnTqSsrIzx48ej0+m47777mD9/fmuUUTRXWT4ANT6hFFeo6wyE1RoKFR1sCSxkKJSwqTaaMJrUlHoNtTL7ab2ID/XldGE5R3JKGJkSdqaKKDoQVxIBAPSwpAjNK6XaaMLbs1NlRhctxFKffF0MLE4XllNRbWy0/omuyfX6pJ6f0oukPjnj9tlao9Hw6KOPUlBQwJ49e/j999/Jzc3lySefbI3yiZZQVgBAuVcwABoNBPk6Tt7uFmzrsTCZJDezUFlOtND4xWDPKPVkezinuFXLJDou2yruDX/1xIX44qf1pNqocEIyjYl6WC8EtQ3Xp8gAHYE6LxQFTuTLWjvCOVfmEwJEBGgJ9FHrk2RCrKvJzUBarZZ+/foxcuRIAgICWrJMoqWZAwu9JhBQeys8PTQOu0QG6PDQQI1JIa9Uei2EynIhCI2PY+4Zrdavw9nSPSycc7XHwsNDYwtUpT6Jetjm7DRcnzQajV2mMalPwrlyF3ss1Pok83bq43ZgUVpaymOPPcaYMWPo0aMH3bt3d7iJdqhcDSyKMAcWteZXAHh5ehBpnsCdpZd5FkJl34Kj0Wga3LeH+ULwSI58cQvnKl0cEw/QI0o9Xx2SwELUwzJnp7ELQZAJt6JxlsDCR9t4fUqVeRb1cnuOxa233sr69euZPn063bp1a/RiQ7QD5jkWeaYAoNxpYAEQE+xLtqGSLH0Fg+LPYPlEu+VqlgywDYU6lC1DoYRzrg41AOgVba5PMrRO1KO8yvVANdWSIlQaPkQ9LPXJtUBV1tqpj9uBxbfffss333zD2LFjW6M8ojVE9gHfULKr/YFywgOcBxbdgnzYCbJMvbCyjYl3pYVZvRDMKa5EX1ZNsJ93I68QXY2rQ6EAelmH1klgIZxzdbIt2PWASaAq6uFefbLMKZRAtTa3h0KFhoYSFiYZXzqUmf+FeVtIN4YAzodCAcSYJ3BnylAoYVbuRgtzoI+3NQnAkVz58hZ1uROo9ox2zAwlhD1FUWxj4l0YutI7xjYHzCgJSoQT7tQna8NHTrHUp1rcDiyefPJJFi9eTFmZZFboaPJLLYvj6Zw+bwksZI6FsHCnhRnsWnFkXLyoRVEU65h4VwLVuBBf/CUzlKhHtVGxS4Xd+PkpMcwPnZcHlTUmThbI9Ytw5JhavfH6lBTuj9bLg4pqE6ekPjlweyjUCy+8wNGjR4mOjiY5ORlvb8fhDtu2bWuxwomWlV+iBhbh9fRYdJPAQtRiaWFubEEzi55Rgfx6OE+6h0UdVUYTirlhz5X6pNFo6BEdyM5TRRzKLrEOZRECbK3L4NrQFU8PDT2jA9iTbuBQdjEp5sm3QoBjfXKl4cPTnLlub4aBg9nFJEt9snI7sJg6dWorFEO0msM/wtfzIGksBaW3AtQ7xyLGvPq2zLEQFmVV6oKKfq4GFtEy7lQ4Z5+62JUvblATAqiBRTGXDOzWWkUTHZClN9XTQ4O3p2tJZHpFB6qBRVYxk/vHtGbxRAdTYU5d7KEBrYsLcvaODmRvhtSn2twOLJYsWdIa5RCtxXAaijOhqoR88/oUjc+xKEdRFMn4JSitVE+2/jrXThWWzFBHZMKtqKXSfCGoceOL25IZSobWidqswzS9Gk+FbWEZF39Qzk+iFvs1LFyuTzFSn5xp8gJ5ooMozlZ/BkRTUGoZCuV8jkW0uceiotqEobzmjBRPtG+llWo9CNC5N8ciQ19BcUV1q5VLdDxN+eK2LLooKYxFbe5MtLXoLfVJ1KM59UkaPhy5HVgYjUb+9re/MXLkSGJiYggLC3O4iXamOBMAk380ReXqhV59PRY+3p7W5zL05WemfKJdKzEHFq72WIT4aa0LLR6V/N7Cjrt1CWwtzGl5pVTVSGYoYePqqtv2LC3Mx3KlPglHzalPR3NLpD7ZcTuwWLp0KcuXL+e6665Dr9ezaNEirrrqKjw8PHj88cdboYiiWUrUHguDdziKAl4emnonbwPEhqi9FumFElgI+x4L1y8Ge1ozQ0mroLCxDKtzpy7FBvsQoPOixiSZoYSjcjfWHLCwr0/HpT4JO1KfWo7bgcUnn3zCO++8w7333ouXlxc33HAD7777LosXL+b3339vjTKK5ijOAiDfQ+1Nig7ywcOj/mEIiWF+AJKOTwBQWtX0VmaZwC3slVp7LFz/4tZoNNbhdYdkuIGwY2lhdmfoikajsc7bOZglDR/CRupTy3E7sMjKymLgwIEABAQEoNfrAbjsssv45ptvWrZ0ovnMgUWWeXG86CDn8yssEkLVwOJUoQQWAkrcnLwN9mtZyIlW2FiHQmndyxli+eKWcfHCXkkTelPBtlCe1Cdhr6n1qZfM26nD7cAiPj6ezEx13H5qaio//PADAFu2bEGna/iiVZxhigKRvSC8B6eqgwDoFuzb4EsSzD0WsuCLAPcnb4PdUCjpsRB2mjKsDtS1UUBd4VYIC0OFWp8CfZpWn6SFWdgrrmheYCH1ycbtwOLKK6/kp59+AmD+/Pk89thj9OzZkxkzZjB79uwWL6BoBo0GZvwb5m8lrSoYsGV+qk+CDIUSdpoz4fZ0Ybn1YlKIptQlsK2NIkOhhL0S64WgdyN7OpIeC+FMsTVQlfrUXG6vY/HMM89Y71933XUkJSXx22+/0bNnT6ZMmdKihRMtJ9O8mnZMcMO9SonWHgtZy0LYj4t3/VQR6q8lIkBLXkkVR3NLGBQf0kqlEx2Ju2uiWFgC1ePmzFBaL8mSLrCms3a3x8JSn04UlFFRbXQrC5DovEoqm1+fyquMbs3R6KzcOkNXV1cze/Zs0tLSrNvOPvtsFi1aJEFFe6Qo1ruW1bRjGhkKFRfii0ajZkjIK6lq1eKJ9q+5w1eklVlYWBIBuPvF3S3Yh0Bz5pW0PMm8IlSWHjB361NEgJYwfy2KAkdkuKYwK27i0Dr7+nQ0V+oTuBlYeHt788UXX7RWWURL27Acnu8B6561po+NDW54KJTWy4Nu5uFSMoFbWMYxB7nZPWwZviLj4oVFUydvazQaesgEblFLUy8EJZOPcKZE6lOLcbtPeerUqfzrX/9qhaKIFld4HEpzqTHWWBe8S47wb/RlMoFbABhNCkVlaq9VqL+7gYWsSCoc6cvUoQZBvm6PwKWXZQK3BBbCrKlj4kEy+Yi6DE2cswOyonttbp/he/bsyRNPPMHGjRsZPnw4/v6OF6p33313ixVONFPhCQDyvbuhKOCv9WxwcTyLpHA/NqcVcExWTu7SDOXVmMyj6UJ8G6839myZoeREK1RF5eYg1c+9ugQygVvUZZlj4e4wTbBr+JChUMKsqXMswLYC9wHpsQCaEFi89957hISEsHXrVrZu3erwnEajkcCiPck/CsApogFICvd3aTJ2aqT6JX5MxjN3aYXm3opAnZfbE2YtLYKnCsopq6rBz83hL6LzKSxVv7hD/JrRwiyBqjDTl6v1Kdi3CfUpSobWCUdFZc2oT+bzk8zZUbn9bW8/cVu0YxUGMJwGYL8xHjhNcoSfSy+1BBZH5Z+kS7MEFiFuDoMCCPPXEu6vJb+0iqM5pQyMD27p4okOxjKsLqQJPRbWzCv5ZVTWGNF5SeaVrs6SXCQ8oCk9YLaU2NLwIQAKStX6FObCqI7aepivmdKL1BTr7ma+62yalbdPURQUu8xDoh3JPaj+DIhhf5H6Z04Ob3x+BUBqlKXHogSTSf6+XZWlhbkpQ1fAfviKtAoKKCyz1Cf3A9XoIB2BPl4YTYoM0RSYTIq14SPc3/2FeS0NHwBHc6Q+dXWKYqtPTQks1BTraj2UzFBNDCzee+89BgwYgI+PDz4+PgwYMIB33323pcsmmiN3v/ozqi/7Mw0A9OkW5NJLE0J98fbUUFFtsk76Fl1PQVnTx8SDXcpZGb7S5VVUGymvVtexaEqPhZp5RSZICpWhohqjudHL3cQSFtLwISwMFTVUG9X61JTAAuzmFco8MPcDi8WLF3PPPfcwZcoUVq9ezerVq5kyZQoLFy5k8eLFrVFG0RS+oZAyAVPC2dYUaH3NE4wa4+XpYe3dOCqtg11Wc7qGAWsKviNyou3yLK2Bnh4aAps4TMBSn+SLW+SX2uZ/NXVYnKXhQyZwC8t3nb/Ws8kLJtpSrEt9cvsM/+abb/LOO+9www03WLddfvnlDBo0iPnz5/PEE0+0aAFFE/WdAn2ncDKvlLLv16H18iDFhVSzFqmRARzOKeFoTgkTekW2YkFFe5VlXq09OqjhtU/q00N6LISZpS5FBerw8Gg8gYQz1vokLcxdXn4z5ldYWC4Ej8j5qcsrKK0EIKwZ9alHlNQnC7d7LKqrqxkxYkSd7cOHD6empqZFCiVazj7zMKje0YF4ebr+506NsvRYSPTdVWWah8HFhjQtsLC0MFsmSIquK9vQvCAV7HospEWwy7NeCDaxNxVsF4KSwlhkG9T6FBng/nwdix5Rcn6ycDuwmD59Om+++Wad7W+//TY33nhjixRKNFN5EZTmA7DleAEAQxJC3HoLS2YoSZ/WdWWaW5m7Bfs26fXhATrC/LUoikyQ7OosPRYxzQosLJmhSqkwz9cQXVNGkbk+BTe/Pp0qLKO8SupTV5ZRZGlEa9p3HdiG1p0sKOvy5yeXhkItWrTIel+j0fDuu+/yww8/cPbZZwOwefNmTp48yYwZM1qnlMI9Oz6B7x+BodPZcvI6AM5KCXPrLSz/JBJYdF22wKLpX949owLYnFbA4ZxiSTnbhWWZWwSbcyEYFagjyMcLQ0UNx3JL6RfrWjIK0fmcLlQvBONDXUuh7ky4v5ZQP28Ky6o5mlvCgDg5P3VVlkA1rhmBRUSAlhA/b4rM9al/bNetTy71WGzfvt162717N8OHDycyMpKjR49y9OhRIiIiGDZsGHv37m3t8gpXHF4DQEVoT/ZlqEOhRia7F1j0iApAo1EnyeWVVLZ4EUX7VlFttP7dmxVYyIrJAsgyD6trzlAo+8xQsqJ715ZeVAY070JQo9HYTeCW+tSVWXosmvNdp9YnGekBLvZYrF27trXLIVpKVSmc2AjAH57DMSl6ksP93G4p9NV6khTmx/H8Mg5lFRPRo+ljD0XHczy/FEVRVyFtzjhm24qk8sXdlR3PVy8EE8Oa3sIM6sJmf54olAncXVx6kaXHoumBBagNH38cL5BMY11chr75Q6FATTCx5Xhhl69PzVogT7RDR9eCsQpCEvnXKXUC9qS+0U16K8tF4UH5Eu9yLC0uqZH+aDRNy+IDMkFSqItPHTMngege6XpmOmd6SQ9Yl6coCifNgWpzhkKBbS6hJCnpuhRFIc2cVj/JxUWE62Ndy6KLN6S5nW62oqKCV199lbVr15KTk4PJZHJ4ftu2bS1WONEEu1YBYOxzOT//kQvABf2aFlj0jgnkh33Z1nUwRNdhmWxt+eJtqtoTJH21TcsRLjqugtIqDBU1aDS4lfLaGcvrT+RLMoCuKttQiaGiBk8PDckRzQssLK8/YQ5URNeTZaiguFKtTy13fura9cntwOKWW27hhx9+4Oqrr2bkyJHNas0ULaysAA5+B8CmwAsoKislIkDH8KTQJr1d7xjpseiq9mboAVtg0FTh/lqCfb3Rl1dzoqCUPjEy4barsaRfjA32bfLiUxaWFsWTBWUoiiLfP13QgSx13mBKhH+TF8ezkPokLMOWksP90Ho1bxBPYrgaqHb1+uR2YPHf//6X//3vf4wdO7Y1yiOaY+cqMFVDzCDePeQHlHL18Hi31q+w19t8UXkoq7hL/5N0NYqisP1UEQBDE0Oa9V4ajYbEMD92p+s5mV8mgUUXtNNclwa2QNaduBBfPDRQUW0ip7iyWZPBRcd0wNyDbhkW1xzxob5oNFBWZSS3pJKoQKlPXc1ec4Kb5jaigVqfPKQ+uT/HIi4ujsDA5v8B3JWcnIxGo3G4PfPMMw777Nq1i3HjxuHj40NCQgLPPffcGS9nmxp5G0xdQcawe1l/SB0Gdd1ZCU1+u+QIf7w9NZRWGa2T5UTnd7qwnNziSrw8NC2SgtG+FUd0PTtPFwEw2M21dJzRenlYJ1h29eEGXdWfxwsB99dmckbn5UmseZ2ek1KfuqStJ9T6NCyxaSM77Om8PK3rPnXl+uR2YPHCCy/w4IMPcuLEidYoT4OeeOIJMjMzrbf58+dbnzMYDFx44YUkJSWxdetWnn/+eR5//HHefvvtM17ONuPpDUNu4LljSSiKOreiOWMGvT09rGPsJQtL1/Hr4TwABsYHN3voCtgyAUlg0fUoitKiF4IAyeEyz6KrMpkUtp5QF30d4WYK9fpY5lkc78IXgl2VoihsO6men4YnNz+wAEgKl/rk9lCoESNGUFFRQffu3fHz88Pb29vh+YKCghYrXG2BgYHExMQ4fe6TTz6hqqqK999/H61WS//+/dmxYwfLly9nzpw5rVamdqGsALx0oPVnT7qer3dmAHDPxJ7Nfute0YEcyCrmQFYx5/dp2iRw0bH8tD8baHo2sdqSwmSCZFe1N8NATnElvt6ezR5WZ5EY7gdHJFDtivZlGigsq8bX25MBLbQAWWKYPxvJ56QEql3OnnQDBaVV+Gk96d9CC24mhfvx29GuXZ/cDixuuOEG0tPTWbZsGdHR0Wd03P0zzzzDk08+SWJiIv/3f//HwoUL8fJSf4VNmzYxfvx4tFpbzv3Jkyfz7LPPUlhYSGhoy0Sj7Y6iwL/uhKzdGKe8xqM/6DApcNmgbi0yjKV3TCDsVOdZiM6vsLSKX4+oPRYT+0a1yHtaeixOyYVgl/PT/hwAzukZ0SK9X2ALVLtyi2BX9e2eTADO7R3Z7Im2FsnSwtxlfb83C4AJvSKbnQjAIjHM3KPahb/v3A4sfvvtNzZt2sTgwYNbozz1uvvuuxk2bBhhYWH89ttvPPzww2RmZrJ8+XIAsrKySElJcXhNdHS09TlngUVlZSWVlbZVpQ0GQyv+Bq1k81tw6Dvw1PH5wUp2nionQOfFY5f1a5G3t61lIXm+u4LVW09RVWOif2yQdfJ+c1nmWJwuLMdoUvD0kCQAXYHJpPD5tlMAXNjElNfOWIYadOUWwa7IaFKsvfEXDXA+cqEpkmRoXZdkNCn8e2c6AJP7t2R9kh56t0P+Pn36UF7eMhN5H3rooToTsmvfDhw4AMCiRYs499xzGTRoEHfccQcvvPACr776qkNg4K6nn36a4OBg6y0hoekTndvE4TXw/SMAnBj2II/8pgCw+LJ+LZYtpY855ezRnBJqjKZG9hYdWWllDW//kgbA9LOTWqw3sluwL14eGqqMJrIMFS3ynqL9W38ol1MF5QT6eHHZoNgWe19pEeyaftqfzamCcoJ9vZu8NpMzCWHqZFtJUNK1/Hwgx1qfWiOw6MpDNd0OLJ555hnuvfde1q1bR35+PgaDweHmjnvvvZf9+/c3eOvevbvT144aNYqamhqOHz8OQExMDNnZ2Q77WB7XNy/j4YcfRq/XW2+nTp1yq/xt6uRmWD0TFCOG3ldzxZ8DMJoUrhwaxzUj4lvsY+JCfPHTelJlNElXcSf3yk+HySupJDHMj6uGtVwd8vTQEB9qyeQjrYJdgcmk8Pz3BwG4bkRCiy6MaOkBKyqrpriiusXeV7RfRpPCSz8eBuCGkYn4ad0ebFGvOHOWsbySKiqqjS32vqL9MpoUlq85BMD1I1v2/JRgHqpZUFpFaWVNi71vR+L2f+dFF10EwMSJEx22W9Y5MBpd/8eMjIwkMjLS3SIAsGPHDjw8PIiKUseBjx49mkcffZTq6mrrhPI1a9bQu3fveudX6HQ6dDpdkz6/TR1bD5/eANWlVCSM47K0aygqr2FwQgh/vXJAi8578fDQ0DM6kJ2nijiYVUyPqObnDhftz9oDObz96zEA/nJp3xYbv2yREObH8fwyTheUQ2qLvrVoh1b+dpx9mQYCdV7cdV6PFn3vAJ0XoX7eFJZVc6qgnH6x3o2/SHRoH2xMY1+mgSAfL24bl9L4C9wQ7OuNn9aTsiojmfqKZq++LNq/Dzamsd9cn24f37JfSEE+3tZFYU8Vds21m9wOLNauXdsa5WjQpk2b2Lx5M+eddx6BgYFs2rSJhQsXctNNN1mDhv/7v/9j6dKl3HLLLTz44IPs2bOHl19+mRdffPGMl7fVbXkXqkspS5jAxVl3cLLYSPdIfz6YeVaLtuRY9I4OUAOL7GIupVuLv79oWxsO53HnJ1tRFLU18MIW7Ba2sLTinC6UXq/ObvOxfJ75Vh3C+sBFvQnz1zbyCvclhPlRWKbnVGEZ/Voom4tonxzrUx/CA1q2MVCj0RAb4suRnBIyisolsOjk/kgr4NnvbPWpNc5PXX1RWLevQidMmNAa5WiQTqdj1apVPP7441RWVpKSksLChQtZtGiRdZ/g4GB++OEH5s6dy/Dhw4mIiGDx4sWdM9XsFa+R7RHJlH0TySlX6BkVwEe3jGqVfxCwTeCWzFCdS7XRxDu/HuOFHw5hNCmc3yeKpZf3b5XPSgg1Z4YqlHHMndmGw3nc/tGfVBlNXNQ/hpvOTmqVz0kI82PXab1kGuvkfjuSx61//5Mak8Jlg7px46jEVvkcS2Ah8yw6t9+O5HHb3/+k2qhw6cDWq08JYb7sTtd32e+7JjVv//rrr7z11lscO3aM1atXExcXx0cffURKSgrnnHNOS5eRYcOG8fvvvze636BBg/j1119b/PPbXEEa7PgHnPcICvCPnUU8vmMS1UaFwfHBrJw1ktBWCirAnHIWWSSvszCZFH46kMMLPxzkgDlYvGpYHMuuHNjiQ6AsLHMs5EKwc6qsMfLW+mO89OMhTAqMSQ3npeuHtFo6cmugKvWpU6o2mlix7igv/XQYo0lhXM8Inr96cKvVp7gQNdlJhgQWnVJVjYm3fznK8jXq+Wlsj3BeuLb16lNXPz+5HVh88cUXTJ8+nRtvvJFt27ZZszLp9XqWLVvG//73vxYvZJelKLDt72rmp6oSqoMSeOzEEFZtUSeZXzqwG89dPQh/XcsPf7JnSTt6PL+Uimpji+WjF2dWtqGC/+zMYNWWUxzJUdMHh/p58/AlfblmeHyrrkljGQp1SoZCdSpGk8K3ezJZvuYQx3LViflXD4/nqakDWvU8Ycnk01VbBDsrk0lhzf5snv3ugLU+XTk0jqevGtiq9Sk2WK1PElh0LiaTwg/7snn++wMcNdcnSyNa656fuvbQX7evSJ966ilWrFjBjBkzWLVqlXX72LFjeeqpp1q0cF1acTZ8PR8Ofw9AacxIbl+rZUP+KTQaeGByH+6Y0P2MLFAYGaizTpY8klPSIgvvidZnMinsyzTw6+E81h/KYXNaAYqakZhAnRc3np3E7eO7t2pvl0WCucci21BJZY2xxRYjEm0jS1/B1zvT+WTzSWu+9ogAHY9d1pfLB8e2+nmpq7cIdjb5JZV8vTODD387bs0+GO6v5ZFL+nLVsLhWr0+xIZJytjPJK6nkvzsz+PumExzLUwOKM1mfrA1pBV2zPrkdWBw8eJDx48fX2R4cHExRUVFLlEns/Rf8dyGUF6B4atmQeCe3HBxJlUlDVKCOF64dzLieTcum1RQajYZe0YFsTivgYFaxBBbtkKIoZOor2JdhYF+mgb0ZerYcL6SgtMphvxFJoVwxJJapQ+MI9Dlz2XTC/LX4entSXm0kvbCc7pGSXayjScsrZf3BHNbsz+a3o/nWIDXUz5uZY1KYOTaZYN8zU6fse8AsGQlFx5KpL2fdwVy+2ZXJb0fzMJnrU5CP2uhx57mpBJ2hc5QlsMgoknV2Oqr0onJ+OZTLd3uy2HAkD6O5QgX5eDF9dBJzxqeeufOTuSHtZEHXPD+5HVjExMRw5MgRkpOTHbZv2LCh3jUnhBvWLIGNLwFQGdGf+41z+Xp/CKAOfXpq6oAz0sJcW+8YNbCQeRZtr9po4mhuiRpEmAOJfZkGisrq5vT313oyOjWccT0jOb9PlPWC7EzTaDQkhPlyKLuEUxJYtHsmk8Lx/FJ2p+vZdqKQ9Ydy66xjc1ZyKFOHxnHl0LhWyUbXkLgQXzQaqKg2kVtSSVRgyywIKlqPvryaP9IK2Hgkj18P51qHplgMjAvmmhHxTBsW3+rDe2uLs+ux6IoXgh2NoihkGyrZdbqIP9IKWH8ol8Pm4b0WgxNCuHJILFePSCDgTNenUPX8VF5tJL+0iogWzmTW3rl9tG+77Tbuuece3n//fTQaDRkZGWzatIn77ruPxx57rDXK2LX0moyy6XX2JM/kxsPnYqjWEOjjxZNXDOCKIa0/xKDeYpnnWRyUwOKMKq6o5kBWsUMQcTC7mKqauquge3po6BkVQL9uQfSLDWJQfAhDE0Pw9mydCdnuSgj141B2SZcdd9peKYrC6cJydqfr2Xm6iN2n9exO11Nc4bi4k7enhhFJYZzbO5JLBnZrsyAVQOvlQbcgHzL0FZwqKJfAop0praxhb4aBXaeL2J2uZ/dpvXVIioWHBgbGhzC5fzSXDuxGUnjbpXmNDtah0aiTfLvihWB7l1dSye7Tenad1rPrdBG70vXkFlc67OOhgaGJoUzoFcmUwbFtmjZY5+VJTJAPmfoKThWUdbn65HZg8dBDD2EymZg4cSJlZWWMHz8enU7Hfffdx/z581ujjF1KTugwnu72IV/tUx+PSQ3nb9cMtnbVthVrZihJOdsqLC0w+zL17MswsNccRJyoZ7XzAJ0XfbsFWoOI/rHB9IgKaNcT67v6uNP2IttQwa7TenafLmKnOYioPWQOQOflQX9zgDo6NZyxPSLOeMtfQ+LD/MjQV3C6sIzhSc4XQRWtr7zKyL5MA7vNF3y7T+s5kltiHSpnLzncj3N6RnBOjwhGd48g2K99LG6o8/IkMkBHTnElGUXlXe5CsD3Rl1WzO13PrvQidp1Sz0/O5r5YGtKGJIQwrmck5/RoP/UJ1Ia0TH0FpwrLGZrYtc5Pbn9LaDQaHn30Ue6//36OHDlCSUkJ/fr1IyBAhjY019oDOSz6bAeFZWqL3IMX9WHWmGQ8PNq+W7ZXlBpYZOgr0JdVt6t/4I6mxmjiWF6pbRiT+aezizuAbsE+1gDC8jMh1K9d1At3WFPOSo/FGVNQWqW2Gp/Wm4OIIrINlXX28/bU0CcmiIHxwQyOD2ZgXAi9ogPwaie9Xc4khPrxR1qBTOA+gyprjBzILDYHEEXsOq3ncE6JdTy7vW7BPgyMC2ZQfDAD40MYGBfcamsttYS4UF9yiis5XVjOoPiQti5Ol1BSWcPedHNPhLlO1R5yCaDRQPcIfwbFhzAoXq1T/boF46ttvw1p8WG+/HG8ayaYaHLzk1arpV+/fi1Zli7PV+tJUXk1/WODePG6IdbhR+1BsJ838aG+nC4sZ2+GnjE9Itq6SB1CaWUNB7IMDkHEgaxiKusZypQa6W8XRATTLzaoXX8ZuyM+1JKCT3osWoOhopo91i9odViTs2PtoYGeUYHWL+hB8SH0jgls171dziSae8BOdsEv7jOhqsbEoexitfXYHJQezCqm2lg3iIgI0KkBqblODYgL7nDD0+JD/dh+skiGaraSimq1Z2vXKbVna9dpPUfr6dlKDPNzaOQYEBd0RpONtISunLnO5cBi9uzZLu33/vvvN7kwXd3Z3cP5YOZZjEmNaLWFyppjUHwwpwvL2XlaAgtnisqq2JthYE+6nj0ZBvam60nLL3V64vTXetK3Vi9Er+iOd3HnDsvaA6e74Im2pZVV1bAvw6D2QphbjmuPYbfoHuFvvuBTW/v6xwad8cnWrcG6loUMrWu2GqOJI7kl5iFyanC6P9PgdC5XqJ+3tS6pPRIhRAfpOvyE54RQqU8tparGxMGsYnal23pLD2UX19uzZWngsPRwhfh1/Ma0xC68dpPL3y4rV64kKSmJoUOHoji7UhIt4tzeUW1dhHoNig/hf7uz2J1e1NZFaXM5hgq7IELPnnRDvTnQo4N0dXohksI63lCm5rL0WOSXVlFaWXPGM790VPbDT3adUifDHsouxsl3NPGhvtYv6UFxwQyIDz5jKTvPNFl0sWmMJoW0vBLzRFh1DPveDD0V1XWDiCAfL/WCLz6YQXFqj4SakavznbukPjWNfVBqGXa5P7OYKmPd+hQRoLUGEIMTOmbPlqu68pxCl7/Z77zzTj799FPS0tKYNWsWN910E2FhYa1ZNtHODDKvX7HrtL6NS3LmWDLm7M3QO/RG1M5IYZEU7kd/82TqAXFq67BMBFQF+3oT5OOFoaKG9KLydjXUr72oNqrDT/akm+dEnNZzIMvgdPhJVKDOYczxwLhgwrtQXbMMNcjUV1BjNLXr+SBtxWRSOFFQZr3g25WuZ2+6ntIqY519A3ReDIgLcmg5Tgzz65RBhDNdeeiKq0wmhbT8UjUzk/n8tDfDQHl13foU7Ottd25Sz1Pdgn26Tn0Ks63m3tXOTy4HFq+//jrLly/nyy+/5P333+fhhx/m0ksv5ZZbbuHCCy/sMpWlKxsQrwYWpwvLySup7HQXzJaT5p50NTOTpSdCX153fQgPDaRGBliDh/6xak/EmVqAp6NKCPNjb4aBUwVlXT6wqDaaOJxdwp50vTkLSsPDTwbGh5jHHAczOCGE6KDO2dLnqqhAHVovD6pqTGTqK9o0/W17YGkEUSfC1p82GMDX25MBcUHWC76B8cGkhPt3uV5Ue9ahmoWylgXY6pMlBfWu03r2pOsprqxbn2oHpYPjQ0gI65w9W66KDvRB6+lBlbHrnZ/cGoug0+m44YYbuOGGGzhx4gQrV67krrvuoqamhr1790pmqE4uyMebPjGBHMgqZtPRfKYMjm3rIjVZtdHEkRz1om5vhrpS9b4Mg9OWPG9PdeXxAbHB9I9Tg4i+3QI7xTj1My0+1NcaWHQlNUYTh3NK2J2ufjnvOq0GEc4m8Qf6eDEg1jaxelB8MPGhXftL2hkPDw3xob4cyy3lZEFZl/riVhSFTH2FdVK1ZUiTs0UydV4e6ro2cWp2pkHxwaRGBuDZhYMIZ2JDfPHQQGVN11t0UVEU0ovK2ZNuMPeWFtVbn3y8PegfG2zt1RoUH0L3iK4dlDpjPT/llXKqsGudn5p8ZeTh4YFGo0FRFIzGuhdjonMa1zOCA1nF/HIot8MEFhXVRvZnGqwBxF5zZiZnLcM+3h707RakBhGxQQyIC6ZndAA6r847qfpMSugCmaEsY453m1v4LD0RzsawB+q86B8XxMA4W0rOrjj/pqkSQv04llvaqQNVRVE4WVCmXvSZz1970/XkO0lP7e2poW+3IFua17gQekYHtJtFMtszb08PugX7kl5UzqmCsk4bWFjqk9rIoX4n7knXU+gkiLCvT4PNc216RrXvNNTtSXyYH8fySjldUA6pbV2aM8etwKKystI6FGrDhg1cdtllvPbaa1x00UV4eEhF6wrG94rknV/TWH8oF6NJaXetXoaKanUYk91wpqO5pU6zUQTqvOhnDh4GmHsiukf4y0mzFXW2CZJlVTUcyCpmf6aB/XZrkjgLIgJ0XuYF59T5NwPjgknu4sNPmsuaGaqT1CfLGjd77C769mUYnA4/8fTQ0Ds60DqUaVBcCL1ipBGkObpH+pNeVM7BrBKGJ3X8OaTqRP1S9mao8yEsgamz4XFeHuaeefOQpkHxwfSOCZT61AwJXXTtJpcDi7vuuotVq1aRkJDA7Nmz+fTTT4mIkJSjXc3IlDBC/LzJKa7k18O5bZbFqsZo4nh+GYeyizmQVcyhrGL2Z9W/UnW4v5b+ccEMsE6s7piLzHV0lgvB+v5O7ZX9yuj7M4vZl2lgf4ahwXTC/eNsmXQGxMkY9tZgm3Db8XrACkqrOJhVzMEsAwezS6zBqbPhcVovD/rEBNLfrie1Twdce6S96xcbxK+H89iX2fESlBSUVnEgy8DBrGIOZBZzILuYw9nFlDkZ3qv18qBvTCD9zQ0cA2KDJShtBV11rR2XA4sVK1aQmJhI9+7dWb9+PevXr3e635dfftlihRPtj87LkyuHxvHBxuO8tyGt1QMLo0kho6icI7klHDQHEAeyijmSW+J0KBNAXIivXWYm9WdnyLPeGfTrpiYAOJRdTEllDQHtLOWsoihkGSo4nF3CkZwSjuSWcCS7hMM5xU6HCoC6OFi/2CD6dgukXzdbz5cEEa3P0gN2It/5Gh7tQUllDUdz1PPXwexi68/6Msv5az3pZz5/WYKIHlEynOlM6NctCID9mcVtXJL66curOZZbwtHcUg5mqcN6D2TVX598vdX6NCA2yNy4FizD484Q2/lJAgunZsyYIRdmAoCZY5L5+PcT/Ho4j3/vSOeKIXHNej9FUcgrqSItr5TjeaUcyyslLa9EfZxfVm8A4evtSa/oAHrHBNIrOpA+MUH0jw0itJOsVN0ZxQT7EBeijmPefrKQcT0jz3gZFEVBX17NifwyThaot2O5pRzJLeFoTgklToadgDr0pHuEvzmIsNwCO+1Y7I6gf6x6Ibg3w9CmgWpFtZFTBWXmc1cpabmlpOWr9+u74AO1B693dBC9YwLoHaNe/MnwuLYzwJxSfXe6nrKqmjZL0FFZY+R0YTnHcks5llui/sxTfzqbW2ORGOZH75hA+sSo34e9YwJJifBvd0OWu4q+5kB1X4ahS63d5NYCeUIAJIX7c/v4VF5be4T7Vu/kdGE5N4xMJMzJBb2iKBjKa8grrSSvuJIMfTnpheWkF1WQXlROemEZGUUVTvNgW2g9PUiO8DMHD7YgIj7UV76AO6DRqeF8vvU03+3JapXAwmhSyCupJEtfQZahgmxDhXVCpiWYcDbG2MLLQ0NSuB89ogLoERVAz6hA630ZetK+JIX7kxjmx8mCMjYczuWiAd1a5XNKKmvIKLKcu9Sb/eMsQ4XTIXEWEQE667nLEkT0jAroMhcaHUX3CH8Swnw5VVDOr4fzmNw/psU/Q1EUSipryDZUkl5UzunCMtILyzldqN4/XVhOTgPBKKiplrtH+tM7OpA+3YKsjWvtrQe4q0sO9yM+1JfTheVsONI69ak9kloommThBb04nl/Kf3dl8vz3B3n++4NEB+kI8VWDiyqjidLKGgpKq6hxtkRwLRqNOoQpJcKf7hH+pET4kxIZQPcIf2JDfKXFpRO5cmgcn289zb+2p3P7+FQSwxtOw6coCuXVRgpKqygsraagrIrC0ir1cZn6M7+kiixDBVn6CnJLKp1O1q8tOkhHYpgfCWF+pIT7W4OHpHB/tF4yTKCjuHhADG/9cowV648xsW90o0M8LPVJX16NobwGfXk1hWVV5JVUkldcRW5JBXnF6uPcErVBxFka6toCdV6kRKrnruRwf7pb7kf4d9rVzzsbjUbD5H4xvLshjRXrj3J+nyiX6lNZlVqf9OXVFJVZflaRW1xJTnElucVqXVIfVzhN7lCbr7en+n0Y6U/3yABSI/3pHhFAcoQfgVKfOgSNRsNF/dX69JYb9cn+/GSoUOtUQWklUUE+nNdG81rdoVGUhtpZuhaDwUBwcDB6vZ6goKC2Lk67pygKX2xL570NaezPNDS4b6DOi/AALbEhvsSG+BJnuYWqj7sF+0hrcBdhMilc89Ymtp4oJEDnxYjkUPy1XigoVFSbKKmoobiyhtLKGkoqayipqKHK2PgXsT0PDUQG6ogJ8iE6yIfYEF8Sw/xIDPMjKdyP+FA/fLVS3zqDbEMF5z6/jvJqI8nhfgyIC0br5UFljYnKaiMV1SaKK2swlFert4pqpyuZNybY19vu3OVjPXfFhfgSH+pHRIBWhgt3Aln6Cs7921oqqk2kRPir9cm80Fl5VQ1lVUbKqowUV1Rbg4mm1KdAnRexIb7Eh6q3uFC1HsWZt4X5S33qDDL15Zz3t3XW+tQ/NgidlyeVNeq5qaLa6PL56YJ+0bwzY8QZ/g1U7lwfS2BhRwKLpisqq+JkQRmG8ho0GjXrhK+3J+EBWsL8tZJtQjhILyrntg//ZF8jAak9racHYf5aQv21hPp5E+qvJcxPa/7pTUywDzHBvsQE+RARoJW0wV3I2gM53L1qe4ND3Grz8tAQ5OtNsK83Qb7eRAboiAzUEhmgIyJQR0SAjki7nzLMpOv4aX82C1btcJrmtz7enhqCzfUp2NebED+tuU7piArS2e4H+hARqJUFVruQptQnTw9bfQry8SLMX8vwpFDmnd+zFUtaPwksmkgCCyHOHJNJYdvJQo7klFBRbUSj0eDr7UmAjxf+Oi8CdF4Emu+H+Hrjp/WUFjxRL31ZNRuP5pGlr6DKaEJnbtzw8fYkQOdlF0R4EeQj9Uk0zFKfMorKqTYqaL088Nd64qv1xE+rnptsQYQ3vt5Sn0T9nJ2ffLw98fH2IEBnOzepgUT7Oz9JYNFEElgIIYQQQghh4871sYwVEEIIIYQQQjSbBBZCCCGEEEKIZpPAQgghhBBCCNFsElgIIYQQQgghmk0CCyGEEEIIIUSzSWAhhBBCCCGEaDZZocWOJfOuweD6ol1CCCGEEEJ0VpbrYldWqJDAwk5xcTEACQkJbVwSIYQQQggh2o/i4mKCg4Mb3EcWyLNjMpnIyMggMDCwzVY8NBgMJCQkcOrUKVmkr4nkGDaPHL/mk2PYfHIMm0+OYfPJMWw+OYbN19bHUFEUiouLiY2NxcOj4VkU0mNhx8PDg/j4+LYuBgBBQUHyD9hMcgybR45f88kxbD45hs0nx7D55Bg2nxzD5mvLY9hYT4WFTN4WQgghhBBCNJsEFkIIIYQQQohmk8CindHpdCxZsgSdTtfWRemw5Bg2jxy/5pNj2HxyDJtPjmHzyTFsPjmGzdeRjqFM3hZCCCGEEEI0m/RYCCGEEEIIIZpNAgshhBBCCCFEs0lgIYQQQgghhGg2CSzakddff53k5GR8fHwYNWoUf/zxR1sXqd365ZdfmDJlCrGxsWg0Gv71r385PK8oCosXL6Zbt274+voyadIkDh8+3DaFbaeefvppzjrrLAIDA4mKimLq1KkcPHjQYZ+Kigrmzp1LeHg4AQEBTJs2jezs7DYqcfvz5ptvMmjQIGtu8dGjR/Ptt99an5fj555nnnkGjUbDggULrNvkGDbs8ccfR6PRONz69OljfV6On2vS09O56aabCA8Px9fXl4EDB/Lnn39an5fvlIYlJyfXqYcajYa5c+cCUg9dYTQaeeyxx0hJScHX15fU1FSefPJJ7KdCd4R6KIFFO/HPf/6TRYsWsWTJErZt28bgwYOZPHkyOTk5bV20dqm0tJTBgwfz+uuvO33+ueee45VXXmHFihVs3rwZf39/Jk+eTEVFxRkuafu1fv165s6dy++//86aNWuorq7mwgsvpLS01LrPwoUL+c9//sPq1atZv349GRkZXHXVVW1Y6vYlPj6eZ555hq1bt/Lnn39y/vnnc8UVV7B3715Ajp87tmzZwltvvcWgQYMctssxbFz//v3JzMy03jZs2GB9To5f4woLCxk7dize3t58++237Nu3jxdeeIHQ0FDrPvKd0rAtW7Y41ME1a9YAcM011wBSD13x7LPP8uabb/Laa6+xf/9+nn32WZ577jleffVV6z4doh4qol0YOXKkMnfuXOtjo9GoxMbGKk8//XQblqpjAJSvvvrK+thkMikxMTHK888/b91WVFSk6HQ65dNPP22DEnYMOTk5CqCsX79eURT1mHl7eyurV6+27rN//34FUDZt2tRWxWz3QkNDlXfffVeOnxuKi4uVnj17KmvWrFEmTJig3HPPPYqiSB10xZIlS5TBgwc7fU6On2sefPBB5Zxzzqn3eflOcd8999yjpKamKiaTSeqhiy699FJl9uzZDtuuuuoq5cYbb1QUpePUQ+mxaAeqqqrYunUrkyZNsm7z8PBg0qRJbNq0qQ1L1jGlpaWRlZXlcDyDg4MZNWqUHM8G6PV6AMLCwgDYunUr1dXVDsexT58+JCYmynF0wmg0smrVKkpLSxk9erQcPzfMnTuXSy+91OFYgdRBVx0+fJjY2Fi6d+/OjTfeyMmTJwE5fq76+uuvGTFiBNdccw1RUVEMHTqUd955x/q8fKe4p6qqio8//pjZs2ej0WikHrpozJgx/PTTTxw6dAiAnTt3smHDBi6++GKg49RDr7YugIC8vDyMRiPR0dEO26Ojozlw4EAblarjysrKAnB6PC3PCUcmk4kFCxYwduxYBgwYAKjHUavVEhIS4rCvHEdHu3fvZvTo0VRUVBAQEMBXX31Fv3792LFjhxw/F6xatYpt27axZcuWOs9JHWzcqFGjWLlyJb179yYzM5OlS5cybtw49uzZI8fPRceOHePNN99k0aJFPPLII2zZsoW7774brVbLzTffLN8pbvrXv/5FUVERM2fOBOT/2FUPPfQQBoOBPn364OnpidFo5K9//Ss33ngj0HGubSSwEEIwd+5c9uzZ4zA2W7imd+/e7NixA71ez+eff87NN9/M+vXr27pYHcKpU6e45557WLNmDT4+Pm1dnA7J0poJMGjQIEaNGkVSUhKfffYZvr6+bViyjsNkMjFixAiWLVsGwNChQ9mzZw8rVqzg5ptvbuPSdTzvvfceF198MbGxsW1dlA7ls88+45NPPuEf//gH/fv3Z8eOHSxYsIDY2NgOVQ9lKFQ7EBERgaenZ50MCdnZ2cTExLRRqTouyzGT4+maefPm8d///pe1a9cSHx9v3R4TE0NVVRVFRUUO+8txdKTVaunRowfDhw/n6aefZvDgwbz88sty/FywdetWcnJyGDZsGF5eXnh5ebF+/XpeeeUVvLy8iI6OlmPoppCQEHr16sWRI0ekDrqoW7du9OvXz2Fb3759rUPK5DvFdSdOnODHH3/k1ltvtW6Teuia+++/n4ceeojrr7+egQMHMn36dBYuXMjTTz8NdJx6KIFFO6DVahk+fDg//fSTdZvJZOKnn35i9OjRbViyjiklJYWYmBiH42kwGNi8ebMcTzuKojBv3jy++uorfv75Z1JSUhyeHz58ON7e3g7H8eDBg5w8eVKOYwNMJhOVlZVy/FwwceJEdu/ezY4dO6y3ESNGcOONN1rvyzF0T0lJCUePHqVbt25SB100duzYOqm2Dx06RFJSEiDfKe744IMPiIqK4tJLL7Vuk3romrKyMjw8HC/LPT09MZlMQAeqh209e1yoVq1apeh0OmXlypXKvn37lDlz5ighISFKVlZWWxetXSouLla2b9+ubN++XQGU5cuXK9u3b1dOnDihKIqiPPPMM0pISIjy73//W9m1a5dyxRVXKCkpKUp5eXkbl7z9uPPOO5Xg4GBl3bp1SmZmpvVWVlZm3eeOO+5QEhMTlZ9//ln5888/ldGjRyujR49uw1K3Lw899JCyfv16JS0tTdm1a5fy0EMPKRqNRvnhhx8URZHj1xT2WaEURY5hY+69915l3bp1SlpamrJx40Zl0qRJSkREhJKTk6Moihw/V/zxxx+Kl5eX8te//lU5fPiw8sknnyh+fn7Kxx9/bN1HvlMaZzQalcTEROXBBx+s85zUw8bdfPPNSlxcnPLf//5XSUtLU7788kslIiJCeeCBB6z7dIR6KIFFO/Lqq68qiYmJilarVUaOHKn8/vvvbV2kdmvt2rUKUOd28803K4qipmV77LHHlOjoaEWn0ykTJ05UDh482LaFbmecHT9A+eCDD6z7lJeXK3fddZcSGhqq+Pn5KVdeeaWSmZnZdoVuZ2bPnq0kJSUpWq1WiYyMVCZOnGgNKhRFjl9T1A4s5Bg27LrrrlO6deumaLVaJS4uTrnuuuuUI0eOWJ+X4+ea//znP8qAAQMUnU6n9OnTR3n77bcdnpfvlMZ9//33CuD0uEg9bJzBYFDuueceJTExUfHx8VG6d++uPProo0plZaV1n45QDzWKYreknxBCCCGEEEI0gcyxEEIIIYQQQjSbBBZCCCGEEEKIZpPAQgghhBBCCNFsElgIIYQQQgghmk0CCyGEEEIIIUSzSWAhhBBCCCGEaDYJLIQQQgghhBDNJoGFEEIIIYQQotkksBBCCCGEEEI0mwQWQggh3DJz5kymTp3aZp8/ffp0li1b1mrvv2/fPuLj4yktLW21zxBCiM5IoyiK0taFEEII0T5oNJoGn1+yZAkLFy5EURRCQkLOTKHs7Ny5k/PPP58TJ04QEBDQap9z9dVXM3jwYB577LFW+wwhhOhsJLAQQghhlZWVZb3/z3/+k8WLF3Pw4EHrtoCAgFa9oG/MrbfeipeXFytWrGjVz/nmm2+47bbbOHnyJF5eXq36WUII0VnIUCghhBBWMTEx1ltwcDAajcZhW0BAQJ2hUOeeey7z589nwYIFhIaGEh0dzTvvvENpaSmzZs0iMDCQHj168O233zp81p49e7j44osJCAggOjqa6dOnk5eXV2/ZjEYjn3/+OVOmTHHYnpyczFNPPcWMGTMICAggKSmJr7/+mtzcXK644goCAgIYNGgQf/75p/U1J06cYMqUKYSGhuLv70///v353//+Z33+ggsuoKCggPXr1zfziAohRNchgYUQQohm+/DDD4mIiOCPP/5g/vz53HnnnVxzzTWMGTOGbdu2ceGFFzJ9+nTKysoAKCoq4vzzz2fo0KH8+eeffPfdd2RnZ3PttdfW+xm7du1Cr9czYsSIOs+9+OKLjB07lu3bt3PppZcyffp0ZsyYwU033cS2bdtITU1lxowZWDrp586dS2VlJb/88gu7d+/m2WefdeiJ0Wq1DBkyhF9//bWFj5QQQnReElgIIYRotsGDB/OXv/yFnj178vDDD+Pj40NERAS33XYbPXv2ZPHixeTn57Nr1y4AXnvtNYYOHcqyZcvo06cPQ4cO5f3332ft2rUcOnTI6WecOHECT09PoqKi6jx3ySWXcPvtt1s/y2AwcNZZZ3HNNdfQq1cvHnzwQfbv3092djYAJ0+eZOzYsQwcOJDu3btz2WWXMX78eIf3jI2N5cSJEy18pIQQovOSwEIIIUSzDRo0yHrf09OT8PBwBg4caN0WHR0NQE5ODqBOwl67dq11zkZAQAB9+vQB4OjRo04/o7y8HJ1O53SCuf3nWz6roc+/++67eeqppxg7dixLliyxBjz2fH19rT0sQgghGieBhRBCiGbz9vZ2eKzRaBy2WYIBk8kEQElJCVOmTGHHjh0Ot8OHD9fpObCIiIigrKyMqqqqBj/f8lkNff6tt97KsWPHmD59Ort372bEiBG8+uqrDu9ZUFBAZGSkawdACCGEBBZCCCHOvGHDhrF3716Sk5Pp0aOHw83f39/pa4YMGQKo60y0hISEBO644w6+/PJL7r33Xt555x2H5/fs2cPQoUNb5LOEEKIrkMBCCCHEGTd37lwKCgq44YYb2LJlC0ePHuX7779n1qxZGI1Gp6+JjIxk2LBhbNiwodmfv2DBAr7//nvS0tLYtm0ba9eupW/fvtbnjx8/Tnp6OpMmTWr2ZwkhRFchgYUQQogzLjY2lo0bN2I0GrnwwgsZOHAgCxYsICQkBA+P+r+abr31Vj755JNmf77RaGTu3Ln07duXiy66iF69evHGG29Yn//000+58MILSUpKavZnCSFEVyEL5AkhhOgwysvL6d27N//85z8ZPXp0q3xGVVUVPXv25B//+Adjx45tlc8QQojOSHoshBBCdBi+vr78/e9/b3AhveY6efIkjzzyiAQVQgjhJumxEEIIIYQQQjSb9FgIIYQQQgghmk0CCyGEEEIIIUSzSWAhhBBCCCGEaDYJLIQQQgghhBDNJoGFEEIIIYQQotkksBBCCCGEEEI0mwQWQgghhBBCiGaTwEIIIYQQQgjRbBJYCCGEEEIIIZrt/wGFhCZzf5wCTQAAAABJRU5ErkJggg==",
      "text/plain": [
       "<Figure size 800x300 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.figure(figsize=(8, 3))\n",
    "plt.plot(times / u.ms, u.math.squeeze(v_intact) / u.mV, label='intact')\n",
    "plt.plot(times / u.ms, u.math.squeeze(v_ablated) / u.mV, label='no $I_K$', linestyle='--')\n",
    "plt.xlabel('Time (ms)')\n",
    "plt.ylabel('Membrane potential (mV)')\n",
    "plt.title('Effect of ablating the potassium current')\n",
    "plt.legend()\n",
    "plt.tight_layout()\n",
    "plt.show()"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "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.13.11"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}
