{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Thalamo-cortical neuron\n", "# Analysis of Ca2+ and voltage-dependent kinetics of the hyperpolarization-activated current" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Imports" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import matplotlib.pyplot as plt\n", "import matplotlib.cm as cm\n", "from matplotlib import ticker\n", "from matplotlib.colors import LogNorm\n", "\n", "from PySONIC.neurons import getPointNeuron" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Functions" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "def plotIhKinetics(Vm, CCa, gatings, ylabel, cmap='viridis', fs=18, lw=2):\n", " \n", " mymap = cm.get_cmap(cmap)\n", " sm = plt.cm.ScalarMappable(cmap=mymap, norm=LogNorm(CCa.min(), CCa.max()))\n", " sm._A = []\n", " fig, ax = plt.subplots(figsize=(5, 3))\n", " for key in ['top', 'right']:\n", " ax.spines[key].set_visible(False)\n", " ax.set_xlabel('$V_m$ (mV)', fontsize=fs)\n", " ax.set_ylabel(ylabel, fontsize=fs)\n", " for c, gating, in zip(CCa, gatings):\n", " ax.plot(Vm, gating, linewidth=lw, c=sm.to_rgba(c))\n", " ax.xaxis.set_major_locator(ticker.MaxNLocator(2))\n", " ax.yaxis.set_major_locator(ticker.MaxNLocator(2))\n", " fig.subplots_adjust(right=0.85)\n", " cbar_ax = fig.add_axes([0.87, 0.1, 0.02, 0.8])\n", " fig.add_axes()\n", " fig.colorbar(sm, cax=cbar_ax)\n", " for item in ax.get_xticklabels() + ax.get_yticklabels() + cbar_ax.get_yticklabels():\n", " item.set_fontsize(fs)\n", " cbar_ax.set_ylabel('$[Ca^{2+}_i]\\ (uM)$', fontsize=fs)\n", " return fig" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Parameters" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "pneuron = getPointNeuron('TC')\n", "Vm = np.linspace(-100, 50, 100) # mV\n", "CCa = np.logspace(np.log10(0.01), np.log10(10.0), 10) # uM\n", "\n", "# rate constants\n", "alpha = pneuron.alphao(Vm)\n", "beta = pneuron.betao(Vm)\n", "\n", "# proportion of regulating factor in unbound state (-)\n", "P0 = pneuron.k2 / (pneuron.k2 + pneuron.k1 * (CCa * 1e-6)**4)\n", "\n", "# Extend to match dimensions (nCa, nV)\n", "alpha = np.tile(alpha, (CCa.size, 1))\n", "beta = np.tile(beta, (CCa.size, 1))\n", "P0 = np.tile(P0, (Vm.size, 1)).T" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Open form" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZ8AAADzCAYAAABZoxsBAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAIABJREFUeJzt3XecVOW5wPHfM7O9sgtL7wiCIgqoWJEoir3FGk001xKvsZeYm6hRE29ijzdeTSA3sWuMvYXYewUUxIIoIIjSWerWOc/945zFdZm2087OzPP1cz6ze857zvvOsu4zbxdVxRhjjMmkgN8FMMYYk38s+BhjjMk4Cz7GGGMyzoKPMcaYjLPgY4wxJuMs+BhjjMk4Cz7GGJOlRORYEXnd73IkwoKPMcZkGREJisjFwL2A+F2eRFjwMcaY7PMH4AjvNStZ8DHGmOxzs6pOBBb6XZBEWfAxxpgso6rf+l2GZFnwMcYYk3EWfIwxxmScBR9jjDEZZ8HHGGMyTESmicgrEa4NEZFHRWSNd9wtInUZLmLaie3nY4wxmSMipwF/BV5V1UkdrnUHZgJFwK1AAXApsAjYVVWbM1rYNCrwuwDGGNOVTflBua5eE4p4feacpo+Bxnanpqrq1I7pRCQI/Bq4Kkp2FwH9gR1U9VPvvneB54FTgGmdLX9XZcHHGGOiWL0mxHv/HhjxerDP/EZV3TnaM0SkBHgXGAPcDewXIekJwCttgQdAVV8QkXnetZwJPtbnY4wxUShKi7ZGPOJUAlQBx6vqKcBWN4pIDTAUt9mto1lA1ACXbazmY4wxUSjgELVvvFpEpgJPqepTEdKsB4arRo1W/bzXpWGufQtUiUi1qq6LVeZsYMHHGGOicGs+kft8gHWqembUZ6g6gBMjq0rvdXOYaw3eazlgwSdbHXjggTp9+nS/i2GM6Tqirgwdo+aTKm3dINEyixXAskZe9vmsWrXK7yIYY7KEAi04EQ+8ZjcROSzJrDZ4r6VhrpV2SJP18rLmY4wx8VIgFH0+ZMxmtzgt9l77hLnWF6hX1U0pyKdLsOBjjDFRKEpLBprdVLVeRBYC48JcHgvMSHshMigvm92MMSZuCqEoB6lrdgN4BJgsIiPbTojIZGBb4MEUPL/LsJpPDM3NLRQVFfpdDGOMTxShJfp4hFQ1uwFcD/wEeFFEbsKdH/QL3Lk/96Yojy7Bgk8Me51+KwClTcqAijJOPGxXDjlyF59LZYzJFAWcDC2BqaorRWQicAtwDe6w68eBS1W1KTOlyAwLPlEsWbySlrIAIsKGcviEZq6Y/jpXPfEalZsdDtt5BOdddCjBYNDvohpj0kSB5hT3UKjq4CjX5gEHpzTDLsj6fKLoP6AHNdtWQe8gVEKoyCFUCK0lAeprC7hnwQJ2O+2PHPOTW/l6yUq/i2uMSRNHJeJBavt88obVfKJoag2xurGBxkDInXtc6cbqQoTAegeaQQmwAIcjrr6Hfg3Cnf9zOrXdK6M/2BiTNRyEZqK2bqSyzydvWPCJoqSwgBm/OJsvVq5m7rcr+PDrb3n3qyUsWbsOqtwOyKADgQ1KQYPwdSFMuXgq+/buxXV/ONnn0htjUsWr4ZgUsuATQ2EwyKjePRnVuyfHjh0NwOI19fz7s/k8+/HnfLJsBaFqobUKgpuUokCA5+tXMPOEG3n0z/9JVbdyn9+BMSYZitCs1q+batbnk4CBtd04Y49deOyMk/jnaSdy6PbbEggIrRVCQx20VAira4Lsf+4dvPnax34X1xiTBHdV60DEwyTGfnJJGtO3NzcdfTBPnfUT9tlmCBqE5mpoqlGaq4Kc//fpvPj8h34X0xiTIFW35hPpwAYcJMSCT4oM61HL1BOP5PbjDqemrJRQqdDQQ2mpDHDZAy/yzFM5tTKGMXnFQSIeeAMOouzlY8Kw4JNi+207jCfOOJldB/VHC4TG7m4A+s1jr/LR7EV+F88Y00lun09BxMMkxoJPGvSqquBvJx3NYaNHQkBorFWaqwKcccPDbN7YEPsBxpguw/p80sN+cmlSGAxy/ZEHcvy4HSAgNNUoDbVBDjv9f/0umjGmE9pGu0Xp8zEJsOCTRgERrj54P47ecTs3ANUqq7sHOe+iO/0umjGmExwNRDywAQcJseCTZuIFoPED+qJBoakW3li/mnmfLvG7aMaYOCgQIhDxwAYcJMSCTwYUFRRw27GH0beqEqcImmsDnHFNTm3NYUzOUoQWDUY8TGIs+GRIbXkZN//QXai2pUJZ1yPIlVf9w+dSGWNiUYWQBiIeJjH2k8ugsf378h+7jQcRmmqEZxctoamp2e9iGWOisJpPeljwybDzJ+1B/6pKtBAaewY556K7/C6SMSaGGH0+JgH2k8uwksICrjvqQABaymFm83o2bbC5P8Z0VVbzSQ8LPj7YeWB/9h4yEALQ1D3A6Rf+3e8iGWMicLfRjjrU2iTAfnI++eWUSaDQWgbzgg2sq9/od5GMMWHEUfOxeT4JsODjk23qunP46JEg0Ng9wIW/us/vIhljIgghEQ9snk9CLPj46KL99kIcCJXC7OYNfhfHGBOGqtDiFEQ8TGIs+PioT3UlBwwfCkBjjyA33WwfnIzpatyFRaNuqWASYMHHZ+dO3hOA1lJ4eNZnPpfGGNORIrQ4wYhHrhKRESJylIj8TETO9L4enqrnW53RZ8N79mBYWRVfbl7Pht5BXn91LnvvM9rvYhlj2smX+TwiMgo4CzgW6NV22ntVL81y4CHgL6r6aaJ55cdPtIv75RH7AdBaIfz3tOd9Lo0xpj1FaNVgxCMXiMgwEXkYmAucBswGrgZ+AhwMHOJ9fY137XRgroj8U0SGJpKn1Xy6gL2GDaIsFGBz0OGbbg6hUIhgMDd+qY3Jdu7abjnft/MJ8BFwKvCoqm6KllhEyoFjgPO8e0s6m6HVfLqAgAhnTZoAQEt1gNv+d7rPJTLGtFGEVicY8cgRx6nqzqp6T6zAA6Cqm1T1LlUdDxyfSIYWfLqI43bdEVQJFcMjMxJuRjXGpEGMeT5ZT1WfyPS9Fny6iJqyUgYEykBgbY+ArfdmTBcRR83HVjhIgAWfLuSSo34AQEuFcPnv/ulzaYwxbWLM88mJFQ5EZK2IvCQiN4jIiSIyIp352YCDLmTfkdsQaFWcQuGtJcv9Lo4xBnfAQS7P52mnGtgHmOR9ryKyEfgQmAXM9F4/VVVNNjMLPjFMeexvOCj9yqsZUFnNsOpahnfrwajaOmpLylKaV1EwyLjansxYv5KNdUFWrVpPjx5VKc3DGNM5iuDk/mg3gFVABfA48AqwPTAeGAvs7aVRoEFEZgMzVfW8RDOz4BNFyHH4ct0aWtVhfv3qra4PrOzG2Lo+7NV3MBP7DaZXWWXSeV527GSO/b8HCJUJ197wOLdc95Okn2mMSZwCrfmxdcI2uHN7zgYmABer6vkiIsAIYBxuMBoH7ATshjvUOiEWfKIIBgLMPPEclm5az9cb1rF4Yz1f1q/m8/rVfLJmBYs31LN4Qz1PLHBHp42r68tRw7bn0CEjqSkpTSjPHfr0ItiihAqFt1dZ05sxXUE+7NujquuBC0XkL8AtwCMi8hJwvqp+AswDHmhLn+jk0jYWfGKoLi6huriE7Wp7fu98q+Pwef0q3lu2hNeWLuTtZUuYtfIbZq38ht+9/zLHDd+BM0fvwoDKbp3KT0QYWdqNj1vXsbEmSFNTM8XFRal8S8aYTlCVLlfzEZEJwJ9xaySzgVNUdX4qnq2qnwEHicihwE3AhyLyZ+A3qrq2XboFyeTTtX6iWaQgEGC72p6cut14/rb/Mcw84efcOvFQ9u47mKZQK/d89gGTHpnGVe+8yLqmxk49+5LjvOV2yoVb/+fZdBTfGNMJjkrEI9NEpAR4DLgBqAGmA3emOh9VfRq33+dXuEvrfC4i/+k1wyXNgk+KlBUWccSw7bhnynH8+8ifcvSw7VHgzk9nst+jf+WpBfFPHJ0wZADSqmgBPP3Jl+krtDEmJgVanUDEwwc/wB3efb+qNgPXAqO9RUFTSlVbVfVGYAdgPnAbcFkqnm3BJw22ranj5omH8PThp7BLr/6satzMua8+xeVvP0dTqDXm/cFAgH7qLpW0tnOtdsaYFGsb7dZVaj7ASGDL/iuqGgIWAikJPiIyREQOF5FficgDIjIX+AJ3gIEAsf+IxcGCTxptV9uThw46kd/uvj9FgSD3fvYhxz57PysbYi6dxLmHuyMbWysCPPSPN9NdVGNMFF1sM7lyoOMSKJuBpOZ+iMjbIrIeN9A8DlyJG+jeB36BW+Oq8WpCSbPgk2Yiwo9HjuWRQ06if0U1c1Yt44R/PcDyzdG3zT5op1HgKE4R3P38jAyV1hjTkWqXa3bbDHQcTlsGbEzyuRNwY8L9wH5AhaqOVdWfquqtqvqqqq5LMo8tLPhkyA49evP4oSczsqaOL9et4fhnH+Cbjesjpi8uKKC6wf3nWVrckqliGmPC6GLNbp/hjnIDQESCwFDaNcUloQz4EfACMF9EHvaa3w4SkV4x7u0UCz4Z1KO0nAcOPIHR3XuxaEM9P33+YTa2NEVMv//IYQA0VQWoX5vshxpjTCIUIeQEIh4+eBnoLiKnikgR8GtgvjdEOhnVuE1rl+DO52kAjgR+BzwNfCMiX4vIkyJylYgcnkxmFnwyrKaklPumHM+w6lrm1a/ioteexYmwTNI5R00CIFQi3Pg/z2SwlMaY9lLd5yMi00TklQjXhojIoyKyxjvuFpG6tuuq2oC7s+jPgdXA/sBxCRWkHVXd4DWt3aKqJ6vqdrgBaW/gQuAeYC1wEG5/0GPJ5GeTTH1QXVzCtP2O5sin7+G5xfO55YM3uHjc3lul61NdSUGT0losvPLV1z6U1BijSkprOCJyGu421K+GudYdt2ZTBFyH+zf6UmCMiOzqDa1GVWcCu6SsUBF4G8u96R1tZSzBXe9tbDLPtpqPT4ZW13LbpMMJiPCn2W/zzrLFYdP18/oV11XlxcKGxnRBqRlqLSJBEbkSmBYl2UVAf2A/Vb1OVa/F3a56R+CUZN5FZ4jIFSJyqIj063hNVRtV9W1VvT2ZPCz4+GhivyGcu+PuAFz2xnQ2tzRvleaMg3YDoLUcZryfktUzjDGdoBCrz6eHiMxod5zZ8RlebWEW7sKd9wBLI2R3AvCKqm6Zla6qL+Cuq3ZCit9aNFcDTwCLRWSFiPxbRK4TkRNEZNtUrHJgwcdnPx+zOyNr6vhqQz03zHp9q+uHT9jBG3It3HLXiz6U0Jg8p27TW6QDWKWqO7c7poZ5SglQBRyvqqcQZqKmiNTgjlqbGeb+WcDOKXtPsQ0Cfgj8NzADdxXrS4H7gE+A9SKS1ARECz4+KwoGuXHvgwmKcOcnM5mx/Pt9O8UFBZQ2uB8y5rXaiDdjMk2BkAYiHnFaDwxX1YeipGlr4gpXK/oWqBKR6vhLnjhVXaKqj6vqFap6sKr2wh3efQ3QBCzADaYJs+DTBYzu3ouzdpiAAr9972U6bhI4plt3ABqq7Z/LmMyL2edTLSJTReSwSE9QVUdVYy1L07Yh2OYw19pWNCjvfPlTQ1W/UNWrcSeg9gQOTOZ59tesizh7zG70KC1n9qpveXbRvO9du/C4HwAQKhWefvp9P4pnTF5zHIl44C7yeaaqPpVkNm1/j6NtUe0kmUfSVPVt4FncfqGEWfDpIsoLi7hwpz0BuH7mazSHQluu7TSkv7vKdRDu/Ne7fhXRmLzk9u1IxCOF2tbcCrcTZWmHNH77EHeuUcJ8Dz6xJlRFua/Om6i1XETWi8grIrJbJsqcLsePGMPQ6lq+2lDP/fM+3HJeRCj3+n0Wauf2BjLGJC/kSMSDOJrd4tQ236JPmGt9gXpv3k3aiciidkvrTBGRHh2S7EiSTYC+Bp92E6p2w51QdRNwOPC8t2xEpPsqgddwZ/XeAVyB21n3koiMTne506UgEOCX4/cB4LY579DY+l0T8Y7d3X/7xirfPy8Yk3di1HxS0uymqvW4WyOMC3N5LO6os0xZAOyLu7TOs8ByEflKRJ4VkdeAnwJvJJOB33/JEp1Q9UtgW+AwVb1KVW/FXQIC3KW/s9b+A7dhVG1PVjVs4vEFH285f+GJkwFoLRGefvI9v4pnTN7J8H4+jwCTRWRk2wkRmYz79+7BVGcWiaruq6q1wDDgeNzKwafAaNzdTacDZyWTh9/Bp9MTqrzJTacAz6jqa+3uW4a7IN7Wk2WyiIhw1uhdAZg69/0t676NHtAbaVUIwt+mW/AxJmNi9/mkqtkN4HpgDfCiiFwkIr8CHsad+3NvCp7fKaq6UFUfVtVfqeqBqjpQVbur6iGqGn5Zljj5FnySmFA1GLeJ7XnvOSIiFQCqeruqRlu6IiscPGRb+pVXsWDdGl5Y/AXg9ftsdj9lfSXW72NMJqkjEQ9SN9oNVV0JTARm486puQB3Y7eDVDXyEvhZyM+aT6ITqoZ7rytE5AagHtggIl+k6JOH7woDQU7b3o29f5n7XS1nbA93HEZjpd8VVmPyS4wVDhJ4ng5W1UkRrs3zJnZWqGpPVT3VC0ppIyL7JXHv5ETu8/OvWKITqrp5r7/FHep3PvAT7zmPR/pBiMiZbWsvrVyZ1n/HlDh+xBiqi0qYuWIpH6z8BoCLT9ofgNZS4dlnwlUYjTGppgrqBCIeOWK6iLzkLSYajJVYRApF5CgReRV3QEKn+fmTS3RCVbH32g3YU1XvVNV7cKuq9cDvwz1IVae2rb1UVxdzJLfvyguLOH7EDgA8MG82ACP79URaFALw9+k238eYTIlR80lln49fxuKuN/ck8K2I3Cci53vBaA8R2VNEDvP6oR7CbZ16GPdD/06JZOjnfj6JTqhqG+f+qKqubTupqvUi8iRwiohUqGrWL4R2wogdmTr3fZ5a+BlX7LovlUXFlDUImwphodMQ+wHGmBTY0rcTyTpV3Wol62yiqnOBA0Rkd+Bs4AjgRLauHAjuOnWPAneoasJLrvgZfBKdUNXWR7QizLUVuD+cCiDrg8/Q6lom9B7Au8uW8MSCTzh55FhGVnRjJvU0VNj+PsZkTIJ9O9nGWzrnba/pbTywHVCH+xNYCcwFPlDVpJf58a3ZLYkJVXNxV1XdPsy1IUAj7g8pJ5w4YgwAD34+B4Dzjp0IQKgUPpj1pW/lMiZvaMzRbjlHVUOq+p7XrXGDqt6oqnep6sxUBB7wf55PpydUebWhJ4FDRWT7dvcNwV0d4QlVDYW7NxsdOGhbqotKmLt6OR+tWsaEUUORkKJB4bb7X/G7eMbkCYly5ESfT8alNPh4c246rgEUTcwJVSIyVEROFpGh7e77Be7ggpdF5Ncicinu5NIG4FepeC9dRUlBAUdv48bYf8yfg4hQ7I0PnLup3seSGZNHnChHCuf55JOUBR8ROQa3JnOJiDwlInvGuifOCVUTcbedndjuvkW468G9iru73uW4q6zuoaoLUvWeuopjh7uj3p5ZOI8WJ8TAgjIANlX4WSpj8oQCKpEPk5BUDjj4L1UdDyAi3XCbxiZGv8WdUAUcHOX6ncCdYc4vAI5NsKxZZVRNHcO7dWd+/WreWLqI06bsxmWvvUxrmbDsm9X07tvd7yIak9MSnUxqIktls1tIRNrGex+NO5jApICIcMTQ7QB4YsGnHLb3GHAULRBumPpvn0tnTB5wJPJhEpLK4HMc8H8i8jpwEnBOCp+d9w4fOgqA5xbPp9kJUeT1+7y7fJmPpTImP4hGPkxiUhl8jgWW4+6t0wQclMJn572Bld0YW9eXza0tvLD4C3o5hQCsK7PffmPSSqPUelK7mVxeSWWfz38Ao1U1JCJzgKeBh1L4/Lx3+NBRfLDyG55Y8CmHjh/FHfM/orVMaGxspqQk4t57xphkRf+Ml/UrHIjIlUk+4m5vIFjcUhl8FgJ7isg7uEszvJPCZxvgkMHb8tv3XuLVpQv4/bFncceNc3AKhTv+9gIXnh1xzIYxJlkpmVbZpV2FG2IT6cRS3F1NF3XmplQGn+OBM3H7fj4gy3cU7Yp6llUwodcA3l62mDe//YqCRmgthWfnzudCvwtnTK5qG2qd+y4EnujkPbWE35MtppQFH1XdANyUqueZ8A4cNIK3ly1m+lfzqW4MsLpUWVnY6nexjMlpkvs1H4BVqvpVZ24QkYTX0Ix7wIG3esFe3r44PxORA0Uk3H47Jo0OGOTupffK1wvYuX9vAJpskVFjTHJ2BxKZt1Hv3dvp1a3jCj4iMh74FHdFgTuA24FngNUi8lcRGdDZjE1i+pRXslNdHxpDrUzYfwQATrHw9DOR1mE1xiRLHIl45AJVfVdVVydwX8i7N9z2N1HFDD4iMgh4HhiBu2PdtcCvgb/jbmHwH8BsETmks5mbxBw0yA0676/5hkCTgsDfn3svxl3GmIRojMOGWicknj6f/8Ld8nqKqj7f/oKICO78nj8Cj4jIkao6PfXFNO1NGTSC3894lReXfEFZo7KxWFhkm8sZkzYx+nyyfqi1H+JpdjsAuKtj4AFQ10PAjsAnwN0iUpPiMpoOBlfVMLKmjg0tzfTs4S4y2lCeG9V/Y7qk6DUfk4B4gk9fIm/sBmxZnfoIoBxbVicjDvSa3vqPrgPczeU+/niJn0UyJieJujWfSIdJTDzBZwNQGiuRqi7B3QDuyGQLZWKbPHAbAOZsWgEtDgSEm+95zudSGZOj8mhLBRHZIxP5xBN8PgPiLcwHwNCYqUzStq/tSZ+ySpZv3kih49b9526wzeWMSYc8q/k8nYlM4gk+jwBHiMiYOJ9XnFyRTDxEZEvtp6zaHTeyyWZdGZMe+dXnk5HqXDzB58/AEuAJERkVI+1EoFMzZE3i2oJPQR93UdHWMmH5cqv9GJNS+dfnk5GQGjP4qGojcAzQDXhfRK4Rkb4d04nI2cBRdH5tIJOg3XoPoLygkG9aNgIOGhRunGYj3Y1JuRyu+YjIQhFZ0HYAVd7XC73v0yKutd1U9QMR2Rd3i4TLgV+KyPvAl0AZMAYYBiwGrktTWU0HxcEC9uk/lGcXzUOCiobg7W+/9btYxuScHN80blK7rwWYDfwg3ZnGvbabqn6AO5/nMtwgsztwMu6W2cNwO6n2UtW1aSiniWDyALfpTSrd79eV5mY7gDG+yuEVDlT1q3bHIiDU/ly68u3Uqtaquhm4AbhBRAYCg3B3uvhUVdekoXwmhh/0H0pAhKYKpWCd0mKbyxmTWpp3Kxx0mQEHYanqYlV9XVXftMDjn5qSUsb37IeD4hQ5aKFw67REFqc1xkSUw30+YdyfiUwSDj6m69hvwDD3iyL3/4Tpn33pY2mMyS1Cfo12U9WfZyIfCz45oK3fxylXFGV1ScjnEhmTQ9RbYifCYRJjwScHDKuuZWBlN5wC0CKludzt9zHGpIgT5TAJseCTA0RkS9ObU+z1+/yfrfNmTKpYzSf1LPjkiLamNy3x+n0++cLP4hiTW7JwwIGIHCsir8eZ7o8i8lMRKehw7Zl0lc+CT47YpVd/KguLcIpBg8oq6/cxJjWybHkdEQmKyMXAvcQYNi0i5wN/wt0w9JfAmyJS2y7J3ukqpwWfHFEUDDKx3xAAQqUOLRXW72NMymRXzecPuPur/SGOtGfj7lJ9GrAdMBN4qV0AStucHws+OeSAgcMBcEoULRBusfk+xqRENtV8gJtVdSKwMI60fVR1NoCqhlT1bOBF4GUR6U4aw6sFnxwyqf9QgiJosaKi/HuezfcxJmnRaj0+1HxE5BgR0TDHnQCq2pkFHleKyJD2J1T1YuBl79jSByQi+6eg+FtY8Mkh1cUlTOg9AMSt/awutX4fY5IldLnRbo/h9tF0PH6WwLNeBE7teFJVLwBeAUranf5d2xciMltE/i4i54nI3iJtq0vGz4JPjtky4bTUoaVcWLtmvc8lMib7daXg4zWPbQxzNCXwuHOIsBOBqp4HDG73/YR2l3+N26y3L+7AhrUiMl9EHoo3Yws+OWZ/b4M5p0RxgnDln57yuUTG5IAcnWSqqs3egtGRri+OcP5pVb1GVY9U1UFAL+DnwIx487bgk2MGVHZjZE0dBECLlXdXLfO7SMZktzQsryMi00TklQjXhojIoyKyxjvuFpG6JN5B2qnqalV9TlWvj/eeTm2pYLLD/gO34bO1KwmVKhsr7POFMclK5ag2ETkNOB14Ncy17rgd/UW4zWEFwKXAGBHZVVXjnj+hqncCdyZYRgEuBg4HioH5wAfALGCWqq5L5Lnt2V+mHLRlyHWpQ2up8vpbn/hcImOyXApGu3mTP68EpkVJdhHQH9hPVa9T1WuBY3A38jwlobIn5lLgeqAvEAJ+hLuX2wvAGhH5sjP9O+FY8MlBo7v3on9FNQRBi+HGf77id5GMyV6xVzjoISIz2h1bbSwnIiW4tYargXuApRFyOwF4RVU/3ZK96gvAPO9apvwUeAvYFrf2A3AW8HugGXcU3JRkMrDgk4NEhIMGjQDAKVW+osHnEhmT5aLXfFap6s7tjqlhnlACVAHHq+opQGvHBCJSAwzFXWWgo1nAzil4J/EaCPxDVUN8V7/7QlUvBw7ALf/IZDKw4JOjDhzsBp9QqUNDJbS2bvW7boyJQxybyVWLyFQROSzKY9YDw1U1WlNVP+81XK3oW6BKRKoTeQ8J2IRbw6HdawmAqr6OOz/oimQysOCTo8bW9aVnaTkUgFMGv7Mh18YkTFQjHsA6VT1TVSP+T6aqjqrG+gTYNlEz3NDntuaL8s6XPiFfAsMAVHWDl/+Adtff47vmuIRY8MlRAREOGrwtAE6Z8q8F8SzzZIzZSuZWtW77exxtGEOmZhY9DxzS7vv3gEPbfd8T6JZMBhZ8ctiBg75reltX1TWX3zUmK0Tv84mn2S0eG7zX0jDXSjukSbdbgGtEpC3fO4BDvDlHvwYuBOYkk4EFnxy2a6/+1BSVQgGEKuHuf2w1rcAYE4cYNZ+YzW5xaltNoE+Ya32BelXdlGQecVHVtar6D1Vt8L5/CLgROAn4Le7w6wuTycOCTw4LBgIcuc12gNv09n9vfOBziYzJQmlY4SBsNqr1uOuljQtzeSydWLomEeI6S0R2i1C+X+COgtsDGKyq7yaTnwWfHHfkUDf4hMoa7yE9AAAVXklEQVQcVpbbKtfGdFYco91S6RFgsohsGcYsIpNx59s8mPLcvu9HwP8Ce0ZKoKpLVfUdVd2YbGYWfHLcmB696V1UDkFoqVXefHee30UyJvuoRj5S1+cD7qoCa4AXReQiEfkV8DDu3J97U/D8aE7Cncx6c7REIvIzEfmNNy8pYRZ8cpyIcOL2OwEQKlOuud92NzWmU2KPdktVnw+quhKYCMwGrgEuAB4HDkpwy4TO2Bn4p6rGakx8HneOzzHJZGbBJw8c4TW9OaXKktIWn0tjTPZJdbObqg5W1UkRrs1T1YNVtUJVe6rqqV5QSrdqYFGsRKq6AHeNt0NjpY3Ggk8eGFxVw7CSbhCA5h7Ko8+873eRjMkqKVjhIBusAbrHmfY1YFQymVnwyRM/HbsLAKEK5abpb/hcGmOyiBKrzydlzW4+mwfsE2fabwk/JDxuFnzyxJHDtiMQcjeYW9U9y7dfNCbDMjjazU+PAweJyISYKSHpze0s+OSJisJiDh3kjt5srVWuve1Jn0tkTHYQMjPPpwv4C/AN8KiIjI2R9gfAF8lkZsEnj5w5zv1AEyp3ePSL+T6XxpgsEa3JLebAsOzhrWZwDO4Cp2+JyA0iMqxjOhE5D3dbhceTyc+CTx4Z3b0X1c2FEID1fZVlK9b6XSRjskKeDDhAVd/D7fdZhLuN9jwR+UBEHhSRR0Tkc9x13xYQYz5QLBZ88sxlEycB0Fqt/Md/3+dvYYzJEjGa3XJlwAEAqvoB7rbdZwMfAWOA44CjcDe7exrYx9tqIWEWfPLMUcNHE2gFLVK+7N7od3GM6foUCGnkIweparOq/llVx+Jun7AzMAGoU9UjVPWbZPOw4JNnSgsKOXnYGMCd83PR79O9XJQx2S9PBhyEpaqrVXWWqr6vqilrqy9I1YNy1QvLPqLVCSEIIiAIARECCCIBAiIIQlCEgPd9gQQokCAFgaD7dSBIUaCAwkABxYFCSoKFFAcKEBFf3tPFe07i7gVz0BLl35u+9qUMxmQTcXI7yojIHsA8VV3dyfuCuDWijzrbDGfBJ4Zr5z7KptbUL6kkCKXBQkoLiqkoKKayoJTqojK6FZZTU1xOXXEVPUuq6F3Sjf5l3akuKktZ3tXFJexR0pe3Wr6hoY/y4BNvc8IRu6fs+cbklO82jctlrwM/Bu7v5H3dvHv3B17qzI0WfGLYr/cONLY2o6j3O6g4qlte274OqYOj7mtIHVqdEK3q0KohWp0QLU6IZqeVJqeFxlALzU4rm0PNbA41s7op9geGqsJShlT0ZHhlb0ZU9mWHbgMZUlFHQBJrOb3l6COY8MAdOKXKte+8acHHmAgEkOh9O9UiMhV4KosHHQjQXUQGdvK+Wu/eTrPgE8Plo49Oy3ND6tDQ2symUBMbWxrZ0NrA+pYG1jZvYk3TRlY1rWd54zq+bahn6ebVrG9pYPbar5i99qstz6gsKGHn7sPYrcdw9qobSV1JVdz59yqrZAdq+UjWsHGActtdL3DOKZPT8VaNyXoSfT7POlU9M1NlSaM/ekdnJVQvtODjk6AEqCgsoaKwhF4l1VHTqiqrmzbw5cblzN+wjE/Wfc2c+sWsaFzHy8s/5uXlHyMIO9YM4oA+Y5jSZ0cqC8NtA/99d570I8bfextaqty+8EPOwYKPMVtRhRzv8wGuTvL+BZ29wYJPFhARepRU0aOkigk9hm85v3TzGt5dNZ+3Vn3OO6vm8+HaRXy4dhG3fvYvDugzhuMH7cGIqshr/3UvKWNy6UBeaF5MQ1+Hq/7nca4678hMvCVjskquj2pT1WSDT6dJ7H2Dcs/OO++sM2akdTv0jNvY2sgbKz7jqaUzeX/1l1vO795jBKcO3YextUPC3tcUamXUtJtxiqFoRYA5559LSUlxpoptTFcRsd+iqqKfTtjp7Ig3vvDm5TNVdee0lCqH2TyfHFFRUMKBfXfif3c5jYf3vogTBu1JSbCQt1d9zs/em8YFM+7k8/VbzwsrDhZw1iB3DcHmHg77X35HpotuTNeXB2u7ZZoFnxw0sLwHF406hCf3uYzTh+1LebCYt1Z9zslv3cbVcx5mVeP676X/xUH7U1UfhAB8PayZB594y6eSG9M1iaMRD5MYCz45rFtRGWcOn8xj+1zCjwbvRaEEeeabWRzz+s3cveA1Wp3QlrQvnnEm0gJaAld89hYtLbbdtjFbRK/55MzCoplkwScPdCsq54KRB/PgXhcwsecoNoeaue3z6Zz81p/4YM1CAOoqKjm1ZjsAWuocdv3NrX4W2ZguQ1SRUOSDHFtYNFMs+OSRAeXduXHcj7l1/Kn0L6tlwcYV/Oy9afxu7qOsa97Mb449lF6ri0Bg7TCHQy+/3e8iG9M1WJ9PylnwyUO7143g/j3P5/Rh+1IoQZ78egbHvXELzyydxZsXnktxvUAA5g7cyCV/fMjv4hrjrzxc1ToTLPjkqZJgIWcOn8x9e57LuJohrG3exNUfPcw5M/7GP398HIEGoBAeKVvEOTc+4HdxjfGVqEY8TGIs+OS5wRU9uWPX0/nNDsdQU1TOrLUL+dnMaZw4fhDBZgcthKerl3DK7+/0u6jG+ETBcSIfJiEWfAwiwiH9xvHQXhdy9IBdUVWeXzubwUNbKCtsggLl1Z4rGH/lTTYKzuQfxfp80sCCj9miuqiMX25/JHfv8XPG1QxhQ6iBih4NdO+xnpLyJlYPaWXUn/7IPU/bPCCTX2KMdjMJ8D34iMgQEXlURNZ4x90iUpeu+0xsI6r6cseup3PTuB8zrKIXwUKlqlsD3Xuup6h/A1eueI0dr76RRUtX+V1UYzLDaj4p5+vCoiLSHXgZKAKu88pzKTBGRHZV1eZU3mfiJyLs3XMUe9RtywvffsRdC1/liw3LqKhqpLyykaaaQg546S8Uf1XEjZOmMGWvMX4X2Zj0UIWQ9e2kmt+rWl8E9Ad2UNVPAUTkXeB54BRgWorv67RvGpagxPeLJ53eU2nr9OGeIAhI29Plu3O4QWLLfyJ4G3xvOR+QAAECBCRIQIIEJUCAYNxbeAclwJS+O3JAnzG8vepzHvrqbd5a+TklpS2UlLbgdN/MJcsf5cK/PUXtN2Vcse9kpuyxQyd/DsZ0cdFrOLmwmVzG+bqqtYh8CSxU1ckdzn8GLFXV/VJ5X5vOrGp9w9zjaNGGuNKmQrR/DQ0TmlS3Pud4D1LvHvdgSwh1NIAQRAgSCBQiUkRBoIgCKaYgUERhoJgi7ygOlFASLKUkUEpJsIzSYBmNIeXDNV/z3KLPaCpUWp0grRoAhNaWAC2NBbCugMr6IibWDOXiow6grjb+je6M8UHET2PVxb11j34nR7xx+sKbbFXrBPhW8xGRGmAo8HCYy7OAQ1J5X6LO6/4hhQltEpt9nC0B67vDARwV99WBUEhwEHaoghPHCCEVWlVoIUCzE6BJAzQ7QRqdIA1OAQ1OkGe/+TMbFxeyobmYDU3FbGgoYfOmEloaKgg0daNboA+Duw1gzID+7LztUCrLY2+EZ0zmKKg1u6Wan81u/bzXpWGufQtUiUi1qq5L0X0JaaUYR1tT8aiYOh/jtq4nSZjrstX1784L0NYCF4hUgO/tpJWemrKjXqBbD5vWecEOCHmBTxFCCg7i9vPy3Xn1zrc9p+3dtdX80O+XWju88r1rsvW1ON5yuFqp8df0z3fh0ilTU/MwG1iQcn4Gn0rvdXOYa23tXOVAxyCS0H0iciZwJsDAgQPjLmR5n4/iTpvNVB0gBIRQpwWlCdVG1GnEoQFHm1BnM6oNOLoZdTai2oDqJlQ34zgbCTkbaHU24DjuObQRoRENNYE2UyAhCkUpECXY9goEcQPflqGXGQh2JvfNqA73+TQBig04SAM/g0/b35pof13C/YsndJ+qTgWmgtvnE08B84lIAPdHW4gES/guxqefqqLaSEtoLas3fsvC5Yv5eu3XbG5ax+aWdYS0CZVGNNCMSCsirQTEIRBwCEgIRAkElAAOAVHEO9z3pQSkrRHxu1qehKsVhtkrOdy5rdIk9e5Numyo3y11D7OaT8r5GXw2eK/hGvhLO6RJxX2mixIRREopDpTSt6YvfWvG+10kY76jCqFQ7HRdhIgcAfw3MAD4HLhAVd/wt1Rb83OS6WLvtU+Ya32BelXdlML7jDEmMVkyyVREhgJ3A2cD3YBbgCdFpNrXgoXhW/BR1XpgITAuzOWxQNix0IneZ4wxifEmmUY6upZBwDRVfVVVHVW9D7cbYqTP5dqK38vrPAJMFpEtPxgRmQxsCzyYhvuMMaZz1B2QE+noSlT1ZVW9pO17EdkNqMBtfutS/A4+1wNrgBdF5CIR+RXu/J2ZwL3gViNF5GSvOhn3fcYYkzJdqOYjIseIiIY57uyQbhvcD+pXqOrajBc0Bl+X11HVlSIyEbdd8hrc4dOPA5eqapOXbCLwd+CnwIJO3GeMMclT7Wr79jxG+OGoW/Y7EZEJwJPA7ap6Q6YK1hm+Lq/jFxFZCXzViVt6ALaEs8k0+73LnFWqemC4C9XBHrpbaeSFU57bdHeXWl5HRA4GHgAuVtW/+l2eSPxeWNQXqtqprRdEZEZX+uUy+cF+77qKrjeqLRIRGQT8AzhVVR/xuzzR+N3nY4wxXZvizvOJdCRARKaJyCsRriWzV9mFuCu83CUiG9sdkxIqaBrlZc3HGGPipYA6qav5iMhpwOnAq2GuJbVXmapeAFyQssKmkQWf+KRodUJjOsV+77oCVTQFKxyISBD4NXBVlGQZ26vMb3k54MAYY+IlItNxB39EUgI0tvt+qreWZPtnlADvAmNwVyDYD/hCVSd1SJfUXmXZxGo+xhgTRaRRcJ1UAlQBx6vqQyKyqGOCTO9V5jcLPsYYk37rgeGqUTcHy+heZX6z0W4diMgwEWmINDpERP5TRD7z0nwkIickk86YjpIc7WS6IG+dtVi7Usa7V1lOsODTjoh0w509XBLh+iXA7cAc3BElS4EHROT4RNIZ01G70U674Y52ugk4HHheRIr8LJtJu0T3OMtK1uzmEZFRwKNEWP3VC0xXAfer6kneuWnAK8ANIvKwqobiTZfed2OyWN6MdjJbyau9yqzmA4jIKcCHQHcg0nIUh+FWee9oO6Hukra3427atEcn0xkTzgnAK22BB0BVXwDmeddM7sqrvcos+LjG4G7FsAPwZoQ0bcuczOpwflaH6/GmM+Z72o12mhnm8izsdyen5dteZRZ8XP+lqqeo6vIoafoBa1W1Y2fgt97rwE6mM6ajuEY7ZbA8JvPyZq+ynO3zEZHeMZJsVNWNALGWrPBUEt8olHjTGdNRvKOdcmKorQnreuAnuHuV3YQ7+OkX5OBeZTkbfPiuphHJtcDlnXhegPCjUNrOOZ1MZ0xHeTXayWwtn/Yqy+Xgc0aM6x908nkbCD8Kpazd9c6kM6ajvBrtlM9UdXCUa/OAgzNXGn/kbPBJwyZKi4FaESnu8Amkr/e6tJPpjOkor0Y7mfxmAw7iNwsQYKcO58d6r+93Mp0x35Nvo51MfrPgE79ncFeuPbfthIgEgLNxt+R+p5PpjAknb0Y7mfyWs81uqaaqq0XkD8BVXjB5CfghsDdwXNuqBfGmMyaCvBntZPKbBZ/OuQbYBPwcOAr4HDhWVTsugR5vOmO+J59GO5n8ZpvJGWOMyTjr8zHGGJNxFnyMMcZknAUfY4wxGWfBxxhjTMZZ8DHGGJNxFnyMMcZknAUfY4wxGWfBxxhjTMZZ8DHGGJNxFnxMXhKRfiKySkSGJPmcl0UkFG3nXBGpEZFmEXnE+/6v3rptxuQtCz7GFyJyjYioiOweJU2BiHwmIhtFpH+Ki/BH4EFVXZjkc+7F/f/o6ChpfggU8t3CoNcAZ4nImCTzNiZrWfAxfpnjvY6OkubnuFsJ/F5Vv05Vxt7CnUcC16XgcQ8DTcAxUdIcD6zF3W4DVV0MPIC7eKgxecmCj/HLR95r2OAjIrXAlbibq6W6iepC4HVVXZLsg1R1HfA0MFFEena87p37AfBPVW1ud+l+YF+r/Zh8ZcHH+OUL3E33ItV8rgJqgYtVtTFVmYrIAOAw3G0K2p9fJCK3icjpIvK5iDSIyPsisquI9BaRh0Rkg4gsFZFrvb2a2twHBHG3z+joGO9ax714XgPWAOek6r0Zk01sSwXjGxGZCfRX1V4dzm8LzAVeVdXJKc7zDGAqsI2qftnu/CLcIBHE7Q8KAL8G6oF1XnnaNgY8ADhVVe/y7i0ClgGzOpZXRF4DBgJDtMP/bCJyHzBRVQek8j0akw2s5mP8NAfoKSJ1Hc7f7L2en4Y898Ld6G9BmGv9gCmqer2q/gG4A+gPzFXVE1R1Ku7AgmbcAASA15z2MDBJRHq0nReRfl5+93UMPJ45QP9kR9wZk40s+Bg/bdXvIyL7AwcDt6vqx2nIcyiwKEIw+FJVP2r3/efe62NtJ1R1E7AC6NPh3ramtyPbnTsOECJvf90WAC34mLxjwcf46XvBR0SCuLWeVbh9PunQHVgf4dryDt+3eq8rOpwPsfX/O68Bi4Fj2507Abcp7tMI+bWVo0eE68bkLAs+xk8daz5nel9frqpr05SnQ+Tf+9YI52N2jHo1qQdwR7DVishgYFci13poV45QrOcbk2ss+BjfqOoyYCUwWkSqgauBD4FpHdOKyGEi8rGIXCEiX4vIWm9k2ngRedMbifaiiJTFyHY5bu0nHe4FCnBH0x2DG1QeiJK+rRwda1zG5LwCvwtg8t5HwHjcOT11wA9V1QmTbhxuf81CYDDwM+AG4HncEWjNuB34hwIPRcnvK2A3EQmqakprHKo6V0TmAIcAvYAXvQAbSduqDYtTWQ5jsoHVfIzfPgKqcUe2Paiqr0dINw6Yqqr3qmorMBsoAc5S1WWqugZYQuwPVC8BZURfWSEZ9wIHAnsSvckNYDfgC2/FA2PyigUf47e2fp8m4BdR0o0Fnmz3/Y7AG17QQUQE2A74JEZ+/8bt99k7odLG9gBQjvt+HouUyJukujvwrzSVw5guzSaZmi7PmzuzEqhtG4ggIn8FVqrqf3nfb4M7EbRSVVtiPO8xoE5V90pvyaOWYX/gOWAnVZ3tVzmM8YvVfEw2GIc7N2dth3MzO3z/UazA47kR2NMLWH75CfC8BR6Tryz4mGwwFpjV9o2IFALb8/3g87000ajqm8BTwGUpLGPcvBUNjsFdvseYvGTNbiYveQuMzgZ2ab/GW4by/juwTlUvyGS+xnQlFnyMMcZknDW7GWOMyTgLPsYYYzLOgo8xxpiMs+BjjDEm4yz4GGOMyTgLPsYYYzLOgo8xxpiM+3/uxZqsLqROkAAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "O = pneuron.k4 / (pneuron.k3 * (1 - P0) + pneuron.k4 * (1 + beta / alpha))\n", "fig = plotIhKinetics(Vm, CCa, O, ylabel='$O_{\\infty}$')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Locked-open form" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAaIAAADzCAYAAADErhKJAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAIABJREFUeJzs3Xd4XNWZwOHfd6erV/fesA3YGBtMD72XNEqAhJBCEtjdFALspofd7CaQhGSzgYQ0AkkooThUB2MwYIrBNjYY3HuRbfU2KjNzv/3jjowQI2mkuSONRud9nvuMNPfec87I4/nmdFFVDMMwDGOwWINdAMMwDGN4M4HIMAzDGFQmEBmGYRiDygQiwzAMY1CZQGQYhmEMKhOIDMMwjEFlApFhGMYQJyKXisjLg12O/jKByDAMY4gSEY+I3Aj8BZDBLk9/mUBkGIYxdP0YuCT+OGSZQGQYhjF0/VxVTwG2D3ZBUmECkWEYxhClqhWDXQY3mEBkGIZhDCoTiAzDMIxBZQKRYRiGMahMIDIMwxgkIvI7EVnWzbnJIvKoiNTEj3tFpHyAizggxOxHZBiGMfBE5PPA74EXVfXULudKgVWAH/gl4AVuAnYAx6pq+4AWNs28g10AwzCMoeKc03K1uiaW8Nyqt9veBVo7PXW3qt7d9ToR8QDfBn7QQ1bfAMYBR6rq+vh9K4AlwDXA7/pT/kxlApFhGEaSqmtivPHPCQnPeUZvblXVBT3dLyJBYAUwB7gXOKObS68AlnUEIQBVfU5ENsbPZVUgMn1EhmEYSVKUiEYTHkkKAgXA5ap6DfChG0WkGJiC0zTX1Wqgx2A3FJkakWEYRpIUsOm2X71QRO4GnlDVJ7q5pgGYrtpj5Bobf9yb4FwFUCAihapan0yZhwITiAzDMJLk1IgS9xEB9ap6XY/3q9qA3Us2+fHHcIJzLfHHXMAEomxx7rnn6uLFiwe7GIZhZI4eV7HuoUbklo4uk54y6i2YDSnDvo+oqqpqsItgGMYQoUAEO+FBvGlORC5KMZvG+GMowblQl2uywrCvERmGYSRLgVj3cy97bZpL0q744+gE58YAdara7EI+GcMEIsMwjCQpSiTNTXOqWici24GjE5yeB6xMawEGwbBvmjMMw0iaQqybA/ea5gAeAc4UkZkdT4jImcBhwAMupJ9RTI3IMAwjSYoQ6X4sg1tNcwC3AZ8BlorIz3DmH92MM7foLy7lkTFMIEpSXf2v8Yfv/MBz+oE3pCR4Xj78s1jxny3AQrFQvIAHxIviB/GDBEFCILmIVYBlFWFZ5Xi8o/F5J+LxTMSygul7wYZhfIgC9gAsz6mqlSJyCnAHcCvOUO5FwE2q2pb+EgwsE4iSVN22lckSSW8m2uWxF21qESFE1CpGrTFYvtkEAsfhD5yIZQXSVkzDGK4UaHexR0NVJ/VwbiNwvmuZZTATiJJUb53Bzw5swHq/ghOv6zh1Hjn0SKffFRGnI84S53cPiiXOo0cUDzYeUXxi40XxWjZ+bHxi44//HLRsgqIELJugOM8F4vf4aAZthtgeiL0BrfdgKzRLETHv4QRzPkEgdDYi/kH5uxlGtrG126a5ZFZWMBIwgShJM4pOpjg0/VBtxY7PJ1MUVLHj09xUbRQbW21sbGxVVG1sYsQ0RkyjxDRG1I4S1QgRu512u52IttEca6XNbqU11kJLtJnmWDPN0QZaYl0nWCteYhR7Whnjb2NSwMNoXztFVhO5NJArMXKpg+gr0PAKkXqLFu8cQvlfxhc4DZEe5+sZhtENG6EdT3en3ewjGlZMIEpSnreAPG/BoOQdsdtpiNRTG6mipq2Syrb9HGjbS0XLbta2VLC2Y9EPRuC3/BxbMIkjc5op1HfItXeSZ9nkx9ZA3ZcJU4SV968Ec69CxAyaNIy+6qFGZPSTCURDgM/yUxoopzRQDnmzPnAuHG1iZ3grmxvfZX3DWva17mJ53SaW10FZYDZnlN/I1EAFsfBfKdYdhKQOmv6TcNOv8BX9FH/wlEF6VYYx9ChCu3ZbIzL6yQSiIS7Hm8esgrnMKpjLxWOvpLJtPytrlvNGzUtUte3nwT1/oshXysfG/SeFoVIq6r5LUWw1IasOrf0Cjb7jyS3+NZYnb7BfimFkPGf1bdOS4DbzF80y5YFRnDf6k3xn9h18ZuK/MDo4nrpINX/a/gvu2/MIoeJfQdnj7NNxKJAbfY3WylOJRTYPdtENI+OpOjWiRAfuTmgdVkwgylIe8TC/5ERunvljLh3/OUKeHNY3rOGnG7/F3rYIY0ctpTLnqzTZFkEaiFVdTHv4ycEutmFkPBtJeBAfrGBGzPWdCURZzhKLk8rO4luzfsbM/Dk0Rxu5c8t/s6zyGUYX3oCW/JVKO4RXYlj1N5pgZBg9cPqIvAkPo/9MIBomCnxFfGnqLZw98qMoyqK99/Hs/scoDM6naMQSKuxSLFGs+huJtDw72MU1jIzU0UeU6DD6z/z1hhFLLC4YczlXT7weQXiq4iGW7F9EwDuCshFPUmEXYYlC3b8SbVs92MU1jIzTMWqumz4io59MIBqGjik5masmfgVBeLLiQV6tWkrAW0pp+eMcsPPwiBKp/TxqZ81OxIbhGluthAdmsEK/mUA0TB1TcjKXT/giAA/vuYedzVsJ+kaRV3ofDbZFgGbC1Z9BNat2JDaMlCgQw0p4YAYr9JsJRMPY8aWncXLZ2cQ0yh+330FTtIH8wOG05N5Mu0Iotp7WhtsHu5iGkTEUIaKehIfRfyYQDXMfHftpJuVMpy5SzV933oWqMrrwc+zznASAN/xH7OiuXlIxjOFBFWJqJTyM/jN/vWHOa3m5dvLXCHlyea9hDW/VvQ7AhLJfcSAWwiNKS80NqA7AJiyGkeFMjSg9TCAyKPKXcPGYTwHw6J4/E44247Vykfz/IKIQsjcSbXlskEtpGJmhhz4io5/MX88A4LjS05iSexiN0Xqe2Hc/ACPzL2cv0wCwG24lCzeGNIw+MTWi9DCByACcOUaXj/8CHvHwavVS9oS3IyKMKPklDbbgI0xb4529J2QYWczZKrzb4dtGP5l1KZK0u7ma16s24bO8+CwPfstLwOMjaPkIevzkeP3keALk+YLkegN4huBeP6NC4zi57ByWVT7NMxUP88WpN5EXmM4O78kU2C8h4T+h+dcjYrYhN4anjhpRN8wOrf1kAlGSNjTs5fb1yb23BCHfF6TIl0tpII/yYAEjg0WMDRUzLqeUKXkjKA3kZ+ROqWeOvJhXq5eyrmE1O5u3MjF3KiOKfkBD1RkUWK20Nf6GYMFXB7uYhjFoYnT7/9bs0NpPJhAlaWxOCZ8Yv5B2O0pUY7TForTbUdrsCC3RdsKxNpqjbTRFW2mOttEQaaEh0sKucFXC9Ap9OcwsGMvc4gkcVTyJucUT8VmD/8+R7yvk5LJzWHrwcZ7Z/zBfnnoLOf5xHPCeSIG9HAn/Ac3/EiLBwS6qYQw4VSFiD/7/02xj/qJJml04jtmF45K6NmrHaIy2UNveTHVbE5WtDexvrWNvuIZd4Sq2Nh6gPhJmRfVmVlQ7+wDlegMcVzads0fP5eTymXitwev8PH3kBSyvepb1DWvY3ryJybkzGFH0feqrzqbQaqW96Q8E8m8YtPIZxmBxFj3NvJaMoc4EojTwWh6K/XkU+/OYkjfyQ+dVlQOt9ayr383btTt5o3oL25oOsnT/OpbuX0eJP5eLxx3DlZNOpMifO+Dlz/MWcEr5uSw5sIilB57gC1NuJNc/kYOeBRTqm9jN96J512dk06JhpJMiROzhNUJORGYAhwMjcGJxJbBOVV3bTdMEokEgIowKFTEqVMSZo44EYG+4hmUH3uXxvavY3nSQe7Yt48Gdr3LZhOO5ZupHyPMObFPYKeXn8PzBJ1hXv4qa9ipK/GUUF95IW+0VBKQWu+01PMETBrRMhpEJhsOcIRGZBXwZuBTo+Dbd8c1T49ccAB4Cfquq61PJL/v/okPE2JwSrpp8Mg+c+FV+v/BLnFh+GC2xdv68/UWuWP5Llh/cMKDlKfAVMbdoIYryStVzABQG5nFQnfdka9MvBrQ8hpEJFCGqnoRHNhCRqSLyMLAO+DywFvgh8BngfOCC+M+3xs99AVgnIn8XkSn9zdfUiDKMiDCneCJ3zL+GdXW7+en6J3ivfg/fWH0vF42dzy2HX4J/gAY1nFJ2DqtrX+W16uc5d9TH8Vl+vDnXoC23EYisRWOViKd8QMpiGJnAWWsuq5uk3wPeAT4LPKqqzT1dLCK5wCeBf4vf26+mG1MjymBHFI3nD8d9ma8edj4By8cTe1dx/Ru/p7qtcUDyn5Q7nbGhiTRHG1lTtwKAkQVXUmn7sERpbbprQMphGJlCEaK2J+GRJS5T1QWqel9vQQhAVZtV9c+qOh+4vL+ZmkCU4TxicdXkk/j9cV9iZLCQt+t2ce1rd7EnXJ32vEWEk8vOBuDlSmf7cK+VQ4vvI875lkVmMVRj2IkhCY9soKr/GIx7TSAaIg4rGMM9x1/PEYXj2d9axw1v/IEDLXVpz3d+yYmEPDnsDG/hQOteAEoLbqBVwU8T2r4q7WUwjEzRS43I7NDaTyYQDSGlgXx+dcznOKJwPBWtddzw5h/T3kzntwLMKTwWgFW1rwJQEDycai0BoKX5D2nN3zAyjY0kPMiSHVpFpFZEnheR20XkU/Hh22llAtEQk+sN8IsFn2V6/mh2hav497f+RtSOpTXP+SUnArCyZvmhpjgNnAeAt305qtG05m8YmUIVIrYn4ZFFCoGPADcCfwXWi0i9iLwoIneIyNUiMltcnEhoAtEQVOAL8b8LPsuIQAFr63Zy5+Zn05rf9LzZFHiLqG4/yM7wVgBK86+h2QYfbdhtL6c1f8PIFIpga+Iji1QBbcADOHOJ/g9nJN084KvAvfHfG0XkFRH531QzNIFoiCoN5PNfR12BRyz+sv1lXjqY0nyyHllicXSxM3l1Ve1yAHL9k6ihY07RPWnL2zAyiQJRtRIeWWQa8FvgE8AtwAuqehJOTWkWcBVwB/BG/PeU1/vKqr/ecHNU8SSun+GMarv1nYepaWtKW14L4s1zq2tfI6ZOU6CELgHAH33TbJpnDBvZvh+Rqjao6teBucAm4BERWQLMUtWNqnq/qn5TVU9X1RJgeqp5Zs9fb5i6etLJLCydRkOkhV9tXJy2fMaFJjMiMIamaAObGtcBUJ5/FfW24CFKrHVZ2vI2jEyhKhlVIxKRhSLylog0i8irIpJyUOigqhtU9TzgEmACsEZE/ldEirtcty3VvEwgGuJEhJtnO6stPLVvNatrtqctn3nFxwHwTv1KAELe0TQwGoC28INpydcwMk2m9BGJsxfLY8DtQDGwGLjH7XxU9UmcRU+/hbO8zyYR+YoZrGB8wPjcUq6Z4kwy/cl7/yBip2cU25GFCwBYV7/q0Og5CZ4JgDeyykxuNbKeAlHbSngMgtNwhoz/TVXbgR8BR8QXLHWVqkZV9afAkcBmnAEMt7iVvglEWeIzk09hfE4p25sO8siuFWnJY1xoEoW+EuojtexpcWpexbmX0WKDjxY08k5a8jWMTJFho+ZmAodWQ1bVGLAdZwBBykRksohcLCLfEpH7RWQdsAU4Dmclbte+8ZpAlCUCHh//dpgzt+fP21+iNRZxPQ8R4YjCowF4p95ZUSHPP50azQegrdk0zxnZr4cJrQMtF2jp8lwYyEklURF5TUQacILOIuB7OEHvTeBmnJpYcbyG5AoTiLLIKSNmMbNgDNVtjTy6Oz21os7Nc+AEp5h/IQDa/kJa8jSMTKGaUU1zYSDU5bkcINXhswtxYsPfgDOAPFWdp6rXquovVfVFVa1PMY8PMIEoi4gIX5x2BgD3bnuJ1li763lMz5tNwAqyt2UnNe2VAOTmXEZUIahVaKzC9TwNI5NkUNPcBuDQ8jsi4gGm0Km5LgU5wJXAc8BmEXk43kR3noh8eNvpFJlAlGVOKp/JrIKx1LQ38ciuN1xP32v5mFUwF3i/VlSScwLVtrNHUnv4MdfzNIxMoQgx20p4DIIXgFIR+ayI+IFvA5tVNdVAVIjT/PZN4H6c5r+PAv8FPAnsE5E9IvK4iPxARC5OMT8TiLJN51rR33YsT8s6dEcUzgfe7yfyiJ9Wz0wAoi1Pu56fYWQSN/uIROR3IrKsm3OTReRREamJH/eKyKGdKFW1BWfH1BuAauAs4LJ+FaQTVW2MN7/doapXq+psnOB0MvB14D6gFjgPp/8o5W+fZofWLHRi+WFMzi1ne3MlLxx4l7NGz3E1/VkFcxGEbU0babfb8FsB/KGLoGUdfnsrqhFEfK7maRiZQBXXaj8i8nmcrbZfTHCuFKfG4wd+gvNZfRMwR0SOjQ/XRlVXAce4UqAexDfJeyV+dJQxiLP+3LxU0zc1oiwkIlw68XgA/r7rddfTz/MWMC40iahG2NrktAKU5J5Lkw0eYmaPIiOLpT58W0Q8IvI94Hc9XPYNYBxwhqr+RFV/hLMl91zgmlReQV+IyHdF5EIRGdv1nKq2quprqnpnqvmYQJSlzhszj1xPgDW1O9jU4P4AgsMKjgRgY6MzdyjkHU09RQC0hvu9UaNhZDSFnvqIykRkZafjuq73x2sRq4Ef4jRx7e0mqyuAZap6aDVjVX0O2Bg/N1B+CPwD2CUiB0XknyLyExG5QkQOc2t1BROIslSuN8AFY505P3/f9Zrr6c/Md5r7NjS8fei5mNfJj/ZXXc/PMDKCOs1ziQ6gSlUXdDruTpBCECgALlfVa0gwKTS+ltsUIFHTwmpggWuvp3cTcVbh/m9gJXAUThPhX4H3gAYReaX725NjAlEWu3SCszbc4n1raYh0nfeWmsm5M/BbASpad1MfqQUglHMRtkLArkDtBlfzM4xMoEBMrYRHkhqA6ar6UA/XdDSDJaotVQAFIlKYfKn7T1V3q+oiVf2uqp6vqiNxhozfirNn0TacwJoSE4iy2MS8co4tnUabHWFJxdu939AHXsvHtDxnJZGNDU7zXEnoZOpsQQSzGreRpXrsIyoUkbtF5KLu7lZVW3vf0jg//hhOcK7jG2Vu38vuDlXdoqo/xJnsOgI4N9U0TSDKcuePcQa0PLPvLdfTPqyjea7RCXI+TwHNMgaAtpYnXM/PMDKBbUvCA2cB0utUNdU3f8fnck+rCNsp5pEyVX0NeBqnHyklJhBluVNHzibo8fF23S72hKtdTXtmgROINja+g63x/xcBZydXT2S1q3kZRiZw+oMk4eGixvhj1+V7Oj/XmODcYFiDM5cpJSYQZbkcb4DTRh4OwDP71ria9sjAGIp8JTRFG6ho2Q1AXugiIgp+Gs1yP0ZWitmS8CCJprkk7Yo/jk5wbgxQF5/Xk3YisqPT8j7niEhZl0vm4kIzYUYFot5mEvdwX3l8hvIBEWkQkWUictxAlHkoOK9T85ybewaJCNPznCC3uek9AIpD86mLL/cTaX3etbwMI1P0UCNypWlOVetwtnM4OsHpeTij1wbKNuB0nOV9ngYOiMhOEXlaRF4CrgWWp5pJxgSiTjOJj8OZSfwz4GJgSXwdpe7uywdewlna4i7guzijTp4XkSPSXe6h4JjSqZQF8tkTruGdut2upj0tfzYAW+KByBIfrdZ4ACKtz7qal2EMtgHcj+gR4EwRmdnxhIicCRwGPOB2Zt1R1dNVtQSYClyO89m8HjgCZ9fWxcCXU80nk5b46ZhJfGTHJC4RWQEswZlJ3N0s5H/H+cc5VVVfit/3IE4kvxlna9thzSMW544+ir/seJl/VqxhTvEE19KenucEoq1N67HVxhIL8Z8Ake14Iutcy8cwMkK8j6gbhSJyN/CECwMWbsP57FoqIj/DmX90M87cor+kmHafqep2nFraw+lIP2NqRPRjJnF8Vu81wFMdQSh+336clWNfTmuJh5CO9eZeOPDu+wMLXFDiL6fYV0Y41vx+P1HOeZ36ifa7lpdhZAK1JeGBe6PmUNVK4BRgLc6cna/hbFJ3nqq2pZp+psmIQJTCTOJJOM1wS+LpiIjkAajqnara01pOw8rMgjGMChZR1dbIuro9rqUrIkzLd+YTbW56F4DC4FHU2R4Aoq1mszwju/SwskI/0tJJqnpqN+c2xieR5qnqCFX9bDxApY2InJHCvWf2915XA1E8EHQdVZGM/s4knh5/PCgitwN1QKOIbHFh5EpWERFOjY+eW3bgXVfT7hiwsKXJqcx6xH+on6i99Z+u5mUYg0kV1LYSHllisYg8H1/o1NPbxSLiE5GPiciLOIMZ+sW1v56IfBKng+2bIvKEiJzYh9v7O5O4KP74nzhj2b+K064aBhZ1F6FF5LqOhQkrK9P6BSOjnD7KCRgvHHjX1dFz07r0EwGI3xm0aPqJjGzTQ43IreHbg2kezvp3jwMVIvJXEflqPDCdICInishFIvINEXkIp6LwMM5n7lH9zdTNwQr/oarzAUSkCOeFnJLkvf2dSRyIPxYBM1S1Np7/E8BW4H9wtrr9gPhihHcDLFiwwL1P5Ax3ZNEESvx57G2pYXNjBTMKxriSbmmgnBJ/GTXtVexr2cW4nEnkhs4j2v4AfmlAY/sRzyhX8jKMwXWoPyiRelX90IrbQ4mqrgPOFpHjgeuBS4BP8eHPZsFZN+9R4C5VfTOVfN2sT8ZEpCMifhxnhEWy+juTuGNS16MdQQgOjcN/HJjf0WdkOKPnPjLSqb284HLz3LQu84mKgvOojfcTxVo/tO+XYQxd2s2RReL7DH0aKMaZUvN5nBHKtwCfw9mMr0RVP5dqEAJ3A9FlwB9E5GXgKuBf+nBvf2cSd/QpHUxw7iBO1DaBqJOOVRaedz0QOQMWtnb0E1lBWi2n66+tdYmreRnGoNEeR81lHVWNqeobqnqPqt6uqj9V1T+r6ipV94bfuhmILgUO4EwobcPZzzwpKcwkXhfP6/AE5yYDrcDw6QRKwoKSKeR5g2xvOsjecI1r6U7Nc+bdbWvacKj/SfzODsZW9B3X8jGMwSfdHFnRRzQo3AxEnwMuUtVlwNU44977os8zieO1pMeBC0Xk8E73TcZZleEfqhrrYzmymtfycFyZM9hweeUG19It9Y+g0FdMc6yJA61ORTUUPIeYQkBrUbvOtbwMY1DZ3RwuziMabtwMRNuBE+PL8VwCvN7H+28DanBmEn9DRL6FMxrj0ExiEZkiIleLyJRO992MM2z7BRH5tojchDORtQX4VkqvKEudWH4YAK9UbnQtTRFhSq7zHWJrsxPgikILqI83WcTa3N8l1jAGnAIqiQ+j39wMRJfjdGD9PJ7uzX25OcmZxKfg7PN+Sqf7duB0pr2Is4Xtd3CWJj9BVbf1/+VkrxPKZyAIq6q3EY66N0m7o3lua5MTiHxWHmEZAUB7i1l3zsgObk5oNRyuDd9W1UachUpTSWMjcH4P5+8B7knw/DacPiojCcX+PA4vHMe6+t28Wb310Ei6VHUNRADqmwOxJWD2JzKyRZYOTBhMWTMd2OibE0e43zw3KjiOHE8udZFqatqdMSKB0Nmogt/ej2pLLykYRuYTTXwY/WcC0TB1UrlTe3m1cqNrqyxYYjE51wlwHbWi4uCJNKpgiaLt7m9XbhgDSsWpESU6zKi5fsukbSCMATQjfzTlgQIOtjW4usrC1LzDeLdhNVubNnBMyckEvGUc1AIKqKe1ZTE58a3EDWPI6v5725BfWUFEvpdiEvfG++37JG2BSESeB14Bft551QMjM4gIJ5TP4B97VvJq5SYXA5EzsXVbp36imG822K+h7StcycMwBpV7u6hkoh/ghNr+dIQpzm6tO/p6Yzqb5l7EWRVhWRrzMFLQMZ/o9arNrqU5LjQZn/g50LaPpkgDAL7A6QD47V2oRl3LyzAG3PAYvv11nAUB+nLMp3/BC0hjjUhVf5iutA13HFM6DQthbd1OmqNt5HoDvd/UC6/lZVLuNDY3vcfW5o3MLTqGopxTaW7+EblWDI1sQPxmB3dj6JLsrhEBVKnqzr7cICJNqWSYthpRfG+iw9KVvpG6Al+Iw4vGE1ObVTXuTbma0mm5H4Ac7wQa1Fm7tt2sO2cYmex4oD+biNXF7+3XAqjpbJrzA++lMX3DBccfap7b5FqaU+Ij57Y1O0PDRYSIZxoAsbblruVjGINBbEl4ZANVXaGq1f24Lxa/N9EuCb1K9/Dt7PjXyWLHlc0A3O0nmpQ7HQuLPeHttMVaAfAETgLAG9vs6qZ8hjGgutsCIns2xhsUKQUiETm1my28O5hPnAw3q3AsBb4Qe8I17G7u8xehhIKeEGNDE7Gx2dHsBLiC0Fm0KfhohdiuXlIwjMwlduIDs+hpv6VaI3oGqBGRLSLykIj8u4icJSIlbhTOSD+PWBxb6jSbudo819FPFG+eyw/MpM72ARBpXeZaPoYx4IbBxngDLdVAlI+zh9CPcPYiugh4DGcPIPfaeoy0Sscw7q7rzlnipd0aD0Ck7XnX8jGMgSTaY43I6KeUhm+rMylkbfz4E4CIWMBMnHHl81ItoJF+C+OBaFXNNqJ2DK/lSTnNjgELO8NbiGkUj3gR/0KIbMOKuLs7rGEMqOyaM9QtETlBVV8diLySrhHFh2OfJCLXiciXRORcEcntep2q2qr6nqrep6rfcLe4RjqMDBYyObeccKydd+rc6b/J9xUyIjCadruN3eEdAOSEziaqEKABjbnTH2UYA20Y1YieHKiMkgpEIjIfWI+zWsJdwJ3AU0C1iPxeRManr4jGQOioFa2o3uJaml3nExUF51FvO2+5qBnGbQxVw6ePaMCqfr0GIhGZCCwBZgBP4/QHfRunKe4gzhbha0XkgjSW00izjkD0RpV7gWhafN25rU3rAfBaOYRlFADtrWajPGMIGl59RAMWXpPpI/oPnEEJ56jqB6bFi4jgbEj3C+AREfmoqi52v5hGuh1dPBmveHivfg/17WEK/Tkpp3lowELzRmy1scQC/9EQ3YdEzJYQxhCVnbUfRGQ7H3x1BSKyDadmpKo6JV15J9M0dzbw565BCJySqepDwFycVRTuFZFil8toDICQ18/c4gnYKCtrtrqSZom/nGJfGS2xZipadgMQDJ2DrRCwq1C7X5OwDWNQZfH9GJynAAAgAElEQVTGeKcCp8WP04Hm+M8dz6dNMoFoDLCypwtUtRK4BMgF/sWFchmDYGFpvJ8oHc1zzU7zXEnweOptQQTsttddy8cwBkyWrqygqjs7HTuAWOfn0pl3MoGoEQj1dpGq7gYeAD6aaqGMwXFsmTOxdUWVe8vwdJ1P5PMU0CxlALS1PuNKHoYxYHruI8q2lRUyZ7ACsAFIdlvNt4C0tSMa6XVYwRgKfTlUtNaxK+zO8OqOjfK2NK0/FNzUd5Rzsn2VK3kYxoAaPqPm/jZQGSUTiB4BLhGROUmml/qmNsagcJb7mQo4tSI3lAdGUeAtoinawMG2fYDTT6QKfns/qi2u5GMYA0EYPqPmVPWGgcormUD0G2A38A8RmdXLtacAaW1LNNKrYzXu11xad05EPtQ8VxQ6iQYVLFG0vcfuR8PILN0MVMiSwQqDptdApKqtwCeBIuBNEblVRMZ0vU5Ergc+BvzD9VIaA2ZhvJ9oVc022m13tvXu3DwHEPCU0EQRAK3hp13JwzAGjN3NYfRbUisrqOpbOMP5KoDvADtE5BURuVdEHhaRTcCvgF3AT9JWWiPtRgQLmZo3ktZYhLW17lRuO0bObW5891A/Ucx7JADavsKVPAxjoJgakfuSXmsuHozmArfgBJzjgauBjwNTcdYlOklVa9NQTmMAHe/yZnmjguPI9xbSEK071E8UCJ2LKgTsvagddiUfwxgQQ2ywgohcKiIvJ3ndL0TkWhHxdjn3VPpK2MdtIFQ1rKq3q+o0YBLwEeBkoFxVL1HVvWkoozHA3N4WQkSYnn84AJsanZW3S3JOO9RPZJtakTFUDKElfkTEIyI3An+hl6HYIvJVnFatfODfgVe67Ct3ctoKSgr7EanqLlV9WVVfUdUaNwtlDK65xRMJWD42N1ZQ1ebO6gcz8joC0ToAAp5SGnHe520t2TLtwhgWhk6N6Mc4Cw38OIlrr8dZxu3zwGxgFfB8p2CU1jlFqW6MZ2ShgMfH/JLJgHvDuGfkHwHAlqb3sNX5+qi++c7J9jdcycMwBsJQqREBP1fVU4DtSVw7WlXXAqhqTFWvB5YCL4hIKWkOtSYQGQl1NM+9VunOMO7SwAhK/OWEY83sbdkBQDDnovi6cwdRu8GVfAwjrbqrDQ1CjUhEPikimuC4B0BVK/qQXKWITO78hKreCLwQPw71GYnIWS4U/wNMIDISOqHc2WH1tapNRO2YK2l21IoO9ROFTqTOthCBWNtLruRhGOkkZNSoucdw+nS6Hl/qR1pLgc92fVJVvwYsA4Kdnv6vjh9EZK2I/ElE/k1EThaR/H7kbQKRkdiE3DLG55TSGG11bdfWrv1EPiuPsOXsT9TWktZBOYbhmkwJRPEmtKYER1s/kvsXupl6o6r/hjM4reP3hZ1Ofxun6e90nEERtSKyWUQe6kvmJhAZ3Tqp3FkR4ZXKja6k1zFyblvzRqLxybLqc97TVvtqV/IwjLTLwgmtqtquqt3Oo1DVhN9GVfVJVb1VVT+qqhOBkcAN9LJjQ1cmEBndOnGE0zznViAq8BUxKjiOdruN7c1Omrk5FxFVCFCLxva7ko9hpI3LS/yIyO9EZFk35yaLyKMiUhM/7hWR8hRKn3aqWq2qz6rqbX25zwQio1vziieR6wmwtekA+1vqXElzdoGz8vZ7DWsAKA4eQ43t9INGWsxyP0bmc2vUnIh8HvhCN+dKcQYJHIfTZPYz4GJgiYj4+5KPqt6jqif1vYTOLtwi8k0ReUlEVojIX0TkRhE5TUQK+5NmIiYQGd3yWd5DexQtr9zgSpqzCuYCsL5hLQAeK0CLx8kj0vKkK3kYRlqlOGouPtH0e8DverjsG8A44AxV/Ymq/ghnzc+5wDX9K3i/3ATchrNBagy4ErgdeA6oEZGtfe0PSsQEIqNHJ5a72zw3JXcmAStIRetuatudPY98wfMA8Mc2oBpxJR/DSIueV1YoE5GVnY7rut4uIkFgNfBD4D6gu9VorgCWqer6Q1mrPgdsjJ8bKNcCrwKH4dTIAL4M/A/QjjOa7pxUMzGByOhRxzDuldVbaYm2p5ye1/IeGsa9Pt48V5J7MU02eIiabSGMzNd9jahKVRd0Ou5OcHcQKAAuV9VrgA8tcS8ixTgbjCbaOXI1sMCNl5GkCcCDqhrj/XrfFlX9DnA2TvlnppqJCURGj8oC+RxROJ42O+raHkVd+4lyfGOpjy/30xJ+2JU8DCMdetkYr1BE7haRi3pIogGYrqo9NWeNjT8mqi1VAAVu9s/0ohmn5kOnxyCAqr6MM//ou6lmYgKR0avTRjnDrl848K4r6c2KB6JNjesODeO2/ccBIO2vupKHYaSLqCY8gHpVvU5Vu108UVVtVe1to6+OSaGJhlN3bGmc2/eS98tWnN0VUNXGeP7jO51/g/eb7PrNBCKjV6eNdALR8soNrmyWV+wvZXRwPG12K9uanUEQeTmXElUIajUaO5ByHoaRFgOz+nbH53JPQyAGaubSEuCCTr+/AVzY6fcREN/lMgUmEBm9GpdTyvT80TRH21hZvdWVNDtGz71b/xYAxaFjqbV9ALSHzSa/Rgbrvo8omaa5ZHQseR9KcC7U5Zp0uwO4VUQ68r0LuCA+p+nbwNeBt1PNxAQiIymnjZwNuNc8d2Sh09/6dv0bqCqW+Ah7nD7PaMtjruRhGOnQQ42o16a5JHWsYjA6wbkxQJ2qNqeYR1JUtVZVH1TVlvjvDwE/Ba4C/hNnSPfXU83HBCIjKR3Ncy8eeM+VRVAn5U6nwFtETXsVu1u2ARDIvdRZjTu2DbXrU87DMFzn8soKCbNQrcNZv+3oBKfn0cflc/oqPon1yyJyXDfluxlnNN0JwCRVTXlnSxOIjKRMyRvJ+JxS6iJh1tTuSDk9SyzmFh0LwJpaZz+i8twLqbUtLFEi4UUp52EYbutl1JybHgHOFJFDQ6NF5Eyc+TwPuJ7bB10J/Bo4sbsLVHWvqr6uqk1uZGgCkZEUEeGMUUcC8M+Kta6kObfIWfB0bbx5zmfl0eyZAUCkxQzjNjKUauLDvT4icFYzqAGWisg3RORbwMM4c4v+4kL6PbkKZ+Lsz3u6SES+JCLfj897Som390sMgGd3bubrLz1J0Osj6PUS8vjI9fnI8fop8AcoCAQp8gcpCeZQGsphZE4eI0N5jMsvIM8XGOziu+K8MUdxz7ZlLN2/jm/OuoiAx5dSelPzZpLnLaCqbT/7WnYxNmcigdBlaPhWAtFNqN2EWHkuld4wXKA91n7qVfVDqyn0KxvVShE5hfhgAZyh3IuAm/q5zUNfLADuUtXeGhyXAJuAffS8XFGvTCBKUjjaTnM0QnO070vQlAZzmFxQzIziMmYVj2Bu+ShmlYzAZ3nSUNL0mZw3gpkFY9nQsJeXKzdwZryG1F+WWMwpPIZXq5eytv4NxuZMpCz/Euqa/otij02k5Qn8uZ9yqfSG4Q43m+FUdVIP5zYC57uXW9IKgR29XaSq20TkOZzh3CYQDYQLJ8/i9PHTaI1GaItFCUcjhKMRmiPtNLa3Ud/eSm1rK7VtYapawhwIN1ERbmRPUz3VrWGqW8OsPPj+ROmAx8uCEWP5yNjJnDZ+CtOLygbx1SXvvDFHsaFhL8/sfSvlQAQwt+hYXq1eypra1zlv1CfxWfk0eaZRzCbaww+aQGRknB4CUaGI3A084cLIucFUA5Qmee1LJNjZta9MIEqS17KcJjh/35rZbFUOhJvYWl/NhtpK3qs+yJrKfWxrqOWVip28UrGT/165jJnF5Vw8ZRafmHY4I3P6tdvugDh79Bz+d+MzvFq1idr2Jor9qTWdTc+fTZ63gANt+9gV3srE3Gn4Qpeh4f8iGF2P2jWIVeJS6Q0jRUpHf1AirjXNDbKNwEdwhmn3poLEw8z7xASiNLNEGJ2bz+jcfE4aM+nQ8zWtYZbv28mLe7ezdPcWNtRWsmFVJXe8tZwLJs3kC4cv4IiyUYNX8G6UBvJZWDqNV6s2saTiHS6beHxK6XnEy4Lik1hW+TSvVy9jYu40yvM/TnXTjynzRGlruo9gwVddKr1hpC4NI+QyzSLgpyKyMImh2a5s1GdGzQ2SkmAOF0+Zxc9OPp83Lr+B35/xcc6dOIOYKou2vceFT9zLV55fxJa66sEu6oecN2YeAE/uXUXv/Zm9W1h6KgCra1+l3W7DZ+UR9joTXu2WB1NO3zDcIqR/HlEG+C3OAIRHRWReL9eeBmxJNUMTiDKA3+PhzAnT+M3pH+WlT17HFw5fQMDj5Zmdmzh70R+5dcVSmiOpb8HgllNHzqbQl8OGhn2sq9+dcnpjQuOZmDOVVruFtXXOnKKCghtoVwhqFXZ7yiuIGIY7uhu67cIXskwRX0XhkziLr74qIreLyNSu14nIv+FsBZHypD8TiDLMuLxCvnPs6bz4iS9y5WHOemx/fG8VZz32B17Ys22QS+cIeHxcMs6psfx95+uupNlRK1pRvQyA4uCxVKmzlmK46U5X8jAMN6S4DcSQoKpv4PQT7QBuBDaKyFsi8oCIPCIim3CGlm+jl/lGyTCBKEONys3nv084h39c+GmOKB3JvuZGrl3yMLeuWEpbLPUVsFP18fELsRCe2/8O1W2pr794dPEJ+MTH5qb3qGo7gIhA8KMA+NtfJv1TJwwjOT00zbm11lxGUNW3cLYmvx54B5gDXAZ8DGfjvieBj8S3h0iJCUQZ7siyUSy68NPcMv8jeMXij++t4hNP/ZU9TYO7FtuYnGJOGjGTqMZYtOfNlNMLeXI4qshZ2urlymcBKC+8jgZb8BKhvfm+lPMwjJQpENPERxZS1XZV/Y2qzsPZ8mEBsBAoV9VLVHWfG/mYQDQEeC2Lr8xZyN8vuJLxeYWsqz7AJU/cx5sH9gxquS6b4IyYe2zXG64shHrqiPMAeK36ecLRZgKeMmo9ziZ6dtPdqGb/cCUj8w2DwQoJqWq1qq5W1TdVtdbNtE0gGkLmlY/hyYuv4eQxk6huDXPl4gd4ePO6QSvPMaVTmZxbzsG2Bp7Ztybl9MblTGZG3hG02a28UvUcACWF36bFhgB1RFueSjkPw0iV2JrwyAYicoKIJDuZtfN9nvi9/ZoEaQLREFMYCPKnsz7JtbPnE7Ftvrn8ae58+3VXhlH3lYjw2amnAvCHrc+7Uis6faSz+eNLlYuJ2hEKg3OokikAtDfekXL6hpGS7jbFy444BPAycE4/7iuK33tMfzI1gWgI8loW3194Bt9feAYC3LbqJX644nnsQQhGZ4+ey8TcMva11PLUvrdSTm9m/hzGBCfQEK1jZe1yAHILbyGiENI9xNreSDkPw+gvASSmCQ+yY9ScAKUiMqEvBzA+fm+/mEA0hF07ez6/OvVi/JaHe9av4ublzxC1B7YfxSMWn596OgB/3PpCyrUiETlUK3p2/yKidoTS0Kkc1BEAtNZ/L7UCG0aKRDXhQfaMmvsFzsZ8fTlWkUK90CzxM8RdOHkmRYEgX1z6GA9vWUdLNMIdp1yI3zNwK3ufNXoOf9r6AtubK3l8z0o+PmFhSukdXXwCzx34B/tb9/JS5T85feSF+PNvItJ0EyF7G5GWxfhC57pUesPoA1XIkv6gbvwwxfv7NdnRBKIscNKYSdx39qVcu+RhntqxkZZohDtP+yhB78D883rE4ovTz+Rba+7nrs1LOH3UkRT5c1JIz8PFY67i7m238eyBRSws/Qgj8i5mZ9P/MYGdROq/jzd4FiJDaxsNIztk8wg5VU01EPWLaZrLEgtGjuP+c6+gOBDi+T3b+MLSRwgP4LJAZ4w8gmNKp1IfCfN/mxannN7sgqM4LP9IWmLNLN7/KCJCcfFPabEhSC2tjb91odSG0UfaYx+R0U8mEGWRI8pG8cB5V1AWymX5vp1cs+RhGtoHZkUCEeHmWRfjEw+P71nJ2tqdKad3ydirEITllUvYE95BYXAuVR5n7pIV/jV2LPMWhDWGgSxfa24wmECUZQ4rLufB865gVE4ebx7Yw5WLH6CmNTwgeU/MK+fTU04B4H/efYzWWN93s+1sbGgiJ5WdhU2M+3b+mojdzqiSn1IT8+IjQkvN5wdl2LoxvGXzPKLBYgJRFppaWMrfz7+KiflFrKs+wGVP38++poYByfuzU05lfE4p25oO8vP1T6ac3kVjPsWIwGj2t+7hyX0PEvCVE8n/ljOcO/YerU0p7VBsGH3XfY0oG4ZvDwoTiLLU+PxC/n7+lcwoKmNLfTUff+ovbKipTHu+QY+P/znqU/gtL4v2vMnTe1ObWxTwBLl64vVYWCyrfJoNDW8zuuBqKqz5AHiafo4dSXk7FMNIimji/qF4H1G2DN8ecCYQZbEROXk8dP6nOHbkOPaHm7j06b/yyr7U+m6SMaNgDN+c5Xwp/PF7i9jUkNq6iBNzp3HOqI8DcM+OX1LRspuxpb+h0s7BKzZt1ZeZ/iJj4Jg+IteZQJTligIh7j37Ms6bOIPGSDufefYh/vjuyrT3rVwybgHnj5lHayzCv678EzuaDqaU3tmjPsbcomNpiYX57dbbCKtNoORuGmyLAE20VH0M1VaXSm8Y3Rhmq28PFBOIhoGg18v/nXox1885jpgqt77xPDe+/HRad30VEb51xMc4rmw6te3N/Mubf2RvuKbf6VlicfXEG5iUO53aSBV3bfkfxDOdaMFttNhCSPfTXPkxnM0lDSN9elhZwegnE4iGCY9lcfP8U/jVRy4i6PHy6NZ3ueAf97C2siJtefotL7fNu4p5xZM42NbAF1b8lnV1/d9a3G/5+eKUbzIiMIaK1t38YtP3iXkX0Jj7DdoVcuytNB88CzvW/4BnGD1TsO3Eh9FvJhANMxdNmcXjF32GmcXl7Gis4xNP/ZWfrHyRlmhqQ627E/T4+dn8zzC/ZArVbY18+Y3fsTiFLSPyvAV8dfr3mZgzlZr2Sn656fvUykLqc2+mxRZy9CCtlWcQaVvl4qswjDjF9BGlgQlEw9CM4jIWXfhprp09n6ja3PXOCs587A88s2NjWvqO8rxBfrXgWj4+/lja7Sjfe/shvr3mfmramvqXnq+AG6Z9h9kFR9Eca+I3W3/Cy/UemvNvp9H2EKQZqfkUjXXfM5vpGa4zKyu4L6MCkYhMFpFHRaQmftwrIuXpum84C3q9fH/hGTx6wdXMKhnB3qYGvvLCP7jw8T+zZNdm17eU8Foebpl9Cf8++xKCHh9L9r/DZcvv4MGdr9LWj4mvAU+QL065iYvGfAqPeHi56ll+vfNZNni+w367DEsgt/UBwgcW0tL8mKuvxRjmTI3IdZIpM9PjuwKuAvzAL3EWZL0J2AEcq6oJe9b7e1+HBQsW6MqVK915EUNU1La5f+NafrX2VQ62NAMwMb+IT8+ax8emHk5psP8LmCayN1zDj99dxIpqZ/5PeaCAKyedxAVj51Hkz+1zenvC2/nbrt+yt8UZmj4qOIbzCyqZJisIWc77u4lyrNzPkpt3LSJmrV+jR93uq1MYHKUnTLgm4bnFm29bpaoL0laqLJZJgehHwC3Akaq6Pv7cmcAS4DpVTTiFvr/3dTCB6H2t0Qh/2biGP727ir3NzkoMHhFOGD2RcyfO4JSxkxifX+RKXqrKiwfX87stS9nc6AyY8ImHk0fM4rSRh3N8+QwKfKGk07PVZmXNcp6qeIi6iDOnqMwrfLRoLzO8FfjiHy1t6qPVO4dAzqcI5VxgVvA2Euk5EI3/TMJzi7fcvgV4AXjCTGrtm0wKRFuB7ap6ZpfnNwB7VfUMN+/rkGwgWldxgIffWoeIHHqXHvpZwBJBEESc5y0ES+KHJXhEsCwLryV4LAuPZeGNHz6PB5/Hwu/xEPB6CXg9+L1egj4vIa+XoM9HyO8j1+/D7/Eg0u+NEJMSs22W7t7K3zauZfm+HUQ79bNMzC/i6BFjmFs2mtklI5hRXEZRIPmA0ZWtNi8f3MBju9/gtarNaHxvLY9YzCoYy5ziiRxeOI7p+aMYl1OK1+o5cETtCGvqVvBy1bPsaN4MQLEnzBl5+5gbbCDPev/9HlUhLCOwfUfi8x9DKHQmHs+4tP99jYzXfSAKjNITxl6d8Nzi7T8zNaJ+yohAJCLFQA1wm6re0uXc34ALVLXQrfs6SzYQPf3uRr7+6NO9XpduXssix+8jL+An1+8nL+AnLxAgPxD/ORggr9PzHdflBnzk+Pzk+L2E4oEt6PXisXruJqxtbeHZXZtZtmcbr1TsTLiad1kwh3H5hYzPK2RUTj4jc/IoC+VQFAhRHAiR7w84h89PwOPt9oP+QEsdz+1fx/LK9ayp3Umsy0ADv+VlTKiYsTkljAoWURbMpyxQQLE/lyJfDnm+IHneILneICGPj4NtFaytW8HbdW+yp2UHHmIcGapiXqiGKb6WDwSlDhEVWjRIuxRiSwl4RuL1jMbnHYvHMwqPZyRezwg8nlJE8kzQyk49BKKResKYqxKeW7zjDhOI+ilTGsvHxh/3JjhXARSISKGq1rt0X5/t3VFD8EAv22B3eftq59/jH1iHnhPnUHn/PrU6Pyfv/24BlvN8FJuG1jYaWl3a3sFWZ48VBez4Y/zzWbTTOZzmtGDAwg4qdgjsgGIHoKo1TFVrmDXJzEnqnI/N++nrB38WKx9vIIo3FMUbiOEJxGj3RdnRXMmO5uTWzFMb1Hb+gB6ZRIGvlR2BUl7wt5Hrb2NSoI4jgzWM94cZ4WmnyBPDJ4pPWoAWYD/E3oMYkKCnMabOS4gpRBFsFWJADMFWsOPPHXqpKvGXLO+PAo7/4zshVz7Q561dHt9/Pv5eSuE7pHb/WZt1IrbF0VNdHM6fAV/es02mBKL8+GOi/Qo6psrnAl0DSr/uE5HrgOsAJkyYkFQBG2tb8cQGvz9B4VBwUuFQgOp4/EAwS/Rc5+AngCXvp0vXD73EH1YSA0+TcygKHlCPol7nESv+uwVYeih/Oh49H8wn8X9rIYoP2n1OEGgEEcXy2Hg8MTwe52fLsrEsxbIUsdS5RhSxiB8az0Oo1xD1rSGIrwS0Rsax2BMl0HFY7YwLNDI1VMdoXzOl3jaKPO3kWjFyrBgBsfGJ4hfFA3jEeUlegQDa7SsxBlfEzX8WBWJmSoDbMiUQdbQP9fSWSfSv36/7VPVu4G5wmuaSKeDZCw/74F71/Xxza4Ib9UM/dP4x/kHa5bau6Wj8c1B7ur6jZnPoHsVGiakStW1sFFuVWPzRRlFVYvH7VTX+cavxb/Pvl0QP5emcI9qlbPGrbZx0bYnnI/F8AI3/DGBLR+qdXlNHevLhv+KHf7ed60SdYChOdUs7ql0ST1egDaVNnLsqgdXS6e/ekaCoE/QkhsQfA94Iub42cn1thDxtBL1Rgt4IPonis2y8EsMbf7QsxcLGE3+0xAmsnngV0Or0Cqx4/oI6b/BOe1O/3z/Z/0/X4VMXcsTUwxWjXUzQ1IhclymBqDH+mKjXO9TlGjfu67Mjpo/hiOlj3EjKMIyhShVivTTRZwgRuQT4b2A8sAn4mqouH9xSJZYpE1p3xR8TfW8ZA9SparOL9xmGYfTPEJjQKiJTgHuB64Ei4A7gcRHpcfDWYMmIQKSqdcB24OgEp+cBCYe19fc+wzCM/lGnjyjRkVkmAr9T1RdV1VbVv+J0U8wc5HIllBGBKO4R4EwROfSHik9MPQx4IA33GYZh9I2Cqp3wyCSq+oKqfrPjdxE5DsjDaaLLOJkUiG7DmRO0VES+ISLfAh7GWb7nL+BUN0Xk6ni1M+n7DMMwXJMhNSIR+aSIaILjni7XTcP5wv5dVa0d8IImIVMGK6CqlSJyCk5b5q04Q7IXATepasekmVOAPwHXAtv6cJ9hGEbqVDNp76HHeH8KS2eHVhEWkYXA48Cdqnr7QBWsrzJiZYXBJCKVwM4+3FIGVKWpOIbRHfO+GzhVqnpuohOFnjI9LnRBwpuebb43o1ZWEJHzgfuBG1X194Ndnp5kTI1osKhqn7aLEJGVmfRmM4YH877LFJk3Qi4REZkIPAh8VlUfGezy9CaT+ogMwzAym+LMI0p09IOI/E5ElnVzLpV91r6Os6rMn0WkqdNxar8KmmbDvkZkGIaRLAXUdqdGJCKfB74AvJjgXCnOlhJ+4Ce8v8/aHBHpdZ81Vf0a8DVXCjoATCDqu7sHuwDGsGTed5lAFU1xZQVxNsH6NvCDHi77BjCOD+6ztgJnn7VrgB73WRtqhv1gBcMwjGSJyGKcgSOJBDm0pC4Ad8fXtex8fxBYAczBWfngDGCLqp7a5bqU9lkbakyNyDAMI0ndjabrgyBQAFyuqg+JyI6uF8T3WZuCMx+yq9VA4mF7Q5gJRIZhGAOnAZiuqtEerhmwfdYyhRk11w0RmSoiLd2NMhGRr4jIhvg174jIFalcZxiJpDhyysgw8XXfegpCkPw+a1nDBKIERKQIZ9ZysJvz3wTuBN7GGZmyF7hfRC7vz3WGkUinkVPH4Yyc+hlwMbBERPyDWTYjrfq7P9uQZZrmuhCRWcCjdLNKbTxI/QD4m6peFX/ud8Ay4HYReVhVY8lel95XYwxxw2rklHHIgO2zlilMjagTEbkGWAOUAt0tiXERTrX4ro4n1Fl6906cDahO6ON1htGdK4BlHUEIQFWfAzbGzxnZadjts2YC0QfNwdk64kjglW6u6VhmZXWX51d3OZ/sdYbxIZ1GTq1KcHo15v2TtYbjPmsmEH3Qf6jqNap6oIdrxgK1qtq1I7Ei/jihj9cZRiJJjZwawPIYA2tY7bOW9X1EIjKql0uaVLUJoLdlM+LySW40S7LXGUYiyY6cypohvMYH3AZ8BmeftZ/hDJy6mSzdZy3rAxHv10C68yPgO31IzyLxaJaO5+w+XmcYiQy7kVPG+4bbPmvDIRB9sZfzb/UxvUYSj2bJ6XS+L9cZRiLDbuTUcKSqk3o4txE4f+BKM3iyPhClYUOoXUCJiAS6fDMZE3/c28frDCORYTdyyojcb+4AAAWLSURBVBi+zGCFvlsNCHBUl+fnxR/f7ON1hvEhw3HklDF8mUDUd0/hrLD7rx1PiIgFXI+z5fjrfbzOMLozrEZOGcNX1jfNuU1V/7+9ewuxqorjOP79NRphhDR5KRxDQwhKKpVE8UJFVlSC5Wg+RQ9dBIOKoIgumBDZheghCkyIYGzEDKGCKNFI8ymUnNEK0TQ18JaTlZST+u9hr0PHreeMOGfcZzy/z8txr71mr3UeDj/X3muv9ZukxcDCFCxrgdnANGBuabWEs61nVkVDzZyyxuUgOjeLgKPAAuA+YBswJyLyy7afbT2z0zTazClrXN4Yz8zMCuVnRGZmVigHkZmZFcpBZGZmhXIQmZlZoRxEZmZWKAeRmZkVykFkZmaFchCZmVmhHERmZlYoB5E1PEkjJB2SNLqX1/la0olquwJLulxSt6RP0vHStI6cWcNyEFnhJC2SFJImV6kzQNJPkv6S1FLjLrwNLI+Inb28ThvZb+r+KnVmAwP5f9HSRcB8STf0sm2zfstBZPWgI32OrVJnAdn2B69GxN5aNZwWFZ0FvFaDy60EjgGtVeo8AHSRbRNCROwG2skWNjVrSA4iqwed6fOMQSSpGXiJbKO4Wt/GegpYHxF7enuhiDgCfA5MlzQsfz6V3Qp8HBHdZac+Am7zqMgalYPI6sF2sk0EK42IFgLNwNMR8U+tGpU0EphJtrVCefkuSe9IeljSNkl/S/pO0kRJV0paIelPSb9KeiXtN1WyDGgi2/YjrzWdy+8ltA44DDxeq+9m1p94GwirC5I2Ai0RMTxXfi2wBfgmIm6vcZuPAEuAMRGxo6x8F1lgNJE9P7oIeB74HTiS+lPa6PAO4KGI+DD97cXAPmBTvr+S1gFXA6Mj98OTtAyYHhEja/kdzfoDj4isXnQAwyQNzZW/lT6f6IM2p5JtXPjzGc6NAO6MiNcjYjHwHtACbImIeRGxhGxSQjdZGAGQbrmtBG6RNKRULmlEam9ZPoSSDqCltzP3zPojB5HVi9OeE0maAdwNvBsRW/ugzWuAXRWCYUdEdJYdb0ufq0oFEXEUOABclfvb0u25WWVlcwFReYvvUhg6iKzhOIisXpwSRJKayEZDh8ieEfWFK4A/Kpzbnzs+nj4P5MpPcPrvaB2wG5hTVjaP7HbdjxXaK/VjSIXzZhcsB5HVi/yI6NH07xcioquP2jxJ5d/A8QrlPT5UTSOsdrKZcM2SRgETqTwaoqwfJ3q6vtmFxkFkdSEi9gEHgbGSBgMvA98D7+frSpopaaukFyXtldSVZrhNkLQhzWhbI2lQD83uJxsV9YU2YADZrLxWsoBpr1K/1I/8SMzsgjeg6A6YlekEJpC9MzQUmB0RJ89QbzzZ852dwCjgMeANYDXZTLZusof/9wIrqrT3CzBJUlNE1HQkEhFbJHUA9wDDgTUpbCsprRaxu5b9MOsPPCKyetIJDCabIbc8ItZXqDceWBIRbRFxHNgMXALMj4h9EXEY2EPP/9FaCwyi+ooOvdEG3AVMofptOYBJwPa00oJZQ3EQWT0pPSc6BjxTpd444NOy4xuBb1MAIUnAdcAPPbT3Jdlzomnn1NuetQOXkn2fVZUqpRdiJwNf9FE/zOqaX2i1fiW9m3MQaC5NYpC0FDgYEc+l4zFkL51eFhH/9nC9VcDQiJjatz2v2ocZwFfATRGxuah+mBXFIyLrb8aTvfvTlSvbmDvu7CmEkjeBKSm8ivIgsNohZI3KQWT9zThgU+lA0kDgek4NolPqVBMRG4DPgGdr2MezllZSaCVbQsisIfnWnDW8tPjpZuDm8jXnzlPbHwBHIuLJ89muWT1xEJmZWaF8a87MzArlIDIzs0I5iMzMrFAOIjMzK5SDyMzMCuUgMjOzQjmIzMysUP8B/0Lub4mceJ8AAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "OL = (1 - O * (1 + beta / alpha))\n", "fig = plotIhKinetics(Vm, CCa, OL, ylabel='$O_{L, \\infty}$')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Global gate activation" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAaIAAADzCAYAAADErhKJAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAIABJREFUeJzsnXeYVOXZuO9n2s72Su+9CYqgYAEVEVuMvSQaa2ISY9TE5EsxGjXJl2gSNb98JkZN7B2NXRFRFFFBAVF6hwUWlu11dtrz++PMIuKWmTln+3tf13vNMuec530ZZZ59uqgqBoPBYDB0FK6OPoDBYDAYejZGERkMBoOhQzGKyGAwGAwdilFEBoPBYOhQjCIyGAwGQ4diFJHBYDAYOhSjiAwGg6EbICLni8iijj5HMhhFZDAYDF0YEXGLyI3A44B09HmSwSgig8Fg6Nr8CTgz9tolMYrIYDAYujZ3qepMYGtHHyRZjCIyGAyGLoyqFnX0GexiFJHBYDAYOhSjiAwGg8HQoRhFZDAYDIYOxSgig8Fg6EBE5AERWdjMtWEi8oKIlMXWoyLSq52P2OaImUdkMBgMHYOIXAU8CLynqscfdC0fWAb4gL8BHuDnwDbgSFUNtuth2xBPRx/AYDAYugonn5CupWWRZq8v+7xhNRA44K37VfX+g+8TETdwE3BrC9v9FBgITFTVtbHnlgDzgcuABxI9f2fFKCKDwWCIk9KyCEvnDW72urvfxoCqTm1Jhoj4gSXAJOBR4MRmbr0IWNiohABU9W0RWR+71m0UkYkRGQwGQ5woSkjDza448QNZwIWqehnwtQdFJBcYjuWaO5jlQIvKrqthLCKDwWCIEwWitBhXzxaR+4FXVPWVZu6pAkaptqi5BsRedzVxrQjIEpFsVa1s7cxdAaOIDAaDIU4si6j5GBFQqapXtyhDNQpEW9kqM/Za18S1+thrOmAUUXfglFNO0TfffLOjj2EwGDoPLXawbsUicorGsElLm7WmzLoMPT5GVFJS0tFHMBgMXQQFQkSbXcRccyJyhs2tqmOvqU1cSz3oni5Pj7eIDAaDIV4UiLRce9mqay5OdsRe+zVxrT9Qoaq1DuzTKTCKyGAwGOJEUULt4JpT1QoR2Qoc3sTlycCnbX6IdqTHu+YMBoMhbhQiLSycc80BPA/MFpGxjW+IyGxgDPC0A/I7DcYiMhgMhjhRhFDLuQxOueYA7gQuBRaIyF+x6o/+B6u26HGH9ugUGEWUAJW7J+KRMNrs/4iy32hXjd0jcsD7EnvWBbhiP7tB3CgeEA+KF8QH+EBSQdIQVwbiykZcuXhc/fB4BuD1jsTlKkCkS46oNxi6JApE26k9p6ruE5GZwN3A7Vip3C8CP1fVhvY5RftgFFEC+KUBb7zf+8noBz3otRUiCkF8hCSTiKsA8YzE45tMmv9E3J4BrQswGAwJoUDQ4YiGqg5t4dp64DRHN+yEGEWUAK81nEN0/y8iCqooUSyDXYlqNKZ/IrH3oygR0ChCBIgS1XDs5zCiIUTCiIZxEcYjileieETxSRSvREiR6FeW3xV7FcUr4CeIn1KIlkJwPQRfg5rfU6c+gq5BuPwnkpF+MS5PU8k3BoMhUaLa4m+Z8XRWMByEUUQJcM7QP7WZ7IhGCEUbaIgGCETqqY/UURuuoS5SQ3GokqpwBRXBMsqDJZQE91ITriTdFaSPp44+nnp6e+rp5w1R4A6S7QrhlyB+3Qz1m6H+fqrohaSeTXrG93G5M1s/kMFg+BpRhCDulm5xMkbUYzCKqJPgFjdudxp+dxrZ3tbvrw1Xsyewi8K6LWyv3cTC2vVUVpUB4CHMhNR6pmaE6e/aS7bUkiH7oP5+wnUPUu+bRWbObbjc3W6+lsHQ5rRiERmSwCiiLkq6J5MRGWMZkWFldqoqRYGdrKlawadlH7CyvpCV9eAin6NyRjMtbRtZkU/IdgXJDL1NpHgBtf5zyci5DZE4NJ/BYEARgtqiRWRIAqOIugkiQv/UQfRPHcSJvc9gV/02Fu57g2Vli1lcsY6llV5O6nMbk1K24a1/lDxXPekNc6nb+xbenHvw+Y/t6L+CwdDpsbpvm/JLpzGfaDdERBiYNoxLhlzDb8bfzZTcYwhpiNf3zOWRPXuJ5L3GHu8F1EWFVKpwlV9JXdXfOvrYBkOnR9WyiJpbOFvQ2mMwiqibk5/Sm0uHXss1I2+iV0pfigKF3L3ht5R7zkEKXqVIe+MS8NfdS03pFbQ8IsVgMESRZhexZAWTMZcYRhH1EMZkHsLPxvwvk7KPpD5SxwNb/sJnVdvp2/cdCuVIIgppocXUllyMtjxvxWDosVgxIk+zy5Ac5pNLgHvWvU4oGsbn8uBzeUhxe/G7vaS6faS6faR7/GR6/WR6/GT70snxpuFxdZ7Apt+dypXDbmD+3hd5rehZni38N4JwVO/H2Fl+K30aniItsoK60stJy38EEfN7isFwICZG1DYYRZQAL+38hNpwYp01crxp9EnNoY8/mwGpeQxKz2d4Rh9GZ/Yjw+tvo5M2j4gwp+/ZeF0pvLjrMZ4pfBCPeDgy/zYKS8P0CT5HangJdeXXk57393Y/n8HQmTFZc22DUUQJcMPY0wmEgwQ1TDASpiEaIhAJE4gEqQs3UBtpoDpUT1WonopgLZWheipCdVSE6lhftftr8gal5XNY7lCm5A1nesEo8lIy2u3vckLv01CN8tLuJ3i68AF6+wcwJO/37CitY0DoNVKD82iofZKU9G+325kMhq5AVFu0iExnhSQwiigBzhw4NaH7IxqlvKGGvYFKigIV7KwtpbCulI3Ve9hcvYfCOuvPr+xahgvhsNyhzO43kVP7Tybdk9JGf4svmdXnG5SHSnh/3zwe2no3Px/zRwbn/5XtezcymA24qm4n6puCyzumzc9iMHQFFIi07JoznRWSwCiiNsQtLgr8WRT4s5jAoK9cC0cjbKguYkXZVpaWbuaT0s0sL9/K8vKt3Lt+HqcNmMylw2bSJzWnTc94Zv9LKKzbytbaDTy87f9xzchf07/gMYr3nkBvdx0NpZfi77MIEV+bnsNg6AooQsi45hzHRN06CI/LzfjsgVw8bAZ/m3o5b826idsmnc/k3KHURhp4bsfHnLfoLv6+/k2qQ/VteA4PVwy7nkxPNhtrVvNByXx8nlxS8u6jNiqkUE595R/abH+DoSuhChF1NbsMyWE+uU5ChtfPqf0n869pV/PEMT/mpL4TaYiGeWzr+1z4wT18tG9Dm+2d7c3jgkFXAfDq7meoCJaRmzqdMp9Vk+etf5poaGub7W8wdBUaLaLmliE5jCLqhIzK7McfDvsWDx91DRNzBlPSUM31yx7mL2teIRxtmxqfSTlHMDF7Kg3Rev6761EABuTdzp5IGm5R6st/hGo7TQQzGDoxEVzNLkNymE+uEzM+eyD3T7uaH46ag1tcPLvjI36y7NGEU8jj5dyBl+NzpfBZxRJWVy7H40rDnXULIYXU6CZC9S+0yb4GQ1fBWERtg1FEnRy3uLhixPE8MO375PrSWVK6kR8ufYDShmrH98r15XNav/MBeGX300Q1Su+MsymS0QCEq+5ENer4vgZDV8EaFe5qdhmSw3xyXYRDcgbx4LQfMCA1j3VVu7nu04epCQcc32dGwRxyvPkUBQpZWbEEESEv54/UR8FPubGKDD2aOCwi0/Q0CYwi6kIMSs/nwenfZ1BaPhuri/jViicdjxl5XF7m9D0LgDf3PE9Uo2T5J7LPNQ6ASPVfjFVk6NFEkGYXpulpUhhF1MXIT8nkb1Mvj7npNvHH1S86nkQwLe948nwF7AnsYkX5RwDk5PyW+iikaBlhYxUZeiiqQijqaXYZksMooi7IwLR87p5yGX63l1d2LePNos8cle9xeZjT9xwA3tzzAlGNkuM/nGKX1WEhVH2XyaAz9EispqctjoEwJIFRRF2U8dkD+dk4yw1955qX2V1X7qj8I/NmkOstoLhhN+urvwAgO/s3BBX8WkI0uNTR/QyGroAihKLuZld3RURGi8jZIvJ9Ebk69vMop+QbRdSFOWPAFI7vM57acAO3fvEcEQdjN27xcEzBiQB8sO8tAHL9R1KsvQCor77Hsb0Mhq5ET6kjEpFxIvI3EdkNrAXmAv8E7ov9vE5EdovIPSIyzs5e3euT62GICL+ecDb5KZl8Vr6Nuds/dlT+UfmzcIuH1VUrKG3Yh4jgTrscVfCHlqPRMkf3Mxg6O4oQVnezqzsgIiNEZC6wCrgKWAncBlwKnAacHvv59ti17wKrROQ5ERmezJ5GEXVxcnzp/HL8mQA8sHkBlcE6x2RneLOYnDMdRVlc8jYAfbIupiTqxSVKoPqfju1lMHQFrF5z0uzqJqwBhgKXA31U9VRVvV1Vn1DVN1X1jdjPt6nqqUAf4EpgeOzZhDGKqBsws/c4puYNpypUz783v+Oo7Bm95gDwcek7hKJBPK406r0zAZD6uWasuKFHoQjhqLvZ1U24QFWnqupjqlrb2s2qWquqj6jqFODCZDY0iqgbICLcMPZ0BOG5HR+zvWafY7KHpI1kYOpQaiM1rKywEhTysq+jLgo+aok0LHRsL4OhK9BKHVGXR1Vfau9njSLqJozO6sc3B04holH+vuFNx+SKCEflzwLg0/LFAGSljKNM+gMQqHnIsb0Mhs5OHBaR6ayQBI4rIrEocFquoXW+P+okUlxe3i9ey+bqPY7JPSx3Gi7crK/6nJpQFQCuVKvOKCW8HNW2acJqMHRGWqkj6hadFUSkXETeEZE/i8i3RGINJ9sIRxWRiJwHPA/8TEReEZFjbMh6QEQWxnnvUhHRJtbcZPfvihSkZPLNgVMAeGzrIsfkZniyGJs1kShRVlRYmXm9Mr5NZVRwEyZc/4ZjexkMnRlVekodUTZwHHAj8ASwVkQqReQ9EblbRC4RkfEi4og/0umeFL+KBawQkRzgZWBmokJE5CqslMD34rhXgHHAi1hK8EC2J7p3V+fbQ4/lhcKlzCtayfdHzaZfaq4jcqfkHsOaqs9YVr6YGb3mkOIpYI8MJZutBGsfw5t2liP7GAydGUWIdp/suJYoATKwvlcXAhOAKcBkYEbsHgXqRWQlsExVr0t2M6ddcxEROSz28zlAQmM9RcQtIrcADyTw2FCsD+wlVX38oOWcWdBFGJCWx+y+E4lolCe3feCY3InZU/G5Uthau4HShmIAPGkXoAopkVVotMaxvQyGzooCYXU1u7oRI4F/AecCvwDeVdVjsSylccDFwN3A0tiff2RnM6c/uQuAf4vIIqyDXhvvgyLiB5ZjFU49BuyK89EJsde1CZyzW3PpcMsIfWnnp1Q4VFeU4vZzSLbl9lte/iEAvTLOozzqwoWa8RCGHkNPmEekqlWq+hPgUGAD8LyIzAfGqep6VX1KVX+mqrNUNQ+w1e7H6U/ufGAvcDPQAJyawLN+IAu4UFUvA8JxPvcVRSQi6Qns2S0ZldmP6QWjCERCvLF7hWNyp+RaIb/GOJHPnU2124phhup6VDjO0ENRlU5nEYnINBFZISK1IvKhkz3gVHVdrGj1TGAw8JmI/D8RyT3ovi129nH6k7sSOENVFwKXADck8GwVMEpVn01wz0OAauAuEakGakRks4hclKCcbsXZA48A4KWdnzjWKXts5kRSXH521W+nLGjVKqWkngeAL7IBVecH9RkMnY2oSrOrvYl5kv4L/BnIBd4EHnZ6H1V9FeuX/l9jtffZICI/dCpZwWlFtBU4RkR8WBo07uZnqhpV1XitoAOZAGQCOVgf0JVYiukpEflOUw/Eusd+KiKf7tvnXPFnZ+LY3mPJ86WzpaaYVZWFjsj0uLyMzZoEwKrK5QAUZJwRy56LEgm868g+BkNnRYFw1NXs6gBOwEoZf1JVg8AfgEPsNiFtClUNq+pfgInARuD/sOJHtnH6k7sQOAK4Kyb7fxyW3xT3A9eq6nmq+l9VfQg4CtgC/FlEvpZTqar3x1pYTO3Vq1c7HLH98bo8nDbgcABeLPzEMbmHZFlxolWVywDwuXOpYgAADXWJGrMGQ9eiMWuus1hEwFhg3f7zWT23tmIlENhGRIaJyDdF5Nci8pSIrAI2AdMBIf4QSos4qohUtVpV/6qq16rqv5O0cBLd8z5Vvfeg9+qxEh76AOPb+gydlTNj7rn5ez6nJuyM22x89mEIwqaaNdRHrEQIl382AJ7QMjNG3NDt6WSD8dKB+oPeqwPS7AgVkY9EpApL6bwI3IKl9D7BMjBOAHJjFpJtuk+ax9cpjr1mdOgpOpAh6QVMzh1KIBJiftHnjsjM8GQxPH0MEY2wrmolADnp51MfBS8BNPSFI/sYDJ0R1U7nmqsDUg96Lw2wW08xDUs/PAmcCGSo6mRVvUJV/6aq76lqpc099pPUJ9fW0/oSOMcAEVkdqz06mLGx14RqmbobZwycCsC83Ssdk9mYxv1FzD2X4RtJmWYCEKh92rF9DIbOSCdzza0D9rffiYUihnOAu84GacC3gbeBjSIyN+aiO1VE+jggfz9xK6L2nNYXL6q6C6vA6nsiknXAWQdhzdJ4V1Wda7rWBTm+z3h8Lg8ryrdREqhyRGajIlpT9RkRDSMiRHzTAZDgQkf2MBg6I4oQibqaXR3Au0C+iFweSxK7CdioqnYVUTaW++1nwFNY7r+zgN8DrwK7RWSniLwsIreKyDftbNbqJ9cR0/paOMvwWI+jA+VeCwwEPhSR60XkN1h+zDA2q327AxkeP9MLRqEo7+xd7YjM3v5+9E7pT32klm21mwBITzufsEKKlqKRvY7sYzB0RpyOEbXUVzOWLPCCiJTF1qMisj/DKhYPPx3ru64UOAmrsYAtYvH+91T1blW9RFXHYymnGcBPsGLw5Vi1ordgpZAnTTy95tYAX2BZGC+0NigpVlB6HnBd7Fm/nQMexEzgIeAKrKw4VPVFETkLK7/9DizNvRCr750T5mmXZ3bfibxfvJa393zBBUOOckTm2KxJFO/bzfrqLxiRMZa8tKMpqfDQyx0mVD8PX8aljuxjMHQmVHHU8mmpr6aI5GNZPD6s7zYP8HNgkogcGUvXRlWXYWUrtymx7/7FsdV4Rj9W/7nJdmTH84m2+7S+mJyhqnr8Qe89rKqiqg8f9P5LqjpNVf2qmquqZxsl9CUzeo8jxeVhZfl2igPOxBfHZk4EYF2VlQThFh8B9zAAggHTjdvQXXEmfTvOvpo/xfL2nKiqd6jqH7B+yT8UuMzO3yIRRORmEfmGiAw4+JqqBlT1I1X9h509WlVEHTGtz+As6Z4Uju41xnLP7VnliMyRGeNxi5sddZupC1sJOu4Ua4CeN7zasW4OBkNnQqG1GFFBY7F8bF19sIwE+mpeBCxU1f19NFX1bWB97Fp7cRvwErBDRIpFZJ6I3CEiF4nIGCe6K3Tn9G3DAczua1kwb+9xJr06xe1nWPpoFGVDjRV7ykr7BgFtTONe78g+BkOnQi33XHMLKGkslo+t+5uQ0mpfzVgvt+HAsiaeXw5Mdezv1DpDsLpw/y/wKXAYlovwCazwS5WILG7+8dZpU0UkFmPacg9DfBzTawwpLg+fV+ygpKHaEZljMq12P+urLOWW6RtNedQqaWio79IDKg2GJlEgoq5mV5zE01ez0Q3WlLVUBGSJSHb8J08eVS1U1RdV9WZVPU1V+2CljN+O1dx6C5ZiTZq2toh8WBrT0MGkeVI4In8kAIv3OWOtjInFidZXW4pIRAh5rGbo0YZ3HNnDYOhctBojyhaR+0XkjOYkxNlXMzP22tQcl8ZOCh02aUBVN6nqbVjFrr2BU+zIaw/XXI8YZ9gVOLaXZZx+UOxMHsegtGGkuTMoDRZT0mClbPtSTwMgJboV1QZH9jEYOhPRqDS7sBqQXq2qdl0Cjd/NLQVbO7yflqp+BLyOFUdKGtuKSESOb8VENFHrTsIxva1mE0tLN9EQCdmW5xIXozMtC6gxey43bTZVUcFFFA0utb2HwdCZsGJB0uxykEb/+cHtew58zxkfu30+w6plShonLKI3gDIR2SQiz4rIL0XkJBHJc0C2wUH6+LMZndmP+kiQ5WXOdD5qdM9tjCUspHr6UkUOAIG6Vx3Zw2DoTESi0uwiDtdcnOyIvfZr4lp/oCKechonEJFtB7T3OVlECg665VBsugnjKWhtjUysmUCHx9YZwG+wtHa8474N7cSM3mPZUF3EB/vWcVSv0a0/0AqjYhbRxpo1qCoignoPh8gCYxEZuiWtWD6Vqvq1lO3E99AKEdmK9Z16MJOxstfaiy3ALOAcYh4uEdkJrMZqKn0MMM/OBrYtotiwpJWq+pCq/lhVj8HKoJiI1ffob3b3MDjHsb0s99wHxescqfUp8PUh25tHbbiaPYGdAKSknowqpER3Y3UgMRi6B+08j+h5YLaINDZwRkRmA2OAdusurKqzVDUPGIHVpOAOrH6jh2AZIW8CP7CzR0IWUaxw6RisGT8CbAcWHWwiqjWUZk1sPWbngAZnGZc9gDxfOkWBCrbU7GVEZl9b8kSEkRnjWFa+mE01a+mXOojc1BlUVQnZLkWDnyIpMxw6vcHQwWirFlG2iNwPvOJAwsKdWH08F4jIX7Hqj/4Hq7bocZuyE0ZVt2JNM5jrtOxEum9PwdKC72F13f4H8BpQKiIPxjped2t21lRSHWzo0l0DXOLimJhVtHjfBkdkjsqwZg82xolS3PnUkAtAoO51R/YwGDoLGpVmF85lzaGq+7D6a67Eqtm5AWtI3anazVJS47KIRGQIMB/IwVI+K7Dy20cCc4ArgXNE5Duq+lobnbXDOevVxympryXV46VPajoDMrIZlJnN0KxcRucUMDa3F/3SM3Gg40WbMr1gFK/sWsbS0k1cOnymbXkjMy1FtKlmLVGN4hIXEc8kiC5Eg0tsyzcYOhNO/x6qqkNbuLYea8pBuyEiJ6rqgiSfnR1rQ5QQ8brmfoWVlHCyqs4/aGMBzgfuAZ4XkbNU9c1ED9LZiaqS5fVRGwpSHw6xrbqCbdUVVo3zAfRKTWdyr/4c1W8wxw8YxtCs3E6nmI7IH4EgfFa+jUAkiN/tsyWvMU5UGSpjT2AX/VMHWXGimoWksAvVAFZ7LYOha6MK2jFzh9qTN0VkEXAX8IaqRlq6WUS8wDewLLajsBoZJES8imgO8MjBSghALT/VsyLyLlbmxKMiMkZVyxM9TGfGJcI7534PVaU6FGRvXTU7a6oorK5gS2UZGypKWFNWzL76Wt7asZG3dmzkNmB4Vi5nj5jAWSMmMCizXTpytEqOL50xWf1ZV7WLFWXbbGfPfTVOtIb+qYPISz2O6mohSxrjRMc6dHqDoWNpxSJyMkbUUUzGUkIvAyUiMh9YCmwGyrDyA/KAUcB0rIy6XOAtrD50CROvIupPK+mCqrpPRM7EGlF7LfC7ZA7U2RERsnwpZPlSGJXz1XR6VWVbVTmfFO9i0a6tLNq9jS1V5fx1xQf8dcUHzBk8imsmTeewXk2VBrQv0wtGsq5qF0tLNzmTxp0xfr8imtnrZFI8BRSTTRYVBOpeJ80oIkO3YH8sqDkcSd/uSFR1FTBHRI4CrgHOBL7F15sTCFbfvBeAf6rqJ8nuGa8iqqbpCt+voKqFIvI01kjZbqmIWkJEGJadx7DsPC4YNZFwNMoHu7fx382reWP7hv2W0vEDhvPbabMYlt1xNb9H5o/k4S3vsaRkkyPyDowTNdYTWXGi902cyNC96Lq5SgkRa9/zkYi4gSlY2dK9sD6BfVhTu1fEsqRtEa8iWgccDdwdx70rsAqfejwel4vjBw7n+IHDuamuhofWLOOxdStYuGsLH764ne9PPJIfH3o0Pre73c82KXcIfreXTTV7KAlUUeC31TyXAl8fsjw5VIUr2Nuwm77+AfhS50Dt+6ToLlSDiNiLRRkMHY7SmkXU7YjFiJbGVpsQb9TteeBMEZkUp8yU5I/UPemdlsEvph7HwnO/x3kjDyEYjfD3lR9x4RtPsbumqt3P43N5ODzXmqi6tHSzbXkiwogMKy18c401xys3dSbV+/vOfW57D4OhcyAtLMda/PQo4lVE9wGFwEsiMq6Ve2diFboamqAgNZ2/zDiNZ0/9Fv3TM1mxbzenv/wIHxa1/0c2rWAUAEtKNzoib/h+RWR190719KVaMwBoCNjqAGIwdB6iLSwH64h6EnEpIlUNYM1KzwE+EZHbRaT/wfeJyDXA2VhjZQ0tcGTfQbz6zcuYOWAY5Q31XD5/Lu8U2rdMEjpDbD7RJ6WbHSnSHZFuKaItNV/OOwp7rNET0YYPbcs3GDocBVSaX4akiDshXlVXYKXpFWE1Nd0mIotF5NFYZ9YNwN+xusbe0San7Wbk+dN4+KTz+M7YyQQjEa5e8F9e39Z+I7aHZ/Qmz5dOSUM122v32ZbXL3UQqe40ykMllAVLAPCmHG+9RrfhQEzTYOhwWhkVbkiChCqzYsroUOAXWArnKOASrOSEEcCrwLHdrYaoLXGJcPv02Vx9yBGENcp1C19h8e72cdOJCFPzRgDwadkW2/Jc4mJYupUKviXmnstKO5H6KHgIoaH2U7IGQ5sRleaXISkSLhFW1TpV/bOqjgSGAscBM4BeqnqmqprRDwkiIvxq6vF8b4KljH747otsqihtl72n5scUUal9RQQwPP2rcaIM73Aq1cpdCTUk1TXEYOhUiDa/DMlhq1eFqu5Q1UWqulhVy5w6VE9ERPjVEcczZ/AoqoINXPX285QH2n6EwtT84QAsK9tC1AHXWWPm3JZaSxGJCCG3pexCgfdtyzcYOhRtwRpydjBejyLuMRAikgJchTWDYi/wpKp+Ld0qNi/j16o6y7FT9hBcItwz83TOf/1JVpcV88vFb3LfrLPatFfdgNQ8+viz2RuoZFP1XkZn2ev6MDhtOB7xsiewi9pwNemeTFwpR0PDGryR9fuLXQ2GLkvLlk+X76wgIrfYFPGoqm5L5IF4u2+nAYuBScSS5YFfi8jNqnrnQbf3wXLXGZIgzevjvllnc+pLDzFvx0bmblrF+aMmttl+jXGi13Yv59OyzbYVkcflZUj6CDbXrGNLzXom5kwl0z+HYOBBfFIPkZ3g6fYTQwzdme6fc3MrlrpN5jcloHSgAAAgAElEQVRGBT4AtiXyULyuuRuwkhT+F0sZnQ4sB/4oIv9IZMOuzJtrNvDexq2s3VNMWW1dm80lGpSZzW3TZwNw25IFFFZXtMk+jTS655yKE+1P4661khOy/BOoiFq/84Qa3nNkD4OhQ+g56ds/AYYluKaQnPKK2zV3PvCMqt4c+/MqEZmHNQb8RyLiVdXvJXOArsSNz71G+ADVner1MCg3h2H5uYzuXcDY3gUcNrAfBRnptvc6Z8QE3t6xiTe2b+DnH7zBU6dc1GYural5liJaUb6VcDSCx2Wv5dDwjLGw98uEBZd4aXANBLYRCryDL/0Su0c2GDoM6f4WEUCJqiaUvisiNcluFq8iGoE1lXU/sUZ3PxaRCuAmEXGr6pXJHqSzE1WFesXtFqJuUDfUh8JsKC5hQ3EJ89Z+GS4blJPN0cMHc9zIYRw1bDBpPm/C+4kIfzh6Dkv2FPLxnkJe27aebwwb2/qDSdAnNYdBafkU1pWyvmo3E3Lsuc6GpY9CEArrthKMNuBzpSC+IyC0DVdolUOnNhgMbcRRQDLdkCtiz65J9MF4XXMBoMlv05iVdDtwuYg8lIDMLkWwIUTe9ghZ28Nk7oiQUaik7Vb8+8BXDt4acDUAUSisqOSZ5V9wzbMvc9Rf7uNXL8/jk+07E3bl5fnT+PkUa4Lq/37yLnWhYBv8zSymxKwiJ+qJ/O40BqQOIUqE7bXW/8+pqXOIKKRQgUZNmZmh6yJRaXZ1B1R1iaomXD+iqpHYs9WJPhuvRbQZawDSvc0c4FYRUeC3wPGJHqIr4Pf7WPTUz7/y3oZ1O5k373OWrdnBzqoaqt1KKM1NxC9EUiDih4AvzAsr1/DCyjVM6Nuba2ZOZ9bo4bjidLNdMGoiT6z/jFWle7nviyX89PAZbfHXY0recF7c+QnLy7Zy2XD7uSbDM8aws34bW2rXMypzArn+qVSUC3luJdLwEZ7Udp1+bDA4g9Ja1lx3GIzX7sSriOYD14tIZnPaTlVviymjW+khEztGjx3I6LEDv/JeXV2A//z7XV77eC2l7ijBTBfhdCGUBqv3FPOjZ19mfJ9e3P6N2Uzs37fVPdwuF7dOO5HzXn+S+1Yt5cLRhzIgw97IhqaYkmd14v6sfJszcaL0sby/b97+vnMeVzp10ps89tJQ/5ZRRIYuSysxoi6fvt0RxOtGewz4B9Zo2GZR1duBHwOP2jxXlyUtzc+1Pz6VN574KUsf/Rm3HD+d/oVR0ncpvkqQCKzZu4/z/v0Ut7/xDvWhUKsyp/YZyDeHjyMYifB/Kz9qk3MX+LMYkl5AfSTImkr7zTGGZ1jNTrfWbiASG3mv3kMBkNAK2/INhg5DW1iGpIi3+/ZGVf2Vqi6P4957VfUK+0frHpxz/tG89dSNvPm7KxhT4iZ9ZxRvDaDwxKcrOeu+x9lW2nrM5IbDjsElwnMbv6CwurJNznp4LE603IE4UbY3lwJfHxqiAXbX7wAgJfUkVMEX3YPV0N1g6FqIWhZRc8uQHN0ysaAz0rtvLnMfvo7/3vwd+u6I4C9WJAzbKio445+P8t7GrS0+Pzw7j7OGjyesUe79vG2sosY07mUOKCL40ipqTOPO9R9DtQouUTRorCJDF6Vn1BEBICJHt8c+RhG1M4OH9Obtp37GtZMmkl4YxV0PQY3y/ade5PXVLXen/vGhR+GOWUU72qDIdXIsTrSyYjuhaNi2vOEH9Z1L8RRQrVZ8K1D/pm35BkNH0MMsolfbYxOjiDqI7119Es/ddDE528J4a6xfpn7y/Ou8sHJ1s88My87jrBETiKjyj88/dvxMBSmZDEvvRSASciRONOKATtyNqesRz3gANLjEtnyDoUPoWTGidjHzjCLqQIaP6Muix35Kv8Io3ipA4FcvvcWC9c1Pav3RpOkI8MLm1ZTU1zp+psP31xPZnxbbK6UvmZ5sasJVFDcUAeD1nwCAL7odVftWl8HQrvS8GFG7qFejiDoYt9vN/CdvZEgR+5XRtc+8zKqivU3ePzw7jxMHjSQYifDYOufjLFMax0I40HdORPaPhdhcsxaAnLRZ1EbBTQQNJVyAbTB0PN3YIhKRrSKypXEBWbGft8b+3CbYUkRiMVhEfE4dqKfy2uM/YXAReOqssSYXP/g0RZVNFyh/95CpADy+7jMC4dbTvxOhsZ7o84odNETsyx6RMQ74MmEhzTOISk0DoKF+nm35BkN7080H4x0PnBBbs4Da2M+N77cJdi2iPGArcKwDZ+nxvPLIdeRtC+NqgABRLnvgGSLRr9v70/oM4pD8PpQG6vjvZmetilxfBiMz+xKMhllVUWhb3pcW0ZeD8sIeK5suElxkW77B0O60bBF16cF4qrr9gLUNiBz4Xlvt64RrrvvlLHYQbrebN+7/EZmFESQC2+uruWv+17+sRYTvTTgCgH+v/tTxcRRT82Ljwx2IE/XzDyLVnU55qITShn0AeFKOB8AX2Yw6MBXWYGg3Wo8RVarq1d2ovY9JVuiJZGSm8cBPziVlXxQUHlyyjGWFX89gO23YGPqmZbCpspSPinY4eoYjHJxP5BIXI9Ib64msOFF26knUR8FDCA1/bcivwdC56cYxoiZ4sj02MYqoE3LoYcM4t+9gvDUKInz3P3MJhL6aYeZ1ublotNUy5/H1nzm6/+TcYbgQVlUWUhdusC1vf5woVk+U4RtBpfoBCAbesi3fYGgvhJ6VNaeqP2qPfRJSRCIy88AFNFbdTmrimsEGt9x8Pv13gYSgzhXldy/O/9o9F42ehFuEt7ZvpLgu6ZlUXyPD62dc9gAiGmVluX238JcJC5ZFJCIE3Zb7LxwwE1sNXYgWEhW6SbJCh5CoRbQQePeA9WLs/b8e8F7jPQab/Pe+H5K612oYOnftWjbt++qIkL7pmcweNJKwRnl6w+eO7v1lnMi+e25g2lB8rhT2NeyhMlQGgDvFGmfhjWxos5HrBkObEG1hGZIiUUXUmNLXuM6Jvf/zA95rvMdgk/TMVC6fMA5PreWiu+r+Z7/2pX3J2MkAPLV+JeEmMuySZer+OJH9hAW3uPd3WdhYbVlFWamn0KDgJQCRNitPMBgcx1hEzpOQIlLV9w5cwAexS581cc3gANdfdzoFhVGIwJ5ogLnLvvjK9WP6D2FoZg5FddW8u9O+0mjk0NwheMXNuqrdVAbrbMsblTkBgI01VgujrJQxlEet8rOg6Ttn6Ep0wWQFETlfRFqtl4jdd4+IXCEinoOuvdZW5zPJCl2A/9x+CSkVlrVz28sLCIa/TFxwifDtMYcB8IyD7jm/28ek3MEo6kga9+iYItpQvQoAETdBtzXeKhx427Z8g6Fd6GItfkTELSI3Ao/TSiq2iFwP/B3IBH4JLBaRvANuaZvx0BhF1CUYMaIvhzekIyEIeeHvb33wlevnjJyAR1y8u3MLe+sSHhffLNPyLUWxpGSTbVkDUoeS6k6nLLiP0oZiADz+E4HGOFEn/FdsMDRF17KI/gScGXttjWuAk1X1KmA8sAx45wBl1GY1RUYRdREe/Pt38ZdYiQsPLl1OVeDLwXIFqenMHjySiCpzNzXfvTtRjiwYCcCSko22Ewpc4mJUhtV5e0PMPZeTesqX9UShdfYOazC0E13JIgLuUtWZWB1wWqOfqq4EUNWIql4DLADeFZF82lDV2lVEZcAwYLEDZzG0gMfj4Yw+g3A1QNQt3PTk61+5fuGoiQA8u+Fzx7LQxmT1J9ubRlGggsK60tYfaIX9caJqSxGl+0ZQsb/vXHcpRDd0a1qyhjrAIhKR80REm1gPA6hqUQLi9onIsAPfUNUb+TIjen/MSEROcuD4+7GliNRiu6rar3o0tMrvfnshGXus+ND8wm2U1dXvvzZzwDD6pmWwvbqCj/fY7xEH4BYXR+RbadxLS+2750ZnfKmIVBURIeSxaowiDSa/xdD5ETpd1tx/sWI6B6/vJyFrAXD5wW+q6g1YZTn+A97+feMPIrJSRB4SketEZIaIZCa6sXHNdTGuOWYK7oCiLuGXj7y8/323y8X5jVbRRueSFqYVWHGij0vst+Lp4x9AlieHqnAFewNW2yKv/2QAUqJbUXW2k7jB0BZ0JkUUc6HVNLGSMQ6uBe5oZp/rgKEH/HnaAZdvwnL9zcJKiigXkY0i8my8GxtF1MW44opZpO+xYkXvFe/6ilV0QUwRvb5tA1VBZ4zUaflWnGhZ2RbC0YgtWSKy3z23ocbKnstLP5maxvlEQWdbFRkMbUI3LWhV1aCqNluroapNNrVU1VdV9XZVPUtVhwB9gB8Bn8a7t1FEXZDrjjsCd0DBJfzi4Zf2vz8oM4ej+g6mIRLmlS1rHdmrb2oOQ9ILqA03sKrSvstvdOYhAKyrsuqhUj39qNRsABrqXmz2OYOhU9AGLX5E5AERWdjMtWEi8oKIlMXWoyLSy8bfoM1R1VJVfUtV74z3mUR7zYmIHCsiV4vI90XkFBFJT/yoBjt859LjySiyrJP39+2mov7LDLoLRje6575o8tlkaEzj/mjfBtuyxmVNAqzC1lA0CEDEa3WHUDOfyNAFcDJrTkSuAr7bzLV8rCSB6Vgus78C3wTmJzqMVFUfVtWk5sbFvvd/JiLvi8gSEXlcRG4UkRNEJDsZmQcTtyISkSnAWuA94J/AP4DXgFIReVBEBjlxIEN83HDi9P1W0U2Pvrr//VOGjCbT62NlSREbyksc2evY3lZ7ng/2rbctK9ubx4DUIQSjDfuH5aWmn0tUISW6B42W2d7DYGhTHMiaixWa3gI80MJtPwUGAieq6h2q+gfgPOBQ4LKkzp4cPwfuBPoDEeDbwJ+Bt4EyEdmcSDyoKeJSRCIyBJgPjAZeB/6AFaB6CCgGrgRWisjpdg5jiJ9vXzyD9GLrV7AFu3dQF7QC/akeL2cMtzLRnnPIKjo8bxipbh8bq4vYU19hW964LGt8xZoqKyaUnzqDsqgLEQiZdj+GzkzrnRUKROTTA9bVB4sQET+wHLgNeAz4+sAxi4uAhaq638+uqm8D62PX2osrgA+BMVgWGcAPgD8CQaxsupPtbBCvRfQrrJTAk1X1DFW9RVX/pKrfBYZgfSgB4HkROcXOgQzxc9WUibiCinqEPz735ZiIxuy5FzavJmQzwQDA5/IwLVbc+sE++4Wn47MsV9zaqpUAeFxp1Lut8oVQ/UvNPmcwdApatohKVHXqAev+JiT4gSzgQlW9DAgffIOI5ALDsbobHMxyYKoDf5N4GQw8o6oRvrT7Nqnqb4A5WOcfa2eDeBXRHOARVf3aUJxYLdGzWObiGuDR2IdoaGN+8IOTSd1n/Ro2d/06ghFL6RxW0I9ROfmUBupYsMOZRqjH9rL+P1vsgHtuaPooUt1pFDfspqRhLwCulDkAeMOrsf5/Nxg6H3EMxssWkftF5IwWxFQBo2Lfm80xIPbalLVUBGQ5FZ+Jg1osy4cDXv0AqroIq/7oZjsbxKuI+tNKKp6q7sPqaZSOlY9uaAfOHzocCUHUK9z76vuAlSZ94WgrKeAZh2qKjulljfv+pHQz9eFgK3e3jFvcjMm0rLZGqygv/WzqouAhiIZW2juswdCGiGqzC6hU1atVtdlWIaoaVdWvWUEH0VgU2lQ6dWPNRnslim0GRgCoanVs/wNzApbypcsuKeJVRNVAams3qWoh8DRwlp1DGeLnV788h9RS61exf3+ygmisvc85Iybgdbl4b9dWimrtN0LNT8lkQvZAgtGwI924D44TpfuGUqFZAARq59qWbzC0Ce3Xfbvxu7mlFIj2qlyaDxwY/18KfOOAP/cGcuxsEK8iWseXY8FbYwWWb9PQTszO6YtEIJQiPP+R9cWe509jzuBRRFWZ61DSQqN7blGx/ThRoyLaWL2ahoiVfh7xHWldDJp2P4ZOTMsxonhcc/HQ+NtjUwZA6kH3tDV3A7eLSOO+/wROj9U03QT8BLDleolXET0PnCkik+KUmZL8kQyJ8sfbLsRXbv1y9IdXF+5//6LR1pf9Mxu/2G8p2WFGLI37/eI1RGyObcj25jEkbSQhDbK22nLFpaWdT1jBr/vQ8E7b5zUY2oJWLKJWXXNx0tjFoF8T1/oDFapaa3OPuFDVclV9RlXrY39+FvgLcDHwO6yU7p/Y2SNeRXQfUAi8JCLjWrl3JrDdzqEMieHxeJgSzYQo1KfCB+ss19kx/YcwICOLnTWVfFhk/z/JqMx+DErLpyxYy4qyeLrKt8yhOZYFtLJ8CQD5acdSEpva2lD7lG35BoPjtEFnhSa3Ua3A6t92eBOXJ5NA+5xkiBWx/kBEpjdzvv/ByqY7Ghiqqkvs7BeXIlLVAFYhVQ7wiYjcLiL9D75PRK4BzgZMDm47839/uhRflfUv4cZYgatLZH//uSfX208AEBFO7GvJe3uPfXdfoyJaXbWCUDSIS7wEPda/u2jg1ZYeNRg6hDiy5pzkeWC2iOxPjRaR2Vj1PE87vttX+TZwL3BMczeo6i5V/VhVa+xuFndnBVVdgdVdtQj4DbBNRBbH/IRzRWQD1pjZHTTTwdXQdqSn+xlZ6QWFitQo64qsKagXjpqEW4S3tm90ZHrr7JgienfvattNUAtS+jAwdSgN0QDrqy3FlpZxMREFvxahkT22z2swOI5q88u5GBFY3QzKgAUi8lMR+TUwF6u26HEH5LfExViFs3e1dFOs1dtv7ZbsJNRrLqaMDgV+gaVwjgIuAc7BSu97FThWVcvtHMqQHPf/76V4aiyr6Af3WiUKfdMzmTN4FGGN8vQG+6ncozL7MjitgPJgLcvLHXTPVSwFID/tBEqjXgAaap+xLd9gcJTWs+acihE1lsTMBFYCtwM3AC8Cp7bDDLipwHPa+pTN+Vg1ROfZ2Szh7tuqWqeqf1bVkVjzKY4DZgC9VPVMVW2uXYWhjenVO5uBJQIKRb4gO0qt3we+M9bqZPDk+pW2Oy2ICLP7WVbRgiLn3HNfVC4jHA3jFh8NnsMAiARebulRg6FDcNo1p6pDVfX4Zq6tV9XTVDVDVXur6uUxBdXWZAPbWrtJVbdg9Zz7Rmv3toTdCa07VHWRqi5WVdOtshPwj5u+hadOQYQf/N2yio7qN5gR2XnsratxpNPCSTH33DsOuOf6+AfQ1z+Q+kgtG6qtGUWpGRcTVfBHCtFIe/ybMxjix4HOCl2BMiA/znvfB1pLYmuRVhWRiJyYrPBYYM3QjowY0ZfeVniIzVJLWW0dIrLfKnp03XLbewzP6MOw9F5Uhur4qMT+aIjDc48CYGmZ1RkiP/VESqNeRKCh5mHb8g0Gx1BaixE55prrYNZjebvioYim08zjJh6L6E0ReUdEviEi7tZuFhGviJwtIu9hdeo2tDN3/fhM3PXWiIjv3WOlQZ8z8hDSPF4+LNrB2rJiW/JFhG8MnALASzvtZ5EekTcTQfii8lNqwzW4XSnUeS2XnQbm0rqb2mBoP9oxa64jeRE4VUSmtXon2B7UF48imozVXfVloEhEnhCR62OK6WgROUZEzohldTyLpR3nYvVIOszuAQ2JM3nyCPL2WF/eq8KVlNfVk+VL4YJRVj3yA6s+sb3Haf0n4xYXi/etZ1+gypasPF8BozMPIawhlpd/CEBO5o9oUEjRcjRk34ozGJxAaJ86ok7Av4DdwAsiMrmVe08ANtnZrFVFpKqrVHUOVj75POAMrJYPLwGLsPyDL2JV2s7BUljTVfVUVV1j53CG5Lnr6jP2D8676q4nAbhqwhRcIry8Za3t/nP5KZnM7D2OiEZ5bZd9RTEtz/ICLCldCECOfwolamWE1lX/w7Z8g8ERWnLLdSPLPdZF4Tys5qsfisifRWTEwfeJyHVY3/sv2tkvkTqij1T1O0Au1ujaq4BfYqVyXwkcAeSp6pWqav9XboMtjpw+er9VtDpcSUV9gEGZOZw6ZDRhjfLwmqbGnCTGWQOPACz3XNRmy5+JOUeQ6k6nsH4ru+q3IyKQamWE+kIfYdVUGwwdTw9JVkBVl2LFibYBNwLrRWSFiDwtIs/HakfvBrbQSr1RaySTvh1R1aWxGeh/VtW/qOojqrpM1ea3kcFR7vnBN3EFFNzClX95AoCrD7FiL0+u/4zqoL1ShCMLRtLHn82u+jKW2Wz543P59ictNFpFvbOuojLqwkOYUO1ztuQbDE7RimuuuyQrAF+pHb0G+AKYBFyA1UFnOFbt6HGx8RBJYyt929C5mXrEKAoaraJoJUUVlRzaqx/T+g6iOhTksXUrbMl3i4szBlhJC89u/8j2eY/KPwGAJaXvEYjUkeLOo8o1HoBw7QMmacHQ8SgQ0eZXN0RVg6p6n6pOxhr5MBWYxpe1o7vt7mEUUTfnn9efh7vOihV9+89WV5BrJ1mWx7++WEqVTavonMHT8Lk8vF+8lm019rLxBqUNZ0TGWALRej4seQeA9KwbCCr4dQ/RhsW25BsMTtBDkhWaRFVLVXW5qn7iZAcdo4i6OYdMGsLgfS5Q2J3SwMothRzbfwhH9hlIZTDAf1bbS78uSMnk9P6TUZTHti6yfd5ZvS3X+nv73iCiYfJTZ7BX+wIQqPqTbfkGg10kqs2u7kAsGzreYtYDn3PHns1s/e6vYhRRD+CJO67EW2N1W7jq/rmICDcePgOAf6/+hPJAfSsSWuaSYTNxIbyx+zP2BiptyRqfdRh9/AOoCJWxvPwjRISUzOuJKKRGNxANrbcl32CwRUtD8bqHHgIrG/rkJJ7LiT17RKIPGkXUA8jPz+LwhgyIQnU6PDLvY6b1HcSM/kOpDgX516qltuQPSs/nhL4TCGuEp7fZc5+5xMWs3tZU4gV7X0FV6ZNxFnuj2QDUVf7BlnyDwQ4CSESbXXSPrDkB8kVkcCILGBR7NmGMIuoh/OevV5ISm+J6x6IPCYbD/CxmFf1nzadsq7Ln7r10mFUH9ELhUkpsFrhOzT2WLE8ORYFCVlYsRcSNK/0qVMEf+phoyH6/PIMhWUS12UX3yZq7B2swXyJrGUnahR775zV0BTweD9dOnMRd278gkiJ864+P8PzNV3HuyEN4ftMqbv14AQ+ddK5Vv5ME47IHMLP3ON4vXss/N87n5onnJn9Wl5eT+53Dc4X/4eXdT3JI9uH0yb6SPbX/op+7lvqKn5He679JyzcYkkYVukksqAVus/n8lkQfMIqoB3H1VSfx6Pc+Z99gYVW0ks82F/LLqcfx1o6NLNy1hfk7NjFnyKik5V835lQ+3LeBV3ct54IhRzEm62tDfOPmqPxZLNo3jz2BXSwqeYsTep+OO/N/CNf+ltTIakKB9/H6ZyYt32BIlu6eHaeqdhVRwhjXXA/jmd9eZg3Pcwnf+fdz5KakcuPkYwG4bekC6kLBpGUPTi/g/MHTUZS7175mq+7HLW6+2f9iAObt+S+14Rr6ZF7EHizlFqr8NaZ+2tDuaKsxIkMSGEXUwxgwsIDT/f2QiBJMFc69/UEuGTuZ8Xm92VVTxe1L37El/6qRs8j2prG8fCtvFdmbCDs+6zBGZxxCfaSWV3Y/iYiQlXsnAQW/FtNQc58t+QZDUvSAXnPtjVFEPZA7b/sWvXYqKKxz1fD0/KXcNeN0fG43T2/4nDe2JZ8ineVN5UejrczPO9e8xN76iqRliQjnDLwUt3j4qPRdVlUuJ8d/JMUuq02Ru+b/mcQFQ7vT3euIOgKjiHoor9/1Q3wVURDh94s/xBeAm6ZaLXZ+uXiere7cZw6cyrG9xlIdDnDbF3NtNUTtlzqIb/S7EICndvyLmlAV/fP/j+JIKm6JUl92Jar2psQaDAnRskXUHdK32x2jiHoomVlp3HrsMbgaFPUKZ/zfY5w9dCyzBg6nMhjg6gUvJB0vEhFuOuRscn3pfFq2hce3fmDrrMf3Po2RGeOoCVfxVOH9eNxZ+HLuJqCQqkXUVd5iS77BEC+izceHYjGi7pK+3a4YRdSDOfec6ZyZ0g8JK+FU4YRb/8Udx5zCkMwcvijdy4/fe4VwNDlrJj8lk98ccg4A926Yx7t7Vyd9Tpe4uHjID0l1p7Gqchmv7H6KvPRZlHitwtfUwHMEah5LWr7BkBAmRuQ4RhH1cP5087eYVJYKUajJhNNvfZAHTzyHnBQ/Cwo3c/NHbxFN8h/YjN7j+P6o2SjKLSuf4fPyHUmfM8/XiyuG3oALN+8Uv8oH++YzMP9OCnUYAJ7q3xMMvJu0fIMhLnpg9+32wCgiA8/+/YcM2gVEoTJTueiPj3HfCWfhc7t5asPn/GzR64SiycVhrhx+AmcOnEpDNMxPlz/CmsqdSZ9zTNZELhr8PQDm7nyIT8s+on/vueyJ5uMShfJrCAZMh25D29JKZwVDEhhFZADg7ft/wrDdLogqNZnwvb/N5feHnUiax8sLm1fzw3deTCpmJCL8YvyZzOg1lqpQPT9Y+gAf7duQ9Dmn5R/Haf3OR1Ge2PFP3i9ZSG7B85RE0/BIBCm/kkDtC0nLNxhaRiEabX4ZksIoIsN+3vzX9Ywt9lk1RunCLc+/w3m+UeSk+Hm7cDNnvPIoa8sSnznkcbm5Y/LFnNZ/MoFIiJ8uf5Snti1OOpvu5L7ncPaA7wDw8u4neXnP6/jzX2BvNAe3KJ6qX1Jb8XtT8GpwHsXEiNoAo4gMX+Gle6/lTPcAK5vOIzyzfgO+VRGGZ+ayubKMM199jAdXfZKwq87jcvPbiedx2bDjiGiUu9e9xvWfPkxxkmMjju99GpcN/TFu8fBh6QL+vvleGjIfZLcOwiWQGniUmuLZRELJx6UMhqYwnRWcp1MpIhEZJiIviEhZbD0qIr3a6jlD09xx84X886STSSmLAkJlKuz6oope1SkEIxF+/8m7nPLiQ7y7c0tCbXxEhB+NOZk7J19CjjeNJaWbuGDR3Ty4aQF14cQnxR6eezQ/Hf07+voHUNywm3s2/YkPguex3XUiQYV03Umk5CSqym4gGq1JWL7B0CTGInIcsdMPzEliEwGXAT7gb7zSHpAAAA+VSURBVFgNWX8ObAOOVNUmAxTJPtfI1KlT9dNP7U0p7a6Ew2G+deP9rMoK8P/bu/vgOOr7juPvz96DZFmykC3bwpaNoYBJ4jixg8FAYEgggUBpQyGBmWYgnQRCk3ZKQhKaJmkoU1JiSsNMO2TG0EmTAUx4KBlCWgdKeAozPJrYhoAdsI2JCUb4+UGWdLff/rF7tny+O91JJ99J933N3Ozd7m93f7JP+t7v9ru/b5iOZuUOk1mYbAyko/fNByZP44p5i/jUnLk0JcqfQ/e9vl3c+MrPefLdVwGYnJ7IRbMXc2H3IjqbJ1XUz/6wn4fevpunen5FSEg6aOLMjm5OTD7E1GBv1MaS7Et/nLZJf08i1V3R8V3DKToFfXtzl506+/KiOy7//ZIXzezEUenVOFZPgegG4Frgg2b2arzubOAR4Eozu62a++V4IBrayy9v5Eu33seWTrCUMIxsa0jYGmJx7JmUbuK8OXP51Jy5nDy9m+Zkqqxjr9i6jn9fs5xX4my6hAJO7Tyej3XN4/SpJ9Cebim7n5v3vc1Db/+MVTuiQn/C+GT7Pk5pXkt7kAGiGfx3B92o6Rza2r5AkOis4F/CNYjSgWjWZUV3XP76Ta8DjwG/8Jtay1dPgegNYL2ZnZ23/jVgk5mdVc39cjwQlW/N2k1cdcu9bD4iJNssTEbYYmRbs1j6QLsE4sRpM1nUNYsPTe3ihI5pzGydRFCk1pGZ8eLWddy78RmefPdVsnGSgRDHtXWxcPLRvL+9m+PajuSoiZ0kg0TJfr7du5Gneh7mhW2/oT/sIyDLyS3vcNrErcxI9JPrhhnspZWB5HEE6YU0N51BOn0iQVBeEHXjVvFA1NRlp878XNEdl6+/2UdEw1AXgUhSB7AVWGJm1+Ztuws438zaq7XfYB6IKpfNZvnnWx7kl+vXsXuSyDaLMGWELSFhc3hQUMpRFlL9oiWToENNTG2eyJyODubOmM4J3dPpap9EW6qJvrCP57b+nsc2/46Xtq5nIG8euYQCZkzoYFbLFKY3tzOtuZ0pTW0ckZ7IEekW2pITaEs105JsQmR4bedKVu54jtd2rqYv7OXI1E5OaunhfeldTE1kCfL+5IQGvaToo5WMOiCYRpDoIpmYSSo1gyCYSjIxnWTQiYLJBEMERTcmlQhE0+3UuDxJIcs3/NAD0TDUS2G8mfFyU4FtfwQmSWo3s/wUq+Hu50YgkUjwvWsu5Hvx6+WPvsR//s+zrB/YQ++EJJkJRrYVLGWEacNShiWgf4LRT4btZFjPHp7b9S6sWQOFJvsOAWslncqQasqSTGVJprOQCnlr7xbe2rulrL5aCGaCsIuJqX7Wp6ewMj2D1lQfnendzGvZyrHpnXQl+5iSyNASGBMZYCLbgG0Qrov6MgDsO/T4AwZZIGMiBLImsojQICRaF5qwQUsj2kbchnhd9JFQ+695D/6IePDzA38nh/s50or/rR2Xlq9dxDfOWVqdg9XBh/fxpl4CUVu83FtgW2+8nAjkB5Rh7SfpSuBKgNmzZ1faV5fn3LMWcO5ZCw5a98zza3nwiVW88mYPWwf2sTuVpb8FMs0QNoEljTAJliDK3QwMC4g+iwbxA9EfpujvTR3438RIJEMSiegRJEKCwAiCaCkZipdBAAqia0UkoJc0vf1piNNXhPF0IkNzYoB0IkNTkKUjuYdjJ2xnVtMupib30ZHsoy3I0BJkaQ5C0gpJYyQFSUFCkADS+8t2+h+pevRCe6HPqsNgQNbvT6u2eglEuTTyUr/Fhf73h7WfmS0FlkL01Vw5HXSVWbzoeBYvOr6strt272X9hh7e2bKDzT072bJrD9v29LKzbx97Mxl6wwH6MyF9liUTZsmakZFFSRNmoOgrtVAhZobFH/azWBzUQiyI2llgkHsu2CeLBjoCZLwlWDW4FvT+gYMRJIwgyKJEFOiSiQFaU320NfXTku6nOTHAhGSGVJAlnciSDrIkg5BkEJJQSEJGICOIn0sQEAdP2f4YTPx8/6ll+9/oKvFWV5k1rBtrLBTZtX1x9Q7mI6Kqq5dAlCt+M6HAtgl5baqxn6sjba0tzJ93FPNr3RHnhmIG2bFT/0rSnwPfB2YBa4GrzWxkdVlGQb3c0Jq7/f3IAttmANvNbE8V93POueEZIze0SjoG+CnwZeAI4IfAg5JKJnDVQl0EIjPbDqwHFhbYvAAomNY23P2cc254LLpGVOxRX44CbjOzJ8wsNLM7iS5VnFDjfh2iLgJR7H7gbEn7/5HiG1PnAnePwn7OOVcZA7Ow6KOemNljZvb13GtJi4FWoq/o6ko9BaIlRPcEPSrpa5L+AbiPaPqeOyAaakr6XDzkLHs/55yrmjoaEUm6WJIVePxXXrtjiT60f9fMth32jg6hXpIVMLMeSWcQfY95PVFK9s+Bb5hZbkbMM4AfA38FrKtgP+ecGzmzeqs79AAHbmMZbCD3RNLJwIPArWZ20+HqWCXqYmaFWpLUA7xZwS6dwHuj1B3nivH33eHznpmdW2hDe6LTFk84v+iOD+/5aV3NrCDpPGAZcI2Z3V7r/hRTNyOiWjGzispFSHqhnt5orjH4+65e1F92XDGSjgJ+BnzezO6vdX9KqadrRM45V9+M6D6iYo9hkHSbpMeLbBtJrbWvEs0s8xNJuwc9zhxWR0dRw4+InHOuXAZYWL0RkaQvAF8EniiwbQpRSYk08AMO1FqbL2nIWmtmdjVwddU6O4o8EFWuSjMnOlcRf9/VAzOsCjMrSEoA3wauK9Hsa0A3B9dae5ao1trlQMlaa2NJwycrOOdcuSQtJ0ocKaaZg+dpXxrPbTn4GM3As8B8opkPzgJeN7Mz89qNqNbaWOIjIuecK1OxbLoKNQOTgEvM7B5JG/IbxLXWjiG6JzLfCqB46t4Y5IHIOecOr53AcWaWKdGmoWqtedZcCZL+RFJvsSwTSX8t6bW4zWpJl46knXP5Rpg15epQPO9bqSAE5ddaGxc8EBUh6Qiiu5abi2z/OnArsIooM2UTsEzSJcNp51y+QVlTi4mypm4G/gx4RFKBguxuHBlujbYxyb+aK0DS+4D/psgstXGQug64y8z+Ml53G/A4cJOk+8wsW2670f1p3BjWMFlT7hANVWvNR0R5JF0O/BaYAhSbEuMComHxj3IrLJp691aiAlSnVtjOuUIuBR7PBSEAM/s/YE28zY1fDVVrzQPRoeYTlY/4IPB0kTa5qVZW5K1fkbe93HbOHWRQ1tSLBTavwN8741qj1VrzQHSob5nZ5Wa2uUSbmcA2M8u/kPjHeDm7wnbO5Ssra+ow9scdfg1Ta60hrhFJ6hqiyW4z2w0w1LQZsTbKy2Ypt51z+crNmhoX6buuoCXAZUS11m4mSpz6JuOw1lpDBCIOjECKuQH4TgXHCyiczZJbF1bYzrl8DZU15Q7VSLXWGiUQXTHE9pcqPN4uCmeztAzaXkk75/I1VNZUIzOzOSW2rQHOO3y9qY2GCESjUBBqIzBZUlPeJ5MZ8XJThe2cy9dQWVOusXmywvCsAAR8OG/9gnj5fIXtnDtIo2VNucbmgWh4fkk0w+7f5lZICoAvE5Udf6bCds4V0jBZU66xNcRXc9VmZlsk3QhcFweWXwMXAacDn83NllBuO+eKaJisKdfYPBAN3/XAHuArwIXAWuAzZpY/bXu57Zw7SCNlTbnG5oXxnHPO1ZRfI3LOOVdTHoicc87VlAci55xzNeWByDnnXE15IHLOOVdTHoicc87VlAci55xzNeWByDnnXE15IHLOOVdTHohcw5M0U9J7ko4e4XEek5QtVRFYUoekfkn3x69vj+eRc65heSByNSfpekkm6ZQSbZKSXpO0W1J3lbtwC3C3ma0f4XHuIPqd+osSbS4CUhyYtPR64CpJ80d4bufGLA9Erh6sipfzSrT5ClH5g38xsz9U68TxpKKfBn5QhcPdB/QBF5docwmwjahECGa2EVhGNLGpcw3JA5GrB6vjZcFAJGky8I9EheKq/TXWV4GnzOytkR7IzHYADwFnSJqWvz1e9zHgXjPrH7TpLuDjPipyjcoDkasHrxMVECw2IroOmAxcY2b7qnVSSbOAC4hKKwxev0HSf0j6oqS1knolPS/pJEldku6RtEvSJkk3xLWmcu4EEkQlP/JdHG/LryX0JLAV+Jtq/WzOjSVeBsLVBUkvAt1mNj1v/VzgZeAJMzu7yue8AlgKHGtmbwxav4EoYCSIrh8FwLeB7cCOuD+5IoefBD5vZj+J900D7wAr8vsr6UlgNnC05f3iSboTOMPMZlXzZ3RuLPARkasXq4Bpkqbmrf+3ePl3o3DOjxIVLVxXYNtM4BwzW2JmNwI/ArqBl83sUjNbSpSU0E8UjACIv3K7DzhTUmduvaSZ8fnuzA9CsVVA90gz95wbizwQuXpxyHUiSZ8AzgNuNbNXRuGcxwAbigSGN8xs9aDXa+PlA7kVZrYHeBc4Mm/f3Ndznx607rOAKF7iOxcMPRC5huOByNWLgwKRpATRaOg9omtEo2EKsLPIts15rzPx8t289VkO/T16EtgIfGbQukuJvq57tcj5cv3oLLLduXHLA5GrF/kjoivj598xs22jdM6Q4r8DmSLrh7yoGo+wlhFlwk2WNAc4ieKjIQb1IzvU8Z0bbzwQubpgZu8APcA8Se3APwG/BW7LbyvpAkmvSPqupD9I2hZnuH1E0tNxRtujklqGOO1molHRaLgDSBJl5V1MFGCWlWif60f+SMy5cS9Z6w44N8hq4CNE9wxNBS4ys7BAu4VE13fWA3OALwE3AY8QZbL1E138/1PgnhLnexNYLClhZlUdiZjZy5JWAecD04FH42BbTG62iI3V7IdzY4GPiFw9WQ20E2XI3W1mTxVptxBYamZ3mFkGWAk0A1eZ2TtmthV4i6E/aP0aaKH0jA4jcQdwLnAapb+WA1gMvB7PtOBcQ/FA5OpJ7jpRH/DNEu0WAA8Oev0h4DdxAEKSgPcDvxvifL8iuk50+rB6O7RlwESin+eBYo3iG2JPAf53lPrhXF3zG1rdmBLfm9MDTM4lMUi6Hegxs2/Fr48luum0zcwGhjjeA8BUM/vo6Pa8ZB8+ATwMfNjMVtaqH87Vio+I3FizkOjen215617Me716qCAU+1fgtDh41cplwCMehFyj8kDkxpoFwIrcC0kp4AMcHIgOalOKmT0N/AK4top9LFs8k8LFRFMIOdeQ/Ks51/DiyU9XAosGzzl3mM79Y2CHmV19OM/rXD3xQOScc66m/Ks555xzNeWByDnnXE15IHLOOVdTHoicc87VlAci55xzNeWByDnnXE15IHLOOVdT/w83TF3Nt+O7jAAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "gHrel = O + 2 * OL\n", "fig = plotIhKinetics(Vm, CCa, gHrel, ylabel='$(O + 2O_L)_{\\infty}$')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "**Ih steady-state activation increases with intracellular Calcium concentration**." ] } ], "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.7.3" } }, "nbformat": 4, "nbformat_minor": 4 }