{ "cells": [ { "cell_type": "code", "execution_count": 1, "id": "70a888c1-7d89-465f-b571-e157ef1ff440", "metadata": {}, "outputs": [], "source": [ "from neuron import h" ] }, { "cell_type": "code", "execution_count": 2, "id": "6dbc08fa-7c57-42a7-a38e-f01c6d9b71a8", "metadata": {}, "outputs": [], "source": [ "from neuron.units import ms, mV, µm\n", "from matplotlib import pyplot as plt" ] }, { "cell_type": "code", "execution_count": 3, "id": "394b9400-2fbb-440e-84de-689bcfff061b", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "1.0" ] }, "execution_count": 3, "metadata": {}, "output_type": "execute_result" } ], "source": [ "h.load_file(\"stdrun.hoc\")" ] }, { "cell_type": "code", "execution_count": null, "id": "8edaf9c3-5530-490f-a106-84a4287a6a36", "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": 4, "id": "74b300e2-a972-4e6d-ace9-f1a3cd8c7076", "metadata": {}, "outputs": [], "source": [ "IAF = h.IntFire1()\n", "IAF.tau = 20 # time constant\n", "IAF.refrac = 5 # refractory period" ] }, { "cell_type": "code", "execution_count": 5, "id": "c29f107e-cd7e-42c4-b936-9a3de81e65d4", "metadata": {}, "outputs": [], "source": [ "stim = h.NetStim()\n", "stim.number = 25\n", "stim.interval = 2\n", "stim.start = 0\n", "stim.noise = 1 # 0 - deterministic, 1 - Poisson process\n", "ncstim = h.NetCon(stim, IAF, 0, 0, 0.3)" ] }, { "cell_type": "code", "execution_count": 6, "id": "f647b209-598b-4935-a4e9-e654e7396811", "metadata": {}, "outputs": [], "source": [ "mvec = h.Vector().record(IAF._ref_m)" ] }, { "cell_type": "code", "execution_count": 7, "id": "d12ec6d4-4671-4f1b-b510-15ad7d056767", "metadata": {}, "outputs": [], "source": [ "#modify advance() procedure in stdrun.hoc\n", "h(\"\"\"\n", "proc advance() {\n", " fadvance()\n", " m1 = IntFire1[0].M()\n", " }\"\"\")\n", "Mvec = h.Vector().record(h._ref_m1) # M vector" ] }, { "cell_type": "code", "execution_count": 8, "id": "646a4e95-c29d-4a22-b587-83c25df022ea", "metadata": {}, "outputs": [], "source": [ "spike_times_vec_IAF = h.Vector()\n", "idvec_IAF = h.Vector()\n", "nc = h.NetCon(IAF, None)\n", "nc.record(spike_times_vec_IAF, idvec_IAF, 1)\n", "\n", "spike_times_vec_input = h.Vector()\n", "idvec_input = h.Vector()\n", "nc = h.NetCon(stim, None)\n", "nc.record(spike_times_vec_input, idvec_input, -1)\n", "h.finitialize()\n", "\n", "soma = h.Section(name = \"soma\")\n", "t = h.Vector().record(h._ref_t) # Time stamp vector" ] }, { "cell_type": "code", "execution_count": 9, "id": "cbec7100-3841-4e80-aaa6-b19747ab1ce6", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "0.0" ] }, "execution_count": 9, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# simulation control\n", "tstop = 100\n", "h.dt = 0.025\n", "h.finitialize()\n", "h.continuerun(tstop * ms)\n" ] }, { "cell_type": "code", "execution_count": null, "id": "02e86aa5-52a4-42ff-a04a-1dab79490795", "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "id": "99ae8b80-bb73-49c7-bee4-c334dffd1ce7", "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": 10, "id": "0b5e60f6-20b7-4608-82c7-c3ad56754d4f", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAiwAAAGwCAYAAACKOz5MAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8hTgPZAAAACXBIWXMAAA9hAAAPYQGoP6dpAABDnUlEQVR4nO3de1xUZcIH8N8wyIAig4rMQKJSmpYpmK6EWtlGIuvHV3rfzFzzwpv2yZU2o3KjTeziLl02tVo33twMfTfT2ox2y8V8SXTdUPJCZhdXC8ULA16C4aKgcN4/jGMTtznDOXNuv+/nMx/lzJkzz3Pm8vzmOc9zjkUQBAFEREREGhagdgGIiIiIOsPAQkRERJrHwEJERESax8BCREREmsfAQkRERJrHwEJERESax8BCREREmheodgHk0NzcjFOnTqFnz56wWCxqF4eIiIi8IAgCampqEB0djYCAjvtQDBFYTp06hZiYGLWLQURERD44fvw4+vXr1+E6hggsPXv2BHC5wmFhYSqXhoiIiLzhdrsRExMjtuMdMURgaTkMFBYWxsBCRESkM94M5+CgWyIiItI8BhYiIiLSPAYWIiIi0jxDjGEhIiLSi6amJly8eFHtYvhNt27dYLVau7wdBhYiIiI/EAQBLpcLVVVVahfF78LDw+F0Ort0rjQGFiIiIj9oCSuRkZHo3r27KU50KggC6uvrUVlZCQCIioryeVsMLERERApramoSw0qfPn3ULo5fhYSEAAAqKysRGRnp8+EhDrolIiJSWMuYle7du6tcEnW01LsrY3cYWIiIiPzEDIeB2iJHvRlYiIiISPMkBZbs7Gz87Gc/Q8+ePREZGYnU1FQcOnSo08e9++67GDp0KIKDgzF8+HBs3rzZ435BEJCVlYWoqCiEhIQgKSkJhw8fllYTIiIiMixJgWX79u1YuHAhdu3aha1bt+LixYuYOHEi6urq2n3Mp59+ihkzZuC+++7D/v37kZqaitTUVBw8eFBc54UXXsArr7yCnJwc7N69Gz169EBycjIuXLjge82IiIjIMCyCIAi+Pvj06dOIjIzE9u3bccstt7S5zvTp01FXV4cPP/xQXHbTTTchPj4eOTk5EAQB0dHReOSRR/Doo48CAKqrq+FwOJCbm4t77rmn03K43W7Y7XZUV1dr+uKH1fUXUdOgzMmC+va0wRbY9RPzyKm+8RLO1TUqsu3w7kEItWlrktuZ2gZcuNikyLaj7SEICFDm2HeF+wIuNjXLvl1rgAXOsGDTHrMn+rELFy6gtLQUsbGxCA4OVrs4ftde/aW03136xq+urgYA9O7du911ioqKkJGR4bEsOTkZeXl5AIDS0lK4XC4kJSWJ99vtdiQkJKCoqKjNwNLQ0ICGhgbxb7fb3ZVq+EVx6Tn8cvUuXGr2OR92KKZ3CD55ZAK6WbUxLOlsbQMmvFiImoZLimw/pJsVHz98C2J6a2PE/YbiMjy+6QvFtn/rtX2x9r/HyL7d5R8fwiufHJF9uy1+mdAfv79zuGLbJyLz8DmwNDc3Y9GiRRg3bhxuuOGGdtdzuVxwOBweyxwOB1wul3h/y7L21vmp7OxsPP30074WXRVfnqrGpWYBARbIHioaLjXj+Lnz+L6+EZE9tZHcS8/UiWHFFih/fc9fbMKRylrNBJbPT1wO79YACwJl7AkRBKCxqRkHTlTJts0fK/mh3IEBFlhlLHezIOBik6BYuYmMQBAEnFeoV7YjId2skno+J0yYgOHDh8NqtWLt2rUICgrCsmXL8Mtf/hLp6en461//CofDgVdffRUpKSmKldvnwLJw4UIcPHgQO3fulLM8XsnMzPTotXG73YiJifF7OXzxi+FR+OMvb5R1m7GZH8H3A3vKio3ogW2PTpB1m//xx5048ENDqzWLbh+MB28fLNv2jlTWIGn5Dtm2154X7hqB/7yxn2zbKzxUiblvfibb9oiM6PzFJlyftcXvz/vVM8noHiSt+V+7di0WL16M4uJibNy4EQsWLMD777+PO++8E0888QRWrFiBWbNmoaysTLFzzfj00zc9PR0ffvghtm3bhn79Ov6SczqdqKio8FhWUVEBp9Mp3t+yrL11fspmsyEsLMzjRkRERMqIi4vDk08+icGDByMzMxPBwcGIiIjA/PnzMXjwYGRlZeHs2bM4cOCAYmWQFLEEQcCDDz6I999/H4WFhYiNje30MYmJiSgoKMCiRYvEZVu3bkViYiIAIDY2Fk6nEwUFBYiPjwdwucdk9+7dWLBggZTiERER6UZINyu+eiZZleeVasSIEeL/rVYr+vTpg+HDr4xPaxnW0XLNICVICiwLFy7E+vXr8cEHH6Bnz57iGBO73S5eK2D27Nm46qqrkJ2dDQB46KGHcOutt+Kll17C5MmTsWHDBuzZswevv/46gMtnv1u0aBGWLVuGwYMHIzY2FkuWLEF0dDRSU1NlrCoREZF2WCwWyYdm1NKtWzePvy0Wi8eyljExzc3yzzhsIWlPvfbaawAuD8D5sTfffBNz584FAJSVlSEg4MqRprFjx2L9+vV48skn8cQTT2Dw4MHIy8vzGKi7ePFi1NXV4f7770dVVRXGjx+P/Px8U079IiIiotYkHxLqTGFhYatl06ZNw7Rp09p9jMViwTPPPINnnnlGSnF0RauDYpVisupCrzXuwmmYiIj8Shsn7SAiIiLqgD4OnhmI2c76aa7aAnp9efVabiJSXltHTo4ePdpqmdI9tuxhISIiIs1jYCEiIiLNY2AxEpONnxRMVGG91pRjeolILgwsREREfmLWmXly1JuBxU+UfItqcbykkp9Jc9VXi7XtnNkGlxN1puUka/X19SqXRB0t9f7pCeik4CwhIiIihVmtVoSHh4unru/evbspgr0gCKivr0dlZSXCw8NhtUq/LEALBhY/M/7b8ydMVmG9fgFZzPZCEamg5YK+Sl5vR6vCw8PbvaCxtxhYiIiI/MBisSAqKgqRkZG4ePGi2sXxm27dunWpZ6UFAwsREZEfWa1WWRpws+GgWyIiItI8BhZShFmn7ukNXyYi0gsGFj9hA25sfHmJiJTFwOJnOp1E4jOTVVe3zPa+JCL9YWAhIiIizWNgMRCzHZUw02EYvdZVr+UmIu1hYCEiIiLNY2AxAL2eXdVnJqqvXquq02ITkYYxsJAieCRAHwS+UkSkEwwsRDJgw09EpCwGFj8zW1e56Q5XERGRIhhYiGTEfEZEpAwGFiIiItI8BhYiIiLSPAYWUgRPGKYPfJ2ISC8YWIiIiEjzGFj8xKy/ZJUcg6qlfap0WZS+2rdSs7k09BIRkc4xsPiZktN8tdSAExERyYmBxQA4k1Y7LHw1AHB6NxHJj4GFdMdMbaGZ6kpE1BHJgWXHjh2YMmUKoqOjYbFYkJeX1+H6c+fOhcViaXUbNmyYuM5TTz3V6v6hQ4dKrgwREREZk+TAUldXh7i4OKxatcqr9V9++WWUl5eLt+PHj6N3796YNm2ax3rDhg3zWG/nzp1Si0Yawmvr6APHPRGRXgRKfUBKSgpSUlK8Xt9ut8Nut4t/5+Xl4fvvv0daWppnQQID4XQ6pRaHiIiITMDvY1jeeOMNJCUlYcCAAR7LDx8+jOjoaFx99dWYOXMmysrK2t1GQ0MD3G63x03rWnoczDYmwSyDL/XeUWGSl4mIdMyvgeXUqVP4xz/+gXnz5nksT0hIQG5uLvLz8/Haa6+htLQUN998M2pqatrcTnZ2tthzY7fbERMT44/iExERkUr8GljWrl2L8PBwpKameixPSUnBtGnTMGLECCQnJ2Pz5s2oqqrCO++80+Z2MjMzUV1dLd6OHz/uh9ITdc4sPUpERP4meQyLrwRBwJo1azBr1iwEBQV1uG54eDiuvfZaHDlypM37bTYbbDabEsUkIiIiDfJbD8v27dtx5MgR3HfffZ2uW1tbi2+//RZRUVF+KBmReXE2FxHpheTAUltbi5KSEpSUlAAASktLUVJSIg6SzczMxOzZs1s97o033kBCQgJuuOGGVvc9+uij2L59O44ePYpPP/0Ud955J6xWK2bMmCG1eKQVfmgHzdTU6rWuSl8DiYjMQ/IhoT179uC2224T/87IyAAAzJkzB7m5uSgvL281w6e6uhrvvfceXn755Ta3eeLECcyYMQNnz55F3759MX78eOzatQt9+/aVWjxT469l8hXH3hCR1kkOLBMmTOjwV1Nubm6rZXa7HfX19e0+ZsOGDVKLoTviLlOgYdByY2OWa+uwI8GTWV53IvIfXkuIdEfLAU1uSl7dm4hITxhYiGTEeEFEpAwGFiIiItI8BhZSBId06APH3hCRXjCwEBERkeYxsJCiOGZUHzirh4i0joHFT67MambDYEQ8Bw4RkbIYWIiIiEjzGFhIt7R42nfFDoFpr6pERH7FwEJkYsxBRKQXDCwG0DIuRksdDloqC/kfB1sTkdwYWIiIAYOINI+BxU9aehzYMHSdmXahmepKRNQRBhYiOfAQGBGRohhYiIiISPMYWIhkxBMDEhEpg4GFFMEzv+oEXyYi0gkGFiIiItI8BhYi4oEsItI8BhY/aTlEYraGwcJ53EREJAMGFtItLQ2/ULosWqqrFDzjMRHJhYGFiIiINI+BxQh+OOrCH7Pq09sRMKVmc+lsNxCRDjCwkCJ4KICIiOTEwEK6Y6aBvP6qqol2KRHpFAMLKYrtIBERyYGBxU94tWYiIiLfMbAQyUDgoB0iIkUxsBAREZHmMbAQmRg7hohILxhYSBFsB4mISE6SA8uOHTswZcoUREdHw2KxIC8vr8P1CwsLYbFYWt1cLpfHeqtWrcLAgQMRHByMhIQEFBcXSy0amYyZegeUHyOjzGhwpU5MR0TmIzmw1NXVIS4uDqtWrZL0uEOHDqG8vFy8RUZGivdt3LgRGRkZWLp0Kfbt24e4uDgkJyejsrJSavE0z2Kyib6cFUVERHIIlPqAlJQUpKSkSH6iyMhIhIeHt3nf8uXLMX/+fKSlpQEAcnJy8NFHH2HNmjV4/PHHJT+X2TATkObwTUlEMvPbGJb4+HhERUXhjjvuwL/+9S9xeWNjI/bu3YukpKQrhQoIQFJSEoqKitrcVkNDA9xut8eNOLVWTdzzRETKUjywREVFIScnB++99x7ee+89xMTEYMKECdi3bx8A4MyZM2hqaoLD4fB4nMPhaDXOpUV2djbsdrt4i4mJUboaRIbEoEVEeiH5kJBUQ4YMwZAhQ8S/x44di2+//RYrVqzA//7v//q0zczMTGRkZIh/u91uhhaNUbK3R8tHG+S+zpHZxjwREbVH8cDSljFjxmDnzp0AgIiICFitVlRUVHisU1FRAafT2ebjbTYbbDab4uUkMgsOjiYirVPlPCwlJSWIiooCAAQFBWHUqFEoKCgQ729ubkZBQQESExPVKB7JiA0hERHJQXIPS21tLY4cOSL+XVpaipKSEvTu3Rv9+/dHZmYmTp48iXXr1gEAVq5cidjYWAwbNgwXLlzAn//8Z3zyySf4+OOPxW1kZGRgzpw5GD16NMaMGYOVK1eirq5OnDVkJGzAiYiIpJMcWPbs2YPbbrtN/LtlLMmcOXOQm5uL8vJylJWVifc3NjbikUcewcmTJ9G9e3eMGDEC//d//+exjenTp+P06dPIysqCy+VCfHw88vPzWw3E1TPO4CEiIvKd5MAyYcKEDhvf3Nxcj78XL16MxYsXd7rd9PR0pKenSy0OkSYwjxIRKYvXEiIyMfb8EZFeMLCQIvzTDGqvsVVqiJL2auod5iEikgsDCxHxbC9EpHkMLAag5ZlHPPGZOfF1JyK5MbD4mZLhgt3vRERkVAwsfsIwQURE5DsGFtIdLR4CUyqPKl1X5mgi0gsGFlIGW0IiIpIRAwuRjLTY++MNua8yTUQkNwYWIiIi0jwGFr8z1y9Z/nAnIiI5MLAQERGR5jGw+AnHoMrPTFPF9VpXnRabiDSIgYVIBnq9iKBOi01EJsTAYgBaPA26wN/WpsaxS0QkNwYWIhnptZ3Wa7mJyDwYWEhRbAiJiEgODCx+xq5yIiIi6RhYSHe0OGaHiIiUxcDiJ5yNQVrEtyUR6QUDC5EM2PATESmLgYUUwR4lIiKSEwMLkYz0etVjnRabiEyEgYWUxZaQiIhkwMDiZ2y+5WOmo056PXOwXi9ZQETaw8BCREREmsfA4idK/kLmURfymUI9IHxLEpHcGFhIEaY7EmC2+hIR+RkDi4GYLiSQbNhLR0Rax8BCijJbO8iGn4hIGQwsfsYGTQYm2od8vxARXSY5sOzYsQNTpkxBdHQ0LBYL8vLyOlx/06ZNuOOOO9C3b1+EhYUhMTERW7Zs8VjnqaeegsVi8bgNHTpUatGIiIjIoCQHlrq6OsTFxWHVqlVerb9jxw7ccccd2Lx5M/bu3YvbbrsNU6ZMwf79+z3WGzZsGMrLy8Xbzp07pRZN0zi+hIiIyHeBUh+QkpKClJQUr9dfuXKlx9+///3v8cEHH+Dvf/87Ro4ceaUggYFwOp1Si0NEXcAcTUR64fcxLM3NzaipqUHv3r09lh8+fBjR0dG4+uqrMXPmTJSVlbW7jYaGBrjdbo8baYvZGkK9nomWiEgv/B5Y/vCHP6C2thZ33323uCwhIQG5ubnIz8/Ha6+9htLSUtx8882oqalpcxvZ2dmw2+3iLSYmxl/FJw0x02E2petqUWgks4leIiJSmF8Dy/r16/H000/jnXfeQWRkpLg8JSUF06ZNw4gRI5CcnIzNmzejqqoK77zzTpvbyczMRHV1tXg7fvy4v6pAEpltlovJqktE5DeSx7D4asOGDZg3bx7effddJCUldbhueHg4rr32Whw5cqTN+202G2w2mxLFVJxSv2SJiIiMzC89LG+//TbS0tLw9ttvY/LkyZ2uX1tbi2+//RZRUVF+KJ3+MQKR1ljM1rVGRIqT3MNSW1vr0fNRWlqKkpIS9O7dG/3790dmZiZOnjyJdevWAbh8GGjOnDl4+eWXkZCQAJfLBQAICQmB3W4HADz66KOYMmUKBgwYgFOnTmHp0qWwWq2YMWOGHHXUBB7LJy0y0zggItI3yT0se/bswciRI8UpyRkZGRg5ciSysrIAAOXl5R4zfF5//XVcunQJCxcuRFRUlHh76KGHxHVOnDiBGTNmYMiQIbj77rvRp08f7Nq1C3379u1q/UyFM1WIiMioJPewTJgwAUIHP8tyc3M9/i4sLOx0mxs2bJBaDNK4jt4jRqT76vIIDhFpHK8lRLpjpraVY0GIiC5jYCFFma65ZcAgIlIEA4ufsT0jIiKSjoGFiIiINI+BxV90PyqTjIgzy4hILxhYSBH+aAbN1Njqtqa6LTgRaQ0DC5EM9N6BxqFVRKR1DCwGoOWpr1ouGymHLzsRyY2Bxc/4PW5sfH2JiJTBwEJERESax8BiIHofR0H+x/cMEekFA4ufsF0gIiLyHQMLKULJX+5mGtDpr6pycDQRaR0DC5EMzHROGCIiNTCwkKL4u52IiOTAwOJn7Ho3Nr68RETKYGAhIiIizWNgId0y1ZRcheqq9D4000tERMpiYPETUzWuANhUERGRnBhYDIDDJqir5H4P8T1JRHJjYCFFmWUQqvl60IiI/IuBhYiIiDSPgcVA+CNffRYeDCEiUgQDC5GJMeQSkV4wsJDumKkXwyxjgIiIOsPA4idmu9YMB6HqC4MREWkdAwsRERFpHgMLKcosh2/YoUREpCwGFj9j17t8zBQS9HpIUeCxQSKSCQMLkYwYSImIlMHAYgRsJMlHSvWAMLgRkdwkB5YdO3ZgypQpiI6OhsViQV5eXqePKSwsxI033gibzYZBgwYhNze31TqrVq3CwIEDERwcjISEBBQXF0stGhERERmU5MBSV1eHuLg4rFq1yqv1S0tLMXnyZNx2220oKSnBokWLMG/ePGzZskVcZ+PGjcjIyMDSpUuxb98+xMXFITk5GZWVlVKLp1lmO5RvsurqnlkGRxORfgVKfUBKSgpSUlK8Xj8nJwexsbF46aWXAADXXXcddu7ciRUrViA5ORkAsHz5csyfPx9paWniYz766COsWbMGjz/+uNQiEhERkcEoPoalqKgISUlJHsuSk5NRVFQEAGhsbMTevXs91gkICEBSUpK4zk81NDTA7XZ73EijMzJM8sNdi7ueiMhIFA8sLpcLDofDY5nD4YDb7cb58+dx5swZNDU1tbmOy+Vqc5vZ2dmw2+3iLSYmRrHyy41d70RERNLpcpZQZmYmqqurxdvx48fVLhL5kZZnoMhdNAZcIqLLJI9hkcrpdKKiosJjWUVFBcLCwhASEgKr1Qqr1drmOk6ns81t2mw22Gw2xcpMRERE2qJ4D0tiYiIKCgo8lm3duhWJiYkAgKCgIIwaNcpjnebmZhQUFIjrkP5wTIe+aLnXiogI8CGw1NbWoqSkBCUlJQAuT1suKSlBWVkZgMuHa2bPni2u/8ADD+C7777D4sWL8c033+BPf/oT3nnnHTz88MPiOhkZGVi9ejXWrl2Lr7/+GgsWLEBdXZ04a8gI2H4TERH5TvIhoT179uC2224T/87IyAAAzJkzB7m5uSgvLxfDCwDExsbio48+wsMPP4yXX34Z/fr1w5///GdxSjMATJ8+HadPn0ZWVhZcLhfi4+ORn5/faiAu6Y+SP9w1OStKIXqtqk6LTUQaJDmwTJgwocOGoq2z2E6YMAH79+/vcLvp6elIT0+XWhzdYde7UbFpJiJSki5nCZEnZiDSHr4riUheDCxEMtJbD5peDzURkfkwsBAREZHmMbCQIgSO6dAVnXUMEZEJMbD4iT+63hkRiIjIqBhYSFF6G9NBRETaxMDiZ2y/u06LIUipHjQt1pWISA0MLEQmxrFGRKQXDCxEMuLVlYmIlMHAQorg+T3kpfjuVChn8X1ARHJhYCEiIiLNY2DxEyXHClg0PDKTh0jMScNvSSLSKQYWIiIi0jwGFj/jL09j4lANIiJlBapdACJSTm3DJbxacBinaxvavL+8+oKfS0RE5BsGFiI5aawHreDrCvzPju86XS88JMgPpSEi8h0Di4EoMYX0XF0j3tp1DHWNTZIed6SyVv7CkGQXLl5+3QZFhuLu0f3aXCemV3dcF9XTn8UiIpKMgUVBZ2sbcKn5coqoa7ikcml8879Fx7Di//7t8+N72KwyloZ8NbBPd9x/yzVqF4OIyGcMLAp5peAwlm/1vaHXitqGiwCAEf3sGDOwt6THWq0WTBvV9q/6rjDTVOmu1pQnbiMio2BgUci+su8BAAEWIOCHqUE9bIG45dq+ahZLspYGb+w1EXg8Zai6haEuME/IIyJjYmBRSEtD/8JdcbhLgV4Gf+EPdO8IGu3KaCkVp9MTkd7xPCzkFS02eBrNCIrwNRCpvY94NWgikgsDi0LEX7aqlqLr1G7w9Earr7dWy0VE5C0GFoX48xCBP3o/2ODpk1o9HHy/EJHcGFgUpsVDKVK0NHh6r4dZteRmvn5EpHcMLApjQ0FaYKap4ERkTAws1CHxFzobPF3iECQiMgpOa1YIG3pjEAQBT//9KxyurOlwvS9Puf1UIt+wp4+I9I6BRSHqjP1Q7ve0WRu8b0/XIvfTo16vHxkWrFxhfCFwDBIRGQMDC3WoZbaTWdu7xkuX698zOBDLUm/ocN2+oTbcdHUffxTLazwkRERGwcCiEJ6/RDn+7C1o6SkL6WbF1Pir/PfELWSqKw9NEpHe+TTodtWqVRg4cCCCg4ORkJCA4uLidtedMGECLBZLq9vkyZPFdebOndvq/kmTJvlSNM24Mp1U3w2FmLt0Xg9f6X1aMIMzERmF5B6WjRs3IiMjAzk5OUhISMDKlSuRnJyMQ4cOITIystX6mzZtQmNjo/j32bNnERcXh2nTpnmsN2nSJLz55pvi3zabTWrRSAFabvD8eVI0tXsofK2pcGX0tyq0/P4hIn2R3MOyfPlyzJ8/H2lpabj++uuRk5OD7t27Y82aNW2u37t3bzidTvG2detWdO/evVVgsdlsHuv16tXLtxpphDjoVuVyyMUo9TArvn5EpHeSAktjYyP27t2LpKSkKxsICEBSUhKKioq82sYbb7yBe+65Bz169PBYXlhYiMjISAwZMgQLFizA2bNn291GQ0MD3G63x83MlGyMzH7xOt0fElLpefV+KJSItEdSYDlz5gyamprgcDg8ljscDrhcrk4fX1xcjIMHD2LevHkeyydNmoR169ahoKAAzz//PLZv346UlBQ0NTW1uZ3s7GzY7XbxFhMTI6UafqH3hu6njFIPqYzSU8YAQUR659dZQm+88QaGDx+OMWPGeCy/5557xP8PHz4cI0aMwDXXXIPCwkLcfvvtrbaTmZmJjIwM8W+326250HLlas36bijMfgI8vQ+e5hgSIjIKST0sERERsFqtqKio8FheUVEBp9PZ4WPr6uqwYcMG3HfffZ0+z9VXX42IiAgcOXKkzfttNhvCwsI8bkTU2pXgTESkb5ICS1BQEEaNGoWCggJxWXNzMwoKCpCYmNjhY9999100NDTg3nvv7fR5Tpw4gbNnzyIqKkpK8bTFIIeExAZP5/XwlVHqr/fyExFJniWUkZGB1atXY+3atfj666+xYMEC1NXVIS0tDQAwe/ZsZGZmtnrcG2+8gdTUVPTp43km0NraWjz22GPYtWsXjh49ioKCAkydOhWDBg1CcnKyj9XSDr23E2Y/pCDofAfovfxERC0kj2GZPn06Tp8+jaysLLhcLsTHxyM/P18ciFtWVoaAAM8cdOjQIezcuRMff/xxq+1ZrVYcOHAAa9euRVVVFaKjozFx4kQ8++yzuj4Xixqza5Rsm/QevHxlmB4WtQtARNRFPg26TU9PR3p6epv3FRYWtlo2ZMiQdn/phYSEYMuWLb4UQ9OMM0vI3L/QzT7omIhIK3w6NT9JYYyGTv/By1fqXu24q0FJ77OciIhaMLAoxCj9EmzwLtNr7Y1yHhkiIgYWhbQcAtN7O6/lMZv+KJtWApuWX4eO6LXcRKQ9DCxEHdB7eyvwRCxEZBAMLArxZzuh5K9/Xkvo8r9s76Xh/iIiuTGwKEztQwlyMUg1JBN0nliMcokIIiIGFoUY5di9UerRVXpt7o0zvZ6IzI6BRSFGGzpg1l/oV04cZ876ExFpBQOLwvTezhnlTK++0nsPE6c1E5FRMLAoRe8t3Q8MUg2f6b3B5yEhIjIKBhaF+bOhUDJbmLa9Y4NPRKQJDCwKMcrsDNNPa/7hX72/jnovPxERA4tCjHYoRUs9DGoMgFXtWkIa2u9ERGpiYFGa3hscDQcvf56aX6+McokIIiIGFoUY7VCKWQ8p6P115KBbIjIKBhaF+PMEqUo+B6c1X/6X52GRhruLiOTGwKIwvTd0gt6PiXSR3k8AeOXV02sNiIguY2BRiMnbecPRa+7k+5CIjIKBRWE6bedEZm/vjDJoVe/lJyIKVLsARmW0sR96P7TlKz0Gtkr3Baz+53eobWjCwZPVAPQfnImIGFgUYpSxH0aoRlV9I/636BhqGy5JfmzZuXoA+pol9dbuMqz+Z6nHsrCQbiqVhohIHgwsCtNTQ9cWvQ86BYCNnx3HS1v/3aVthNr081G5cLEJADB6QC/ccm1fhHSz4r9G9VO5VEREXaOfb2GNq2u4hLVFR1FVfxEAUFnTAMDP1xJSsDdEz0eEWnpWhkWHYdygCMmPD7BYMDU+Wu5iKablbTBqQC/8+vbBqpbF6ARBwJv/OopjZ+vULgqR4qwBAciacr1qz8/AIpOPvijHC/mHWi3voaNf5m0xwqGtlir8bGBvPPGL69QtjB8Y4TXTi6/La/DMh1+pXQwivwgKZGAxhLoffsVf07cHkq5zAAD69QpBXD+7msXqsg8PlAPQ1iEhqWXR89lqu7TftfSiGVR94+XPvT2kG+YkDlC5NETKsgaoO7GYgUUmLT9qh0XbkWmQX/GVNRfE/0eGBatYkrbpN4b4RhAEr2ZrXTnLsvqJxei9PS2169MjCBkTh6haFiKj43lYZGLEr+WLTVdqlTzMqWJJusZs19Mx4ntRqwyex4g0hYFFZmo0iko/Z1BgAKwBJmntDUTNgKaF3h2/Mll1idTAwCIT8YyoKpdDTkbpzr8yNdtIr077DPKy6YIRP/dEWsXAQu3y5xWnlWTWBlzvr5semPStRaQKBhaZGfEU9kapklHq0ZmWWVFmqa8WGPFzT6Q1PgWWVatWYeDAgQgODkZCQgKKi4vbXTc3NxcWi8XjFhzsOeNEEARkZWUhKioKISEhSEpKwuHDh30pmmqM0hthRGIDrnI5/MWsPUpq4OeeyH8kB5aNGzciIyMDS5cuxb59+xAXF4fk5GRUVla2+5iwsDCUl5eLt2PHjnnc/8ILL+CVV15BTk4Odu/ejR49eiA5ORkXLlxoZ4vao+dzfbRHS9Nju8R4L41XdP+66YARP/dEWiU5sCxfvhzz589HWloarr/+euTk5KB79+5Ys2ZNu4+xWCxwOp3izeFwiPcJgoCVK1fiySefxNSpUzFixAisW7cOp06dQl5enk+VUpUB2wij9HYbpR6kPXxvESlPUmBpbGzE3r17kZSUdGUDAQFISkpCUVFRu4+rra3FgAEDEBMTg6lTp+LLL78U7ystLYXL5fLYpt1uR0JCQrvbbGhogNvt9ripTQu9EXL/2jPKr0dj1EI6NqJ+oIHPPZFZSAosZ86cQVNTk0cPCQA4HA64XK42HzNkyBCsWbMGH3zwAf7yl7+gubkZY8eOxYkTJwBAfJyUbWZnZ8Nut4u3mJgYKdVQhBEbRaMdnzfLwEhOtSUiI1J8llBiYiJmz56N+Ph43Hrrrdi0aRP69u2L//mf//F5m5mZmaiurhZvx48fl7HEXWPENlFrDb3U4ui5Afdl3xsxPGuVeI4fPb65iHRGUmCJiIiA1WpFRUWFx/KKigo4nd6dur1bt24YOXIkjhw5AgDi46Rs02azISwszOOmNiPOzNB6lbw9sZ1RXhvJ9dBAK2qQXd8uo7y3iPRAUmAJCgrCqFGjUFBQIC5rbm5GQUEBEhMTvdpGU1MTvvjiC0RFRQEAYmNj4XQ6Pbbpdruxe/dur7epBUacOqvnnok2GaYiHWMjSkRGJPlqzRkZGZgzZw5Gjx6NMWPGYOXKlairq0NaWhoAYPbs2bjqqquQnZ0NAHjmmWdw0003YdCgQaiqqsKLL76IY8eOYd68eQAud3kvWrQIy5Ytw+DBgxEbG4slS5YgOjoaqamp8tXUT9T5Uavwk+q8oTfdqfk1EJ410LnjF1dO0meSChOpSHJgmT59Ok6fPo2srCy4XC7Ex8cjPz9fHDRbVlaGgIArHTfff/895s+fD5fLhV69emHUqFH49NNPcf3114vrLF68GHV1dbj//vtRVVWF8ePHIz8/v9UJ5rTMiL9qDVglU2Ebqjwjfu6JtEpyYAGA9PR0pKent3lfYWGhx98rVqzAihUrOtyexWLBM888g2eeecaX4miKkX7FG2WWkFgPvVfES2xE/c8kby0iVfFaQjIzYqOo9+5uo5xPRiojhWet4iwhIv9hYJGJtzNW9MVYdTJLm2KsV03bjPm5J9ImBhaZGPGwg1HqZJR6eMts9dUC7msi5TGwyM5431zGqxGRPMw2A41ITQwsMtFCx7DcvdNaqJOczNOoqD+tmYhIbgwsMjFiN/yVOum7UmYdZ6Dzl00fDPi5J9IqBhbqlNa+i1vKo8Mz1Uv24yJ7W18t5TMtlUUJZp2BRqQGBhaZaOHsonIzypfxlXEG5mCUnjE94Z4mUh4Di0zUPCSk1HMa8TAXkZwEnoiFyG8YWMgL+v4yNlujYpSeMT0w+iEvIi1hYJGJEac3GuXL2GwNOHvG/I+7mkh5DCxyEVqu2qpyOWR05Uq0KhdEJgapBmmIyTrviFTFwEKd0vt3sdl6HIzY26dVZp0yT6QGBhaZGHEmilG+iw1SDdIwI33uibSKgUUmRp5KapQqmaXHwWw9Smq6ckiIO5tIaQws1Cm9N/RG6SnylhHPCURExMAiEy3MRJH9WkLqV0lW/BFMchN7s9QtBpEpMLDIxGiNO2CkWUIm63HgISE/MuAHn0ijGFhkZsRGQmtVEscLeNlW6DlM/vj95O2MFC1VVws9j/5gxM89kdYwsMjEiFNJ9dzQt8VsjYqR3otadeWQEPc1kdIYWGSiZuOu1FelUWZAGHkGV1sEDZzE0CS72iT9R0TawMAiM7N8URPRj/BzT6Q4BhaZGHEqqVHO4mmWcRQtzFVbdXGWEJH/MLDIxYCthFGuk2KQ3OU1s9VXTWYLw0RqYmCRmd4b97YYpU5GqYe3zDJmRwu4q4mUF6h2AbTuTG0DNn52HPWNlzpcb3fpOQDabSR2fXcW/zx8WtJjXNUNCpVGGWVn6/HevhO41NzssfyLk9UAjDuT46f1/rrcDYCHKbqi4VIT1u8uw5najj8DhytqARj3vUWkJQwsnXjzX6VYte1br9fvHmRVsDS++9Vb+3CurtGnx/YI0sfb5MWPD+Hvn59q9/4eNm2+Nl3VXr2NWl9/2PZNJZ7++1der899TaQ8fbREKqq5cLlnJS4mHDf2D+9w3VBbIGYmDPBDqaSruXARAHD36H7oYfP+ZbfAgskjopQqlqxa6njz4AgMigz1uM8e0g1T469So1iKa6vevboHYfKIaDWLpWvuHz73/XqF4I7rHR2uGxhgwbTRMf4oFpGpMbB0omUA463X9kXGHdeqW5hOdDQAsKUeGXcMgdMe7KcS+VdLHafGX4W7RvVTtzB+ZNZ6K+qHfXqtoyeWThmmblmICAAH3XaKswC0y2yvjR5ryxlLRCQXBhYv6X1InVGmKAPtvxZCJ/frkTeDOY1Yb7UZ8bxKRHrnU2BZtWoVBg4ciODgYCQkJKC4uLjddVevXo2bb74ZvXr1Qq9evZCUlNRq/blz58JisXjcJk2a5EvRZKeHX4jehBDxdO0Kl0VNWjglvRq0WG+9z5rRw+eeyGwkB5aNGzciIyMDS5cuxb59+xAXF4fk5GRUVla2uX5hYSFmzJiBbdu2oaioCDExMZg4cSJOnjzpsd6kSZNQXl4u3t5++23faqQQLTUGROQf/NwTaYfkwLJ8+XLMnz8faWlpuP7665GTk4Pu3btjzZo1ba7/1ltv4Ve/+hXi4+MxdOhQ/PnPf0ZzczMKCgo81rPZbHA6neKtV69evtVIZka5CrP4g1Hf1fCKWRsZs9ZbCVc6WLhTibRCUmBpbGzE3r17kZSUdGUDAQFISkpCUVGRV9uor6/HxYsX0bt3b4/lhYWFiIyMxJAhQ7BgwQKcPXu23W00NDTA7XZ73IiIiMi4JAWWM2fOoKmpCQ6H53kJHA4HXC6XV9v4zW9+g+joaI/QM2nSJKxbtw4FBQV4/vnnsX37dqSkpKCpqanNbWRnZ8Nut4u3mBjlzoEgXtxM5z+0rlykTecV6YAZ6tgWs9ZbSUb53BMZiV/Pw/Lcc89hw4YNKCwsRHDwlXOB3HPPPeL/hw8fjhEjRuCaa65BYWEhbr/99lbbyczMREZGhvi32+1WMLRw9B0REZHaJPWwREREwGq1oqKiwmN5RUUFnE5nh4/9wx/+gOeeew4ff/wxRowY0eG6V199NSIiInDkyJE277fZbAgLC/O4Kc0oP7SM/ItRnIpq4Dq2xaz1VhKnNRNpj6TAEhQUhFGjRnkMmG0ZQJuYmNju41544QU8++yzyM/Px+jRozt9nhMnTuDs2bOIilL/lPBGmN4oCOYYQmiE18oXZq23krhPibRH8iyhjIwMrF69GmvXrsXXX3+NBQsWoK6uDmlpaQCA2bNnIzMzU1z/+eefx5IlS7BmzRoMHDgQLpcLLpcLtbWXr3JaW1uLxx57DLt27cLRo0dRUFCAqVOnYtCgQUhOTpapml2nh1+v/JIlkpcePvdEZiF5DMv06dNx+vRpZGVlweVyIT4+Hvn5+eJA3LKyMgQEXMlBr732GhobG3HXXXd5bGfp0qV46qmnYLVaceDAAaxduxZVVVWIjo7GxIkT8eyzz8Jms3Wxel13ZfCdfr+5fhxk9FyPzhjhtfKFWeutJKOczoDISHwadJueno709PQ27yssLPT4++jRox1uKyQkBFu2bPGlGGRyZutR0mN9dVhkItIoXkuoE0a4wN6Pa2CE34vtdSQYcqCkN5ddMGK91abByx0QmR0DSyf08KtWSre1Gb6AzVDHtmip3loqCxEZAwOLl/T8BSzoIXXJwCTVbMWs9VaSka5uTmQUDCydMEJb4HlIyLjfwGYdKGnWeiuJIZBIexhYvMTGgMh8+Lkn0g4Glk4Y4ZoigtFG3bbHAK+VT8xabwUJVy7QREQawcBiMmZo1ExQxTaZtd5EZA4MLJ0wwpRRI0zN9oZZ6vlTZq23kq6MCyIirWBg6YzB2gIzfAGboRepLWatNxGZAwOLl/TcGJhlxoNZ6vlTZq23kni5AyLtYWDphNHaAiN9Af/0tTH6NanbO/Sj5XrrNUzptNhEhsbA4iVOb9QS87wWBsqXusTdT6QdDCyd0MNZYjtr1Dyu1qxsUVQlmPT6L1qst5bK4gs9fO6JzIaBxUt6/wJuYZR6dMQEVWyTWeutJDN8Xoj0goGlE0b4nWWWaa/mqGVrZq23PzCvEGkHA0snjNYzbIaxOEYaWCyFWetNRObAwOIlPTcGRgtd7TFLPX/KrPVWEqc1E2kPA0snjNYWmOH71wRVbJNZ601E5sDA4iU9NwZGC13tEU+nrucXywdmrbeSjHBJDiKjYWDphBGmNxqhDl4xSz1/yqz1VhB3KZH2MLB4ySi/Xo1Sj46YoY5tMWu9FcV9SqQZDCyd0NNVW9v7VWiWH4tmqedPmbXeSrryudfDJ5/IHBhYOmOw1sBIX8DtBTQj1fHHOjtMoc16G+wDRESqYWDxkp6nNxrteHx7L4XR6gl417NnxHqr7cq0ZnXLQURXMLB0Qg9niZXynWqKL2Az1LEtGqq3Nnt7iEjPGFi8pOuG3iwXPzTpVFSz1ltJ3KdE2sPA0gkjdLfroZdIDkZ4rXxh1norifuUSHsYWDohHstWtxiy0fNYHG+ZoY5tMWu9lcRdSqQdDCze0vE3l1l+LZqlnj9l1nr7A8fiEGkHA0snjHY4xQxfv2aoY1vMWm8iMgcGFi/puTH4ceTScUdRp8x6TR2z1ltJLZez4D4l0g6fAsuqVaswcOBABAcHIyEhAcXFxR2u/+6772Lo0KEIDg7G8OHDsXnzZo/7BUFAVlYWoqKiEBISgqSkJBw+fNiXosmO3e1ERETqkxxYNm7ciIyMDCxduhT79u1DXFwckpOTUVlZ2eb6n376KWbMmIH77rsP+/fvR2pqKlJTU3Hw4EFxnRdeeAGvvPIKcnJysHv3bvTo0QPJycm4cOGC7zWTmZ5/af344odGHpgp/irWdX+YdGatt5J44jgi7QmU+oDly5dj/vz5SEtLAwDk5OTgo48+wpo1a/D444+3Wv/ll1/GpEmT8NhjjwEAnn32WWzduhV//OMfkZOTA0EQsHLlSjz55JOYOnUqAGDdunVwOBzIy8vDPffc05X6dYkgCDjx/XnVnl+qCvcF9Pq+W6vlVfUXVSiN8r6vb8SJ7+vFvy82NatYGuWdrDoPW2Dr3xharvelZsHjNdIL9wVjfmaI9ExSYGlsbMTevXuRmZkpLgsICEBSUhKKiorafExRUREyMjI8liUnJyMvLw8AUFpaCpfLhaSkJPF+u92OhIQEFBUVtRlYGhoa0NDQIP7tdrulVMNrF5sEfFV+edt6+PU6b90etYvgVy9uOYQXtxxqtdyov4pvf2l7h/drsd5V9Rcx/vltahejCzS4U4lMSlJgOXPmDJqamuBwODyWOxwOfPPNN20+xuVytbm+y+US729Z1t46P5WdnY2nn35aStF9ZgsMQESoDWNie/vl+XwxJS4aa4uOdjreJnmY0z8FUljyMCd2fXcWjZda9yz06xWCG6LtKpRKGaG2QNw+NBI7j5zpcD2t1Ts2ogfi+tnxjatG7aL4LCTIip8PjVS7GET0A8mHhLQgMzPTo9fG7XYjJiZG9ucJCgzAoWUpsm9Xbpm/uA6Zv7hO7WL4zV2j+uGuUf3ULoZfWCwWvDH3Z2oXQ7KgwAB8kD5e7WIQkYFIGnQbEREBq9WKiooKj+UVFRVwOtv+9e50Ojtcv+VfKdu02WwICwvzuBEREZFxSQosQUFBGDVqFAoKCsRlzc3NKCgoQGJiYpuPSUxM9FgfALZu3SquHxsbC6fT6bGO2+3G7t27290mERERmYvkQ0IZGRmYM2cORo8ejTFjxmDlypWoq6sTZw3Nnj0bV111FbKzswEADz30EG699Va89NJLmDx5MjZs2IA9e/bg9ddfB3C5y3vRokVYtmwZBg8ejNjYWCxZsgTR0dFITU2Vr6ZERESkW5IDy/Tp03H69GlkZWXB5XIhPj4e+fn54qDZsrIyBARc6bgZO3Ys1q9fjyeffBJPPPEEBg8ejLy8PNxwww3iOosXL0ZdXR3uv/9+VFVVYfz48cjPz0dwcLAMVSQiIiK9swiC/s/l6na7YbfbUV1dzfEsREREOiGl/ea1hIiIiEjzGFiIiIhI8xhYiIiISPMYWIiIiEjzGFiIiIhI8xhYiIiISPMYWIiIiEjzGFiIiIhI8xhYiIiISPMkn5pfi1pO1ut2u1UuCREREXmrpd325qT7hggsNTU1AICYmBiVS0JERERS1dTUwG63d7iOIa4l1NzcjFOnTqFnz56wWCyybtvtdiMmJgbHjx/ndYoUxP3sH9zP/sN97R/cz/6h1H4WBAE1NTWIjo72uHByWwzRwxIQEIB+/fop+hxhYWH8MPgB97N/cD/7D/e1f3A/+4cS+7mznpUWHHRLREREmsfAQkRERJrHwNIJm82GpUuXwmazqV0UQ+N+9g/uZ//hvvYP7mf/0MJ+NsSgWyIiIjI29rAQERGR5jGwEBERkeYxsBAREZHmMbAQERGR5jGwdGLVqlUYOHAggoODkZCQgOLiYrWLpFvZ2dn42c9+hp49eyIyMhKpqak4dOiQxzoXLlzAwoUL0adPH4SGhuK//uu/UFFRoVKJjeG5556DxWLBokWLxGXcz/I5efIk7r33XvTp0wchISEYPnw49uzZI94vCAKysrIQFRWFkJAQJCUl4fDhwyqWWH+ampqwZMkSxMbGIiQkBNdccw2effZZj+vPcD9Lt2PHDkyZMgXR0dGwWCzIy8vzuN+bfXru3DnMnDkTYWFhCA8Px3333Yfa2lplCixQuzZs2CAEBQUJa9asEb788kth/vz5Qnh4uFBRUaF20XQpOTlZePPNN4WDBw8KJSUlwi9+8Quhf//+Qm1trbjOAw88IMTExAgFBQXCnj17hJtuukkYO3asiqXWt+LiYmHgwIHCiBEjhIceekhczv0sj3PnzgkDBgwQ5s6dK+zevVv47rvvhC1btghHjhwR13nuuecEu90u5OXlCZ9//rnwH//xH0JsbKxw/vx5FUuuL7/73e+EPn36CB9++KFQWloqvPvuu0JoaKjw8ssvi+twP0u3efNm4be//a2wadMmAYDw/vvve9zvzT6dNGmSEBcXJ+zatUv45z//KQwaNEiYMWOGIuVlYOnAmDFjhIULF4p/NzU1CdHR0UJ2draKpTKOyspKAYCwfft2QRAEoaqqSujWrZvw7rvviut8/fXXAgChqKhIrWLqVk1NjTB48GBh69atwq233ioGFu5n+fzmN78Rxo8f3+79zc3NgtPpFF588UVxWVVVlWCz2YS3337bH0U0hMmTJwv//d//7bHsP//zP4WZM2cKgsD9LIefBhZv9ulXX30lABA+++wzcZ1//OMfgsViEU6ePCl7GXlIqB2NjY3Yu3cvkpKSxGUBAQFISkpCUVGRiiUzjurqagBA7969AQB79+7FxYsXPfb50KFD0b9/f+5zHyxcuBCTJ0/22J8A97Oc/va3v2H06NGYNm0aIiMjMXLkSKxevVq8v7S0FC6Xy2Nf2+12JCQkcF9LMHbsWBQUFODf//43AODzzz/Hzp07kZKSAoD7WQne7NOioiKEh4dj9OjR4jpJSUkICAjA7t27ZS+TIS5+qIQzZ86gqakJDofDY7nD4cA333yjUqmMo7m5GYsWLcK4ceNwww03AABcLheCgoIQHh7usa7D4YDL5VKhlPq1YcMG7Nu3D5999lmr+7if5fPdd9/htddeQ0ZGBp544gl89tln+PWvf42goCDMmTNH3J9tfY9wX3vv8ccfh9vtxtChQ2G1WtHU1ITf/e53mDlzJgBwPyvAm33qcrkQGRnpcX9gYCB69+6tyH5nYCFVLFy4EAcPHsTOnTvVLorhHD9+HA899BC2bt2K4OBgtYtjaM3NzRg9ejR+//vfAwBGjhyJgwcPIicnB3PmzFG5dMbxzjvv4K233sL69esxbNgwlJSUYNGiRYiOjuZ+NhEeEmpHREQErFZrq5kTFRUVcDqdKpXKGNLT0/Hhhx9i27Zt6Nevn7jc6XSisbERVVVVHutzn0uzd+9eVFZW4sYbb0RgYCACAwOxfft2vPLKKwgMDITD4eB+lklUVBSuv/56j2XXXXcdysrKAEDcn/we6ZrHHnsMjz/+OO655x4MHz4cs2bNwsMPP4zs7GwA3M9K8GafOp1OVFZWetx/6dIlnDt3TpH9zsDSjqCgIIwaNQoFBQXisubmZhQUFCAxMVHFkumXIAhIT0/H+++/j08++QSxsbEe948aNQrdunXz2OeHDh1CWVkZ97kEt99+O7744guUlJSIt9GjR2PmzJni/7mf5TFu3LhWU/P//e9/Y8CAAQCA2NhYOJ1Oj33tdruxe/du7msJ6uvrERDg2VxZrVY0NzcD4H5Wgjf7NDExEVVVVdi7d6+4zieffILm5mYkJCTIXyjZh/EayIYNGwSbzSbk5uYKX331lXD//fcL4eHhgsvlUrtourRgwQLBbrcLhYWFQnl5uXirr68X13nggQeE/v37C5988omwZ88eITExUUhMTFSx1Mbw41lCgsD9LJfi4mIhMDBQ+N3vficcPnxYeOutt4Tu3bsLf/nLX8R1nnvuOSE8PFz44IMPhAMHDghTp07ldFuJ5syZI1x11VXitOZNmzYJERERwuLFi8V1uJ+lq6mpEfbv3y/s379fACAsX75c2L9/v3Ds2DFBELzbp5MmTRJGjhwp7N69W9i5c6cwePBgTmtWy6uvvir0799fCAoKEsaMGSPs2rVL7SLpFoA2b2+++aa4zvnz54Vf/epXQq9evYTu3bsLd955p1BeXq5eoQ3ip4GF+1k+f//734UbbrhBsNlswtChQ4XXX3/d4/7m5mZhyZIlgsPhEGw2m3D77bcLhw4dUqm0+uR2u4WHHnpI6N+/vxAcHCxcffXVwm9/+1uhoaFBXIf7Wbpt27a1+Z08Z84cQRC826dnz54VZsyYIYSGhgphYWFCWlqaUFNTo0h5LYLwo1MFEhEREWkQx7AQERGR5jGwEBERkeYxsBAREZHmMbAQERGR5jGwEBERkeYxsBAREZHmMbAQERGR5jGwEBERkeYxsBBRl8ydOxepqamqPf+sWbPEqyUr4auvvkK/fv1QV1en2HMQUed4plsiapfFYunw/qVLl+Lhhx+GIAgIDw/3T6F+5PPPP8fPf/5zHDt2DKGhoYo9z1133YW4uDgsWbJEsecgoo4xsBBRu1wul/j/jRs3Iisry+PqxKGhoYoGhc7MmzcPgYGByMnJUfR5PvroI8yfPx9lZWUIDAxU9LmIqG08JERE7XI6neLNbrfDYrF4LAsNDW11SGjChAl48MEHsWjRIvTq1QsOhwOrV69GXV0d0tLS0LNnTwwaNAj/+Mc/PJ7r4MGDSElJQWhoKBwOB2bNmoUzZ860W7ampib89a9/xZQpUzyWDxw4EMuWLcPs2bMRGhqKAQMG4G9/+xtOnz6NqVOnIjQ0FCNGjMCePXvExxw7dgxTpkxBr1690KNHDwwbNgybN28W77/jjjtw7tw5bN++vYt7lIh8xcBCRLJbu3YtIiIiUFxcjAcffBALFizAtGnTMHbsWOzbtw8TJ07ErFmzUF9fDwCoqqrCz3/+c4wcORJ79uxBfn4+KioqcPfdd7f7HAcOHEB1dTVGjx7d6r4VK1Zg3Lhx2L9/PyZPnoxZs2Zh9uzZuPfee7Fv3z5cc801mD17Nlo6mBcuXIiGhgbs2LEDX3zxBZ5//nmPnqOgoCDEx8fjn//8p8x7ioi8psg1oInIcN58803Bbre3Wj5nzhxh6tSp4t+33nqrMH78ePHvS5cuCT169BBmzZolLisvLxcACEVFRYIgCMKzzz4rTJw40WO7x48fFwC0upx9i/fff1+wWq1Cc3Ozx/IBAwYI9957b6vnWrJkibisqKhIACCUl5cLgiAIw4cPF5566qkO63/nnXcKc+fO7XAdIlIOe1iISHYjRowQ/2+1WtGnTx8MHz5cXOZwOAAAlZWVAC4Pnt22bZs4JiY0NBRDhw4FAHz77bdtPsf58+dhs9naHBj84+dvea6Onv/Xv/41li1bhnHjxmHp0qU4cOBAq22GhISIPUJE5H8MLEQku27dunn8bbFYPJa1hIzm5mYAQG1tLaZMmYKSkhKP2+HDh3HLLbe0+RwRERGor69HY2Njh8/f8lwdPf+8efPw3XffYdasWfjiiy8wevRovPrqqx7bPHfuHPr27evdDiAi2TGwEJHqbrzxRnz55ZcYOHAgBg0a5HHr0aNHm4+Jj48HcPk8KXKIiYnBAw88gE2bNuGRRx7B6tWrPe4/ePAgRo4cKctzEZF0DCxEpLqFCxfi3LlzmDFjBj777DN8++232LJlC9LS0tDU1NTmY/r27Ysbb7wRO3fu7PLzL1q0CFu2bEFpaSn27duHbdu24brrrhPvP3r0KE6ePImkpKQuPxcR+YaBhYhUFx0djX/9619oamrCxIkTMXz4cCxatAjh4eEICGj/a2revHl46623uvz8TU1NWLhwIa677jpMmjQJ1157Lf70pz+J97/99tuYOHEiBgwY0OXnIiLf8MRxRKRb58+fx5AhQ7Bx40YkJiYq8hyNjY0YPHgw1q9fj3HjxinyHETUOfawEJFuhYSEYN26dR2eYK6rysrK8MQTTzCsEKmMPSxERESkeexhISIiIs1jYCEiIiLNY2AhIiIizWNgISIiIs1jYCEiIiLNY2AhIiIizWNgISIiIs1jYCEiIiLNY2AhIiIizft/d1I5yVly0E0AAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAi8AAAGwCAYAAABhDIVPAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8hTgPZAAAACXBIWXMAAA9hAAAPYQGoP6dpAABLl0lEQVR4nO3dd3hUVf4/8PdMyqQnpBeSEGqoIYBgAAUVKSKCu8u6ijQFFxdUxF1XXMV11R+6+1VsKKKrWECwALIWEKmCoRM6kWYSQhqEZNJIm/v7AzIhkDbJ3HLufb+eJw9kcmfmc+5Mct5zzrn3miRJkkBEREQkCLPaBRARERE5guGFiIiIhMLwQkREREJheCEiIiKhMLwQERGRUBheiIiISCgML0RERCQUV7ULcDabzYZz587B19cXJpNJ7XKIiIioGSRJQlFRESIjI2E2Nz62orvwcu7cOURHR6tdBhEREbVARkYG2rZt2+g2ugsvvr6+AC433s/PT+VqiIiIqDmsViuio6Pt/XhjdBdeaqaK/Pz8GF6IiIgE05wlH1ywS0REREJheCEiIiKhMLwQERGRUHS35oWIiEgk1dXVqKysVLsMRbi7uzd5GHRzMLwQERGpQJIkZGdno6CgQO1SFGM2mxEXFwd3d/dWPQ7DCxERkQpqgktoaCi8vLx0f2LVmpPIZmVlISYmplXtZXghIiJSWHV1tT24BAUFqV2OYkJCQnDu3DlUVVXBzc2txY/DBbtEREQKq1nj4uXlpXIlyqqZLqqurm7V4zC8EBERqUTvU0XXclZ7GV6IiIhIKLKGl/nz5+OGG26Ar68vQkNDMW7cOKSmpjZ5vy+//BLx8fHw8PBAz5498f3338tZJhEREQlE1vCyZcsWzJw5Ezt27MD69etRWVmJ4cOHo6SkpMH7/PLLL7j33nvx4IMPYv/+/Rg3bhzGjRuHw4cPy1kqERERCULW8LJ27VpMmTIF3bt3R0JCApYsWYL09HTs3bu3wfu88cYbGDlyJP72t7+ha9eueOGFF9CnTx+8/fbbcpaquMLSSpy9WIqswjJIkqR2ObKzXrrc3swCY7Q3r6gcZy+WIq+oXO1Smk2SJJwrKMPZi6UoLDPGCbOIyHFTpkyByWTCjBkzrvvZzJkzYTKZMGXKFFlrUPRQ6cLCQgBAYGBgg9skJydjzpw5dW4bMWIEVq9eXe/25eXlKC+v7SCsVmvrC5VZ8qkLmPjfnaiyXe7EJ94YixfG9VC5Kvkcz7birre3o6LKBgAY3TMCCyf0Ubkq+SzdmYZ/rKodKXxhbHdMTGqnXkHN9MQXB7ByfyYAwM3FhK9mDERCdIC6RRGRJkVHR2P58uVYsGABPD09AQCXLl3CsmXLEBMTI/vzKxZebDYbZs+ejUGDBqFHj4Y76uzsbISFhdW5LSwsDNnZ2fVuP3/+fDz//PNOrVVuR84V2oMLABw8W6BeMQo4nlVkDy4AcEDn7T2YUVj3+7OFDWypLSlXvS6V1RKOZ1sZXogUJEkSyipbdwhxS3m6uTh0JFCfPn1w6tQprFy5EhMmTAAArFy5EjExMYiLi5OrTDvFwsvMmTNx+PBhbNu2zamPO3fu3DojNVarFdHR0U59Drn4ebjCeqlK7TIU4+vhiiK2V/NErZtIdGWV1eg2b50qz330XyPg5e5YJHjggQfw0Ucf2cPLhx9+iKlTp2Lz5s0yVFiXIodKz5o1C99++y02bdqEtm3bNrpteHg4cnJy6tyWk5OD8PDwere3WCzw8/Or80VERETyuv/++7Ft2zakpaUhLS0N27dvx/3336/Ic8s68iJJEh555BGsWrUKmzdvbtZQUlJSEjZs2IDZs2fbb1u/fj2SkpJkrJSIiEhdnm4uOPqvEao9t6NCQkIwevRoLFmyBJIkYfTo0QgODpahuuvJGl5mzpyJZcuW4ZtvvoGvr6993Yq/v799gc+kSZMQFRWF+fPnAwAee+wxDBkyBK+++ipGjx6N5cuXY8+ePVi8eLGcpZKMJOj/6KKrCdteQcsm0guTyeTw1I3aHnjgAcyaNQsAsHDhQsWeV9Zpo3fffReFhYUYOnQoIiIi7F8rVqywb5Oeno6srCz79wMHDsSyZcuwePFiJCQk4KuvvsLq1asbXeRLREREyhs5ciQqKipQWVmJESOUGzWSfdqoKfUt7Bk/fjzGjx8vQ0XaYrhrWqhdgMJEba+odROR8lxcXHDs2DH7/5Ui1vgUERERaYoaB8owvBAREVGzLVmypNGfN3RSWWfiVaVVYICz49fB9opB0LKJyIAYXoiIiEgoDC8qqlmva5RPvDULlEUdmXCUvb0q1+Eoo71ORCQehhciIiKVNOeoXD1xVnsZXlRQcxIzox2SarAjw4Vtr6h1E4nEzc0NAFBaWqpyJcqqqKgA0PrDqnm0EcnOYB8shJsmqmG0T4BEanJxcUFAQAByc3MBAF5eXro/95fNZkNeXh68vLzg6tq6+MHwQkR16PvPJ5F21FxwuCbAGIHZbEZMTEyrgxrDi4r0nrKvZazWGq+9ROQYk8mEiIgIhIaGorKyUu1yFOHu7g6zufUrVhheiIiIVOTi4qLoqfX1gAt2VcClBaRFfFsSkSgYXkh2RusUGU6JiOTF8KIirokgLTLaWiwiEg/DCynGaJ2i0dpLRKQUhhciIiISCsOLCq5dEsE1Evomyut7bZ2ClE1EBsTwoiKjzCoY7cytkuDdvkHelkQkMIYXVRmrmzBWa8Vtr1FCNRGJi+GFFGO0TtFo7SUiUgrDCxEREQmF4UUFBlsCIvgKEOMQfa0OERkHwwuRswmfATjfRUTaxvCiIq6JICIichzDCynIaGnNaO0lIlIGw4sKuLaAtMhoa7GISFwMLyQ/dopEROREDC8aYLSRGL2fcfe6yz8I9vpyLRYRaR3Di4qM1kcYrVMUvb06z5hEJDCGF1KM4H25w0Rtr6h1E5FxMLyooOYTreifzElfONJCRKJgeCHZibbmw+gYqolI6xheVGTiAL0u6X1BMhGR2hheSDH8RE9ERM7A8EIkE4Y1IiJ5MLwQERGRUBheSHZcAiIWrsUiIq1jeFERpxX0iVmNiEheDC+kGH6iJyIiZ2B4UcG1h9IabVrFKM21hzWjNJiISCEML0QEoDZUm+yZi6mLiLSJ4UVFRplEqekCjbLGR/TLPwhaNhEZCMOLikyi9m5EREQqYnghxTCqERGRMzC8qMBoC3SNSrSwxrclEYlC1vCydetWjBkzBpGRkTCZTFi9enWj22/evBkmk+m6r+zsbDnLJJkxrImF05lEpHWyhpeSkhIkJCRg4cKFDt0vNTUVWVlZ9q/Q0FCZKiRyPmY1IiJ5ucr54KNGjcKoUaMcvl9oaCgCAgKcXxAREREJT5NrXnr37o2IiAjcfvvt2L59e6PblpeXw2q11vkibeJ0BBEROYOmwktERAQWLVqEr7/+Gl9//TWio6MxdOhQ7Nu3r8H7zJ8/H/7+/vav6OhoBStuGU4rGINoYY1rk4hIFLJOGzmqS5cu6NKli/37gQMH4tSpU1iwYAE+/fTTeu8zd+5czJkzx/691WoVIsBcTe+dxrVnatV7e69lsOYSEclOU+GlPv3798e2bdsa/LnFYoHFYlGwIucR7IM5NdO1164iIiLn0tS0UX1SUlIQERGhdhlEhmG/thEzGBFplKwjL8XFxTh58qT9+zNnziAlJQWBgYGIiYnB3LlzkZmZiU8++QQA8PrrryMuLg7du3fHpUuX8MEHH2Djxo348ccf5SxTcUbtFDjSREREziBreNmzZw9uueUW+/c1a1MmT56MJUuWICsrC+np6fafV1RU4IknnkBmZia8vLzQq1cv/PTTT3UeQ0/YmeubaK9vzdok0eomIuORNbwMHTq00fn/JUuW1Pn+ySefxJNPPilnSaQCo440ERGRPDS/5kXPTMJd/Yaag1mNiEheDC9EREQkFIYXFVx73hOj4FoKIiJyBoYXIpmIFtZq1iZxOpOItI7hhWRnzHEmIiKSC8OLikT7ZE7NxLRGRCQrhhcNMFpfZ7Q1P7xcABGRczG8qMCofZnR1lKI1t6atyVHBIlI6xheSDFG6xRFbW9N2QbN2EQkAIYXFQnatznOqENNREQkC4YXFZlE/WhOjTLamh4iIqUxvKiAXRsREVHLMbyQYjjOpG32k9RxRJCINI7hhUgmjABERPJgeFGRUTo3TpMREZEzMbwQORkPriIikhfDixoM2rtxLQURETkDw4sGGO308QZrrkDTZpcrZcQkIq1jeCHZGS2s1OBIExGRPBhe1GSwvs1gzRW3vcIWTkRGwfBC5GS6GWnSTUOISG8YXlRg2C6Bn+g1zX6SOnXLICJqEsOLithJEBEROY7hhWRntKOp7JhOiYhkwfCiIh6NQkRE5DiGFxUYdSDCKCTjrmoiIlIEwwsphuNM2lYTuTgiSERax/BCREREQmF4IdkZdRKF4xdERPJgeFGRUTs3o4UZrnEiInIuhhcVGHVBp1HWUogeVozxKhGRyBheSDFG6xRFC2s15+MRrGwiMiCGFyKZiJ4BBB9AIiIdY3hRkVE+4Yo+jUJERNrC8KIik/Cfzak+zGpERPJieFEBRyJIyxiqiUjrGF5IMUaZJhMVMzURiYLhhUgmDGtERPJgeFGRUTo3fqInIiJnYnghcjKuaSIikhfDiwqu7duM1tkZrr1qF9BMNa+LUUYEiUhcDC+kGKMdxWK09hIRKYXhhRRjtE/0RmsvEZFSGF5IdpLR5omIiEhWDC8qYF+ud/p4gfk+JSKtkjW8bN26FWPGjEFkZCRMJhNWr17d5H02b96MPn36wGKxoGPHjliyZImcJRLRFbVXleZ8FxFpm6zhpaSkBAkJCVi4cGGztj9z5gxGjx6NW265BSkpKZg9ezamTZuGdevWyVmmathJEBEROc5VzgcfNWoURo0a1eztFy1ahLi4OLz66qsAgK5du2Lbtm1YsGABRowYIVeZqmF0ISIicpym1rwkJydj2LBhdW4bMWIEkpOTG7xPeXk5rFZrnS8iIiLSL02Fl+zsbISFhdW5LSwsDFarFWVlZfXeZ/78+fD397d/RUdHK1Fqq0g6WdBJ9RN9oStHBIlI6zQVXlpi7ty5KCwstH9lZGSoXRI1gGt8tE3wzEVEBiLrmhdHhYeHIycnp85tOTk58PPzg6enZ733sVgssFgsSpRHREREGqCpkZekpCRs2LChzm3r169HUlKSShXJq2YgQu/TSNdPo+i7vTVqRpp4kj4iIueSNbwUFxcjJSUFKSkpAC4fCp2SkoL09HQAl6d8Jk2aZN9+xowZOH36NJ588kkcP34c77zzDr744gs8/vjjcpZJRFfh7B4RaZ2s4WXPnj1ITExEYmIiAGDOnDlITEzEvHnzAABZWVn2IAMAcXFx+O6777B+/XokJCTg1VdfxQcffKC/w6QN+kHcKH1izcsrXHt5VWkiEoSsa16GDh3a6JB5fWfPHTp0KPbv3y9jVaQWdopEROQMmlrzQkTawbU6RKRVDC8qMspIhN4XJDfEKK8vEZHSGF5UULsmgr0baQ/fl0SkdQwvRE4m6nSLmFUTkRExvJBiOI1CRETOwPBCREREQmF4UZFRRiIEnUVpNS2+vjabQV8MItIVhhcVXLsmwmidu+Haq3YBVxw9Z0XiC+vx3DeHGw0xWgxdRERXY3ihZjt7sRTzfziGtAslapeiaVoJK9dKyShAYVklPk5Ow/P/O1JPiNZq5UREdTG8KGDNgXP465cHkF9SoXYprbJsZzre23Iaf3wvuUUBhofgquvq8+18nJyGV9amMrAQkZAYXhTwzqaT+GrvWdz/wU4UlNYGGNG68rLKagBAjrUc972/E2cvljp0f6NNR2gtrNXklDZebgCARVtO4a2NJ6/bTltVExFdj+FFAZXVNgDA0SwrJv53F6xlVSpX1DJXf0jPLCjDve/vQFZhWdP3k7EmLdNaWKt5HfrHBeKZ0V0BAK+t/xUf/HxavaKIiFqA4UUBNZ2G2QQcyizEij0ZqtbTWn/o2xaxQV7IyC/Dfe/vRK71ktolUXNcSZ8mmDDtpvb46/DOAIAXvzuGT3ek1W53JXUZNXQSkfYxvCjhSi/w3Jju8Pd0U7cWJwj388Cy6TciKsATZ86X4L4PduJ8cbnaZWmGVpeR2C9LcWVEaNatnfCXoR0AAM+uPoySimp1CiMichDDiwJqOo1ukX747MEB8PVwBQC4uYi1+2sWd5pMQFSAJz6ffiMi/D1wMrcY972/A3lFjQcYjc2iGE59oepvI7pg2uA45YshImoFsXpPHejZ1h+fT78R9/SLxh/6tlW7HIfUXlDyspggLyybfiNCfS34NacY9yxORnahvqeQfj6Rh5O5RWqX0SpXr8UxmUz4x+iueGBQbYBxMzNmEpG2MbwowD5iceX7HlH+eOUPvRAb5K1eUS0gXTvvACAu2Btf/DkJkf4eOJ1Xgj++l3zdUUhanUZx1NmLpZj4310Y/eY2bErNbXJ7rUUA6ao1L1czmUx49s6uePbObnjk1o4I9/dQozwiomZjeFFAPX2+kGrOE3JtM9oFe2PFn5MQHeiJ9PxS3PPeDl2eyK6gtBIAUF5lw0Of7MEPh7JUrsgx9gxZz/vQZDLhwcFxeGJ4FyVLIiJqEYYXBUiN9RoCqWlHfSEsOtALX/w5CXHB3sgsKMM97+3AqbxiZQuU2dUjSJXVEmYu24ev9p69fjsFa3KE/fVTtwwiolZjeFGAfcSigV5Dq53dtWrXvNTfkAh/T6x46EZ0CvVBtvUS7nlvB45lWWs3EPwQ3JrXMczPgnv6RcMmAX/98gA+Tf6tqTtqQu0IIOMLEYmN4UUBevnE29jIS41QPw8sf+hGdI3ww/nictzzXjJSMi4qU6BCXEwmzP9dT0wZ2A4A8Ow3R7Boy6nrN9RYSLh27RURkagYXhRQ2+kbo9sI8rFg+UM3ol9sG1gvVWHdkRwA4neaV7+OZrMJz43phlm3dAQAvPzDccz//lidqzVrtb1NvQ2N8j4lInExvJADmv/J3d/TDZ8+OAC3dAmRtyQFXTv7YzKZ8NcRXfDUqHgAwHtbT+OvXx5A1ZXLQRARkTwYXhQk+ufZ5kwbXc3T3QWLJ/XDuN6RAIC2bTxlqkwZV5+k72ozhnTAf/7QCy5mE1buz8Qvpy6oUF3T9DJ9SUTkqnYBRtBQpyealkx/ubmYseCe3ph5S0dUSxK+PSjW4cVXa+yQ9/H9ohHsY8Fflu6zX31ba2oXjjfv9dPL+XmISH848qKApo7SEYXUwsNmTCYTOoX5wkXw9FY7clF/O26JD8XS6QMQ4HX5+lVau44VR16ISC848qKAhqZbROvLHZ020p+mR9D6xLTBd4/ehH1pF5FXVI4tv+YpVFvTdHK6ISIijrwooaUjFlqjlxGklmruNEpUgCfGJETC1UVb+6mpkSMiIlEwvCiAIxb6IurL2NTJEomIRMHwoiDRP/EaPYTp5Qy1YldPRMTwogh9TBo1fGFGoxB9wSuPHiIivWB4UUBTIxaSKL2Kk0ZehGnvNaQWphetrXkSfOCIiIjhRRn6WGtg+AW7V/4VtfW11zYStQVERJcxvChAL0d5tPZke8KHNwdP0qe15jZ3zZLW6iYiuhbDiwIaOzOrSLQ1+aE80df86OV9SETE8KKA2uF6sRnt6tjXETy91S41MujrR0S6wfCiIKP2+Xqhl5GL5tYveFYjIh1jeFFAbScgdq8n+oJVZxF17ZLo015ERDUYXhTQ0JHBonUi+rk6dsvGFEQ8Sd/J3CI8viIFv5w6z/O8EJFu8MKMCtBNp3/lX5GbsfXXPDy2fD/GJUbh7yPj4eHm0uz7au18Lc3x5d6zWLU/E6tTMhEb6AVA/PchERFHXmRwobgc//rfUWw/eR6APjp9AFedpE7cluw4fQEXSyvx0fbfMOatbTicWdjs+4q4YLmy6nLRkgT8dqEUgLjTXlpXVW3Di98exdKdabDZxAu6RCJheJHBd4ey8OH2M5jwwU48981hVNtqRl7E7jS2nsgDIPYn96u7lBO5xbj7ne1YuOmk/TVqzn1Fan7NaFG/2DYI8HIDAHhbOOAqh4OZhfhg2xn8Y9Vh3PfBDmTkl6pdEpFuMbzIoLzSZv//x8lpKK2oBtBwpyfCZ7SLJRUoulQFoPWdt5rtrRk9+V1iFEb1CEdltYT/rEvFPe8lI/1C451Na9fKqKHmufvHBWLtYzfj2Tu74cHBceoVpGMVVbW/9ztO52PUGz9jxe50YS+HQaRlDC8yqPm0Gx/ui1Bfi/12s8BDFsXlVfb/39g+SMVKWqfmtWnj7Y53JvTB/41PgI/FFXvSLmLUG1sbHfIX+VBpkwkI9/fAg4PjEHLVe5KcpyajhPpa0C+2DYrLq/D3rw/hwY/3INd6Sd3iiHSG4UUGNX/EukX44cfHb8a9/aMxJiESbdt4qltYK9S0ydPNBZ3CfFv4KBro9a+6tqLJZMIf+rbFD4/dhP7tAlFSUW0f8v/tfEnD921uMzSUchxZ56KhsoVSE4z9Pd2w4s9JmDsqHu4uZmw8novbF2zFNymZHIUhchJFwsvChQvRrl07eHh4YMCAAdi1a1eD2y5ZsgQmk6nOl4eHhxJlOp8JCPByx/zf9cJb9ybCbBa/VxC9Y6tv9CQ60AufP3Qj5t3ZDZ5uLthxOh8j39iK97eerrMWpvY8KeLsBL0c6SaEq8Kti9mEPw/pgG8fHYweUX4oLKvEY8tTMHXJbpy9yLUwRK0le3hZsWIF5syZg+eeew779u1DQkICRowYgdzc3Abv4+fnh6ysLPtXWlqa3GU6lR6vviziYcL1qe3M6742LmYTHhgch3Wzb8agjkG4VGnDS98fw+/e/QWp2UV1thUpCOjjVRPL1b/3ncN8seovgzDn9s5wdzFjc2oehi/Yig+3nWnWInEiqp/s4eW1117D9OnTMXXqVHTr1g2LFi2Cl5cXPvzwwwbvYzKZEB4ebv8KCwuTu0ynEvFkZk2Rrppu0YOG2hET5IXPHhyAl3/XE74WVxzIKMCdb/2M135MxaUrC7FF2gd6e920rKE1UW4uZjx6Wyd8f2V6srSiGv/69ih+9+4vOJZlVbxOIj2QNbxUVFRg7969GDZsWO0Tms0YNmwYkpOTG7xfcXExYmNjER0djbFjx+LIkSMNblteXg6r1VrnS216GaW4Wu0fZrG7Qal2WKxBJpMJf+ofg/VzhmBY1zBUVkt4c+NJPPHFgZoNZK/TWSSHF+pcdV+uz3BIU7urY6gPlj90I166u4c9GI95axteWXscpRVVjd+ZiOqQNbycP38e1dXV142chIWFITs7u977dOnSBR9++CG++eYbfPbZZ7DZbBg4cCDOnj1b7/bz58+Hv7+//Ss6Otrp7XCUHj/t6ubK2Ff+bc6UXri/B96f1BfvTOiDMD8LyiobP+Rdi/T4XtSq5nxoMZtNmDAgFj89MQQju4ejyibh3c2ncNurW/DdwSwGRqJm0tzRRklJSZg0aRJ69+6NIUOGYOXKlQgJCcF7771X7/Zz585FYWGh/SsjI0PhihvW1IddkUYx7H9SxSm5Xo5O6ZlMJtzRMwI/zRmCBwbFwWwCuka09Ggr9Qj0VhOWI2dgDvPzwKKJfbF4Yl+0beOJrMJLmLlsH+7/706czC1q8v5ERifrqTaDg4Ph4uKCnJycOrfn5OQgPDy8WY/h5uaGxMREnDx5st6fWywWWCzaPG+Fnhbs1hC9RS2d0vP1cMO8Md3w+O2d4O0uzhlq9bh4XKtacgbm4d3DcXPnELy7+RTe3XIK209ewMjXf8YDg+Pw6G2d4MOzIRPVS9aRF3d3d/Tt2xcbNmyw32az2bBhwwYkJSU16zGqq6tx6NAhREREyFWm0+nx8FQRr+tTn9ZOo/h6uAl1yLseF49rVUt/7z3cXPD47Z3x0+OX11hV2SQs3noat/zfZizflc6jkojqIfu00Zw5c/D+++/j448/xrFjx/Dwww+jpKQEU6dOBQBMmjQJc+fOtW//r3/9Cz/++CNOnz6Nffv24f7770daWhqmTZsmd6lOo88OQ1+BTC/taBo7PqW09gzMMUFe+GByP3w05QbEBnkhr6gcT608hDve+Blbfs1zWp1EeiD7mOQ999yDvLw8zJs3D9nZ2ejduzfWrl1rX8Sbnp4Os7k2Q128eBHTp09HdnY22rRpg759++KXX35Bt27d5C7VaWq7i2b+FROgf3HmOkJ1r/Wj/InmtHBtI8NkNQ1o7XvrlvhQDOwYhE+T0/DWxpNIzSnC5A934ebOIXj6jnjEh/s5qVIicSkyoTpr1izMmjWr3p9t3ry5zvcLFizAggULFKhKPno8YEDEKyo3xigjL/ocBdQoJ/7eW1xdMO2m9vhD37Z4a+NJfJL8G7b+modtJ/Iwvm80Zt/eCRH+4l5uhKi1NHe0kR5IOptiAZyz5kUL+0PJEKaB5to58rppqW6RyPF7H+Dljmfv7Iaf5gzB6J4RsEnAij0ZGPKfzXjh26M4X1zuvCcjEgjDi4z02AmI3iaptQsTBKPHEyZqlZxTdLFB3lg4oQ++fjgJ/dsFoqLKhv9uO4ObXtmEf689joLSChmelUi7GF5koMeher2MJtVeXNEY9Phe1ColgnHf2ECs+PON+OSB/kho64+yymq8s/kUbnplE97ccAJFlyple24iLWF4kYEez61Ru45H7DYZrTPX43tRq5SakjSZTLi5cwhWzxyExRP7Ij7cF0XlVXht/a8Y/MomvP7TrygsZYghfWN4kYOuz/Oibh2tZbRJFD0uHtcqpc/vZDKZMLx7OL5/9Ca8dW8i2od4o7CsEq//dAKDXtmIl384jrwirokhfWJ4kYHejswB9LN2onZdgp5enYbpZbpPBGr93pvNJoxJiMT6x4fgrXsTER/ui+LyKizacgqDX9mIf645gnMFZQpXRSQvhhcZNPfTrkgdin7OF2LMztxgzTUklysh5ofHbsIHk/ohIToA5VU2LPnlNwz5zyb89csDOJ5tVbtMIqfghTNkUPtpV39dhl6apJNmNE0n030i0MolNEwmE4Z1C8NtXUOx/eQFvL3pBHaczsdXe8/iq71ncVOnYDwwOA5DOoUIdakLoqsxvJBDRJ9u0cvanebigl0laetINpPJhMGdgjG4UzD2pV/Ef7edwQ+HsvDzifP4+cR5dAz1wYOD43B3YhQ83FzULpfIIQwvMnC0gxRhNYkzO31JxVWkanw6VnO9kB4vEqpVWg7GfWLaoM99bZCRX4qPf/kNy3dn4GRuMeauPIT/rEvFff1jcO+AGEQF8Ky9JAaueZGBHj/t6uX8KHpZeNxcxmqtukT4vY8O9MIzd3ZD8txb8czorogK8ER+SQXe3nQSN72yEdM+3oPNqbmw8UrWpHEceZGBlj+BtZRTLg/gpFpaQ8nXRguvf0sGudResyEqSaDDDH093DDtpvaYMrAd1h/Nwac70vDLqQv46VgOfjqWg5hAL9w3IAZ/7BeNQG93tcslug7Diwz0MkpxNb18DhPh07Ez1Z701RjtVZOIv/euLmaM6hmBUT0jcDK3GEt3puGrvWeRnl+Kl384jtfW/4o7eoRjfL9oJLUP4gJf0gyGFznopae/iprrVJxJJ81wGLsc+Yn+3uoY6oPnxnTHkyPi8b8D5/DpjjQcyizE6pRzWJ1yDlEBnvh937YY37ctogO91C6XDI7hxQE7T1/A57vScXPnEIzoHg5vS+O7T6sfdv+z7jhKyqtxZ68I9I1t49gVhzXapmt9dzALG4/n4vZuYbglPgQW18tHU+j9pG3fHczC5tTL7R7aJZQLdp0gv6QC/157HLFB3hjbOxKRTSxqFX1fe7q74I83ROOPN0TjQEYBvtiTgTUHziGzoAxvbjiBNzecQFL7IIzv1xajekTA051HKpHyGF4c8N9tZ/Dj0RysTjkHD7dDGN4tHOMSI3FTpxC4udSufdbyUP354nIs3HQKALDkl98QHeiJsQlRGJcYiY6hvg3eT7SLMb+18QSOZxfh631n4e/phtG9InB3YlTteU/ULU82Ne3+cu9ZtPFys78v9dpeJWw8novluzMAAP9edxz92wXi7sQojOoZAX9PN/t2epySTIgOQEJ0AJ69sxvWHcnGl3vOYvup80g+fQHJpy9g3jdHMKJ7OO7qHYlBHYLg6sJjQEgZDC8OqKy2AQC83V1QUlGNNQfOYc2Bcwj0dsedvSIwLjEKidEBtZ921Sy2AVXVtWPb3u4uyMgvw9ubTuLtTSfRI8oP43pHYUxCJML8POrcT7TT6l/9WhWWVWLZznQs25luD1+ihDBH1bTb080FF6+6OF9LgrTo0yDOUrNPLa5mlFfZsPNMPnaeyce8b47glvgQ3J0YpftRLg83F4ztHYWxvaOQWVCGr6+c8C49vxRf7zuLr/edRZC3O+7oGYGxvSPRJ6YN18eQrBheHFDzt/yfd3VH5zBfrNqfiW8PnsP54gp8kpyGT5LTEB3oCY8rUxRa7Odrpk3cXEzY88ztWH8sB9/sz8SWX/NwONOKw5lW/L/vj6F/XCBG94zAiB7hCPX1gGin1a95rT6YfANskoRV+zPxw6EslFRUA4BuT8pV0+7/TumHymoJq/dn4tecIgzqGKxqXSKrCXE3dw7BP+/qjjUp57B6fyZSc4qw7kgO1h3Jga+HK+KCvQGI8zvSUlEBnnj0tk6YdUtH7E2/iDUp5/DdoSxcKKnApzvS8OmONEQFeOLOhAjclRCJbhF+mhyFJrExvDjg6sOFa4ZTnxndFdtOnsc3Keew7kg2MvJrL4DW1CiFGr/OV4+geLq74K6ESNyVEIn8kgp8d/AcVu3PxL70Auw4nY8dp/Mxb80R3NAuEF3CfFWruUWutNPFbEJSXBAGdQzGC2N74KdjOTiWZcVdCZHq1ieXK+12NZsxsEMghnQOUbceHbj63EBRAZ54eGgHPDy0A45lWbF6fya+STmHbOslHDxbCECc0cnWMptNuKFdIG5oF4h5Y7ph+8nzWHPgHH48koPMgjK8t+U03ttyGjGBXhjZIxwjuochMZojMuQcDC8OqG8U3dXFjKFdQjG0SyhKK6qw8XguvjuYhYNnC3FzJ+192m1oJiDQ2x0Tk9phYlI7ZOSX4ofDWfjuUDYOZBRg15l87DqTD0Cb63jqU187Pd1dMCYhEmP0GlygywPdVNfQRUm7Rviha4Qf/j4yHrt/y8d3h7Kw7cR5jOoZrniNanO76u/gpcpqbDyeizUp57ApNRfp+aVYvPU0Fm89jVBfC4Z3D8OI7uG4sX1QnbWCRI5geGmBhrpvL3dX3NkrEnf2EqBzbCSDRAd64aGbO+Chmzvg7MVSrD2cje8OZWF/egE6hfooV6MTCJK1nM6o7ZZTQ/vUbDZhQPsgDGgfpGxBGuXh5oI7ekbgjp4RKK2owpbUPKw9ko2Nx3KRW1SOz3ak47Md6fD3dMNt8aG4tWsobuoUUmfxM1FTGF4cINeCPCXPoeLoYuK2bbww7ab2mHZTexSUVsDLvfVvGSVaq6VF00oufNVSu/VCj0cRKcXL3dV+EryKKht+OXUe645kY/3RHJwvrsDK/ZlYuT8TLmYT+sa2wS1dQnFrfCg6h/kIM8pL6mB4aQGRf6dac3r8AC9xThMu2qHdzmLUdstKx0cRKcndtXZq6cVxEvamXcT6o9nYlJqHk7nF9unpV9YeR1SAJ4Z2CcEtXUIxsGOQUz40kb7wHeEA0Q4XbowabVD0Ss720Q71Xis19rEz2i3+u9u5GAidz8VsQv+4QPSPC8Q/RgMZ+aXYlJqLjcdzkXzqAjILyrB0ZzqW7kyHm4sJiTFtMLhjMAZ1DEKvtgFcK0MML47Qw9lZ9XjRyPro4bVqCaO2W056+tCiVdGBXpiU1A6TktqhrKIayafPY9PxPGxKzcXZi2X2UZnX1gM+Flfc2D4QgzoGY1DHYHQK5RSTETG8OEAPJ+0S8eJxLdHQESJ6Z9R2y0kv1/UShae7C26ND8Ot8WGQJAnp+aXYdvI8tp88j19OXUBBaSV+OpaLn47lAgBCfC0YEBeIAXGBuCEuEJ1DfXk4tgEwvDhAD3/Drj5XjZ4ZpZ3XMmq75aSBGUjDMplMiA3yRmyQNyYMiEW1TcLRc1ZsP3U5zOw6k4+8onJ8ezAL3x7MAgD4e7rhhnZtcEO7y9NSPaL8Oc2kQwwvLcCOgch4+FuvPhezCT3b+qNnW3/MGNIBlyqrsS/9InafuYhdv13AvrQCFJbVHZnxdHNBYkwA+sW2Qe+YAPSOboNAb3EOPqD6Mbw4QA9TLrWHfRqDUdp5LWe0W+Ip7wBwNEvLPNxcMLBDMAZ2CAbQCZXVNhw5Z8XuK9ef2pOWj4LSSvxy6gJ+OXXBfr/YIC/0jg6wf3WL9LNfeZ7EwPDiAD0sdpUMsihCzxfJa4xR2y0nowV+kbm5mO2BZPrN7WGzSTiZd/kw7P3pBUjJuIhTeSVIu1CKtAul+CblHADA3cWMbpF+6B0dgJ5R/uge5YeOIT68SraGMbw4wNknq1KjgzHKH2KjnljMqO2WEwOhuMxmEzqH+aJzmC/uvzEWAFBYWokDZwuQklH7lV9SYf9/DXdXM7qG+6JbpD96RPmhe6Q/4sN9dXtRV9EwvDhCB6PoRhkC18MoWUsYtd1EzeXv5YabO4fg5isXLZUkCRn5ZdifcREpGQU4cs6KY+esKCqvwoGzhThw5YKbwOU1Nx1DfNA90g/xEZdDUZdwX4T7eej+b6rWMLw4QK5zaCibiTTwKVKBBmtpvYaSlWip3XphkJlWwzKZTIgJ8kJMkBfG9o4CANhslw/RPnLOisPnCnHknBVHMgtxoaQCqTlFSM0pAvbXPoafhyu6hNeGmc5hvugS5os2XBgsG4YXB+jhj5ge2tAcRh2BMGq75VT7oYU71SjMZhPaBXujXbA3RveKAHB5hCbHWo7DmZfDzK9XQsyZ8yWwXqrC7t8uYvdvF+s8ToivBZ3DfNAhxAftg73RPsQH7UO8EenvyXPRtBLDSwvwb1jLGG23Cfs+EbVumXG3GJvJZEK4vwfC/T0wrFuY/fbyqmqcziu5HGayi+yhJiO/DHlF5cgrKsf2kxfqPJaHmxntgrzRIcQHccHeaB9SG2z8PHh17eZgeHFA7YC8uH/Gaq/TIm4bmsOoC1eN2m45SUZZ5U4tYnF1QdcIP3SN8Ktze0l5FU7kFuPXnCKczivB6bxinD5fgrQLJbhUacPx7CIczy667vGCvN0RHeiFmEAvxAZ52f8fE+iFcD8PjthcwfDiAD0cdcBpI30zarvlxEBILeFtcbUftn21qmobzl4sw+nzxTidV4JTVwWbvKJyXCipwIUrRz9dy93FjLaBnvYwExPohbZtvBAV4InIAA8Eervr/oNpDYYXB+hhKaRxLtxnlHZey6jtlo8eLgtC2uHqYravp7k1vu7Pii5VIj2/FOkXSi//e+UrI78UZy+WoaLadmUUp6Tex7a4mhF5JchE+HsiMsATUQEeiAzwvPK9B7zc9dHt66MVCtHDqEXtH2KRW9E0o14J2KjtlpNxAj+pzdfDDd0j/dE90v+6n1VV25BVeAkZV4WatPxSZF4sw7mCMuQVl6O8yoYz50tw5nz94QYA2ni5IdzfE2F+FoT6WhDm54FQXwtC/Tzs/w/xtWj+elAMLw7Qw3oRo0wr1L5WqpahOKO2W056+NBC4nN1MSM68PIamIH1/LyiyoYc6yVkFlwOM1mFtf+//HUJxeVVuFhaiYullTiW1fBzmUxAoJc7Qq+EmTC/yyEnxNeCEB8Lgn0vB5/YIG/Z2tsUhhdH1Kx5UbmM1tDD9ZmaQ9LBa9USzmw3p0vqYiAkLXN3rQ03DbFeqrQHmzxrOXKsl5BbVPtv7pV/q2ySfe1NQyEnKsAT25+6VabWNI3hxQF6+lSrhzY0h1HaeS2jtltOnIoj0fl5uMEv3A3x4X4NbmOzSbhYWoEcazlyiy4h98q/2dZLyCsqx/niCpwvLkekv6eClV+P4aUFnNcxKP/H0Cifpg3SzOsYtd1y0sNRhkTNZTabEORjQZCPBd3QcMhRm7ZX5GiMnhZD6qENjTHKwuRrGbXdcjLKOjEikTC8OECu68YoORqihT/ESjRXS5+WlX19tdNuveBoFpH2KBJeFi5ciHbt2sHDwwMDBgzArl27Gt3+yy+/RHx8PDw8PNCzZ098//33SpTZJD2cadMwC3av/Kv3dl7LqO2WE0eziLRH9vCyYsUKzJkzB8899xz27duHhIQEjBgxArm5ufVu/8svv+Dee+/Fgw8+iP3792PcuHEYN24cDh8+LHepTdLDIZO1Iy/Kt0LRp1SxnTVUeWYntFvvU4qO4nleiLRH9gW7r732GqZPn46pU6cCABYtWoTvvvsOH374IZ566qnrtn/jjTcwcuRI/O1vfwMAvPDCC1i/fj3efvttLFq0SO5yG1RZbcPZi6UAtH2eF5tNQpb1kn364Fq5ReUKVyQvSZKQVXgJtmvaa9PpodIVVTbkFl1q8OdabXdxeRUKSivULqNFrGVVALS3T4mMTNbwUlFRgb1792Lu3Ln228xmM4YNG4bk5OR675OcnIw5c+bUuW3EiBFYvXp1vduXl5ejvLy2Q7Zara0vvB4Z+aWwXtL+H7EZn+3Fj0dzmtxOw/nLIU+vOozPd6U3+HO9tBMAqm0SRr6+FacbOXtmDS21OyO/FLcv2IJLlTa1S2kVLe1TIqOTNbycP38e1dXVCAsLq3N7WFgYjh8/Xu99srOz690+Ozu73u3nz5+P559/3jkFN8JkMsHiakZUG0/0anv9qZu14sDZAgCXL+DV0B9bkwm4s1ekckXJqObiZW4uJpivaXCPKH9EBah7LgJnKr5UZQ8u7q7mBkO01tr9a06RPbhYXMU8RsDL3QW3xoeqXQYRXSH8eV7mzp1bZ6TGarUiOjra6c8TF+yN1BdHOf1x5bJq5sB6r4+hVx9OuQE3dQpRuwzFHH1+BFw1fu2RayVEB+CbmYPULoOIdEDW8BIcHAwXFxfk5NSdxsjJyUF4eHi99wkPD3doe4vFAovF4pyCiYiISPNk/ejm7u6Ovn37YsOGDfbbbDYbNmzYgKSkpHrvk5SUVGd7AFi/fn2D21NdRjmDbo2GFibrkVznGWr4+Zz0OMZ5iYhIIbJPG82ZMweTJ09Gv3790L9/f7z++usoKSmxH300adIkREVFYf78+QCAxx57DEOGDMGrr76K0aNHY/ny5dizZw8WL14sd6lEREQkANnDyz333IO8vDzMmzcP2dnZ6N27N9auXWtflJueng6zuXYAaODAgVi2bBmeeeYZPP300+jUqRNWr16NHj16yF2q4uQ8esFo5+owXHsFPPRFvIqJSKsUWbA7a9YszJo1q96fbd68+brbxo8fj/Hjx8tclXYoPR2gNiNN9VxmtPYSEclLrMMVqEnsJvVL1MwnaNlEpGEML6QY403tqF1By4haNxEZB8OLThmtAzJce9UuoAWM9hoRkXwYXoiIiEgoDC86I+q6iJYyUntFbarxFmgTkdwYXoiIiEgoDC86ZbT1BQZrrpCvr4AlE5FGMbwQERGRUBhedMdY6wuMdII/pdeOOOvpjPMKEZFSGF6IiIhIKAwvKpJzDYDRTghnuOYKuOhFxJqJSJsYXjTAaEeSGqy5hnt9iYjkxvCiM0brKI3UXlGbaqTXiIiUwfBCijHarIGoU3diVk1ERsLwolMMCqQ1fIWIyFkYXnSGI/SkPXxXEpFzMbyQ0IzULXLtCBHRZQwvOsUhetIao01lEpF8GF5IF4zUMRqprURE9WF40RmlTyFP1BS+JYnI2RheSGhGCmtKX8fJSNeNIiKxMLzoFKcWSGt4ODsROQvDi4p4rRfnMdKeNFJbiYjqw/CiAc6c+RBhoN9AMz0AxHhN5GT09hOR8zG8kNAM1TEaqrFERA1jeNEtTi6oTtCXQLbZTEH3BxFpD8ML6YKR1g8Zqa1ERPVheNEZo60nIe3je5KInI3hhcRmoI7RQE0lImoUw4tOcWaBtIZvSSJyFoYX0gUjhTUDNZWIqF4MLzpjpNPlkxh4mQEicjaGFxKakbpFpXMpczARaRXDi05xaoG0xkhTe0QkL4YXFRn1b7kc0wha3pfOnspjCCAio2N40RmO9JPWcPqJiJyN4YUUI8eIgZYXKDu7uUotfDVpehyLiIjhRbd4CnnSGoYiInIWhhe90e5AhKyMlNVECwEGfUsSkYwYXoiIiEgoDC86JdZn85Yz0qd6DS/vaRYjjY4RkbwYXoiIiEgoDC86I/iH81Yw0Md6wZqq5SPCiEhMDC9EREQkFIYXnTLK+gIjfagXvalGeU8SkfxkDS/5+fmYMGEC/Pz8EBAQgAcffBDFxcWN3mfo0KEwmUx1vmbMmCFnmURERCQQVzkffMKECcjKysL69etRWVmJqVOn4qGHHsKyZcsavd/06dPxr3/9y/69l5eXnGWqpuaTqDPXBIiwvkCOErX8qd7ZzdVwU4mIFCFbeDl27BjWrl2L3bt3o1+/fgCAt956C3fccQf+7//+D5GRkQ3e18vLC+Hh4XKVRkRERAKTbdooOTkZAQEB9uACAMOGDYPZbMbOnTsbve/SpUsRHByMHj16YO7cuSgtLW1w2/Lyclit1jpfpM2zsMpxyQKlrvfTEs5ur1KjanKNYmnxPUlEYpJt5CU7OxuhoaF1n8zVFYGBgcjOzm7wfvfddx9iY2MRGRmJgwcP4u9//ztSU1OxcuXKerefP38+nn/+eafWTkRERNrlcHh56qmn8MorrzS6zbFjx1pc0EMPPWT/f8+ePREREYHbbrsNp06dQocOHa7bfu7cuZgzZ479e6vViujo6BY/v+i0Ow4hLyN9ptfy+p76CLAMi4gE43B4eeKJJzBlypRGt2nfvj3Cw8ORm5tb5/aqqirk5+c7tJ5lwIABAICTJ0/WG14sFgssFkuzH4+IiIjE5nB4CQkJQUhISJPbJSUloaCgAHv37kXfvn0BABs3boTNZrMHkuZISUkBAERERDhaqqGJ9um8pYz0qV70thrlPUlE8pNtwW7Xrl0xcuRITJ8+Hbt27cL27dsxa9Ys/OlPf7IfaZSZmYn4+Hjs2rULAHDq1Cm88MIL2Lt3L3777TesWbMGkyZNws0334xevXrJVSoREREJRNaT1C1duhTx8fG47bbbcMcdd2Dw4MFYvHix/eeVlZVITU21H03k7u6On376CcOHD0d8fDyeeOIJ/P73v8f//vc/OcvUFdE/nbeUHEcyaZVoR+1o+YgwIhKTrCepCwwMbPSEdO3atatz+Gd0dDS2bNkiZ0lEREQkOF7biIRm1JEmJYhwtmYiMiaGFw1wZhchwhC99it0LqNnAKO3n4icj+FFRaKtXdAyI+1JAy3vISKqF8MLKcZofa6o7WU4IiKtY3jRKXZA+iP69IuRjggjInkxvOiM6B0c6Q/fk0TkbAwvpAtG+lBvoKYSEdWL4UWnOERPWsN3JBE5C8MLCc1I5yIR4TB4IiIlMLzoDLs30hq+J4nI2RheSBeMdM4cTgkSkdExvOgUuzdqLWfPyDFzEZGzMLyQ0Iw0JWGg5T1ERI1ieNEAp3ZKInRwItToRAZr7nWMtKiaiJTB8KIiDqM7j5H2pYGaSkRUL4YXndJiZ67FmuQkbnvlKVzY3UFEmsPwQkIz0oyEgZpKRNQohhed4YnMSGv4jiQiZ2N4IRIN51+IyOAYXnTKSCdtIzHw5HpE5CwMLyQ0I02T8ZBjIqLLGF50hv0baQ7fk0TkZAwvpAtGmpEwUFOJiOrF8KJTRurMSR7OHjDhW5KInIXhRQOcuW5DhBF6p7ZXgAY7a62KAE2tl5HWJRGRMhheiIiISCgMLzqlxSF6OQ/f1uKh4XJN3cl9yLF8dcvzuERkPAwvREREJBSGF50x2rlAjNRaUV9aUesmIu1ieCEiIiKhMLzolcHWFxhpPYW4bRW2cCLSGIYXIiIiEgrDi84YbXmBsdZTiNlYMasmIi1jeCEiIiKhMLzolBbPeyIncdeBOE7UphrpNSIieTG8kOCMMymh9BSZsabkiEgkDC8a4MxOQoQOR4QayXn4ehORszG8qIjD6ERERI5jeNEpLQYjOWvS4hofuWqS/dpGgj0uERkPwwsJzUhTEgZqKhFRoxheiEhWEmMXETkZwwvpghanyeRioKYSEdWL4UWn2MGR1hgpYBKRvBheSGhGmpAw0voeIqLGMLzoiMTejTSIb0sicjbZwstLL72EgQMHwsvLCwEBAc26jyRJmDdvHiIiIuDp6Ylhw4bhxIkTcpVIOmKkGQlOvxCR0ckWXioqKjB+/Hg8/PDDzb7Pv//9b7z55ptYtGgRdu7cCW9vb4wYMQKXLl2Sq0zdkvtcIESO0uK5eIhITK5yPfDzzz8PAFiyZEmztpckCa+//jqeeeYZjB07FgDwySefICwsDKtXr8af/vQnuUpVXbVNwtmLpa1+HFGG522Sc9oLXN53Wnepstop7c2xKhviC8sqnVJ3QWmFE6ohIqolW3hx1JkzZ5CdnY1hw4bZb/P398eAAQOQnJzcYHgpLy9HeXm5/Xur1Sp7rc52oaQCg1/ZpHYZirFJMFR7d/92Ucj2frj9DD7cfkbtMoiIrqOZ8JKdnQ0ACAsLq3N7WFiY/Wf1mT9/vn2URzQdQnzQM8ofv+YUOfVx+8cFoo2Xm1Mf0xlCfCwY2CEIe9MuOvVx4yP8EBvk7dTHdIa+sW0QHeiJXGt50xs7YExCpFMf71q3dQ3F94eyUFpR7bTHdHc1Y3j3sKY3JCJqBofCy1NPPYVXXnml0W2OHTuG+Pj4VhXliLlz52LOnDn2761WK6KjoxV7/tbwcHPB/x4ZrHYZijGbTVg2/Ua1y1BMdKAXfn7yVrXLcNit8WHYP2+42mUQETXIofDyxBNPYMqUKY1u0759+xYVEh4eDgDIyclBRESE/facnBz07t27wftZLBZYLJYWPScRERGJx6HwEhISgpCQEFkKiYuLQ3h4ODZs2GAPK1arFTt37nToiCUiIiLSN9kOlU5PT0dKSgrS09NRXV2NlJQUpKSkoLi42L5NfHw8Vq1aBeDyob2zZ8/Giy++iDVr1uDQoUOYNGkSIiMjMW7cOLnKJCIiIsHItmB33rx5+Pjjj+3fJyYmAgA2bdqEoUOHAgBSU1NRWFho3+bJJ59ESUkJHnroIRQUFGDw4MFYu3YtPDw85CqTiIiIBGOSdHZOeavVCn9/fxQWFsLPz0/tcoiIiKgZHOm/eW0jIiIiEgrDCxEREQmF4YWIiIiEwvBCREREQmF4ISIiIqEwvBAREZFQGF6IiIhIKAwvREREJBSGFyIiIhKKbJcHUEvNCYOtVqvKlRAREVFz1fTbzTnxv+7CS1FREQAgOjpa5UqIiIjIUUVFRfD39290G91d28hms+HcuXPw9fWFyWRy6mNbrVZER0cjIyOD102SEfezMriflcH9rBzua2XItZ8lSUJRUREiIyNhNje+qkV3Iy9msxlt27aV9Tn8/Pz4i6EA7mdlcD8rg/tZOdzXypBjPzc14lKDC3aJiIhIKAwvREREJBSGFwdYLBY899xzsFgsapeia9zPyuB+Vgb3s3K4r5Whhf2suwW7REREpG8ceSEiIiKhMLwQERGRUBheiIiISCgML0RERCQUhpdmWrhwIdq1awcPDw8MGDAAu3btUrskoc2fPx833HADfH19ERoainHjxiE1NbXONpcuXcLMmTMRFBQEHx8f/P73v0dOTo5KFevDyy+/DJPJhNmzZ9tv4352nszMTNx///0ICgqCp6cnevbsiT179th/LkkS5s2bh4iICHh6emLYsGE4ceKEihWLp7q6Gs8++yzi4uLg6emJDh064IUXXqhzPRzuZ8dt3boVY8aMQWRkJEwmE1avXl3n583Zp/n5+ZgwYQL8/PwQEBCABx98EMXFxfIULFGTli9fLrm7u0sffvihdOTIEWn69OlSQECAlJOTo3ZpwhoxYoT00UcfSYcPH5ZSUlKkO+64Q4qJiZGKi4vt28yYMUOKjo6WNmzYIO3Zs0e68cYbpYEDB6pYtdh27doltWvXTurVq5f02GOP2W/nfnaO/Px8KTY2VpoyZYq0c+dO6fTp09K6deukkydP2rd5+eWXJX9/f2n16tXSgQMHpLvuukuKi4uTysrKVKxcLC+99JIUFBQkffvtt9KZM2ekL7/8UvLx8ZHeeOMN+zbcz477/vvvpX/84x/SypUrJQDSqlWr6vy8Oft05MiRUkJCgrRjxw7p559/ljp27Cjde++9stTL8NIM/fv3l2bOnGn/vrq6WoqMjJTmz5+vYlX6kpubKwGQtmzZIkmSJBUUFEhubm7Sl19+ad/m2LFjEgApOTlZrTKFVVRUJHXq1Elav369NGTIEHt44X52nr///e/S4MGDG/y5zWaTwsPDpf/85z/22woKCiSLxSJ9/vnnSpSoC6NHj5YeeOCBOrf97ne/kyZMmCBJEvezM1wbXpqzT48ePSoBkHbv3m3f5ocffpBMJpOUmZnp9Bo5bdSEiooK7N27F8OGDbPfZjabMWzYMCQnJ6tYmb4UFhYCAAIDAwEAe/fuRWVlZZ39Hh8fj5iYGO73Fpg5cyZGjx5dZ38C3M/OtGbNGvTr1w/jx49HaGgoEhMT8f7779t/fubMGWRnZ9fZ1/7+/hgwYAD3tQMGDhyIDRs24NdffwUAHDhwANu2bcOoUaMAcD/LoTn7NDk5GQEBAejXr599m2HDhsFsNmPnzp1Or0l3F2Z0tvPnz6O6uhphYWF1bg8LC8Px48dVqkpfbDYbZs+ejUGDBqFHjx4AgOzsbLi7uyMgIKDOtmFhYcjOzlahSnEtX74c+/btw+7du6/7Gfez85w+fRrvvvsu5syZg6effhq7d+/Go48+Cnd3d0yePNm+P+v7W8J93XxPPfUUrFYr4uPj4eLigurqarz00kuYMGECAHA/y6A5+zQ7OxuhoaF1fu7q6orAwEBZ9jvDC6lu5syZOHz4MLZt26Z2KbqTkZGBxx57DOvXr4eHh4fa5eiazWZDv3798P/+3/8DACQmJuLw4cNYtGgRJk+erHJ1+vHFF19g6dKlWLZsGbp3746UlBTMnj0bkZGR3M8GwmmjJgQHB8PFxeW6oy9ycnIQHh6uUlX6MWvWLHz77bfYtGkT2rZta789PDwcFRUVKCgoqLM997tj9u7di9zcXPTp0weurq5wdXXFli1b8Oabb8LV1RVhYWHcz04SERGBbt261bmta9euSE9PBwD7/uTfktb529/+hqeeegp/+tOf0LNnT0ycOBGPP/445s+fD4D7WQ7N2afh4eHIzc2t8/Oqqirk5+fLst8ZXprg7u6Ovn37YsOGDfbbbDYbNmzYgKSkJBUrE5skSZg1axZWrVqFjRs3Ii4urs7P+/btCzc3tzr7PTU1Fenp6dzvDrjttttw6NAhpKSk2L/69euHCRMm2P/P/ewcgwYNuu5w/19//RWxsbEAgLi4OISHh9fZ11arFTt37uS+dkBpaSnM5rpdl4uLC2w2GwDuZzk0Z58mJSWhoKAAe/futW+zceNG2Gw2DBgwwPlFOX0JsA4tX75cslgs0pIlS6SjR49KDz30kBQQECBlZ2erXZqwHn74Ycnf31/avHmzlJWVZf8qLS21bzNjxgwpJiZG2rhxo7Rnzx4pKSlJSkpKUrFqfbj6aCNJ4n52ll27dkmurq7SSy+9JJ04cUJaunSp5OXlJX322Wf2bV5++WUpICBA+uabb6SDBw9KY8eO5SG8Dpo8ebIUFRVlP1R65cqVUnBwsPTkk0/at+F+dlxRUZG0f/9+af/+/RIA6bXXXpP2798vpaWlSZLUvH06cuRIKTExUdq5c6e0bds2qVOnTjxUWm1vvfWWFBMTI7m7u0v9+/eXduzYoXZJQgNQ79dHH31k36asrEz6y1/+IrVp00by8vKS7r77bikrK0u9onXi2vDC/ew8//vf/6QePXpIFotFio+PlxYvXlzn5zabTXr22WelsLAwyWKxSLfddpuUmpqqUrVislqt0mOPPSbFxMRIHh4eUvv27aV//OMfUnl5uX0b7mfHbdq0qd6/yZMnT5YkqXn79MKFC9K9994r+fj4SH5+ftLUqVOloqIiWeo1SdJVpyUkIiIi0jiueSEiIiKhMLwQERGRUBheiIiISCgML0RERCQUhhciIiISCsMLERERCYXhhYiIiITC8EJERERCYXghIqeYMmUKxo0bp9rzT5w40X5FZzkcPXoUbdu2RUlJiWzPQUTNwzPsElGTTCZToz9/7rnn8Pjjj0OSJAQEBChT1FUOHDiAW2+9FWlpafDx8ZHtef7whz8gISEBzz77rGzPQURNY3ghoiZlZ2fb/79ixQrMmzevzhWUfXx8ZA0NTZk2bRpcXV2xaNEiWZ/nu+++w/Tp05Geng5XV1dZn4uIGsZpIyJqUnh4uP3L398fJpOpzm0+Pj7XTRsNHToUjzzyCGbPno02bdogLCwM77//PkpKSjB16lT4+vqiY8eO+OGHH+o81+HDhzFq1Cj4+PggLCwMEydOxPnz5xusrbq6Gl999RXGjBlT5/Z27drhxRdfxKRJk+Dj44PY2FisWbMGeXl5GDt2LHx8fNCrVy/s2bPHfp+0tDSMGTMGbdq0gbe3N7p3747vv//e/vPbb78d+fn52LJlSyv3KBG1BsMLEcnm448/RnBwMHbt2oVHHnkEDz/8MMaPH4+BAwdi3759GD58OCZOnIjS0lIAQEFBAW699VYkJiZiz549WLt2LXJycvDHP/6xwec4ePAgCgsL0a9fv+t+tmDBAgwaNAj79+/H6NGjMXHiREyaNAn3338/9u3bhw4dOmDSpEmoGYCeOXMmysvLsXXrVhw6dAivvPJKnREld3d39O7dGz///LOT9xQROUSWa1UTkW599NFHkr+//3W3T548WRo7dqz9+yFDhkiDBw+2f19VVSV5e3tLEydOtN+WlZUlAZCSk5MlSZKkF154QRo+fHidx83IyJAASKmpqfXWs2rVKsnFxUWy2Wx1bo+NjZXuv//+657r2Weftd+WnJwsAZCysrIkSZKknj17Sv/85z8bbf/dd98tTZkypdFtiEheHHkhItn06tXL/n8XFxcEBQWhZ8+e9tvCwsIAALm5uQAuL7zdtGmTfQ2Nj48P4uPjAQCnTp2q9znKyspgsVjqXVR89fPXPFdjz//oo4/ixRdfxKBBg/Dcc8/h4MGD1z2mp6enfaSIiNTB8EJEsnFzc6vzvclkqnNbTeCw2WwAgOLiYowZMwYpKSl1vk6cOIGbb7653ucIDg5GaWkpKioqGn3+mudq7PmnTZuG06dPY+LEiTh06BD69euHt956q85j5ufnIyQkpHk7gIhkwfBCRJrRp08fHDlyBO3atUPHjh3rfHl7e9d7n969ewO4fB4WZ4iOjsaMGTOwcuVKPPHEE3j//ffr/Pzw4cNITEx0ynMRUcswvBCRZsycORP5+fm49957sXv3bpw6dQrr1q3D1KlTUV1dXe99QkJC0KdPH2zbtq3Vzz979mysW7cOZ86cwb59+7Bp0yZ07drV/vPffvsNmZmZGDZsWKufi4hajuGFiDQjMjIS27dvR3V1NYYPH46ePXti9uzZCAgIgNnc8J+radOmYenSpa1+/urqasycORNdu3bFyJEj0blzZ7zzzjv2n3/++ecYPnw4YmNjW/1cRNRyPEkdEQmvrKwMXbp0wYoVK5CUlCTLc1RUVKBTp05YtmwZBg0aJMtzEFHzcOSFiITn6emJTz75pNGT2bVWeno6nn76aQYXIg3gyAsREREJhSMvREREJBSGFyIiIhIKwwsREREJheGFiIiIhMLwQkREREJheCEiIiKhMLwQERGRUBheiIiISCgML0RERCSU/w+P1wr+4e6GGwAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plt.figure()\n", "plt.plot(t, mvec, label = \"m\")\n", "plt.legend(loc = 1)\n", "plt.xlabel('Time (ms)')\n", "plt.savefig('Fig1.jpg')\n", "#plt.show(0)\n", "\n", "plt.figure()\n", "plt.plot(t, Mvec, label = \"M\")\n", "plt.legend(loc = 1)\n", "plt.xlabel('Time (ms)')\n", "plt.savefig('Fig2.jpg')\n", "#plt.show(0)" ] }, { "cell_type": "code", "execution_count": 11, "id": "1b58c0a8-8892-4461-a56e-6e4bcb1d8efc", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjwAAAGwCAYAAACtlb+kAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8hTgPZAAAACXBIWXMAAA9hAAAPYQGoP6dpAAAp1ElEQVR4nO3de3QUZZ7G8aeTkIYQ0uGWCxAgCmu4E0Aw4BFmjFzkcMyMyyILBFjAw2yiRNxxxF1AxZkMOogziCAOyLjDTRxBRcWNgcCi4U4EFPGGhIV0ADFpEzBAUvsHhx5aQkiGVHd48/2cU+ek3nqr3l/qjfRjdXWXw7IsSwAAAAYLCnQBAAAAdiPwAAAA4xF4AACA8Qg8AADAeAQeAABgPAIPAAAwHoEHAAAYLyTQBdS2iooKnThxQk2aNJHD4Qh0OQAAoBosy9IPP/ygVq1aKSio9q/HGBd4Tpw4obi4uECXAQAA/gHHjh1TmzZtav24xgWeJk2aSLp0wiIiIgJcDQAAqA6Px6O4uDjv63htMy7wXH4bKyIigsADAMBNxq7bUbhpGQAAGI/AAwAAjEfgAQAAxjPuHh4AAGpTeXm5Lly4EOgyjBAaGmrLR86rg8ADAEAlLMuS2+1WUVFRoEsxRlBQkOLj4xUaGur3sQk8AABU4nLYiYqKUlhYGF9me4MufzFwQUGB2rZt6/fzSeABAOAnysvLvWGnefPmgS7HGC1bttSJEyd08eJFNWjQwK9jc9MyAAA/cfmenbCwsABXYpbLb2WVl5f7fWwCDwAA18DbWLUrkOeTwAMAAIxna+DJzMzU7bffriZNmigqKkopKSk6fPhwlfssX75cDofDZ2nYsKGdZQIAAMPZGni2bNmitLQ0bd++XVlZWbpw4YIGDx6s0tLSKveLiIhQQUGBdzl69KidZQIAAMPZGng2btyoCRMmqEuXLurRo4eWL1+u/Px87dmzp8r9HA6HYmJivEt0dLSdZQIAYIwJEyYoJSXFpy03N1fBwcEaPnz4Vf2//fbbq95ZcTgcGjt2rJ8q9g+/fiy9uLhYktSsWbMq+5WUlKhdu3aqqKhQr1699Lvf/U5dunSptG9ZWZnKysq86x6Pp/YKBgDAAEuXLtVDDz2kpUuX6sSJE2rVqtVVfT788EOf19pGjRr5s0Tb+S3wVFRUKCMjQwMGDFDXrl2v2e+2227TsmXL1L17dxUXF+sPf/iD+vfvr08//VRt2rS5qn9mZqaeeuop2+ouLZXCwy/9XFIiNW5s21B1GufBrHNg0u8C+Nt17sqodTf632dJSYnWrFmj3bt3y+12a/ny5XriiSeu6te8eXPFxMTc2GB1mN8CT1pamg4ePKht27ZV2S8pKUlJSUne9f79+6tTp056+eWXNWfOnKv6z5gxQ9OnT/euezwexcXF1V7hAABc4fL/LPiLZd3Y/q+//roSEhJ02223aezYscrIyNCMGTPq3Ufu/fKx9PT0dG3YsEGbN2+u9CpNVRo0aKDExER99dVXlW53Op2KiIjwWQAAwCVLly713o8zdOhQFRcXa8uWLVf169+/v8LDw73Lvn37/F2qrWy9wmNZlh566CGtW7dOOTk5io+Pr/ExysvLdeDAAd177702VAgAQM2UlAS6guo7fPiwdu7cqXXr1kmSQkJCNGrUKC1dulSDBg3y6btmzRp16tTJu27auyW2Bp60tDStXLlSb731lpo0aSK32y1Jcrlc3puhUlNT1bp1a2VmZkqSnn76ad1xxx3q0KGDioqK9Nxzz+no0aOaPHmynaUCAFAtN9M9b0uXLtXFixd9blK2LEtOp1MvvviiXC6Xtz0uLk4dOnQIRJl+YWvgWbRokSRdlSJfffVVTZgwQZKUn5+voKC/v7P2/fffa8qUKXK73WratKl69+6tjz/+WJ07d7azVAAAjHLx4kW99tprmjdvngYPHuyzLSUlRatWrdLUqVMDVJ3/2f6W1vXk5OT4rM+fP1/z58+3qSIAAOqHDRs26Pvvv9ekSZN8ruRI0v3336+lS5fWq8DDs7QAADDQ0qVLlZycfFXYkS4Fnt27d2v//v0BqCww/PrFgwAAwF7Lly+/bp++ffv6vAtTnXdkbnZc4QEAAMYj8AAAAOMReAAAgPEIPAAAwHgEHgAAYDwCDwAAMB6BBwAAGI/AAwAAjEfgAQAAxiPwAABgkAkTJiglJcWvYy5fvlyRkZF+HbOmCDwAAMB4BB4AAAw1aNAgPfzww3rsscfUrFkzxcTE6Mknn/Tp43A4tGjRIg0bNkyNGjXSLbfcojfeeMO7PScnRw6HQ0VFRd62vLw8ORwOffvtt8rJydHEiRNVXFwsh8Mhh8Nx1Rh1AYEHAICaKC3173KD/vKXv6hx48basWOHnn32WT399NPKysry6TNz5kzdf//9+uSTTzRmzBg98MADOnToULWO379/f73wwguKiIhQQUGBCgoK9B//8R83XHdt42npAADURHi4f8e7wSeZd+/eXbNnz5YkdezYUS+++KKys7N1zz33ePuMHDlSkydPliTNmTNHWVlZWrBggV566aXrHj80NFQul0sOh0MxMTE3VKuduMIDAIDBunfv7rMeGxurkydP+rQlJSVdtV7dKzw3C67wAABQEyUlga6gRho0aOCz7nA4VFFRUe39g4IuXRuxrrjSdOHChdopzo8IPAAA1ETjxoGuoNZt375dqampPuuJiYmSpJYtW0qSCgoK1LRpU0mXblq+UmhoqMrLy/1T7D+It7QAAKjn1q5dq2XLlumLL77Q7NmztXPnTqWnp0uSOnTooLi4OD355JP68ssv9e6772revHk++7dv314lJSXKzs7W6dOndfbs2UD8GlUi8AAAUM899dRTWr16tbp3767XXntNq1atUufOnSVdekts1apV+vzzz9W9e3fNnTtXzzzzjM/+/fv319SpUzVq1Ci1bNlSzz77bCB+jSo5LOsGb/+uYzwej1wul4qLixUREXHDxyst/fsN+SUlRl7JrBbOg1nnwKTfBbDDjz/+qCNHjig+Pl4NGzYMdDm2cjgcWrdunV++nbmq81rbr98/xRUeAABgPAIPAAAwHp/SAgCgHjPszpZr4goPAAAwHoEHAIBrqC9XP/wlkOeTwAMAwE9c/nbiuvh9Mjez8+fPS5KCg4P9Pjb38AAA8BPBwcGKjIz0PnMqLCxMDocjwFXd3CoqKnTq1CmFhYUpJMT/8YPAAwBAJS4/+funD9rEPy4oKEht27YNSHgk8AAAUAmHw6HY2FhFRUXdlA/LrItCQ0O9DyP1NwIPAABVCA4ODsg9J6hd3LQMAACMR+ABAADGI/AAAADjEXgAAIDxCDwAAMB4BB4AAGA8Ag8AADAegQcAABiPwAMAAIxH4AEAAMYj8AAAAOMReAAAgPEIPAAAwHi2Bp7MzEzdfvvtatKkiaKiopSSkqLDhw9fd7+1a9cqISFBDRs2VLdu3fTee+/ZWSYAADCcrYFny5YtSktL0/bt25WVlaULFy5o8ODBKi0tveY+H3/8sUaPHq1JkyZp3759SklJUUpKig4ePGhnqQAAwGAOy7Isfw126tQpRUVFacuWLbrrrrsq7TNq1CiVlpZqw4YN3rY77rhDPXv21OLFi687hsfjkcvlUnFxsSIiIm645tJSKTz80s8lJVLjxjd8yJsS58Gsc2DS7wLADLX9+v1Tfr2Hp7i4WJLUrFmza/bJzc1VcnKyT9uQIUOUm5tbaf+ysjJ5PB6fBQAA4Ep+CzwVFRXKyMjQgAED1LVr12v2c7vdio6O9mmLjo6W2+2utH9mZqZcLpd3iYuLq9W6AQDAzc9vgSctLU0HDx7U6tWra/W4M2bMUHFxsXc5duxYrR4fAADc/EL8MUh6ero2bNigrVu3qk2bNlX2jYmJUWFhoU9bYWGhYmJiKu3vdDrldDprrVYAAGAeW6/wWJal9PR0rVu3Tps2bVJ8fPx190lKSlJ2drZPW1ZWlpKSkuwqEwAAGM7WKzxpaWlauXKl3nrrLTVp0sR7H47L5VKjRo0kSampqWrdurUyMzMlSdOmTdPAgQM1b948DR8+XKtXr9bu3bu1ZMkSO0sFAAAGs/UKz6JFi1RcXKxBgwYpNjbWu6xZs8bbJz8/XwUFBd71/v37a+XKlVqyZIl69OihN954Q+vXr6/yRmcAAICq2HqFpzpf8ZOTk3NV28iRIzVy5EgbKgIAAPURz9ICAADGI/AAAADjEXgAAIDxCDwAAMB4BB4AAGA8Ag8AADAegQcAABiPwAMAAIxH4AEAAMYj8AAAAOMReAAAgPEIPAAAwHgEHgAAYDwCDwAAMB6BBwAAGI/AAwAAjEfgAQAAxiPwAAAA4xF4AACA8Qg8AADAeAQeAABgPAIPAAAwHoEHAAAYj8ADAACMR+ABAADGI/AAAADjEXgAAIDxCDwAAMB4BB4AAGA8Ag8AADAegQcAABiPwAMAAIxH4AEAAMYj8AAAAOMReAAAgPEIPAAAwHgEHgAAYDwCDwAAMB6BBwAAGI/AAwAAjEfgAQAAxiPwAAAA4xF4AACA8Qg8AADAeAQeAABgPFsDz9atWzVixAi1atVKDodD69evr7J/Tk6OHA7HVYvb7bazTAAAYDhbA09paal69OihhQsX1mi/w4cPq6CgwLtERUXZVCEAAKgPQuw8+LBhwzRs2LAa7xcVFaXIyMjaLwgAANRLdfIenp49eyo2Nlb33HOPPvrooyr7lpWVyePx+CwAAABXqlOBJzY2VosXL9bf/vY3/e1vf1NcXJwGDRqkvXv3XnOfzMxMuVwu7xIXF1erNTVuLFnWpaVx41o99E2F82DWOTDpdwGA6nBYlmX5ZSCHQ+vWrVNKSkqN9hs4cKDatm2r//7v/650e1lZmcrKyrzrHo9HcXFxKi4uVkRExI2UDAAA/MTj8cjlctn2+m3rPTy1oW/fvtq2bds1tzudTjmdTj9WBAAAbjZ16i2tyuTl5Sk2NjbQZQAAgJuYrVd4SkpK9NVXX3nXjxw5ory8PDVr1kxt27bVjBkzdPz4cb322muSpBdeeEHx8fHq0qWLfvzxR/35z3/Wpk2b9D//8z92lgkAAAxna+DZvXu3fvazn3nXp0+fLkkaP368li9froKCAuXn53u3nz9/Xo8++qiOHz+usLAwde/eXR9++KHPMQAAAGrKbzct+4vdNz0BAIDaZ/frd52/hwcAAOBGEXgAAIDxCDwAAMB4BB4AAGA8Ag8AADAegQcAABiPwAMAAIxH4AEAAMYj8AAAAOMReAAAgPEIPAAAwHgEHgAAYDwCDwAAMB6BBwAAGI/AAwAAjEfgAQAAxiPwAAAA4xF4AACA8Qg8AADAeAQeAABgPAIPAAAwHoEHAAAYj8ADAACMR+ABAADGI/AAAADjEXgAAIDxCDwAAMB4BB4AAGA8Ag8AADAegQcAABiPwAMAAIxH4AEAAMYj8AAAAOMReAAAgPEIPAAAwHgEHgAAYDwCDwAAMB6BBwAAGI/AAwAAjEfgAQAAxiPwAAAA4xF4AACA8Qg8AADAeAQeAABgPAIPAAAwnq2BZ+vWrRoxYoRatWolh8Oh9evXX3efnJwc9erVS06nUx06dNDy5cvtLBEAANQDtgae0tJS9ejRQwsXLqxW/yNHjmj48OH62c9+pry8PGVkZGjy5Mn64IMP7CwTAAAYLsTOgw8bNkzDhg2rdv/FixcrPj5e8+bNkyR16tRJ27Zt0/z58zVkyBC7ygQAAIarU/fw5ObmKjk52adtyJAhys3NveY+ZWVl8ng8PgsAAMCV6lTgcbvdio6O9mmLjo6Wx+PRuXPnKt0nMzNTLpfLu8TFxfmjVAAAcBOpU4HnHzFjxgwVFxd7l2PHjgW6JAAAUMfYeg9PTcXExKiwsNCnrbCwUBEREWrUqFGl+zidTjmdTn+UBwAAblJ16gpPUlKSsrOzfdqysrKUlJQUoIoAAIAJbA08JSUlysvLU15enqRLHzvPy8tTfn6+pEtvR6Wmpnr7T506Vd98840ee+wxff7553rppZf0+uuv65FHHrGzTAAAYDhbA8/u3buVmJioxMRESdL06dOVmJioWbNmSZIKCgq84UeS4uPj9e677yorK0s9evTQvHnz9Oc//5mPpAMAgBvisCzLCnQRtcnj8cjlcqm4uFgRERGBLgcAAFSD3a/fdeoeHgAAADsQeAAAgPEIPAAAwHgEHgAAYDwCDwAAMB6BBwAAGI/AAwAAjEfgAQAAxiPwAAAA4xF4AACA8Qg8AADAeAQeAABgPAIPAAAwHoEHAAAYj8ADAACMR+ABAADGI/AAAADjEXgAAIDxCDwAAMB4BB4AAGA8Ag8AADAegQcAABiPwAMAAIxH4AEAAMYj8AAAAOMReAAAgPEIPAAAwHgEHgAAYDwCDwAAMB6BBwAAGI/AAwAAjEfgAQAAxiPwAAAA4xF4AACA8Qg8AADAeAQeAABgPAIPAAAwHoEHAAAYj8ADAACMR+ABAADGI/AAAADjEXgAAIDxCDwAAMB4BB4AAGA8Ag8AADCeXwLPwoUL1b59ezVs2FD9+vXTzp07r9l3+fLlcjgcPkvDhg39USYAADCU7YFnzZo1mj59umbPnq29e/eqR48eGjJkiE6ePHnNfSIiIlRQUOBdjh49aneZAADAYLYHnueff15TpkzRxIkT1blzZy1evFhhYWFatmzZNfdxOByKiYnxLtHR0XaXCQAADGZr4Dl//rz27Nmj5OTkvw8YFKTk5GTl5uZec7+SkhK1a9dOcXFxuu+++/Tpp59es29ZWZk8Ho/PAgAAcCVbA8/p06dVXl5+1RWa6Ohoud3uSve57bbbtGzZMr311lv661//qoqKCvXv31//93//V2n/zMxMuVwu7xIXF3dpg8sllZZWXWBpqeRwXFpKS69e/2mfK9uvdQy73Mg41/sd7By/Ns5PVcew+/jX6udwSCdP+q7XpIbams/K/k6vrOtax/bX361JanLOOL9AnVPnPqWVlJSk1NRU9ezZUwMHDtSbb76pli1b6uWXX660/4wZM1RcXOxdjh075ueKAQBAXRdi58FbtGih4OBgFRYW+rQXFhYqJiamWsdo0KCBEhMT9dVXX1W63el0yul03nCtAADAXLZe4QkNDVXv3r2VnZ3tbauoqFB2draSkpKqdYzy8nIdOHBAsbGxdpUJAAAMZ+sVHkmaPn26xo8frz59+qhv37564YUXVFpaqokTJ0qSUlNT1bp1a2VmZkqSnn76ad1xxx3q0KGDioqK9Nxzz+no0aOaPHmy3aUCAABD2R54Ro0apVOnTmnWrFlyu93q2bOnNm7c6L2ROT8/X0FBf7/Q9P3332vKlClyu91q2rSpevfurY8//lidO3e2u1QAAGAo2wOPJKWnpys9Pb3SbTk5OT7r8+fP1/z58/1QFQAAqC/q3Ke0AAAAahuBBwAAGI/AAwAAjEfgAQAAxiPwAAAA4xF4AACA8Qg8AADAeAQeAABgPAIPAAAwHoEHAAAYj8ADAACMR+ABAADGI/AAAADjEXgAAIDxCDwAAMB4BB4AAGA8Ag8AADAegQcAABiPwAMAAIxH4AEAAMYj8AAAAOMReAAAgPEIPAAAwHgEHgAAYDwCDwAAMB6BBwAAGI/AAwAAjEfgAQAAxiPwAAAA4xF4AACA8Qg8AADAeAQeAABgPAIPAAAwHoEHAAAYj8ADAACMR+ABAADGI/AAAADjEXgAAIDxCDwAAMB4BB4AAGA8Ag8AADAegQcAABiPwAMAAIxH4AEAAMYj8AAAAOP5JfAsXLhQ7du3V8OGDdWvXz/t3Lmzyv5r165VQkKCGjZsqG7duum9997zR5kAAMBQtgeeNWvWaPr06Zo9e7b27t2rHj16aMiQITp58mSl/T/++GONHj1akyZN0r59+5SSkqKUlBQdPHjQ7lIBAIChbA88zz//vKZMmaKJEyeqc+fOWrx4scLCwrRs2bJK+//xj3/U0KFD9etf/1qdOnXSnDlz1KtXL7344ot2lwoAAAxla+A5f/689uzZo+Tk5L8PGBSk5ORk5ebmVrpPbm6uT39JGjJkyDX7l5WVyePx+CwAAABXsjXwnD59WuXl5YqOjvZpj46OltvtrnQft9tdo/6ZmZlyuVzeJS4urnaKBwAAxrjpP6U1Y8YMFRcXe5djx44FuiQAAFDHhNh58BYtWig4OFiFhYU+7YWFhYqJial0n5iYmBr1dzqdcjqdtVMwAAAwkq1XeEJDQ9W7d29lZ2d72yoqKpSdna2kpKRK90lKSvLpL0lZWVnX7A8AAHA9tl7hkaTp06dr/Pjx6tOnj/r27asXXnhBpaWlmjhxoiQpNTVVrVu3VmZmpiRp2rRpGjhwoObNm6fhw4dr9erV2r17t5YsWWJ3qQAAwFC2B55Ro0bp1KlTmjVrltxut3r27KmNGzd6b0zOz89XUNDfLzT1799fK1eu1H/913/piSeeUMeOHbV+/Xp17drV7lIBAIChbA88kpSenq709PRKt+Xk5FzVNnLkSI0cOdLmqgAAQH1x039KCwAA4HoIPAAAwHgEHgAAYDwCDwAAMB6BBwAAGI/AAwAAjEfgAQAAxiPwAAAA4xF4AACA8Qg8AADAeAQeAABgPAIPAAAwHoEHAAAYj8ADAACMR+ABAADGI/AAAADjEXgAAIDxCDwAAMB4BB4AAGA8Ag8AADAegQcAABiPwAMAAIxH4AEAAMYj8AAAAOMReAAAgPEIPAAAwHgEHgAAYDwCDwAAMB6BBwAAGI/AAwAAjEfgAQAAxiPwAAAA4xF4AACA8Qg8AADAeAQeAABgPAIPAAAwHoEHAAAYj8ADAACMR+ABAADGI/AAAADjEXgAAIDxCDwAAMB4BB4AAGA8Ag8AADAegQcAABjPtsBz5swZjRkzRhEREYqMjNSkSZNUUlJS5T6DBg2Sw+HwWaZOnWpXiQAAoJ4IsevAY8aMUUFBgbKysnThwgVNnDhRDz74oFauXFnlflOmTNHTTz/tXQ8LC7OrRAAAUE/YEngOHTqkjRs3ateuXerTp48kacGCBbr33nv1hz/8Qa1atbrmvmFhYYqJibGjLAAAUE/Z8pZWbm6uIiMjvWFHkpKTkxUUFKQdO3ZUue+KFSvUokULde3aVTNmzNDZs2er7F9WViaPx+OzAAAAXMmWKzxut1tRUVG+A4WEqFmzZnK73dfc71//9V/Vrl07tWrVSvv379dvfvMbHT58WG+++eY198nMzNRTTz119YbiYqlx46oLbdxYsizftp+uV9anJttry42MUxs1/qPHsHtsf/5u1fl7qe0xq7vvT9uud2x//d2apCbnjPML1Dk1CjyPP/645s6dW2WfQ4cO/cPFPPjgg96fu3XrptjYWN199936+uuvdeutt1a6z4wZMzR9+nTvusfjUVxc3D9cAwAAME+NAs+jjz6qCRMmVNnnlltuUUxMjE6ePOnTfvHiRZ05c6ZG9+f069dPkvTVV19dM/A4nU45nc5qHxMAANQ/NQo8LVu2VMuWLa/bLykpSUVFRdqzZ4969+4tSdq0aZMqKiq8IaY68vLyJEmxsbE1KRMAAMCHLTctd+rUSUOHDtWUKVO0c+dOffTRR0pPT9cDDzzg/YTW8ePHlZCQoJ07d0qSvv76a82ZM0d79uzRt99+q7ffflupqam666671L17dzvKBAAA9YRtXzy4YsUKJSQk6O6779a9996rO++8U0uWLPFuv3Dhgg4fPuz9FFZoaKg+/PBDDR48WAkJCXr00Ud1//3365133rGrRAAAUE84LMusjxJ4PB65XC4VFxcrIiIi0OUAAIBqsPv1m2dpAQAA4xF4AACA8Qg8AADAeAQeAABgPAIPAAAwHoEHAAAYj8ADAACMR+ABAADGI/AAAADj1ejhoTeDy18c7fF4AlwJAACorsuv23Y9AMK4wPPdd99JkuLi4gJcCQAAqKnvvvtOLper1o9rXOBp1qyZJCk/P9+WE4bq83g8iouL07Fjx3iuWR3AfNQdzEXdwVzUHcXFxWrbtq33dby2GRd4goIu3Zbkcrn4460jIiIimIs6hPmoO5iLuoO5qDsuv47X+nFtOSoAAEAdQuABAADGMy7wOJ1OzZ49W06nM9Cl1HvMRd3CfNQdzEXdwVzUHXbPhcOy6/NfAAAAdYRxV3gAAAB+isADAACMR+ABAADGI/AAAADjGRd4Fi5cqPbt26thw4bq16+fdu7cGeiSjJeZmanbb79dTZo0UVRUlFJSUnT48GGfPj/++KPS0tLUvHlzhYeH6/7771dhYWGAKq4/fv/738vhcCgjI8Pbxlz4z/HjxzV27Fg1b95cjRo1Urdu3bR7927vdsuyNGvWLMXGxqpRo0ZKTk7Wl19+GcCKzVReXq6ZM2cqPj5ejRo10q233qo5c+b4PLOJubDP1q1bNWLECLVq1UoOh0Pr16/32V6dc3/mzBmNGTNGERERioyM1KRJk1RSUlKjOowKPGvWrNH06dM1e/Zs7d27Vz169NCQIUN08uTJQJdmtC1btigtLU3bt29XVlaWLly4oMGDB6u0tNTb55FHHtE777yjtWvXasuWLTpx4oR++ctfBrBq8+3atUsvv/yyunfv7tPOXPjH999/rwEDBqhBgwZ6//339dlnn2nevHlq2rSpt8+zzz6rP/3pT1q8eLF27Nihxo0ba8iQIfrxxx8DWLl55s6dq0WLFunFF1/UoUOHNHfuXD377LNasGCBtw9zYZ/S0lL16NFDCxcurHR7dc79mDFj9OmnnyorK0sbNmzQ1q1b9eCDD9asEMsgffv2tdLS0rzr5eXlVqtWrazMzMwAVlX/nDx50pJkbdmyxbIsyyoqKrIaNGhgrV271tvn0KFDliQrNzc3UGUa7YcffrA6duxoZWVlWQMHDrSmTZtmWRZz4U+/+c1vrDvvvPOa2ysqKqyYmBjrueee87YVFRVZTqfTWrVqlT9KrDeGDx9u/du//ZtP2y9/+UtrzJgxlmUxF/4kyVq3bp13vTrn/rPPPrMkWbt27fL2ef/99y2Hw2EdP3682mMbc4Xn/Pnz2rNnj5KTk71tQUFBSk5OVm5ubgArq3+Ki4sl/f1Brnv27NGFCxd85iYhIUFt27ZlbmySlpam4cOH+5xzibnwp7ffflt9+vTRyJEjFRUVpcTERL3yyive7UeOHJHb7faZC5fLpX79+jEXtax///7Kzs7WF198IUn65JNPtG3bNg0bNkwScxFI1Tn3ubm5ioyMVJ8+fbx9kpOTFRQUpB07dlR7LGMeHnr69GmVl5crOjrapz06Olqff/55gKqqfyoqKpSRkaEBAwaoa9eukiS3263Q0FBFRkb69I2Ojpbb7Q5AlWZbvXq19u7dq127dl21jbnwn2+++UaLFi3S9OnT9cQTT2jXrl16+OGHFRoaqvHjx3vPd2X/ZjEXtevxxx+Xx+NRQkKCgoODVV5ert/+9rcaM2aMJDEXAVSdc+92uxUVFeWzPSQkRM2aNavR/BgTeFA3pKWl6eDBg9q2bVugS6mXjh07pmnTpikrK0sNGzYMdDn1WkVFhfr06aPf/e53kqTExEQdPHhQixcv1vjx4wNcXf3y+uuva8WKFVq5cqW6dOmivLw8ZWRkqFWrVsxFPWLMW1otWrRQcHDwVZ82KSwsVExMTICqql/S09O1YcMGbd68WW3atPG2x8TE6Pz58yoqKvLpz9zUvj179ujkyZPq1auXQkJCFBISoi1btuhPf/qTQkJCFB0dzVz4SWxsrDp37uzT1qlTJ+Xn50uS93zzb5b9fv3rX+vxxx/XAw88oG7dumncuHF65JFHlJmZKYm5CKTqnPuYmJirPnx08eJFnTlzpkbzY0zgCQ0NVe/evZWdne1tq6ioUHZ2tpKSkgJYmfksy1J6errWrVunTZs2KT4+3md779691aBBA5+5OXz4sPLz85mbWnb33XfrwIEDysvL8y59+vTRmDFjvD8zF/4xYMCAq76e4YsvvlC7du0kSfHx8YqJifGZC4/Hox07djAXtezs2bMKCvJ9uQsODlZFRYUk5iKQqnPuk5KSVFRUpD179nj7bNq0SRUVFerXr1/1B7vhW67rkNWrV1tOp9Navny59dlnn1kPPvigFRkZabnd7kCXZrRf/epXlsvlsnJycqyCggLvcvbsWW+fqVOnWm3btrU2bdpk7d6920pKSrKSkpICWHX9ceWntCyLufCXnTt3WiEhIdZvf/tb68svv7RWrFhhhYWFWX/961+9fX7/+99bkZGR1ltvvWXt37/fuu+++6z4+Hjr3LlzAazcPOPHj7dat25tbdiwwTpy5Ij15ptvWi1atLAee+wxbx/mwj4//PCDtW/fPmvfvn2WJOv555+39u3bZx09etSyrOqd+6FDh1qJiYnWjh07rG3btlkdO3a0Ro8eXaM6jAo8lmVZCxYssNq2bWuFhoZaffv2tbZv3x7okownqdLl1Vdf9fY5d+6c9e///u9W06ZNrbCwMOsXv/iFVVBQELii65GfBh7mwn/eeecdq2vXrpbT6bQSEhKsJUuW+GyvqKiwZs6caUVHR1tOp9O6++67rcOHDweoWnN5PB5r2rRpVtu2ba2GDRtat9xyi/Wf//mfVllZmbcPc2GfzZs3V/oaMX78eMuyqnfuv/vuO2v06NFWeHi4FRERYU2cONH64YcfalSHw7Ku+KpJAAAAAxlzDw8AAMC1EHgAAIDxCDwAAMB4BB4AAGA8Ag8AADAegQcAABiPwAMAAIxH4AEAAMYj8ACoFRMmTFBKSkrAxh83bpz3yeR2+Oyzz9SmTRuVlpbaNgYA+/BNywCuy+FwVLl99uzZeuSRR2RZliIjI/1T1BU++eQT/fznP9fRo0cVHh5u2zj//M//rB49emjmzJm2jQHAHgQeANfldru9P69Zs0azZs3yeRJ4eHi4rUHjeiZPnqyQkBAtXrzY1nHeffddTZkyRfn5+QoJCbF1LAC1i7e0AFxXTEyMd3G5XHI4HD5t4eHhV72lNWjQID300EPKyMhQ06ZNFR0drVdeeUWlpaWaOHGimjRpog4dOuj999/3GevgwYMaNmyYwsPDFR0drXHjxun06dPXrK28vFxvvPGGRowY4dPevn17PfPMM0pNTVV4eLjatWunt99+W6dOndJ9992n8PBwde/eXbt37/buc/ToUY0YMUJNmzZV48aN1aVLF7333nve7ffcc4/OnDmjLVu23OAZBeBvBB4AtvnLX/6iFi1aaOfOnXrooYf0q1/9SiNHjlT//v21d+9eDR48WOPGjdPZs2clSUVFRfr5z3+uxMRE7d69Wxs3blRhYaH+5V/+5Zpj7N+/X8XFxerTp89V2+bPn68BAwZo3759Gj58uMaNG6fU1FSNHTtWe/fu1a233qrU1FRdvtCdlpamsrIybd26VQcOHNDcuXN9rlyFhoaqZ8+e+t///d9aPlMAbHejj30HUL+8+uqrlsvluqp9/Pjx1n333eddHzhwoHXnnXd61y9evGg1btzYGjdunLetoKDAkmTl5uZalmVZc+bMsQYPHuxz3GPHjlmSrMOHD1daz7p166zg4GCroqLCp71du3bW2LFjrxpr5syZ3rbc3FxLklVQUGBZlmV169bNevLJJ6v8/X/xi19YEyZMqLIPgLqHKzwAbNO9e3fvz8HBwWrevLm6devmbYuOjpYknTx5UtKlm483b97svScoPDxcCQkJkqSvv/660jHOnTsnp9NZ6Y3VV45/eayqxn/44Yf1zDPPaMCAAZo9e7b2799/1TEbNWrkvSIF4OZB4AFgmwYNGvisOxwOn7bLIaWiokKSVFJSohEjRigvL89n+fLLL3XXXXdVOkaLFi109uxZnT9/vsrxL49V1fiTJ0/WN998o3HjxunAgQPq06ePFixY4HPMM2fOqGXLltU7AQDqDAIPgDqjV69e+vTTT9W+fXt16NDBZ2ncuHGl+/Ts2VPSpe/JqQ1xcXGaOnWq3nzzTT366KN65ZVXfLYfPHhQiYmJtTIWAP8h8ACoM9LS0nTmzBmNHj1au3bt0tdff60PPvhAEydOVHl5eaX7tGzZUr169dK2bdtuePyMjAx98MEHOnLkiPbu3avNmzerU6dO3u3ffvutjh8/ruTk5BseC4B/EXgA1BmtWrXSRx99pPLycg0ePFjdunVTRkaGIiMjFRR07X+uJk+erBUrVtzw+OXl5UpLS1OnTp00dOhQ/dM//ZNeeukl7/ZVq1Zp8ODBateu3Q2PBcC/+OJBADe9c+fO6bbbbtOaNWuUlJRkyxjnz59Xx44dtXLlSg0YMMCWMQDYhys8AG56jRo10muvvVblFxTeqPz8fD3xxBOEHeAmxRUeAABgPK7wAAAA4xF4AACA8Qg8AADAeAQeAABgPAIPAAAwHoEHAAAYj8ADAACMR+ABAADGI/AAAADj/T9TFLYFRDCK0AAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig, ax = plt.subplots()\n", "for times, id in zip(spike_times_vec_input, idvec_input):\n", " Input = ax.vlines(times, id+0.5,id+1.5,'Red')\n", "for times, id in zip(spike_times_vec_IAF, idvec_IAF):\n", " IAF_cell = ax.vlines(times, id+0.5,id+1.5,'Blue')\n", "ax.legend([IAF_cell, Input], ['IAF', 'Input'])\n", "ax.set_xlabel('Time (ms)')\n", "ax.set_xlim([0, tstop])\n", "plt.savefig('Fig3.jpg')\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": 12, "id": "5f742ab3-d5e7-49e7-b658-3b56f4378088", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "0.025\t0.05\t0.075\t0.1\t0.125\t\n", "0.15\t0.175\t0.2\t0.225\t0.25\t\n", "\n" ] }, { "data": { "text/plain": [ "10" ] }, "execution_count": 12, "metadata": {}, "output_type": "execute_result" } ], "source": [ "t.printf(1,10)" ] }, { "cell_type": "code", "execution_count": null, "id": "7647ff70-804f-4199-8a1b-a6f9d33457d2", "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.11.5" } }, "nbformat": 4, "nbformat_minor": 5 }