{
 "cells": [
  {
   "metadata": {},
   "cell_type": "markdown",
   "source": [
    "# 低阈值与高阈值钙通道的门控特性对比\n",
    "\n",
    "本实例将展示如何使用 `braincell` 库分析低阈值 T 型钙通道和高阈值钙通道的门控变量稳态特性。通过对比它们的门控曲线，你将理解不同钙通道在神经元电活动中的功能差异。\n",
    "\n",
    "## 准备工作\n",
    "\n",
    "首先确保已安装必要的库（`braincell`、`brainunit`、`braintools`、`matplotlib`），并导入所需模块："
   ],
   "id": "61603c8025087c57"
  },
  {
   "metadata": {
    "ExecuteTime": {
     "end_time": "2025-10-13T09:23:58.484973Z",
     "start_time": "2025-10-13T09:23:58.481944Z"
    }
   },
   "cell_type": "code",
   "source": [
    "import braintools\n",
    "import brainunit as u\n",
    "import matplotlib.pyplot as plt\n",
    "import braincell"
   ],
   "id": "2761ba6d900de51a",
   "outputs": [],
   "execution_count": 9
  },
  {
   "metadata": {},
   "cell_type": "markdown",
   "source": [
    "## 代码详解\n",
    "\n",
    "### 创建钙通道模型\n",
    "\n",
    "`braincell` 内置了多种基于经典文献的离子通道模型，命名格式为 `通道类型_文献标识`，如 `ICaT_HP1992` 表示基于 Huguenard & Prince 1992 年研究的 T 型钙通道。\n",
    "如果你想了解更多的离子通道模型，可以查阅我们的[离子通道模型库](https://braincell.readthedocs.io/latest/apis/braincell.channel.html)。"
   ],
   "id": "da624d2453aa1324"
  },
  {
   "metadata": {
    "ExecuteTime": {
     "end_time": "2025-10-13T09:23:58.590117Z",
     "start_time": "2025-10-13T09:23:58.587337Z"
    }
   },
   "cell_type": "code",
   "source": [
    "# 创建低阈值 T 型钙通道（ICaT）\n",
    "cat = braincell.channel.ICaT_HP1992(1)\n",
    "\n",
    "# 创建高阈值钙通道（ICaHT）\n",
    "caht = braincell.channel.ICaHT_HM1992(1)"
   ],
   "id": "32d47cc46172bcb8",
   "outputs": [],
   "execution_count": 10
  },
  {
   "metadata": {},
   "cell_type": "markdown",
   "source": [
    "\n",
    "### 生成电压序列\n",
    "\n",
    "为了分析通道门控特性随膜电位的变化，我们生成一个从 -100 mV 到 0 mV 的连续电压序列："
   ],
   "id": "d649aeccd48d6e05"
  },
  {
   "metadata": {
    "ExecuteTime": {
     "end_time": "2025-10-13T09:23:58.700512Z",
     "start_time": "2025-10-13T09:23:58.697465Z"
    }
   },
   "cell_type": "code",
   "source": [
    "# 生成电压序列\n",
    "vs = u.math.arange(-100 * u.mV, 0 * u.mV, 0.1 * u.mV)"
   ],
   "id": "43f7945b1fe9233f",
   "outputs": [],
   "execution_count": 11
  },
  {
   "metadata": {},
   "cell_type": "markdown",
   "source": [
    "### 计算门控变量稳态值并绘图\n",
    "\n",
    "离子通道的门控变量的稳态值描述了在某一膜电位下，通道门控处于开放状态的概率。我们通过通道的 `f_q_inf` 和 `f_p_inf` 方法计算这些值，并可视化对比："
   ],
   "id": "88ec4ea16b8d4d4a"
  },
  {
   "metadata": {
    "ExecuteTime": {
     "end_time": "2025-10-13T09:23:59.281336Z",
     "start_time": "2025-10-13T09:23:58.808713Z"
    }
   },
   "cell_type": "code",
   "source": [
    "# 创建图形和子图布局\n",
    "fig, gs = braintools.visualize.get_figure(1, 2, 3., 4.5)\n",
    "\n",
    "# 计算低阈值通道的激活门稳态值和失活门稳态值\n",
    "q_inf = cat.f_q_inf(vs)\n",
    "p_inf = cat.f_p_inf(vs)\n",
    "\n",
    "# 添加子图 1\n",
    "fig.add_subplot(gs[0, 0])\n",
    "plt.plot(vs / u.mV, q_inf, label='q_inf（激活门）')  # 横轴转换为 mV 便于阅读\n",
    "plt.plot(vs / u.mV, p_inf, label='p_inf（失活门）')\n",
    "plt.legend()  # 显示图例\n",
    "plt.fill_between([-80, -60], 1., alpha=0.2)  # 标记低阈值通道的典型激活范围\n",
    "plt.title('低阈值钙通道（ICaT）')\n",
    "plt.xlabel('膜电位（mV）')\n",
    "\n",
    "# 计算高阈值通道的 q_inf 和 p_inf\n",
    "q_inf = caht.f_q_inf(vs)\n",
    "p_inf = caht.f_p_inf(vs)\n",
    "\n",
    "# 添加子图 2\n",
    "fig.add_subplot(gs[0, 1])\n",
    "plt.plot(vs / u.mV, q_inf, label='q_inf（激活门）')\n",
    "plt.plot(vs / u.mV, p_inf, label='p_inf（失活门）')\n",
    "plt.fill_between([-60, -40], 1., alpha=0.2)  # 标记高阈值通道的典型激活范围\n",
    "plt.legend()\n",
    "plt.xlabel('膜电位（mV）')\n",
    "plt.title('高阈值钙通道（ICaHT）')\n",
    "\n",
    "# 显示图形\n",
    "plt.show()"
   ],
   "id": "c80b60877a874795",
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<Figure size 900x300 with 2 Axes>"
      ],
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA48AAAE3CAYAAAAKb3Q+AAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjcsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvTLEjVAAAAAlwSFlzAAAPYQAAD2EBqD+naQAAcEdJREFUeJzt3Qd8FNXaBvBnd7PpvRMSCL3X0BFEQRCVq9iwgQj2frkWsFD87O3aewG7V7GCNBEQ6b13kpDee8/ufr9zJhsSkpDdtNny/O8dd3Z2NvtmgH3nzJzzHo3JZDKBiIiIiIiI6Dy053uRiIiIiIiISGDjkYiIiIiIiBrFxiMRERERERE1io1HIiIiIiIiahQbj0RERERERNQoNh6JiIiIiIioUWw8EhERERERUaPYeCQiIiIiIqJGsfFIREREREREjWLjkYiIiIiIiBrFxiMRSYsXL4ZGo8HOnTvrvLZ3717ccsstiIqKgpubGwIDAzFhwgR8/vnnMBgMFn/GwoUL5Wc0towbN66FfzsiIqK2y3dm0dHRuOKKK+p9bf369TKOH3/8UT63JD+KRbyPSC0uqn0ykUoOHTqEQYMGwdXVtd7Xy8vLsWfPnkb3OXLkCEpLSx1ivy5duqAhn3zyCe6++26EhYVh+vTp6NatGwoKCrB27VrMnj0bKSkpeOKJJ2CJq6++Gl27dq1+XlhYiHvuuQdTp06Vr5mJz3rvvfcwZ84cuLjU/zUVHByMuLg4iz6XiMgZMd+pl++a4ssvv6z1/IsvvsCaNWvqbO/VqxfCw8NlDq1PZWWlzKGzZs1qtVjJebHxSE7HZDJh2LBh+Oeff+p9fcSIERbv4yj7NWTr1q0ykY4cORJ//PEHfHx8ql97+OGH5VXbgwcPwlL9+/eXi1lmZqZsPIpt4kpvTceOHcMjjzyCZ599ts7PEScJNRuhRERUF/Nd3f3aKt81xbl5UMQkGo/nbjc3EHNzc+u9wDp37lwYjcZWjZWcF7utElGDFi1aJLvIfP3117USqdmQIUMwc+bM6uevvvoqRo0ahaCgIHh4eCAmJqa6Ow4REZGtYr4jsgzvPBJRvYqLi2VXnbFjx6JDhw4WvefNN9/Ev/71L9x8882yi9B3332H6667DsuWLcPll1/e6jETERHZWr6rqKiQPW3OlZeX12K/A1FbYeORiOp18uRJmfD69etn8XuOHz8ur8Ca3X///Rg8eDBef/11Nh6JiMgp893q1asREhLSojETqYWNRyKqV35+vnysr/tOQ2om0pycHFmZbsyYMfj2229bJUYiIiJbz3fDhw+vd/z+vn375Nh+InvCxiMR1cvX11c+ikpzlhLddUSCFKXOy8rKqreLcSRERETOmO9EdXAx3ce5GqomTmTLWDCHiOolqpmKxHbgwAGL9t+4caMc/+Hu7i5LhItqdaJK3E033XTeCndERERqYr4jshwveRBRvTw9PXHxxRfjr7/+QkJCgpww+XyWLl0qE+mqVavkxMpmYmJlIiIiW8V8R2Q53nkkogYtWLBAXkUVkyXXNxnxrl27sGTJErmu0+lkdx0x7sMsLi4Ov/zyS5vGTEREZC3mOyLLsPFIRA0Sc1i9++67sotOz549MW/ePHz22WeyRPnUqVPlhMzJyclyX1FdTpQ7v/TSS/HBBx/gmWeekUUCRHcgIiIiW8Z8R2QZdlslovO66667MHToULz22mv44osvkJGRAW9vb1mSXHTRueWWW+R+osvPp59+ihdffBEPP/wwOnXqhJdeeklejd2/f7/avwYREdF5Md8RNY6NRyKSZs6cKZf6iMT59ddfN/ozZs2aJZdzLVy4sN79RQU6FhcgIiJHzXeiQdmQcePGnTcHvvPOO3IhsiXstkpERERERESN4p1Hckpbt26Fv79/va+ZB8pbso8j7WeLXn311QavuoquREREdH7Md3X3cxSi9059SktLeceSWo3GxD5jRERERERE1Ah2WyUiIiIiIqJGsfFIREREREREjWLjkYiIiIiIiByjYI7RaJQTs/r4+ECj0agdDhERORAx9L+goAARERHQam3zmirzIBER2UIutIvGo0iYUVFRaodBREQOLCEhAZGRkbBFzINERGQLudAuGo/iSqv5l/H19VU7HCIiciD5+fmyYWbONbaIeZCIiGwhF9pF49HcRUckTCZNIiJqDbbcHZR5kIiIbCEX2ubgDiIiIiIiIrIpbDwSERERERFRo9h4JCIiIiIiopZvPP7999+YMmWKLOMq+sT+8ssvjb5n/fr1GDx4MNzc3NC1a1csXrzY2o8lIiKyCcyDRETkrKxuPBYVFWHAgAF49913Ldo/NjYWl19+OS666CLs3bsXDz/8MG6//XasWrWqKfESERGpinmQiIicldXVVidPniwXS33wwQfo1KkTXnvtNfm8V69e+Oeff/Df//4XkyZNsvbjiYiIVGXPebCk3ICC0grliUb8XwNzYT1NjSp7yrp5u9zxvPvJfWq8B/Vsb/Dn23CVWyIiauOpOrZs2YIJEybU2iaSpbjy2pCysjK51Jx3hIiIGmEyAYZyoKIYqCgFKkuUx4qSs+uVpco+xkrlUS5V68YKwFBRta1q3VjPc6MBMBmrHg1nn8t1Yz3bztm/1nur9pH7isWkPKLq0SR/sQZeMwEdRwK3/g5bZkt5cPmBFDzywz7YqvM2MqsasVoNoNdqodNp4KLVQKcVj9qqRw1cdGKbtsZrVY9yfy3cXLRw1+vgodfBXa+smxePqucerjp4u7nAz0MPXw+9fBSLeI2ISBUiX1YUAeUixxcpuV2uFwMdRwO6tpmBsdU/JTU1FWFhYbW2ieciEZaUlMDDw6POe1544QUsWrSotUOjRhxIzFM7BLvQL9JP7RDI3omGVGkeUJylLCW5QFkBUJZf9VhwzvMa20XykI3Dqkai0tpyHqLha+NsKQ+KdphofIm/JaLtbWvk9QHzSt1Xq9dKIS4itD1XFy183fUI8NQjzNe9anGrfmzn54HoIC/4eerhbHjOYBmeMzg5owEoyVFyfVHm2bxfnFmV+/OB0oZyf0HDOf6xWMAzsE1+hbZpolpp3rx5mDNnTvVzkWCjoqJUjYmIyOpGTVE6kJ8CFNRYCtOA4uyqZFH1KBJJSzf6NFpA7wm4uAN6j9qPOlflCqV8dAW05nW9smj1De8jHrU65efLR93ZR7lNW8+2+vY371f1KJ+bF3GHSQzJ19R4Lu9B1X3dxQ2OqLXy4DUxkXI5l8lkqm6vKQ1L5Ym5kWmq+vtZs01n3n6+95mfWLKf/Axz47GezzVvNxpNqDSaYDAa5WOlQaybt4lHo/JoqL2f2FZhMKGs0oDSCiNKKwzVS0lF7W3ieUFpJfJKKpBfUiEfjSagvNKIzMIyuZxIL2zwOIvGZXSwFzoFeaFrmDf6Rvihb3s/BHq5Nu0PjohsX0UpkBMH5CcC+clK/s9PqlpPVs4BWiTfi+4XnoCrp5LX9V5Kb6I20uqNx/DwcKSlpdXaJp77+vrWe7VVENXoxEJEZLNEF868RCVR5MRWPYolHihIVRqOsoulFdz8lCuHHv6Amy/g5lPjseZi3uatJJDqxqFYdwdcPJRGIMeS2QR7yIOim2jtvy78u1OTaOwWlimNSbFkF5UjPb8MaQWlymN+qVwSc0qQXlCGnOIK5JzJxZ4zubV+Tnt/D/Rt74thnYIwsnMQeob7QCtuBROR/XQdzY0D0g4DWSeB7NNVS6zSULS0YejuD3gGAV7ByqNYPAIAd98G8n+NdZHzVczvrd54HDlyJP74449a29asWSO3ExHZPDGeIPM4kHEUSD+iPIolN0EZq3c+4o6aTzjg00559I0AvEMBT3OyCKydNESDjxwO86D9E41rH3e9XCIDzr9vUVkl4rKKEJdZjNjMQhxNLcCh5HzEZhYhKbdELqsOpVXfoRzVNRgTe4fh4p6h8ucTkY2oLANS9ilL2kEg7ZDSaBTjDRsiGnl+UUq+l0t7wLedsi7OBbxCAI/ANhuf2BqsjrywsBAnT56sVYJclB4PDAxEhw4dZFebpKQkfPHFF/L1u+++G++88w4ee+wxzJo1C3/99Rf+97//Yfny5S37mxARNZcYd5i8F0jeDSTvAVL2K3cTG7qSqHMDAqLPLoGdAP8OVUkiQrmiKLpokkNhHqTz8XJzQZ8IP7nUJKrcikbk3oRcbDmVhR1x2fIO5fL9KXJx1WkxumsQpg6OxKQ+YXBz4XcHUZsqzADiNwEJ24HE7UqjURSMqy/3h/RQlsAuQGDns4u4KOzgvX6sbjzu3LlTzlVlZh6Tceutt8pJj1NSUnDmzJnq10V5cpEg//3vf+PNN99EZGQkPvnkE1Wm6VhzOA3HUvPleA8xqJ2InJgYRCW6msT9A8RvBpJ2AVkn6t9XXCUM7QWE9Kx67AEEdQW8w5Uxe+RU7DkPknrEXcURnYPkcveFXVBhMGJfQi7WHUvHioOpOJ1RhHXHMuQi7khOHRSJmaOi0SHIU+3QiRyTKDIn8v/pdcCp9UDagbr7iJ5C7WOA8H5AWB8grK/SSLTjO4fNpTGZR63bMFEowM/PD3l5eXKMSFNN+3ALtsVmywsCF3QNxnVDomRXEZberh8rp1mGldPsSO4Z4OTaqgbjJmXw+rnEncOIQUDEYCBiIBDaW+lm4uBXEp1ZS+UYZ4+Rmk6cip1ML8Rv+5Lxw85EpOaXyu1iOOSUARG4Z1wX9Ay37T93njNYhucMKhMVTY+vAo78ppwPiGmsahKNw6jhQNQwZQno5DT5P9/CPOPiTF/MorEoWsrbY7Ox8USmXMTVvekjOmL6yGiE+LBID5HDVTxN2AacWAUcXw1kHKn9uqge2n4IED1aSRai0Si6mhIRtfGYym5hPvjPxB54eEJ3/H08A59vjpOPv+5NlotoRD5+aQ9EBvBOJJHVdxiPLAP2fwecXl+7MqkYk9j5IqDLRUCnCwHvEDUjtQtOdefRLD6rCD/uSsTSXYlIziutnrvp6kHtce+4ruwiUoVXES3Dq4g22GCM3QAc/Ak4ugworVHtUEzvEDkM6DwOiL4AiByiVC0jp2YPd/XsIUZqeQeT8vD++lP442CK7GkvzlVuv6ATHhzfzeZ6TfGcwTI8Z2gj4h9M4k5g71fK+YCYK9EspBfQa4qyiO6oTnJnsaXyjFM2Hs3EnE+rDqXio79PywHsgl6nwQ1DO+CBi7si1NcdzoyJwDJMBDZAfI2JO4z7vlO6ooi5E2uOV+x2CdBtItDl4jabRJfshz00zOwhRmo9h5Lz8OyyI9hyWvlu6xTshZev7Y+h0bbzfcZzBsvwnKGVVZYDh34Gtr4HpOw9u92vAzDwJqDftUBwNzUjtFlsPFpBHIKd8Tl4a+0J2ZVVcNdrcf9FXXHH2M5OW/GMicAyTAQqKsoC9n0L7P4CyDx2druY+qL3VUDfq4EOI1nxlOy+YWYPMVLrEucqqw+nYf6vB5GWXyZvltwxpjMendQDep36hbt4zmAZnjO0ktJ8YMfHwLaPgMJUZZuY/1icCwy6Geh4AQvcNYKNxyYS5bNfXnW0emLfzsFeWHRlH4zp5nx9oJkILMNEoALRFUVcVTzy+9ky2npPoM9U5api9FinroRGjtcws4cYqW3klVTgueWH8b+difL50OgAvHPTYISp3FuK5wyW4TlDCysrBLZ/BGx+CyjJUbaJSujD7gBibgO8gtSO0G6w8dgM4pCIimfPLj+CjIIyue36IZGYP6UPvN2c54SUicAyTARtxGgEjq8ANr8NnNlydnu7gUDMrUDfawF3nlSTYzbM7CFGalsrD6bi0R/2oaCsEsHebvhoRgwGdwhQLR6eM1iG5wwt2D1V3Gnc+NrZoSpB3YCxjwB9rgZcXNWO0O6w2mozq55dObA9LuoZitdXH8eSLXHyCt/W09n477QBiOloO2MMiBye0QAc+AH4+xUgq2pidq0e6HcdMPwuZToNIiInc2nfcPQI98E9X+3C0dQC3PTxVrx942Bc0jtM7dCIWteJNcDKeWfnZhbzLl44V+l5xGEqrY6df8/D112Phf/qg2/vGIH2/h44k12M6z7YgvfWn5R3J4mole80HlwKvDcC+PkupeHo7gdc8G/g4QPA1PfZcCQipyYK5yy9ZxQu6hGC0goj7vpyJ/63I0HtsIhab67mr68Hvr5WaTiKOZinvAnctwMYMI0NxzbCxqMFRnQOwoqHx8ipPIwm4OWVx3D/N3tQVFZjnhgiajli4t4PxwA/zgIyjwPu/sD4BcC/DwETFgK+7dSOkIjIJni5ueDjGUMwbUiUPEd5/Kf9+N9ONiDJwS4mb/8YeG+kMm+z6H006gHggV1AzEzWOGhjPNpW3IV87foBGNwxAIt+P4TlB1JwKqMQn9w6hBP2ErWUrFPAqieVsY2Cmy8w8j5gxD3KXUciIqrDRafFi9f0k5Xil2yJx+NL90On0eCamEi1QyNqnuxY4Nf7gPhNynNRQf1fb3O6DRWx8WjlWMhbRnREz3Af3P3VbjnG4Jr3N2PJrGHoGc4CBkRNVl4EbHgJ2PIeYKwAtC7AsDuBsY9yXkYiIgvPUcRQG3H38cut8Xhs6X4EebtiXI9QtUMjahoxdOX3h4GyfEDvpfQ8Gno7p9xQGY9+EwyJDsTvD4xG9zBvOdeSGAe5PTZb7bCI7NPp9UpXlE1vKg3HLuOBezYDl77AhiMRkZUNyGeu7COH2RiMJtz39W4cTGIVVLIz5cXAbw8qQ1dEwzFqBHDvFmD4nWw42gD+CTRROz8P/HDXKAzpGICC0krM+GwbNp/KVDssIvsh5mP65T7giyuB3HjANxK48TvglqVASA+1oyMistsG5IvX9MforkEoKjdg1uIdSM8vVTssIsuL4nx6CbB7ifjbrPRAmrkcCOiodmRUhY3HZvDz1OOr24djXFWVs9mLd2Lb6aq5ZoioYac3KHcb936lJAfRRfW+rUCPyeLMR+3oiIjsmquLFu/fEiN7SKUXlMkifxUGo9phEZ1f/Bbgo4uAtIOAVygw4xfg4qdYEMfGsPHYTO56HT64JQZju4egpMKA2xbvwK54dmElqpehAvhzoXK3sSAFCOoK3LYCuOwVwM1H7eiIiByq0J84P/Fxc8H2uGy8vPKo2iERNWzPV8CSKUBxJhDeH7hzHdB5nNpRUT3YeGyhBuRH02MwplswissNmL1kJ06mF6odFpFtyT4NfDoR+Oe/AEzA4FuBu/4GOo5UOzIiIofUOcQbr1w3QK5/vDEWaw6nqR0SUW1i3vT1LykVVUXdg95XArNWAn6sFGyr2Hhs0QbkEAyM8kducQVmfr4d6QUcY0AknfgT+GgckLxbmbPx+i+Af70FuHqpHRkRkUO7tG84br+gk1yf99N+ZBWWqR0S0dn5G1fOA9Y/rzwf8whw7WKeG9g4Nh5bkIerDp/eOgTRQZ5IzCmRg9SLyyvVDotI3SuKG18Dvr4WKM0DIocC92xSriwSEVGbeGRSD/QI80FmYTnm/XQAJvHdTKQmQyXw673AtveV55e+BIx/mtVU7QD/hFpYkLebnPcxyMsVB5Py8diP+/klTc47d+MPtwJrnznbTVVUTGNXFCKiNu8d9d9pA6HXabD6cBqW7k5SOyRy9obj0tnAvm8BjQ6Y+iEw4m61oyILsfHYCjoGeeGD6TFw0WqwbH8KPtkYq3ZIRG2rMB1YfDlw+FdAqweueEPppuripnZkREROqXeELx6e0F2uP7f8MLKLytUOiZyR0QD8cjdw+BdA5wpM+woYcIPaUZEV2HhsJUOjAzF/Sm+5/sKKI9h0knNAkpPIPKnM0ZS8B/AIBGYuA4bcpnZURERO786xndEz3Ac5xRV4ccURtcMhZxzj+NsDwIEfAK0LcN0SoOdlakdFVmLjsRVNH9ER1wyOhNEEPPTdHmQUcJA6Obgz25SGY04cEBANzF4DdBihdlRERARAr9Piual95fr/diZieyynFqM2IoZwrZwL7P1a6ap6zadsONopNh5bkUajkV/S4iqfGKT+yA/7YBQtSSJHdHKtMn9jSTYQMRiY/ScQ3FXtqIiIqIaYjoG4cViUXJ//60EYeF5CbWHTm8D2D5X1qR8Afa5SOyJqIjYe22CQ+ls3DoKbixYbjmfg881xaodE1PKOrQS+vQGoLAG6TVS6qnqHqB0VERHV47FJPeHr7oKjqQVYuitR7XDI0e3/AfhzgbI+6Xmg//VqR0TNwMZjG+ge5oOnrlDGP7604igOJ+erHRJRyznyO/D9LYChHOh5BTDta87RRERkwwK8XPHg+G5y/dXVxzitGLWe2L+BX+5R1kfcB4y8T+2IqJnYeGwjtwzvgAm9wlBuMOLRH/ehwmBUOySi5jv4E/C/WwFjBdD3GuC6xYCLq9pRERFRI6aP7IgOgZ5ILyjDx3+zKjy1AlH/4H8zlHOEPlOBic+qHRG1ADYe23D84/NX94Wfhx6HkvPx0d+n1Q6JqHmOrQCW3g6YDMCAG4GrPwZ0erWjIiIiC7i56PD4pT3l+od/n0JmIYv6UQvP9fzdzUBJDhAxCLjqfUDLZocj4J9iGwr1ccf8qu6rb/55AifTC9QOiahpTm9Q7jiKhmP/G4Ar3wO0OrWjIiIiK1zWLxwDIv1QXG7Ax7yoTS1ZWfXX+4C0g4BXqDKcRe+hdlTUQth4bGNXD26PC7uHyO6rj/24n9VXyf4k7gS+vREwlCljHK98l1cTiYjstFfUwxO6y/UvtsTz7iO1jE1vAId+BrR6YNqXgF97tSOiFsQzPlW6r/aDl6sOu8/k4oddCWqHRGS5tMPAV9cAFUVA53HAtZ8BOhe1oyIioiYa1yMEA6L8UVJhwIcbTqkdDtm7uE3A2meU9cte5lzPDqhJjcd3330X0dHRcHd3x/Dhw7F9+/bz7v/GG2+gR48e8PDwQFRUFP7973+jtLQUzqq9v0f1lb6XVh5DXnGF2iERNS4/Bfj6WqA0F4gcpnRDcXFTOyoiVTAPkmPdfVQqr365NR4ZBbz7SE1UnF1VC8EIDLgJGDJL7YjIFhqP33//PebMmYMFCxZg9+7dGDBgACZNmoT09PR69//mm28wd+5cuf+RI0fw6aefyp/xxBNPwJnNHB2NbqHeyC4qx2trjqkdDtH5lRUC31wP5CcBQd2Am74H3LzVjopIFcyD5GjGdQ/BwCh/lFYYsYTzUVNTxzmKKTkKkpXzhMteUTsispXG4+uvv4477rgDt912G3r37o0PPvgAnp6e+Oyzz+rdf/PmzRg9ejRuuukmeZV24sSJuPHGGxu9Suvo9DotFl3ZR65/tTUeh5Lz1A6JqH5GA7B0NpC6H/AMBm7+AfAMVDsqItUwD5Ij3n28+8LO1XcfOe8jWW3bB8DxlYDOTRnSwgvMDsuqxmN5eTl27dqFCRMmnP0BWq18vmXLlnrfM2rUKPkec5I8ffo0/vjjD1x22WUNfk5ZWRny8/NrLY5oVJdgXNG/HUTNnGd+PwyTuGpDZEvE38mVc5WE4OIO3PgdENhJ7aiIVMM8SI7qkt7h6BjkibySCvywM1HtcMiepOwHVj+trE96DmjXX+2IyFYaj5mZmTAYDAgLC6u1XTxPTU2t9z3iSuszzzyDCy64AHq9Hl26dMG4cePO213nhRdegJ+fX/Uixoc4qicu6wVXFy22xWZj3bH6uzwRqWbnZ8D2j5T1qR8CUUPVjohIVcyD5Kh0Wg1uv0C5OPjJP6dhYDV4skRludJd1VgB9LgcGHq72hGRvVdbXb9+PZ5//nm89957cmzITz/9hOXLl+P//u//GnzPvHnzkJeXV70kJDhuRdIIfw/cNjparr+44ii/rMl2nNkGrHhcWR8/H+hzldoREdkl5kGyF9fGRCHAU4+E7BKsOlT/xRCiWv5+WZnP0TMImPKm6AOtdkTUyqyqsR8cHAydToe0tLRa28Xz8PDwet/z9NNPY/r06bj9duVKRL9+/VBUVIQ777wTTz75pOzucy43Nze5OIt7x3XF9zsScDytEEt3JeL6obzCTDZQWfV/05Urib2vBC6Yo3ZERDaBeZAcmYerDtNHdMRbf53ER3+fxmX92qkdEtmypF3AxteV9ctfB7xD1I6IbO3Oo6urK2JiYrB27drqbUajUT4fOXJkve8pLi6ukxhF4hU4xk/h56HH/Rd1leui8mpJuUHtkMjZu6D8bwZQmAaE9AKufI9XEomqMA+So5s+Mhp6nQZ7E3JxMInF/KgBFaXAL/cCJgPQ9xr2TnIiVndbFeXJP/74YyxZskSWHL/nnnvkFVRRdU6YMWOG7G5jNmXKFLz//vv47rvvEBsbizVr1sirsGK7OXmS+LLuiMgAD6Tll+GzTbFqh0PObOXjQOJ2wN0PuOFrVkwjOgfzIDmyEB83TO7brroaPFG9NrwEZBwFvEKBy15VOxqy1W6rwrRp05CRkYH58+fL4gADBw7EypUrq4sHnDlzptYV1qeeekqWgBaPSUlJCAkJkQnzueeea9nfxM65uejwn4nd8e/v9+HjjacxY2RH+Ljr1Q6LnM2BH5UiOdAAV38CBHVROyIim8M8SI7ulhEd8du+ZPy6NxnzLusle0gRVUs7DGx+S1m/4nVO3+VkNCY76DMjSpSLanOiaICvry8clSiWc8l/N+B0RhEendQD91V1ZVXLgUR2V7FEv0g/OISsU8CHY4HyQmDsY8DFT6odEVGbsIccYw8xkuMQp4aT3vhb1mJYMKU3bhvd+BRNPGdwknMGoxH4/FIgYRvQ8wqlhxI5BEvzTKtXWyXrymQ/NL6bXBd3HwtKK9QOiZxFZRnw421Kw7HjaODCqiqrRETkdMSdcnH3Ufh62xmOzaWz9nyhNBxdvYHJL6kdDamAjUcbc0X/CHQO8UJucQW+2MKxBtRG1swHUvYBHoHANZ8AOqt7tBMRkQOZOqg9PF11OJleiK2ns9UOh2xBYbpyviBc9CTgF6l2RKQCNh5t8O7jgxfz7iO1oaPLgW0fKOtTPwB8I9SOiIiIVCbqLlw1qL1c/99OzjNKAFY9CZTmAeH9gWF3qh0NqYSNRxs0ZUAEOgfz7iO1gYJU4Nf7lfWR9wPdJ6kdERER2YjrhyjzTq84mMKL2c4ufjNw4H9KQb0pb7CHkhNj49FG7z7ef7FSLOfzTbEoreC8j9QKxBiW3x4ESrKB8H7A+AVqR0RERDZkQKQfuoZ6o7TCiOX7U9QOh9RiNAArqmohxNwKtI9ROyJSERuPNnz3sb2/BzILy/HT7iS1wyFHtHsJcGIVoHMFpn4EuLiqHREREdlY4ZxrY5RxbT/uSlQ7HFLL3q+B1P2Amy9w0VNqR0MqY+PRRul1Wsy+oFP12EcxjQdRi8mOBVY+oayPnw+E9VY7IiIistHCOVoNsDM+B7GZRWqHQ22tNB9Y+4yyLiqxe4eoHRGpjI1HGzZtaJScmFd8Wa85nKp2OORI3U9+vhuoKFKm5Rhxr9oRERGRjQrzdcfY7kqDYSnvPjqfv18BijKAoK4skkMSG482zMvNBTNGKvMsvb/hNOdZopax+W0gYasyR9NV7wNandoRERGRDbsuRimcs3R3IntCOZOsU8DW95X1Sc9zeAtJbDzauFtHRcPVRYt9CbnYHst5lqiZMo4D655T1i99EQhQLk4QERE1ZHyvUPi6uyAlr5TnIs5EzOlorAC6TgC6TVQ7GrIRbDzauGBvN1xXNVhdjH0kajKjEfj9QcBQriSCQbeoHREREdkBd70Ok/u2k+u/7UtWOxxqC2e2AUeXARotMPE5UT1J7YjIRrDxaAdmVRXOWXs0HWeyitUOh+zVzk+BM1sAvRdwxX+ZCIiIyKoq8OY5HysMRrXDodYkhkn9uVBZH3gzENpT7YjIhrDxaAe6hHjLweri3/IXW+LUDofsUV7i2UQwYQHg30HtiIiIyI6M7BIke0PlFlfgnxOZaodDren4KuDMZsDFHRg3T+1oyMaw8WgnbhsVLR+/35mAorJKtcMheyKuOiybA5QXApHDgKG3qx0RERHZGZ1Wgyv6s+uqU1RkX7tIWR9+F+DXXu2IyMaw8WgnLuweguggTxSUVuKnPUlqh0P25MCPwIlVgM4V+NfbrK5KRETN6rq6+lAqSisMaodDrWH/90D6YcDdD7jg32pHQzaIjUc7odVqZOVVYcnmOE7bQZYpzgZWzlXWxz7KcQtERNRkgzv4o72/B4rKDfjraLra4VBLqygF1j2vrF8wB/AIUDsiskFsPNqRa2Mi4eWqw8n0Qmw6maV2OGQP1j4DFGcCIT2B0Q+rHQ0REdkxjUZTfffxd3ZddTy7PgfyEgCfCKXLKlE92Hi0Iz7uelw3RJmo9/NNsWqHQ7YucRewa7GyfvnrnNyXiIia7V9VjUdRAb6QNRgcR3kxsPF1Zf3CxwC9h9oRkY1i49HOzBipTOr+17F0JOZw2g46z4D35WKsggnofwMQPVrtiIiIyAH0aueDzsFeKK80Yh27rjrWdF5F6Uo1ds4DTefBxqOd6RzijVFdgmQBze93JKgdDtmqnZ8BKfsANz9g4v+pHQ0RETlQ19VL+4bL9ZWHUtUOh1pCeRHwzxvK+tjHAJ1e7YjIhrHxaIduGq7M0Scaj5yol+ooSAPWVjUYxz8NeIeqHRERETkQc+NR3Hlk1VUHsOMTpT5CQDQw4Aa1oyEbx8ajHZrYOxzB3q5ILyjD2iPsMkLnWPM0UJYHtBsIDJmldjRERORg+rX3Q4SfO4rLDdh4IlPtcKg5ygqBTW8q67zrSBZg49EOubpocW2MUjjn2+1n1A6HbEncP8ocTdAAV7zOOR2JiKhVuq5OMnddPciuq3Zt+0dAcRYQ2BnoP03taMgOsPFop24YqjQe/z6RgYRsFs6hqiI5Kx5X1ofcBrSPUTsiIiJyUJf2URqPfx5J4xAae1WaD2x+S1m/8HFA56J2RGQH2Hi0U9HBXriga7AsnPPdDt59JAC7lwBpBwF3f+Dip9WOhoiIHNiQ6EAEebkir6QC205nqx0ONfWuY0kOENQV6Hut2tGQnWDj0QEK5/xvZyKv+jk78eVvLpJz0ROAZ6DaERERkQPTaTWY2CdMrq88lKJ2ONSUCqtb36sx1pF3HckybDzasUt6hyHY2w0ZsnBOmtrhkJo2vAyUZAMhPVkkh4iI2sSlfdvJx1WH0mAUXaHIfuz+Qhnr6N8R6HuN2tGQHWHj0Y7pdVpcNyRSrn+znXM+Oq2MY0rXE+HSF1gpjYiI2sTIzkHwcXeRF7GPphaoHQ5ZqrIc2FQ11vGCh3nXkazCxqOdmzZEKZyz8UQGUvJK1A6H1LDqCcBYCfS4DOhysdrREBGRE1V/H99TmUt46+kstcMhS+37FihIBrzDgYE3qx0N2Rk2Hh2gcM6w6EBZOOen3Ulqh0Nt7fhq4OSfgFYPTHxW7WiIiMjJTOitjHvcHsuiOXbBUAn8819lfdQDgIub2hGRMzQe3333XURHR8Pd3R3Dhw/H9u3bz7t/bm4u7rvvPrRr1w5ubm7o3r07/vjjj6bGTOcwd139YWcCTBxz4FzdTlbNU9ZH3AMEdVE7IiKnwTxIpBjbPQR6nQZJuSVIzmUPKJt3+BcgJxbwCABiZqodDdkhqzs5f//995gzZw4++OADmTDfeOMNTJo0CceOHUNoqNJ1oaby8nJccskl8rUff/wR7du3R3x8PPz9/Vvqd3B6l/VrhwW/HUJcVjF2xudgaDQrbToFMc4x6yTgFQKMfVTtaOySuNhSWVkJg8GgdijUyvR6PXQ6XYv8LOZBorN83fUY3ikI/5zMlHcfrxrUXu2QqCFGI7DxdWV9xL2AmzfzoBPR6XRwcXGBRqNp28bj66+/jjvuuAO33XabfC6S5/Lly/HZZ59h7ty5dfYX27Ozs7F582aZvAVxtZZajpebCy7v1w4/7EqUdx/ZeHSSqTn+fkVZHz8fcPdVOyK7I07oU1JSUFxcrHYo1AZEsoyMjIS3t3ezfxbzIFFt43uFysbjttgsNh5t2YlVQPohwNUHGHYH86AT8vT0lD1gXF1d26bxKP6S7dq1C/PmzTvb71WrxYQJE7Bly5Z63/Pbb79h5MiRsrvOr7/+ipCQENx00014/PHHG7wKXFZWJhez/Px8a8J0StcNiZKNx+X7U7BgSh/ZoCQHJq4cluYCob052L0JjEYjYmNj5XdQRESE/BJt7pU4sl3iynpGRgYSExPRrVu3Zt2BZB4kqmtCrzAs+v0wDqfko7C0Et7uPAexOWJY09+vKutDZ8Po5ofYEyeYB50oD5aXl8tcKM5/RC4UuasprPrXnZmZKW9rh4Upg6PNxPOjR4/W+57Tp0/jr7/+ws033yzHd5w8eRL33nsvKioqsGDBgnrf88ILL2DRokXWhOb0hkYHIDrIU3Zd/eNAimxMkoPKTQC2faisT1gEaFumK54zEV+gogEZFRUlr8KR4xMNtri4OJl7mtN4ZB4kqisq0BMdAj1xJlsMn8nGuB51u2+TymL/BpJ2Ai7uwMj7mAedkIeHh+z9IoZNiD9/MWbfJqutir+YYpzHRx99hJiYGEybNg1PPvmk7ObTEHFFNy8vr3pJSOAcho0RV4uujakqnLMrUe1wqDWtew4wlAHRY4Bul6gdjV1r6lU3sj9qXlFnHiRnICq/CzviWHXVJm2umtdx0C2A99nGPfOgc9G2wJ+3VXceg4OD5RXbtLS0WtvF8/Dw8HrfI/rVnluooFevXkhNTZWt3vr63IpKdGIh61wTE4nX1hyXA9bjs4rQMchL7ZCopaUeAPZ9p6xfskicEasdEZFTYR4kqt/wToH4cXcidsXnoNJghIuOjRKbkXZImdZLo5V3HYmaw6p/2SLBiauma9eurXVFVTwX4znqM3r0aNlFR+xndvz48WYP1qS62vl5YEy3ELn+I+8+OqY/F4qe60CfqUD7GLWjIXI6zINE9esW5gM/Dz2Kyg04lMIxujZl8zvKY68pQGBntaMhO2f1ZSFRnvzjjz/GkiVLcOTIEdxzzz0oKiqqrjo3Y8aMWoUExOuiytxDDz0kk6WoSPf888/LwgHU8q6r6rq6dFciDEbO+ehQTm9QrhxqXYCLn1Y7GrJDosKnmFbCGps2bUK/fv3knbOrrrqq1WKzJ8yDRHXptBoM6Rgg10UPKLIR+cnAgR+U9VEPqR2NTWAubB6ry2GJsRqiUs/8+fNll5uBAwdi5cqV1cUDzpw5U6s/rRiIu2rVKvz73/9G//795fxWIoGKKnPU8i7pHSav/CXnlWLLqSxc0C1Y7ZCoJYg7FmvmK+tDZgNBXdSOiOzQjh074OXlZXVDSXzPr1ixQk5zMXXqVFmp7Vyi1LvYZ+vWrXjuuefq3FET84hNnz4dDz/8MPr06VPvlBmim+a2bdvwwAMPYMOGDXXGZpSWluLDD5ViUXfddVedwf7izt6FF16It99+G62JeZCofsM6BWLt0XQ57vH2Czqxeqct2Po+YKwAOo4GItljSWAubJ4m1VK+//775VKf9evX19kmuvKIg0itz12vwxX92+HrbWfw055ENh4dxaGfgJS9ytxMFz6mdjRkxxVHrXXq1Cncfffdco5EQcwJtnfv3jr7zZw5U1YPLSgowGOPPSafn5sbRANLlAsXP6u+XDFixAj5KBpmYnqLc+dCXLhwIUpKSuT6DTfcIJ/XJKqp1jfPYmtgHiSqa1BUAFy0GqTklSIxp0RWYSUVleYDuxYr66MeUDsam8Fc2DwczeyArh6sTNC78mAqissr1Q6HmquyDFj7jLI++iHAixcEWoP4Ihf/Xtp6EZ9rDdE9UnSLFFcrxZi51157DePGjZNXMa3tqiPuCnzyySfyCqoo1S7mfRKJypx8xOtZWVmYNWuWXF+8uOokhIioHh6uOvSP9Jfr29h1VX27lwBl+UBwd6DbJJvNg8yF9oWzuDqgwR0C0DHIE/FZxVh1KBVTBylXSchO7fwMyI0HvMOBkfeqHY3DKqkwoPf8VW3+uYefmQRPV8u/ih999FHZjUVMNi+mf3jiiSewe/du2Z2mKcRcgi+//DJeeeUV2cVFzEUo5oASXS3FldUePXrgmWeekV01/fz8zju9BBGR6Lq6+0yOnO/RPIUYqcBQoXRZNd91tGCKBrXyoMBcaD9459EBiasiUwcpdx9/2p2kdjjUHKV5wIaXlfWL5gGunH7FmRUWFuLTTz/Fq6++ivHjx8vB+6JoixhD0VSiS82NN96Irl27yiIu4jO2b98up5UQU0+I7xORKMW6mGCYiOh8hlYVzTmSko+C0gq1w3FeB38C8pMAr1Cg3/VwJMyF6uKdRwclGo9v/HkCm05mIi2/FGG+tQfTkp3Y9CZQkq10ORl4i9rRODQPvU5e+VTjc60ZcyHmBRw+fHj1tsDAQHlFtKlEARczUUDA19cX6enpTf55ROTcQn3dER3kibisYjnn47geZyekpzYiuoBuriqWMvwuQO9u03nQ/NmWYi5UFxuPDqpjkJcsmb0zPge/7k3CnWNZndMuy2tveU9Zn7AQ0PGfa2sSVxWt6TLjKETZ8XOPQ835CImIrDU0OlA2HnfEsfGoitPrgLQDgN4LGDLL4rc5ax4UmAstx26rDmxqVeEcdl21U+ueBypLgKgRQI/L1I6GbECXLl1kghMlvM1ycnLk3IFERLbUeBR2ncnmnNNq2PSW8jh4OuCp/Fk4EuZCdbHx6MCu6BcBV50WR1MLcDg5X+1wyBrpR4G9XyvrlzwjLoGpHRHZAFFVbvbs2bJQwF9//YWDBw/KcRrnzgFFRKSm7mE+8HF3QVGZQY59pDaUekC586jRASMcs8gec6G6nPPetJPw89RjfK9QrDiYip/3JKJ3RG+1QyJL/bkQMBmBXlOADmf79BOJSnBiIP+UKVPg4+OD//znP8jLy1M7LCKiajqtRg6dWXcsAzvistG3vZ/aITkP81jHPlcBAR3hqJgL1cPGoxMUzhGNx1/2JuPxS3vCRcerMjYvbhNwfIVy1XD8ArWjIRu84vrll1/KxWz58uUWvVfMV1VTffNq5ebmnvc5EZGlXVfNjcfbRndSOxznkJcIHFx6dnoOB8ZcqB62JBycGKge4KlHRkEZNp3KUjscaoz4AlszX1mPuRUI7qZ2RERERFYb1CFA3oFMyClBSl6J2uE4BzGvo7ESiB4DRAxSOxpyULzz6OBcXbSYMiACX2yJx0+7E3Fh9xC1Q6LzOfwrkLRTqZB24Vy1oyE7cuDAAXkltiGie09L6NWrF4YMGVLva2LuKzFZs5gj65133qnzunlMioilvp8RHBxcXQzh2muvrfczJk1SysgvW7ZMLg29TkTq8nZzQZ92vtiflCerrv5rgHPPjdcm80LvWqKsj3oQzoq5sPVzocZU371aG5Ofny8n5hR9mcW8K2SdPWdyMPW9zXDXa7HzqUvkF7olDiSy77gl+kW20FgOQwXw7nAg+xRw4ePARU+0zM+lOkpLSxEbG4tOnTrB3d0x5kAtKSlBUlLDlZXFxMfO7Hx/5vaQY+whRnJu9Z0z/LInCZ9uisXAKH/835V9VYnLYc8Z6psXWvRcCukF3Lul0UJ7jpgHBebC1s+FvPPoBMSXdudgL5zOLMKKAym4bkiU2iFRfXYtVhqOXiEOP1aBWp640unsSZGIbG/co2g8HkzKQ3F5pdPOIdjqKsuBrR8o6+L8wYkrtDMXtj6OeXQCYqJTUThH+HkP53y0SWUFwIaXlHVx19HNR+2IiIiImqV9gAfa+bmj0mjC3gQWHGk1okhOQTLgHQ70q7+bI1FLYePRSVxV1XjccjoLybkcuG6TpbWLMoDALkDMTLWjISIiarG7j4KoukqtQIw+M0/PMfwuwMVN7YjIwbHx6CSiAj0xrFOg/I75ZS/vPtqUgjRgc9WA6vHzAZ1e7YiIiIhaxLCqxuPOuBwYbb/Mhv059ReQfghw9QaGzFI7GnICbDw6kWsGV3Vd3Z1U75w2pJINLwIVRUD7IUDvK9WOhoiIqMX0jvCFh16H3JIKnExvmUqXVIP5ruPgGYCHv9rRkBNg49GJTO7XDm4uWpxIL8TBpHy1wyEh88TZ0tqXPOPUg9yJiMjx6HVaDO6gNGq2s+tqy0rZD5xeB2h0wIh71I6GnAQbj07E112PCb3D5PrS3Ylqh0PC2kWAyQB0nwxEj1Y7GnJw0dHReOONN6x6z6ZNm9CvXz/o9XpcddVVrRYbETkujntsJVuqhrz0mQr4d1A7GrvAPNh8rJnshF1Xl+9Pwe/7kvHk5b3kFUFSyZltwJHfAY0WmLBQ7WjICezYsQNeXl5WvWfOnDkYOHAgVqxYISdenjp1qpwj6lzFxcVyn61bt+K5556Dq6trrdcrKysxffp0PPzww+jTp0+9kzi7ublh27ZteOCBB7BhwwY5ifK581N9+OGHcv2uu+6qM0eV0WjEhRdeiLffrurGRUQ2IaZjAES/mtMZRcgqLEOQN4u6NFteolJlVRh1v9rR2A3mweZj49HJjO0WgmBvN2QWlmH9sQxcUnUnktqYGHP65wJlfeDNQGhPtSMiJxASEmL1e06dOoW7774bkZGR8nlKSgr27t1bZ7+ZM2eioqICBQUFeOyxx+TzmtavX4+VK1fK8dbiZ4nn5xoxYoR8zMjIwG+//SavENe0cOFCOQG0cMMNN8jnNcXFxWHu3LlW/45E1Lr8PV3RPcwHx9IKsCMuB5f2DVc7JPu37QPAWAlEjwEiBqkdjd1gHmw+3nZyMi46La4aGCHXl+5i11XVHPsDOLMFcPEALnpC7WjI3KAvL2r7xcriVePGjcP9998vFz8/PwQHB+Ppp5+2qAjWud11xBywn3zyibyK6unpiW7duslkZU5A4vWsrCzMmjVLri9evLgJB5aICBjaiV1XW0xpHrCz6vt41IP2nwetzIXMg+rinUcndPXgSHzyTyzWHk1DbnG5vCJIbchQCaypuus48j7AV2nMk8oqioHnVfizeCIZcLWuC82SJUswe/ZsbN++HTt37sSdd96JDh064I477rD64xctWoSXX34Zr7zyiuzmcvPNNyM+Ph5RUVHy6mqPHj3wzDPPYNq0aTJJf/DBB1Z/BhHRsOgAfLU1HnsTc1FWaYCbi07tkOyXKLRXXgCE9AS6TrD/PNiEXMg8qB7eeXTSstm92vmiwmCSYx+pje1eAmSdADyDgNEPqR0N2SGR0P773//KhCaSnBgbIZ43hehWc+ONN6Jr1654/vnnUVhYKJOxTqdDeHi4vNIqkqVY9/DwaPHfhYicQ3SQF4K9XVFeacSBxDy1w7FfleXA1veV9ZH3A+eMiXMWzIPq4Z1HJy6c8+zyfPy4OwnTR9buT02tqKwAWP+Csn7hXMDdV+2IyEzvqVz5VONzrSTGRIhkZjZy5Ei89tprMBgMMtlZo3///tXrooiAr68v0tPTrY6JiOh8xHeWqLq64mCqnLJjSFUFVrLSoZ+BgmTAOwzof71j5EHzZ1uBeVA9bDw6qSsHtscLK45iX0KunLS3a2jdik/UCja/AxRlAIGdgZjaA6lJZSIJWdl91BGI0uM1iWQsqrUREbU0c+NRFM0R49NqnvyTBcSYvs1VVTSH3wW4tHDVWuZBiXnw/JzzXjchxMcN47orFad+4pyPbaMg9eyX/vgFgAvHmlLTiDLeNYmy4GKQv7VXW4mI2lL/SD+4umhlxfe4rGK1w7E/p9cDaQcAvRcQcxucGfOgeth4dPLCOcLPe5JgMFpX8ZGaYP2LQEUREDkU6H2l2tGQHTtz5oycd+rYsWP49ttv5QD/hx7i+Fkism2iSM6ASD+5vpNVV623+S3lcfB0wNO5u/0yD6qH3Vad2PheofB1d0FKXim2nMrCBd2C1Q7JcWUcA3Z/oaxf8n9K1xCiJpoxY4ac52nYsGHyKqtImKLSHBGRPXRdFd1WxZQd1w2JUjsc+5F6EDj1F6DRAiPugbNjHlQPG49OzF2vw5QBEfh62xks3Z3IxmNr+nMRYDIAPa8AOo5UOxpygPEZYp6q99+vqrhnITFnVU31zYmVm5t73udERM1tPAKncDS1AHklFfDzqD3ejBqw5R3lUfRcCmChQ+ZBO+u2+u6778pJNt3d3TF8+HBZztYS3333nRyEetVVVzXlY6kVXBOjdF1deTAVhWWVaofjmOI3A8eWAxqdMtaRiOwe8yBR0wR7u6FzsBfEKfuu+By1w7EPuQnAgR+U9VEPqB0NOTmr7zx+//33so+xmCBTJEzR6p80aZLscxwaGnrelv4jjzyCMWPGNDdmakGDovzRKdgLsZlFWHEghV1IWpq4orX6aWU95lYgpLvaEZGD2rhxIyZPntzg62LeqpbQq1cvDBkypN7XxPxXIg+IebLeeafqKvk5c2lptVoZS30/IzhY6f3QpUsXXHvttfV+hsg3wrJly+TS0OutiXmQqPl3H09nFsmuqxf3bPjfDNW462isBDqNBdrHqB2NzWIebJs8qDHVd7/2PESiHDp0aPUBEaVsxUSdYnLOuXPn1vseMefK2LFjMWvWLPkHK27//vLLLxZ/Zn5+vpycMy8vT869Qi3rnb9O4NXVxzGicyC+u/Nsl0pO4muZflWD/xucj+mHmUpltAf3AD5hbRkaNaC0tBSxsbHo1KmTvHPkCMTYj6SkpAZfF5MfO7Pz/Zlbm2OYB4nqsuac4VhqAR75cR88XXX4evZwuOicp37jec8Z6lOUCfy3L1BZAkz/BehyUYvEwTzonEpbIBdadeexvLwcu3btwrx586q3idbzhAkTsGXLlgbf98wzz8hW+OzZs2XSbExZWZlcav4y1HqmDo7Ea2uOY+vpbCRkFyMq0PpJy6kelWXAnwuV9dEPsuFIrUpc7WRibH3Mg0TN1y3MW451FGMeD6fko3+kv9oh2a5tHyoNx3YDgc7j1I7GpjEPtg2rLvVkZmbKq6dhYbVPgsXz1NTUet/zzz//4NNPP8XHH39s8ee88MILsuVrXsQVXWo97f09MLJzkFz/ZU/DV2zISts+AHLiAJ92wMj71Y6G6mFlxwuyYy31Z808SNR8Wo0GQzoGyHXRdZUaUFYAbP9QWR8zp1UqtTMPOhdTC/x5t2o/gYKCAkyfPl0mTHMfXkuIK7rilql5SUhIaM0wqcacj6LqKr9IWkBhOrDhFWV9/HzAzVvtiOicKm1CcTEnqXamO4ZCW08gzTxIdL6qq6LxyKI5Ddr5OVCaBwR1A3pOadEfzTzonIqr/rzNf/5NYVW3VZH4ROJNS0urtV08Dw8Pr7P/qVOnZIGAKVPO/oUXY0PkB7u4yOICYkDoudzc3ORCbWdy33As+PUg4rKKsS02GyOq7kRSE617DigvULqZ9L9B7WjoHOJ7zN/fH+np6fK5p6enrIBJjknknYyMDPnnLHJPczAPErWMQR384aLVICm3BEk5JWgf4KF2SLY39GXLu8r66IdE//gW/fHMg87FZDLJhqP48xZ/7s25kGpVFnV1dUVMTAzWrl1bXWZcJEHx/P7763bL69mzJw4cOFBr21NPPSWvxL755pvshmNDvNxc8K+BEfh2ewK+35HAxmNzpB4Adn+hrF/6Yot/4VPLMJ/omxMnOTYxLrFDhw7NPjliHiRqGZ6uLujb3g97E3KxIz4b7QPaqx2Sbdn3LVCYCvhEAP2ntcpHMA86H39//3ovdFrD6kuwojz5rbfeKsvLDhs2TJYoLyoqwm233SZfnzFjBtq3by/Ha4gqPn379q0TtHDudlLftKEdZOPxjwMpWDilj9rh2CfR5XflPMBkBPpMBTqerV5LtkU0Itq1ayeLmFRUVKgdDrUy0egTDciWwDxI1DKGRgcojce4bFw1kI3HakYDsOlNZX3U/YCLa6t8DPOgc9Hr9S0ydMPqxuO0adNk95/58+fL4gADBw7EypUrq4sHnDlzpsUSNLWtAZF+6Bnug6OpBfh1XxIGRSmD2ckKx/4A4jYCOjdgwiK1oyELiC/Sth4HR/aNeZCoZYhxjx9vjMWh5HwUlVXKXlAE4PCvQPZpwCMAGHxrq38c8yC16jyPauD8Vm3n802xWPT7YfRq54uXru7H/u/WzNlUWQ68N1z5wr9gDjBhgdqhEZGD5Bh7iJGcW1Pnhr7n611IzCnB45f2xAVdLS8q5bDzPIox0R+OAdIOAhfOBS46Oy0QkS3kGV4apVqmDmoPVxctjqTk41RGkdrh2JftHykNR69QpaQ2ERERWVZ1NZZTdkjHlisNR1cfYPhdakdDVAcbj1SLv6errLwqrD5c/5xl1NDUHC/XmJrDR+2IiIiI7KbxuDM+GwajzXeGa12iM+CGl5T14XcCnsqxIbIlbDxSHdOGKtX/NhzPQGmFQe1w7MOaBUBZHtBuADDwJrWjISIisgu9wn3g5aZDfmkljqbmw6kdW6FUbHf1BkbWrd5MZAvYeKQ6RnQKQscgTxSXG7DpZKba4di+M9uAfd8o65e9Bmg56JyIiMgSLjpt9d3Hzaey4Nx3HV9U1ofdwbuOZLPYeKQ6tFoNrh+i3H1cfbj2RNhUTzntP/6jrA+6BYgaqnZEREREdsVcKGfzqUwYbb+OY+s4vgpI2QfoPXnXkWwaG49Ur+tiIqHVAIdT8hGfxcI5DQk8+pXSxcTdj1NzEBERNYGYGsxDr0NmYTlOpBXCqcc6Dr0d8HL8qrNkv9h4pHqF+rpjROcgub78QIra4dgkXUkWwne+ojy5+Gl+2RMRETWBqPJu7rq66ZQTDpc5+SeQvBtw8QBGPah2NETnxcYjNejyfu3k4/pjGSgur1Q7HJsTvuNF6MrzgfD+wJBZaodDRERkt0Z3VS5Yi1oLdjAFecsRv+u655X1obMB7xC1IyI6LzYeqUH92vshKtATJRUG/HU0Xe1wbIpn2i4EHv9eeXI5i+QQERE1x+AOAXBz0SK9oMy55pk+uly56yjGOvKuI9kBNh6pQRqNBpdXzfn4x4EU57oSeD7GCrT/Z65cze52HRA1TO2IiIiI7Jq7XochVV1X/3GWSu+i6N5fzyrrI+4BfMLUjoioUWw80nld1DNUDmJPyCnB/qQ8tcOxCSEHPoJ7zjFUugUgdfiTaodDRETkcFVXneKC9YEfgYwjStG9UQ+oHQ2RRdh4pPPydHWRDUhh+X4WznHNi0Po7jfkesqI+TC4cx4mIiKilhDTIUAWz0nJK0VspoN3Xa0sB9ZXjXUc/RDgEaB2REQWYeORGnVZVdfVbbFZyCwsg9MymdB+0zxoDWUoiLgAuV2vVjsiIiIih+HhqpMNSKfourrnSyAnDvAKBYbfrXY0RBZj45Ea1THISxbPMZqAFQdT4az8Ty6Fd/ImGHVuSL7geTEoVO2QiIiIHMqYbkrX1Q3HMxy362pFCbDhZWV97KOAq5faERFZjI1HsmrajpUHU1BWaYCz0ZVmo93W/5Pr6YMfRrlvtNohERERORwx36OotSCqrh5JLYBD2vYhUJgK+HUAYm5VOxoiq7DxSBYZ0TkIoT5uyC+tdMppOyK2LIRLWQ5KAnoio9+daodDRETksFVXR3ZR5nxcf8wBzzeKMoGNrynrF80DXNzUjojIKmw8kkV0Wg2uHBgh13/dmwyjo3YlqYdv7Ar4n/oFJo0WSWNeBrR6tUMiIiJyWBf1UAr1/XMiExUGIxzK+heAsnyg3QCg/w1qR0NkNTYeyWITeoXBy1WHpNwS7IjLhjPQlWSh/aYn5HpG/3tREjpQ7ZCIiIgcmqizEOjpioKySuw5kwNH4ZZzAtj5ufJk4nOAlqfhZH/4t5asmrbj0r7K2Mef9yTBGURsfhoupVkoDeiB9MEPqR0OERGRU/R2MhfOWXcsA44ifPvzgMkA9LgM6DRG7XCImoSNR7LKlP7t5Jf6oeR8HE9z0IHsVfxO/w7/2GUwaXRIuPB1mHQcl0BERNQWxlV1Xd0em43i8krYO6+kf+CbsBbQugCXPKN2OERNxsYjWSXI2w0Xdgtx+LuPLsUZiNj0lFxPH3g/SoP7qR0SERGR0+gS4oXIAA+UG4zYfCoLds1oQLttzyrrQ2YDwd3Ujoioydh4JKtdNUgpnLP5VCZS8krgcEwmtP9nrlJdNbA3MgY+oHZERERETkWj0VTffVx7JA32LPDYt/DIPgyDqy9w4eNqh0PULGw8ktU6BXtjcIcAGE3Aj7sS4WgCj3wF3zNrYNS6IlF2V3VVOyQiIiKnM75nKLQa4GByPpJySux2nuiwnS/J9bTBcwAvZRoSInvFxiM1ybShUfJRzPmYXlAKR+GWcwzttiljEVKHzUVpUG+1QyIiInJKwd5u8mK1sMZO7z6G73gRLmV5KAnshazeM9QOh6jZ2HikJundzhf9I/1QaTRh6W7HGPuoqSxF1F8PQGsoQ0HkOGT1maV2SERERE5tYp9w+bj2aBoq7WzOR4/03Qg89p1cTx71f0qxHCI7x8YjNdkNQ5S7j2sOpyKrsAyOUELbI+coKtyDkTj2NUDDfx5ERERqGtoxAP6eeuQWV2BnvB3N+Wg0oH1V4b2cbteiOHyY2hERtQieHVOT9W3vJ+9AVhhM+MnOK6/6xK9G8OHFcj3pwldR6alUlCUiIiL1uOi0cuyjsPpwKuxF4NGv4ZF1UBbJSR06T+1wiFoMG4/UrEpo5rGPKw+lIruoHPZInx+PqA1z5Hpm39koiLpY7ZCIiIioyiW9lK6ru+JzkFFg+z2d9IVJCN/xglxPi3mEF6TJobDxSM0yKMofPcN9UF5pxPc7E2CP4xw7rr0buvJ8FIXG8OogERGRjWkf4IF+7f1klfcVB1Ng+9N9PQFdRZE8r8jqNV3tiIhaFBuP1Oy7jzNGRsv1VYdSkZxrX6W0IzY/DY+sQ6h0D8KZ8e9xWg4iIiIbdEX/dtXnGuKCta3yP/kzfBLXyem+ksa8DGh1aodE1KLYeKRmE1cDYzoGwGA04ettZ2AvAo59h8Dj38Ok0eLMRW+j0ktJTERERGRbhncKQoiPG/JLK/H3iQzYIl1JJtptXSjX0wc/jLKAbmqHRGQbjcd3330X0dHRcHd3x/Dhw7F9+/YG9/34448xZswYBAQEyGXChAnn3Z/s04wRHeWj+EI/lVEIW+eRvkfedTSPRyhqf4HaIRGRHWEeJGpbOq0Gl/VVLvL+vj8ZJpMJNsVkkucVLmW5KAnsjYz+d6kdEZFtNB6///57zJkzBwsWLMDu3bsxYMAATJo0Cenp6fXuv379etx4441Yt24dtmzZgqioKEycOBFJSfZdnZNq6xzijbHdlAHhX2yJgy3TFyaj45rb5XyO+R0uQcaAe9UOiYjsCPMgkTom9g6Dq06L0xlFOJySD1vrruofuxwmjQ6JY18BtHq1QyJqFRqTlZduxBXWoUOH4p133pHPjUajTIQPPPAA5s6d2+j7DQaDvPIq3j9jxox69ykrK5OLWX5+vvyMvLw8+Pr6WhMuNcOBxDyr9k/JK8E9X++W3VcXTukju7LaGm1FETr/fg08sg+jJKAnTk/5CUZX72b9zH6Rfi0WHxG1PZFj/Pz8LM4xzINEzT9naKq3/jqBNYfTMLpLEOZO7gVboC9IQLefLoWuogBpMf9B+qCHGtyX5wxk77nQqjuP5eXl2LVrl+xyU/0DtFr5XFxNtURxcTEqKioQGBjY4D4vvPCCDN68iIRJtq+dnweu6Kd0Kfl442lUGGxsQLvJiMj1D8uGY4V7MOInftbshiMRORfmQSJ1/at/hHzcfCoLSTk2UKTPaEDU+odlw1FUV00fcJ/aERG1Kqsaj5mZmfKKaVhYWK3t4nlqqmUTtz7++OOIiIiolXjPNW/ePNnqNS8JCfY3BYSzunFYB/h76JGUW4Jl+5NhS8J2vgy/+FWyAtqZSz5ChU+k2iERkZ1hHiRSV3SwF4ZGB0B0m1u6J1HtcBCy/314pe2AQe+NhHFvAFoXtUMicpxqqy+++CK+++47/Pzzz7LIQEPc3Nzk7dKaC9kHLzcXzBipFM/5dnsCcorKYQuCDn2O0H3vyfWkMS+hOGyI2iERkRNiHiRqvutjlDvx646mI7PwbPfutuaRvhthu16X68kjF6HCVzn/IXJkVjUeg4ODodPpkJaWVmu7eB4eHn7e97766qsyaa5evRr9+/dvWrRkF8b3CkPXUG+UVBiwxAaK5/id+g3ttiils8VYhNxu16gdEhHZKeZBIvX1bOeLvhG+qDSa8PMedQpP6Uqy0GHtPdCYKpHb6QrkdrtWlTiIbLrx6OrqipiYGKxdu7Z6mygUIJ6PHDmywfe9/PLL+L//+z+sXLkSQ4bwjo+j02o0uGtsZ7m+9mg6Dia1zSD6+nglbUTkhn9DAxOyet+K9IEPqhYLEdk/5kEi23DdEOXu46pDqcgrqVBhnOODcC1KQZlfZ9mjCRpN28ZAZC/dVkV5cjFn1ZIlS3DkyBHcc889KCoqwm233SZfF5XjxFgNs5deeglPP/00PvvsMzknlhgTIpbCQtufC5Carme4Lyb1VsYEvbPuJMor2754jkfGPnT8805ojRXI7XQ5kkcs5Jc7ETUb8yCR+gZF+aNriDfKKo34aXfbjn0M3fMGfJI2wujigfjxH8Lo6tOmn09kV43HadOmya438+fPx8CBA7F37155JdVcPODMmTNISUmp3v/999+X1emuvfZatGvXrnoRP4Mc28zRnRDo6SqL53y340ybfrZ75gFEr7wFuooiFEaMQqIcxK5r0xiIyDExDxKpT6PR4KbhHeT6sv0pyGqjsY8+8asRuuctuZ40+gWUBfZok88lstt5Hu1hDi6ynTmbtpzKxPMrjkKrAd6YNhCdglt/agz3zIPotOJGuJTlybLZcZd+2apTcnDOJiL7Zg85xh5iJOfWVvM81iROYR//6QCOpOTj0j7huO+irq36ee5Zh9H596uhqyxGVq8ZSB79rNU/g+cM5FTzPBJZa2SXYIzqEgSjCXh9zfFWn/tRfLF3WnGTbDgWhw5G3KVfcC5HIiIiB737eGtVhfc1R9KQnNt68z66FKej4+pZsuFYGDEaySMXtNpnEdkyNh6p1d19YRf4eegRl1WML7bEt9rnuGfuR6c/xB3HXBSHDESsbDhyHAIREZGj6hPhh5iOATAYTfiilSq8aypL0XHNHXAtSpYFcuLHvw9o9a3yWUS2jo1HanUBnq548GKlK8kve5OwNyG3xT/DK3kTOi+fBpeyHBSHDECs7KrKrl1ERESO7taR0XJ4zKZTWdif2MLnGMZKdPjrXnhm7EGlmx/iJn4Oo5t/y34GkR1h45HaxLBOQZjcV5kD7b9rjrdoWW3f2D8QvfJWpThOu1GInfwNjG4cU0BEROQMOgV74dK+7eT6xxtPy7uQLcJkROTGx+B75k8YdW6Iv+QTlPt1apmfTWSn2HikNjNrdCe09/dAdnE5Xl19rEW+3AOOfi2vCGqN5ciLnoy4SYvZVZWIiMjJ3DysA3zcXOQQmRUHz1Y7bjKTCeHbn0fAiR9h0uhw5uL3UBw+vCVCJbJrbDxSm3HX6zBvck+4uWhl19WvtzVj/KPRgPBtzyLyn3nQmIzI7nGj/GI3ubi3ZMhERERkB3w99JheVTznq63xzZu6w2RC2K5XEXLgI/k0cczLKOh4SUuFSmTX2HikNtUxyAsPXtxNrv+wKxFbTmdZ/TO05QXouGZ29Zd62qCHkHTBi5zHkYiIyIlN7B2O7mHeKCo34L31p+RUHk1qOO58GaF735ZPk0csQG7361o+WCI7xcYjtbmx3UPwrwER1eMfT2cUWvxefX48uvw+Fb4Jf8nxB2cuegfpMf8R9bpbMWIiIiKydTqtRl6gdtFqsD0uG3+fyLS+q+qOFxC67135NHnEQmT1nd06wRLZKTYeSRW3jYrGgEg/lFQYsGjZYWQUNN69xDduFbr9cjncc46jwjMUp6/4EXld/tUm8RIREZF99HC6fkiUXP/w71PILiq37I3GSkRsmoeQ/R/Ip8kjn0FW31mtGSqRXWLjkVThotNi7uRe6BDoKb/YF/1+CEVllfXuqzGUo93WRej45x3QleejOGQQTl65DCUhA9o8biIiIrJt18ZEonOwFwpKK/H6mmMwNtJ9Vczj2GHt3Qg6+g1M0CBp9PPI6jOzzeIlsidsPJJqvN1csGBKbwR46hGfXSzvQJaUG2rt45ofh87LrkXwwU/l84x+d+L0FT+g0kuZ9oOIiIioJr1Oi0cn9ZAF+vYl5mHprsQG99WVZKHTipvgF79aGQ4z/gNk97qlTeMlsidsPJKqQn3csWBKH3i56nAkJR/PLj+MskqDnFsp8PASdPtpEjwz9ioT817yCVKHPwWTzlXtsImIiMiGRQZ44u6xXeT6V9vicTApr84+7pkH0PXXK+CVthMGV1/ETv4K+Z0mqxAtkf1g45FU1yXEG4v+1Rceeh32J+Xho9/Wo+Pym9B+89PQVpagsN1InLxqBQo6TlQ7VCIiIrIT43uFYlz3EIhppV9YcQRp+aXVr/md+hVdfr8GroVJKPPthFNTfuY8jkQWcLFkJ6LW1iPcBwsv74bTy/+LB7J+gLemFAadO9KGzUNW71sBDa9zEBERkeU0Gg3uu6grEnNKcDKjUPZuevmqHui85yUEH14s9ymIvAhnLnoLRjc/tcMlsgs8Iyeb4Jm6DVduuxHzdF/KhuMOY3fMcH0dJ6JvZsORiIiImsRdr8OTl/eS9RVcs48h9LvJ1Q3H9AH3Im7iZ2w4ElmBdx5JVfqCBITtfBUBp36WzyvdAnCw1xzcta87snMMOPbDPjx5WS90DfVWO1QiIiKyQ8GeLvi81270PPgq3A0VyNP6I3P8GyjueLHaoRHZHTYeSRW60hyE7n0bgYe/gNZYLktjZ/e8EWlDHofOPQAvdivBot8PIym3BI//tB8PXdwNY7uHqB02ERER2RG37GOI/OdxeKbvBjTABuMA/Kf0bgw8GYn7o0zQaTVqh0hkV9h4pDalK81F0OHPEXzgE+gqCuS2wohRSBn2BEqD+1fv187PA69eNwCvrj6GXfE5eGX1MRxPK8Cto6JlCW4iIiKihmgqSxC6910E738fWmMFDHpvpA6diyP6S5G9+jj+PJKOskoj/j2hO88riKzAxiO1CZfidAQf+BiBR7+CrqJIbisJ7I3UYfNQ2H6sGNVe7zyQT1/eG19ujcfS3Yn4dV+yrMb6n0u6o2OQlwq/BREREdk0kxH+J39G+I6XoC9OlZvyOk5E8qj/Q6VXO4wG8JhGKy9ObzyRKeeXFnNCerrylJjIEvyXQq3KPesIgg4vhv/Jn6A1lFU3GjMG3oe8Tpc3WgxHdCeZOSoavdv54M21JxCbWYQ5/9uHW0Z0wL8GtGd3EyIiIgJMJngnbUTYzpfgmXlAbir3jkTK8KeQHz251kXq0V2D4abX4oU/jmJnfA4e+3E/nrqiN8J93VX8BYjsg8ZkMplg4/Lz8+Hn54e8vDz4+vqqHY7TOJBYd0Jdixgr4Be3EkGHl8ArdXv15uLQwUgfeD8KosbXe6exMTlF5Xhj7XHsPpMrn3cK9sJ947rKaT7U1C+SVdqI7Jk95Bh7iJGcW5PPGZrLZIJP4jqE7n4Tnhl75CbRRTV94API6nMbTC4NNwjFcJjnlh9BdnE5fNxd8OjEHhjUIaBVw+U5A9kqS/MMG4/UMonAZIJ71kEEnFgK/1O/wqU0S9ms0SEv+lJk9Z6J4vBhTWo01v4YE1YfTsOSzXEoKKsUY99xcc9Q3DS8A0J91LliyERAZN/sIcfYQ4zk3Nq68agxlMMvdjmCDn5SfafRqHNHdq+bkT7gPhg8gi36OVmFZXj2jyM4mV4on185IAIzRkbD1aV1xkHynIFslaV5ht1WqVlc8+LgG/cHAk7+BPec49XbKzxCkN3zJmT3vBmVXuEtOuHvpD7hGNE5CJ9visXao+ly+ftEBi7vF4FrYyLh56Fvsc8jIiIi2+FSlIrAo18j8Og30JdkyG1GFw9k9ZqBzH53otLTusrsQd5uePHqfvj0n1isOJgq6yvsTcjF/Rd1Rc92vFBDdC7eeSTrriLKO4wH4Be3Gr7xK2s1GI06N+R3nIjcbtegQBTB0bb+tYljqQVYsiUOB5KUWF11WozvFYqrBrZHhL8H2gKvIhLZN3vIMfYQIzm31rzzKCqn+satQsDJpXJco8ZklNsrPMOQ1Wu6vFBt8Ahq9udsj83GW3+dQF5JhXx+Sa8wzBjZEf6ermgpPGcgW8Vuq9RiicClOANeyZvgnfyP/NJ2LUqp3sekcUFhuxHI63yFLIBjdGv7L0XxV1iMg/xqW3x1txNRR2d4pyBM7BOGQVEBrVpYh4mAyL7ZQ46xhxjJubV041FTUQyfpL/hG78KvnErqyu1C0VhQ5HVZ6YcFgNty/Y2Eg1HMTRmzZE0+dxdr8WU/hHyorRvC/Rs4jkD2So2Hqnp8pOBhO3IPLIRXkn/wCPnaK2XRfeQgshxyO84CfkdLobRzR+2QPxVPpiUh6V7kuTckGbB3q4Y3ysMY7oGo0Ogp+z62pKYCIjsmz3kGHuIkZxbSzQe9QWJ8JYNxtXwTt5UXaVdKPPpgNyuV8ul3C8are1oSj4+3Hi6+qK0h16Hy/qF47K+7RDajKqsPGcgW8XGI1mmNB9IPwwk7QYStwMJO4D8xDq7lQT1RWHEaBS2H4Oi8GHnrV5mC+KzimRhnXVH02VhHbPIAA+M6hKMEZ0C0SXUG9oWaEgyERDZN3vIMfYQIzm3pjQexfhFr5St8E7eDK+UzXArOFPrddFgLOg4EXnRk1EcNqTZRfesJU6Rt8Vm45vtZ+RUYYLoyDSsUyAm9g7HoCh/uOisK6zDcwayVWw8Um2V5UBuPJB+BEg7CKQdAlIPKNvOJeZeDOuDLP8BKGo3QjYaW2IsgRrKK43YejoLG45nYPeZHFQaz/51F2W5+0f6yy///pF+cn6nptyVZCIgsm/2kGPsIUZybo01HrXl+fDI2A/PzH3y0SNjH1yLkmvtIyq0l4QMQH7UeFlDoSyge5s3GOtjNJnkeMhl+5Oxr8bv6ePmglFdgnBBtxD0ifCF3oKGJM8ZyFax8eiMyouB/CSlQZh1Gsg6CWSfArJOKduqBpjX4RMBtOsPRA4FooYBEYMBN2/15mxqJcXlldgRl4NNJzNlJbWSCkOt133dXdA9zEfOGykexTyS/h76RhuUTARE9s0ecow9xEjOzXzOoC0vhFvuSbjlnpBF9ZTHE3AtTKjzHhM0KA3qg8KIUSgSS9hQGF3Vnbu5MQnZxVhxMAUbT2Yit1gprCO4uWjRr70fBnXwl48dAr3qrbfAcwayVWw8Otpdw+JMoCijaslUxiWKhmKeWBKVrqYlZ8f51UvvCQR3A8L6yTuLCO8LhPUFPAPr3d3RGo81VRqMOJFeKBuRexJycSKtoNZdyZp3J8U4SbF0DPREuJ8Hwnzd5JyS5jmgmAiI7Js95Bh7iJGcgDhlLCsACtOA3DPKhWnxmBOP4vRY6AsSoC/NbPDt5d5RKA7pL+8ulgSLpa/NNxYbYjCacDA5D38fz8COuGzk1GhImhuTXUO95cXoLiHe6BDoIavAD4mu/5yLyKEbj++++y5eeeUVpKamYsCAAXj77bcxbNiwBvf/4Ycf8PTTTyMuLg7dunXDSy+9hMsuu6zFfxm7+MItywdK85SxhuKx+nmNpSRbaSCaG4tim6XEl7BfJBDYGQjqUrV0BQK7AD7hVnX/cOTG47kqDEaczijCsbQCHE8rkI3JlLxSnO8fR6CXK8J83dEt1BshPm7K4l31WLWIO5fWjocgorbVlBzDPEgOc3G6NLfq3KPqUTwXF6OLs5RGYmF61WPVemVpoz9WzPVcFtANpf7daz0a3B2z4SROpeOziuXwGHFBWkwjdm7vJkGcgUUFespGpajBIBqT7f3PPgZ5u1rU9ZWoNViaZ6yeiO/777/HnDlz8MEHH2D48OF44403MGnSJBw7dgyhoaF19t+8eTNuvPFGvPDCC7jiiivwzTff4KqrrsLu3bvRt29ftHm3TvGlZ6gADOVVi3m9vm0NrZcBFaVARXHVUgKUFymPcik+57Wq9fM2RRqh0QFewYBXiPLoHQ74tVcair6RZ9fdeResKcSXteiuKhazskoDEnNKZEI4k10su6qk5ZciraAUpRVGZBeVy+VISv55f7aXq06W9/bz0MPXXQ9fDxf5XKx7uupkBTcP8Whe1+vg7qqDZ9V2d71OxqfXaaoea6+35jQkRORgeZBsl9FY+3zDaMX5iTwnKapxzlFU9VhSY71qEecr5ovV8tykCcwXqgM6Av4dAP+OiDcEo9wnCuU+kTZThb2tiOEt0cFecrl6cKQcIynOH46nFsiL0qKIX0JOCQrLKuX5hFgaIobQBHm7IcBTj0AvNwR5ucLPUw8vVxd4ueng7eYCTzcXeLvpqra5yHMFMc+1uNvpal50Wl68plZh9Z1HkSiHDh2Kd955Rz43Go2IiorCAw88gLlz59bZf9q0aSgqKsKyZcuqt40YMQIDBw6Uibc+ZWVlcqnZEhaf0ewrrl9dA5z8E6oS8xGJBp67r/LoVvVYc/EIqNFQrFrc/QFt234JONOdR2uIfzL5pZVKQzK/VH5JZxSUKUth1WNBGbKLy+UN59Ym2o4iQbhWNSrFul6rgVYsGo1sXIobzjqN8lzZjlrrZ1+r2n7Oa2dvWCsr4rl5k/k1jfifpu42845n95dbz9lP2V5zP+Wtmnr2q/GzVWYDdRxsIAZ1A+gU7Ik7x3Zp1s+w9q6eXefBhO3Anq9qX8ys/qKqeqz1vWWqf59632dq2vus/llo+Rgae5+oGSCeGw2ASSzGqnWx3aA0/Mzr1a8ZzvOec99ftY9a3MznH+LRX1mXF6rDAO/Qqseqda9QwNWzzo/gOUPj5w5ijKSLi0b2dErOLalaSpGUW4LU/FLZFbYliRxubki66ZUGpsjzMudrNfJRV+NcQdmGWttqvV4jQZ97HnDuOcC5+f/c3F8n759zLnHuz7a1fKlpsZ/VMj/oict6wsddb3t3HsvLy7Fr1y7MmzeveptWq8WECROwZcuWet8jtosrtDWJK7S//PJLg58jrs4uWrQILU7nWrUizrbdlOc6vdKgM6/XejzPut5D+fIU4wjFuv7cdfG6V9Xzqm3iy1hMcaH+2Z5FOJav+eMqRSMzv6QC+aUVcuLh/JLKGusVsltLaYUBxeUGlIil6nlJ1bbSqm2VBhMqjEZUGEx1kot4KqrKioXI2QyLDmx249Gp8qAopLZ7Scv/XGp5ovJ5nfOPBs5JxDmNOM+Q5x1V5yPV5yg11s3nJfJitT/g4a9cxNbqmh0uzxksJ6YMO5fI7eK8IKuqV1N2URmyiyrkoziXEHcti6oWZd1QvS7OG8oNynlAzVMEsS56SokFpWenLSPHM+eS7vBpo1n0rGo8ZmZmwmAwICwsrNZ28fzo0doTyZuJ8SD17S+2N0Qk5ZqJ1nzFtdmuW6J8QbbAlyRRY8QdQDEuUiwtyWg825CsqDTWWq80iuRhkmM4RbcZkTjE1U6RlMS6sk15Li+I1/OaXJfbzK8pmUg8VF+Hr75BYaqxfvaFmvuZOzeI/9bat+pKrJny889+Vt331H7NnthBXbJa7Cncdv4ebfp5dp8Hw/sDFz9V9aTGhczqi5rndB+od1tT32fpz0IrxtDE94mGnHkR5xByvep8Qt4a0dV9TT7XWv6azq12I5HnKk5F3NkL8HKVS3MvXJsbkmWVtR/FuYFB5Hmjkt+VddSzTXk0nwMYxI1ycd5QlfPPzc3V5wfn5HuR+2qdDzSQ98/N+W2Z79siP5va4PcQw6DaitVjHtuCm5ubXFqcS8uexBOpQXQrcdPq4Cb+9bbCPxMicuA8KKpsi4WIHPbCtVg8ecpLrcSqQXTBwcHQ6XRIS0urtV08Dw8Pr/c9Yrs1+xMREdkq5kEiInJmVjUeXV1dERMTg7Vr11ZvE4UCxPORI0fW+x6xveb+wpo1axrcn4iIyFYxDxIRkTOzutuqGINx6623YsiQIXJOK1GiXFSRu+222+TrM2bMQPv27eVgf+Ghhx7ChRdeiNdeew2XX345vvvuO+zcuRMfffRRy/82RERErYx5kIiInJXVjUdRcjwjIwPz58+Xg/1FqfGVK1dWFwM4c+aMrDxnNmrUKDmn1VNPPYUnnnhCTo4sKsxxbisiIrJHzINEROSsrJ7nUQ3WzsFFRETkSDnGHmIkIiL7ZWmeadtZ54mIiIiIiMgusfFIRERERERE9jnP47nMPWvF7VQiIqKWZM4ttjyKg3mQiIhsIRfaReOxoKBAPkZFRakdChEROSiRa8R4D1vEPEhERLaQC+2iYI6YQys5ORk+Pj7QaDTNalGLxJuQkMCCA43gsbIMj5PleKwsx2PVtsdJpEGRLCMiImpVSXXEPCjw75dleJwsx2NlGR4ny/FYtf1xsjQX2sWdR/ELREZGttjPEweXfxEtw2NlGR4ny/FYWY7Hqu2Ok63ecWytPCjw75dleJwsx2NlGR4ny/FYte1xsiQX2uYlViIiIiIiIrIpbDwSERERERFRo5yq8ejm5oYFCxbIRzo/HivL8DhZjsfKcjxWluFxahoeN8vwOFmOx8oyPE6W47Gy3eNkFwVziIiIiIiISF1OdeeRiIiIiIiImoaNRyIiIiIiImoUG49ERERERETUKDYeiYiIiIiIyHkbj8899xxGjRoFT09P+Pv717vPmTNncPnll8t9QkND8eijj6KysrLWPuvXr8fgwYNlFaOuXbti8eLFcHTHjx/HlVdeieDgYDnh6AUXXIB169ZZfeycwfLlyzF8+HB4eHggICAAV111Va3XeZxqKysrw8CBA6HRaLB3795ar+3fvx9jxoyBu7s7oqKi8PLLL8OZxMXFYfbs2ejUqZP8+9SlSxdZQa28vLzWfs5+nGp69913ER0dLY+F+He4fft2tUOyKcyDTcc8aB3mQssxD54fc6Ed5EGTg5o/f77p9ddfN82ZM8fk5+dX5/XKykpT3759TRMmTDDt2bPH9Mcff5iCg4NN8+bNq97n9OnTJk9PT/kzDh8+bHr77bdNOp3OtHLlSpMj69atm+myyy4z7du3z3T8+HHTvffeK49DSkqKxcfOGfz444+mgIAA0/vvv286duyY6dChQ6bvv/+++nUep7oefPBB0+TJk0WFZ3lMzPLy8kxhYWGmm2++2XTw4EHTt99+a/Lw8DB9+OGHJmexYsUK08yZM02rVq0ynTp1yvTrr7+aQkNDTf/5z3+q9+FxOuu7774zubq6mj777DP5b++OO+4w+fv7m9LS0tQOzWYwDzYd86DlmAutwzx4fsyFtp8HHbbxaPb555/XmzTFl5dWqzWlpqZWbxNffL6+vqaysjL5/LHHHjP16dOn1vumTZtmmjRpkslRZWRkyC+0v//+u3pbfn6+3LZmzRqLj52jq6ioMLVv3970ySefNLgPj1Pd49GzZ0/5BXdu0nzvvffkyUfN4/L444+bevToYXJmL7/8sqlTp07Vz3mczho2bJjpvvvuq35uMBhMERERphdeeEHVuGwR86B1mActx1xoHebBpmEutK086LDdVhuzZcsW9OvXD2FhYdXbJk2ahPz8fBw6dKh6nwkTJtR6n9hHbHdUQUFB6NGjB7744gsUFRXJbiUffvih7GYSExNj8bFzdLt370ZSUhK0Wi0GDRqEdu3aYfLkyTh48GD1PjxOZ6WlpeGOO+7Al19+KbstnUscq7Fjx8LV1bXWsTp27BhycnLgrPLy8hAYGFj9nMdJIbov7dq1q9b3s/i3KJ478vdzS2MerB/zoOWYCy3HPNh0zIW2lQedtvGYmppa64tMMD8Xr51vH/GFV1JSAkck+uD/+eef2LNnD3x8fGQf6tdffx0rV66U4xgsPXaO7vTp0/Jx4cKFeOqpp7Bs2TJ5fMaNG4fs7Gz5Go+TQvRwmDlzJu6++24MGTKk3n14rOo6efIk3n77bdx1113V23icFJmZmTAYDPUeC2c6Ds3FPFg/5kHLMRdahnmw6ZgLbS8P2lXjce7cufJL/XzL0aNH1Q7Tro+d+IK777775BXWjRs3yoG3YuD7lClTkJKSAkdn6XEyGo1y/yeffBLXXHONvBr9+eefy9d/+OEHOANLj5X40i8oKMC8efPgjJryvSWu5F966aW47rrr5JVqIjPmwaZjHrQcc6FlmActx1zoOFxgR/7zn//IKzfn07lzZ4t+Vnh4eJ2KRKJLgfk186N5W819ROU1UQHKEY/dX3/9Ja8citv+4vcU3nvvPaxZswZLliyR//gtOXb2ytLjZD6B6N27d/V2UYlQvCaqygmOfJys/TslulCI41OTuPp68803y79XDf1bc4RjZe33VnJyMi666CJZJfOjjz6qtZ8jHydriAqYOp2u3mPh6MeBebDpmActx1xoGeZByzEXOk4etKvGY0hIiFxawsiRI2UZ8/T0dHl1URCJQSQK85eg2OePP/6o9T6xj9hubyw9dsXFxdX9pmsSz81XGC05dvbK0uMkrq6KJCD614sS7kJFRYUsMd2xY0eHP07WHKu33noLzz77bK2EIMYmfP/997KstPlYiSvX4hjq9frqYyXGHZm7idkra763xFVWkSzNV+/P/XfoyMfJGmKcizhGa9eurZ4SQHw/ief3338/HBnzYNMxD1qOudAyzIOWYy50oDxoclDx8fGyitWiRYtM3t7ecl0sBQUFtUpHT5w40bR3715ZdjwkJKTeEuWPPvqo6ciRI6Z3333X4UuUiypzQUFBpquvvloeF1F2+5FHHjHp9Xr53NJj5wweeughWWVOlJM+evSoafbs2bKcdHZ2tnydx6l+sbGxdarM5ebmyrLb06dPl2W3Rflp8W/PmcpuJyYmmrp27WoaP368XBdTApgXMx6ns8Tv7ubmZlq8eLGcQuLOO++UJcprVnR0dsyDTcM8aB3mQusxDzaMudD286DDNh5vvfVW+Q/z3GXdunXV+8TFxcm5dsTcMGLOITGHjCg7XZPYf+DAgXIelc6dO8uS545ux44d8ks+MDDQ5OPjYxoxYoQsL12TJcfO0ZWXl8vfWyRJcZzEHFbiS6wmHifLkqYg5lO74IIL5BehOBF58cUXTc5EfLfU95117jU+Zz9ONYk5Bzt06CC/n0XJ8q1bt6odkk1hHmw65kHLMRdaj3mwYcyFtp8HNeI/rXtvk4iIiIiIiOydXVVbJSIiIiIiInWw8UhERERERESNYuORiIiIiIiIGsXGIxERERERETWKjUciIiIiIiJqFBuPRERERERE1Cg2HomIiIiIiKhRbDwSERERERFRo9h4JCIiIiIioka5NL4LETXFhg0bcNddd8Hd3b3WdqPRiAsvvBDbt29HWVlZnfcVFhbi0KFDeOONN/Dll1/CxaX2P9Py8nI8+eSTGDFiBCZPngxPT886P6NTp074+eefmxT3lClTUFFRgZUrV9Z5bePGjRg7diz27duHBQsWIDY2ts4+xcXFWLFiBbp06dKkzyciIsfAPMg8SI6HjUeiVlJSUoIbbrgBCxcurLU9Li4Oc+fOhUajwd69e+u8b9y4cTCZTMjJycE777wjn9e0ePFiFBQUyMQ2atQo+fxcIqE21ezZs3HNNdcgMTERkZGRtV77/PPPMWTIEPTv3x8pKSn1xj9z5kwZGxEROTfmQSLHw26rRA5s/fr1MjmvWrUKgwYNgoeHBy6++GKkp6fLq6K9evWCr68vbrrpJnmlVLjiiisQEhJSJxmLK8E//PCDTKpERET2gHmQqGWx8UjkBMRVX3H1dvPmzUhISMD1118vuwN98803WL58OVavXo23335b7iu6B82YMUMmTXHl10wkTIPBgBtvvFHF34SIiMh6zINELYONRyIn8Oyzz2L06NHyqqu4YirGobz//vvy+ZgxY3Dttddi3bp11fvPmjULp06dkvvV7KojuvH4+fmp9FsQERE1DfMgUctg45HICYixGWZhYWGyuEDnzp1rbRNdeMx69uwpx5F89tln8vnJkydlkQB21SEiInvEPEjUMth4JHICer2+el2M/aj53LxNVL+rSSTIpUuXyqIE4mqrqBonquMRERHZG+ZBopbBxiMR1UuMB9FqtXI8yBdffCG78IjkSkRE5AyYB4nq4lQdRFQvb29vTJs2DfPmzUN+fr4sPU5EROQsmAeJ6uKdRyJqkOiyI+bZmjRpEiIiItQOh4iIqE0xDxLVxsYjkQMzT7Ts7+9fvU1cOc3Nza1Twry+iY5Hjhwp3y/KmBMREdkb5kGilsXGIxERERERETWKYx6JWomYB2rZsmVyOZfo/iKueg4ZMqTe94oB+pGRkXjkkUfqff2JJ56Ah4cHDh48WO/P6NevH1pbr169GoxfxEZERM6NeZDI8WhM4l48ERERERER0Xmw2yoRERERERE1io1HIiIiIiIiahQbj0RERERERNQoNh6JiIiIiIioUWw8EhERERERUaPYeCQiIiIiIqJGsfFIREREREREjWLjkYiIiIiIiNCY/wfLQxZdJQQ7CgAAAABJRU5ErkJggg=="
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "execution_count": 12
  },
  {
   "metadata": {},
   "cell_type": "markdown",
   "source": [
    "\n",
    "## 结果解读\n",
    "\n",
    "运行代码后，你将看到两张门控特性曲线图。\n",
    "\n",
    "分析两张门控特性曲线图，核心结论如下：\n",
    "\n",
    "低阈值钙通道（：\n",
    "   - `q_inf`在膜电位约 -80 ~ -60 mV 时显著上升，表明此范围内通道易激活。\n",
    "   - 与填充区域吻合，说明该通道在神经元静息电位附近即可被激活。\n",
    "\n",
    "高阈值钙通道）：\n",
    "   - `q_inf` 在膜电位约 -60 ~ -40 mV 时才显著上升，激活阈值更高。\n",
    "   - 填充区域对应其典型激活范围，需更强去极化才能激活。\n",
    "\n",
    "## 扩展练习\n",
    "\n",
    "- 尝试修改通道的电导密度参数，观察门控曲线是否变化。\n",
    "- 对比其他离子通道的门控特性，理解不同通道在动作电位中的协同作用。\n",
    "\n",
    "通过这些工具，你可以快速构建符合电生理特性的通道模型，为复杂神经动力学模拟奠定基础。"
   ],
   "id": "5bb1953fccb96af6"
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 2
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython2",
   "version": "2.7.6"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}
