{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "## Regulation of ion transport from within ion transit pathways\n", "-------------" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import matplotlib.pyplot as plt\n", "import seaborn as sns\n", "from scipy import linalg\n", "from scipy.integrate import odeint\n", "sns.set()\n", "%matplotlib inline" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Detail" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ " " ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "Ni = 100*1e-3\n", "Ci = 2*1e-6\n", "F_3n = Ni**2.5/(Ni**2.5 + (17*1e-3)**2.5)\n", "K_con1 = 0.1*1e-6\n", "K_con2 = 20\n", "K_coff1 = 0.05*1e-6\n", "K_coff2 = 0.3\n", "K_12 = 0.2\n", "K_43 = 0.1\n", "K_inact1 = 0.2\n", "K_inact2 = 25*K_inact1" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "k = {\"Ni\": Ni,\"Ci\":Ci, \"F_3n\": F_3n,\\\n", " \"K_con1\": K_con1,\"K_con2\":K_con2,\\\n", " \"K_coff1\":K_coff1,\"K_coff2\": K_coff2,\\\n", " \"K_12\": K_12,\"K_43\":K_43,\"K_inact1\": K_inact1,\\\n", " \"K_inact2\": K_inact2}" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "{'Ni': 0.1,\n", " 'Ci': 2e-06,\n", " 'F_3n': 0.9882245384948203,\n", " 'K_con1': 1e-07,\n", " 'K_con2': 20,\n", " 'K_coff1': 5e-08,\n", " 'K_coff2': 0.3,\n", " 'K_12': 0.2,\n", " 'K_43': 0.1,\n", " 'K_inact1': 0.2,\n", " 'K_inact2': 5.0}" ] }, "execution_count": 6, "metadata": {}, "output_type": "execute_result" } ], "source": [ "k" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "#### Differential Equations" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "$F_{3n} = \\frac{Ni^{2.5}}{Ni^{2.5} + 17^{2.5}}$\n", "\n", "$\\frac{dF_{1}}{dt} = F_4C_iK_{con1} - F_{1}K_{coff1} + F_{2}F_{3n} K_{inact1} - F_1K_{12}$\n", "\n", "$\\frac{dF_{2}}{dt} = F_3 C_i K_{con2} - F_{2} K_{con2} + F_{1} K_{12} - F_{2} F_{3n} K_{inact1}$\n", "\n", "$\\frac{dF_{3}}{dt} = F_2 K_{coff2} + F_{4} K_{43} - F_{3}C_{i}K_{con2} - F_{3}F_{3n}K_{inact2}$\n", "\n", "$F_{4} = 1 - F_{1} -F_{2} - F_{3}$" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "#### Reduced Equation" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "$\\frac{dF_{1}}{dt} = C_iK_{con1} - F_{1}(K_{coff1}+C_iK_{con1} + K_{12}) + F_{2}(F_{3n} K_{inact1}- C_iK_{con1}) - F_{3}C_iK_{con1}$\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "$\\frac{dF_{2}}{dt} = F_{1} K_{12} - F_{2} (K_{con2}+F_{3n} K_{inact1})+ F_3 C_i K_{con2}$" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "$\\frac{dF_{3}}{dt} = K_{43} - F_{1}K_{43} + F_2 (K_{coff2} - K_{43}) - F_{3}(C_{i}K_{con2} + F_{3n}K_{inact2} + K_{43}) $" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Calculation" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "# function that returns dy/dt\n", "def model(y,t,k):\n", " \n", " '''y[0] = F1,\n", " y[1] = F2,\n", " y[2] = F3,\n", " '''\n", " \n", " dy1dt = k[\"Ci\"]*k[\"K_con1\"] \\\n", " - y[0]*(k[\"K_coff1\"]+k[\"Ci\"]*k[\"K_con1\"]+k[\"K_12\"])\\\n", " + y[1]*(k[\"F_3n\"]*k[\"K_inact1\"]-k[\"Ci\"]*k[\"K_con1\"])\\\n", " - y[2]*k[\"Ci\"]*k[\"K_con1\"]\n", " \n", " dy2dt = y[0]*k[\"K_12\"] \\\n", " -y[1]*(k[\"K_con2\"]+k[\"F_3n\"]*k[\"K_inact1\"])\\\n", " +y[2]*k[\"Ci\"]*k[\"K_con2\"]\n", " \n", " dy3dt = k[\"K_43\"]\\\n", " -y[0]*k[\"K_43\"]\\\n", " +y[1]*(k[\"K_coff2\"]-k[\"K_43\"])\\\n", " -y[2]*(k[\"Ci\"]*k[\"K_con2\"]+k[\"F_3n\"]*k[\"K_inact2\"]+k[\"K_43\"])\n", " \n", " dydt = [dy1dt,dy2dt,dy3dt]\n", " return dydt" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [], "source": [ "N = 1\n", "T = 50\n", "y0 = [1,0,0]\n", "t = np.linspace(0,T)\n", "\n", "'''solve ODEs'''\n", "y = odeint(model,y0,t,args=(k,))\n", "\n", "y1 = np.empty_like(t)\n", "y2 = np.empty_like(t)\n", "y3 = np.empty_like(t)\n", "y4 = np.empty_like(t)\n", "I = np.empty_like(t)\n", "for i in range(len(t)):\n", " y1[i] = y[i][0]\n", " y2[i] = y[i][1]\n", " y3[i] = y[i][2]\n", " y4[i] = 1 - y1[i] - y2[i] - y3[i] \n", " I[i] = y2[i]*F_3n" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "(5.014955054667541e-05,\n", " 5.407977907554697e-07,\n", " 0.019835680478156214,\n", " 0.9801136291735063,\n", " 5.344296471883424e-07)" ] }, "execution_count": 9, "metadata": {}, "output_type": "execute_result" } ], "source": [ "y1[T-1],y2[T-1],y3[T-1],y4[T-1],I[T-1]" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAA4AAAAFFCAYAAAC9joN5AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJzs3Xd8HPWB/vHP7K5WvTdLcm9jGeNuY4wNtgHTQocADgGSS0hySS6FS7lL7siRhIT75bgUQi4VCMTUhBoMxtgEg3HvbdybJKtbve3u/P5YSZaNuzUa7ep5v1567U7ZmQcYJ3r8nZ2vYds2IiIiIiIiEv08bgcQERERERGRnqECKCIiIiIi0keoAIqIiIiIiPQRKoAiIiIiIiJ9hAqgiIiIiIhIH6ECKCIiIiIi0keoAIqIiIiIiPQRKoAiIiIiIiJ9hAqgiIiIiIhIH6ECKCIiIiIi0keoAIqIiIiIiPQRPrcDnKNYYApQAgRdziIiIiIiItLTvEAesApoOdMPRWoBnAIsdTuEiIiIiIiIy2YCH5zpzpFaAEsAqqsbCIVst7McIzMzicrKerdjSJTS9SVO0zUmTtL1JU7S9SVO6o3Xl8djkJ6eCO3d6ExFagEMAoRCdq8rgECvzCTRQ9eXOE3XmDhJ15c4SdeXOKkXX19n9ZU4PQRGRERERESkj1ABFBERERER6SNUAEVERERERPoIFUAREREREZE+QgVQRERERESkj1ABFBERERER6SMcnwbCNM0UYBnwCcuy9h23bTzwByAFeB/4omVZAacziYiIiIiI9EWOjgCapnkR4VnpR55kl2eAr1iWNRIwgM87mcdpIdumsbnN7RgiIiIiIiIn5PQtoJ8HvgwUH7/BNM1BQLxlWcvbVz0J3O5wHke9/P4ePvWfb3GgtM7tKCIiIiIiIh/j6C2glmV9DsA0zRNtzgdKuiyXAP3P5viZmUnnnM0JsbExBIIhVu+sYNKYfLfjSJTKzk52O4JEOV1j4iRdX+IkXV/ipGi5vhz/DuApeAC7y7IBhM7mAJWV9YRC9ul37CEXDErjb8D764q4afpgPB7D7UgSZbKzkykv1wizOEfXmDhJ15c4SdeXOKk3Xl8ej3FOA2JuPgX0EJDXZbkfJ7hVNJIM7pdMXmYitQ2tWAeq3Y4jIiIiIiJyDNdGAC3L2m+aZrNpmpdYlvUh8GlggVt5uoNhGMycUMALi3awcnsZhYMz3I4kIiIiEci2bbDtris6l+3O9Xb4Xqrj9+vY1rF4/P6c+DN2x8aP7dN1vX3cpq53YtkneHuCbMfs9vF1JzzmMTd82SfYDxqbE2mpajhxnhMsfnz7me9vH7/vqZxo39N+/DQ7nOL8p492nnfPnc0/e3d+Fs47+vloYSAQ616AbtTjBdA0zTeB/7QsazXwKeD37VNFrAV+2dN5utul48MFcPX2Mj515Uh8Xk21KCIi5862bQgGsUMhCAWxgyEIhbA733d5DbVvC4a3h/cLnfAVu+O9fYJ9bbCPfW+HQuFf3kKhcKbQiddh213ed6wPl5eaWB9NTS3hc9qhzvWd+7Xva9uho4Wn43jnsQz20SJ0/H6dJSrU/rY9b2dZ6jje0X06ttl0HKt9G3QpaR37dSkrHytwx64/71+O+7h9bgeQqHbI62XI//tffCkpbkc5bz1SAC3LGtzl/bVd3m8ApvZEhp4yKC+FgqxEiioa2LqvmrHDMt2OJCIi7exgEDsQwG5rC78G2l/bjr4PtbWF37cFsIMBCASxgwHsztf2n2AwXMy6LB/dHn5PMNi+PthZ4sLvA+GS1rE+GGwvdEeX6dg3ikpBjdsBIpFhHPu+fdk4yfpjP2McXTzZduPYdUbHis7dj1sOH+zYTV2PjXGCtydY1+UzHz/n8Z85xfYu63xeD8HgcX9ejsl2/DFOsOL4/T+2u3Git6c+5pkc9zTnOqftp/zoeT6n4pQfdy63m5IHFOBNTHQ7Rrdw8yEwUWtqYQ4vL93Lym2lKoAiIidh23a4aLW0EmptwW5pIdTaSqilBbvztYVQaxt2Wyt2Wxuh1vBr+Kc1XNa6rOu6PdTWZd/20hexZcrrxfB4wOMJv7YvG15v+zoveNu3eY7b91Svhufo54yj6/EYXdYZ4XUd7zv2M4yjxzKMo+u6vvd4wr9oGuFjJqfEU1/fEl7u2M/T/hmj4/OE19NeajzhwtFxbuDo8eHoOTrKTPv5wqs8nes69/EYncfuLEZG13VdjtNZkNrP32W5M2vn++O2dRyiy3LXfU5a4IzjtskZ640P6ZDoEU3XlwqgA6YW5vLy0r2s3VFOWyBIjM/rdiQRkfNi23a4VDU3t/80EWpuDpe2Lssf+2kJv9qtLYRaWttfjxa8Hi9khoERE4Ph87X/xGDEtL92rDtme5cfrw/D6+1cpuN956sPw+c9us3r7fzM0WVvuKD52pc9Hug4rtdzzOfwGOHXKCoE0fQLlIhIpFIBdEBuRgKD+iWz/3AdG3dXMcnMdjuSiPRxdigULmRNjYQaGwk2dnltaiLUFH5/JNRGY1UNwfb9Ovdpbgp/p6ubGT4fht+PJzYWwx+Lx+/HiD3u1e/HExODEdP1fQyG34/hiwnv077u6Hs/hj8mvG+XcofXGzVlSkRE5FyoADrkosJc9h+uY+W2UhVAEelWdiBAsKGBYH09wYb68Gt9HaH6jvfh5c596usINTae92ib4fNhxMXhiYvDE9v+2vkTf9zy0R8jNg5PbCwef+zHy51Xd0iIiIj0JBVAh0wZlcMLS3axYVcFza0B4vz6Vy0iJxdqaSFQW0OwpoZATQ3B2tqjy7U1BGtrCNbVEayvJ9TUdE7n8MTF4UlIwBOfgDchof19fJf3CaTlZtAQ9OBNSMQTH17vjY/HEx8fHkETERGRiKb/N3dIZmocwwtS2VVUw/pdFUwb3c/tSCLigmBjI4HqKgJVVQSqqwnUdhS89pLXXvjsluYzP6hh4E1MwpuUhCcp/OpNSmpfl4w3KbF9XTKexCS8yUl4ExLPaLRN39ESERGJbiqADppamMOuohpWbi1TARSJQqGWFgLVVbRVdRS8KtqqKjvfB6qqCDWfWbEzfD68qan4UlPxpqTiS0kNL6ek4k1JaX9NxpuYhCchof0piSIiIiJnRwXQQVNG5fDsuzvZtKeSxuY2EuJi3I4kImch1NZKW3kFbWWltJWV0VpeRqC94LVVVxGqrz/tMQy/H19GBjHpmfjS0/GlpbUXvJRjCp8nPl4PJxERERHHqQA6KDUpllED09m2v5q1OyqYMTbP7UgicpxQcxOtZWW0lZXRVl5Ga3vZaysvI1BdfeoHp3i9xKRn4MsI/8RkZOJL73ifgS89A09iooqdiIiI9BoqgA6bWpjDtv3VrNxWqgIo4hI7EKC1rJTWoiJaD5ccLXllZQTrak/+QY+HmKxsYnJyiMnJwZ+dgy8zK1zuMjLwJqfoVkwRERGJKCqADptk5vDMwh1s3VdNbWMrKQl+tyOJRC07FKKtvIyWoiJai4toLTpES3ExrYdLIBg84WcMny9c8LJziMnJxd/lfUxGhp58KSIiIlFFv9k4LCk+hguGZLBxdyVrrHJmTyhwO5JIxLNDIQJVlbQUF9FaVNT52lpSjN3WdsLPxGRn488vwJ+Xjz8nt31ULxdfWppG8URERKTPUAHsAVMLc9i4u5KVW0tVAEXOkh0K0VZ6mOa9e2net4fmfftoKSo66bQJvvQM/AUFxOYXhF8L+uPPy8cTG9vDyUVERER6HxXAHjBhRDY+r8WOg0eormshPVm/iIqciG3bBKqraN67p73w7aVl/74TTnzuTUnBn99e8PILiC0owJ+fjzch0YXkIiIiIpFBBbAHxMf6GDcskzU7ylm1vYy5Uwa4HUmkVwjW14dH9fbuDZe+fXsJ1n78oSy+9AziBg8hbsgQYgcPIW7AQLzJyS4kFhEREYlsKoA9ZEphDmt2lLNyW6kKoPRJtm3TWlxE4/ZtNO/aSfO+vbSVl39sP09CInFDwmUvbvBQ4gYPwZeW5kJiERERkeijAthDxg3LIjbGy57iWsqPNJGdFu92JBFH2bZNW+lhGrdvo3H7Npqs7QTr6o7Zx/D7iR04iLgh4aIXN3gIMTk5mjdPRERExCEqgD0k1u9l/IgsVmwtZeW2Uq67eLDbkUS6lW3btFWU09Re+Bq3bydYc+SYfbxpaSSYhcSbJvFDhuLPL8Dwel1KLCIiItL3qAD2oKmFOazYWsqqbWUqgBIV2ior20f3woUvUFV5zHZvcjLxZiEJo0aRMGo0Mbm5Gt0TERERcZEKYA8aMyST+FgfB8rqKalsIC9TTyuUyBJqbaVx6xYaNq6ncds22srLjtnuSUgkwRxFfGEhCWYh/vx8FT4RERGRXkQFsAfF+DxMGpnNB5tKWLmtjBtnDHE7kshpBevradi0gfp1a2nYvAm7tbVzmyc+nvgRI0kYVUj8qEJi+w/QpOoiIiIivZgKYA+bOjqnvQCWcsMlgzU6Ir1SW2Ul9evXUr9uLU07LAiFOrfFDhpM0vgJJI65kNiBg/QdPhEREZEIogLYwwoHpZMUH0NJZSMHy+oZmKu5zMR9tm3TeuhQZ+lrObD/6EaPh4TC0SROmEjSuAnEZGa6F1REREREzosKYA/zejxMGZXDknVFrNxWpgIorrFDIZp27Qzf2rluLW0VR+fkM2JjSRxzIUnjJ5I4dhzeRH1fVURERCQaqAC6YGphRwEs5dbLhuo2UOkxtm3TvGc3tR9+QP3aNQTrj87L501OIXH8eJImTCShcDSeGL+LSUVERETECSqALhjRP420JD8VNc3sKallWH6q25EkygVqaqhdvozaD5bSWlLcuT4mJ5ekCRNIGj+JuGHD9AAXERERkSinAugCj8dgyqhc3ll9kJVby1QAxRF2IEDDpo3UfLiUho0bOh/k4k1OIWX6dFIuvgR/QX+NQIuIiIj0ISqALpk6OidcALeXcsec4Xg8+iVcukdLcRG1Hy6ldtkygnW14ZUeD4njJ5B6yUwSLxyL4dMffREREZG+SL8FumRoXgpZqXFU1DSz89ARzIHpbkeSCBZsaqJu1QpqP1hK857dnev9efmkzJhJyrSL8aWmuZhQRERERHoDFUCXGIbB1MJc3ly+n5XbylQA5azZoRBNOyxqPlxK/ZrVnRO0e+LiSJ46jZQZM4kboocMiYiIiMhRKoAumlqYw5vL97NqexnzrhyBVw/gkDMQammhZun7HHl3IW3lR6duiB9VSOolM0maOAlPbKyLCUVERESkt1IBdNGAnCTyMhMoqWxk2/5qxgzRBNtycsH6eg68u4Ci1/9OqL4eAF9GBinTZ5ByyQz82TkuJxQRERGR3k4F0EUdt4G++sFeVm4tUwGUE2qrqqR64dvUvP9e522ecUOHknHNdSSOm6CpG0RERETkjKkAumxqYQ6vfrCXNTvK+fRVJjE+/TIvYS3FRVS/9Sa1K5ZDMAhA2sQJJF1+FfEjTX23T0RERETOmgqgy/IyExmQk8TBsno2761kwohstyOJy5p276Jqwd9pWL8uvMIwSJ46jfSrr2HApDGUl9e5G1BEREREIpYKYC8wtTCHg2X1rNxWpgLYR9m2TePmTVQt+DtNOywAjJgYUmbMJH3u1fp+n4iIiIh0CxXAXmBqYS5//cce1u0sp6U1SKzf63Yk6SF2MEjd6pVULXiT1kMHAfDEx5M2+3LSLr8SX2qqywlFREREJJqoAPYC2WnxDM1PYU9xLRt2VzC1MNftSOIwOxik5sOlVL/5d9oqwlM5eFPTSJ97FamXzsIbH+9yQhERERGJRiqAvcTUwlz2FNeyYmupCmCUa9y2lbLn5tNadAiAmNx+ZFx9DcnTpuOJiXE5nYiIiIhEMxXAXmLKqBxeWLyLjbsrqa5rIT1ZE3lHm9bSw5S/+Hznw118WVlk33I7SZOnaCoHEREREekRjhZA0zTnAd8HYoCfW5b16+O2TwR+C/iBg8DdlmUdcTJTb5WeHMvEkVmstsr5x/oibpo51O1I0k2CjY1U/f01qhe9A8EgRmwcmdd9grQr5+KJ8bsdT0RERET6EMeGHUzTLAB+DMwAxgP3m6Y5+rjdfgH8p2VZ4wAL+Fen8kSCORP7A/CP9cUEgiGX08j5skMhjvxjCfu+9x2q334LQiFSZsxkyI9/Ssa1n1D5ExEREZEe5+QI4BXAYsuyqgBM03wJuA14qMs+XiCl/X0CUOVgnl7PHJhGQVYiRRUNrLHKuWi0vgsYqRq3baXs+Wc7n+wZP2Ik2XfMI27wYHeDiYiIiEif5mQBzAdKuiyXAFOP2+ebwELTNH8ONAAXnc0JMjOTziugU7Kzk8/5szdeNozH/7qRpZtK+MRlw7sxlfSEppIS9j3xZ6pWrAQgNiebwffdQ+b0izEMo1vOcT7Xl8iZ0DUmTtL1JU7S9SVOipbry8kC6AHsLssG0Hlfo2ma8cAfgSssy1ppmuY3gT8D153pCSor6wmF7NPv2IOys5MpL68758+PGZRGfKyXrXurWLO5mIG50XGhRbvw9/xep3rRwvbv+cWSce0nSJ97FXaMn4qK+m45z/leXyKno2tMnKTrS5yk60uc1BuvL4/HOKcBMScL4CFgZpflfkBxl+UxQJNlWSvbl38L/NDBPBEhzu/jkjF5LFpziMVrD3HfNYVuR5JTsEMhapa+T+UrfyVYF/4fhZTpM8i65TZ8aWkupxMREREROZaTBXAR8APTNLMJ3955K3B/l+27gAGmaZqWZVnAjcAqB/NEjNkTC1i05hDLt5Ry++zhJMZpbrjeqHnvHkr//AQtB8Pf84sbPoKcO+cRN3iIy8lERERERE7MsQJoWVaRaZrfA5YQnubhD+23er5J+Mmfq03TvA94wTRNAygDPuNUnkiSl5nIBYPT2bKvmg82lnDV1IFuR5Iu7GCQqjffoPL1VyEUwpeRSfbtd4Tn8+um7/mJiIiIiDjB0XkALcuaD8w/bt21Xd4vABY4mSFSzZnUny37qlmytogrpwzAo2LRK7SWlXH4D7+lec9uANKuvIqsm2/F49eUDiIiIiLS+zlaAOXcjRuWRWZKHGVHmti8p4qxwzLdjtSn2bZN7dL3KXt+PnZLC770dPp99vMkFB4/taWIiIiISO+lAthLeTwGsycW8NJ7u1m89pAKoIsCdbWUPvUEDevXAZA8ZSo5d9+LNzHR5WQiIiIiImdHBbAXmzk2j1eW7mXT7krKqhvJSU9wO1KfU79xA6VP/pFgbS2e+HhyPvVpki/qvjn9RERERER6ksftAHJyyQl+LirMwQaWrCtyO06fEmppofSZP1P8y/8lWFtLvDmKQT/4ESnTpqv8iYiIiEjEUgHs5eZM6g/A0g0ltLQFXU7TNzTv3cP+hx6k5r3F4PWSdfsd9H/g28Rk6jZcEREREYlsugW0lxuSl8KQvBT2ltSyYmspl47LdztS1LKDQaoW/D08vUMwiD+/gLzPf4HYAZqGQ0RERESig0YAI8DlkwoAWLzmELZtu5wmOrWWlXHwv39C5St/g2CQtCuvYuB/PKjyJyIiIiJRRSOAEWDKqByeX7yLA2X17CqqYUT/NLcjRQ3btqn9cCllz87HbmnGl55O7mc+R+LoC9yOJiIiIiLS7VQAI0CMz8ul4/L5+0f7Wby2SAWwm4Samzn8xB+oX7MaaJ/e4VP34E1KcjmZiIiIiIgzdAtohJg1vgDDgNXbyzhS3+J2nIgXOFLNwf/+CfVrVuOJj6ff5+6n3/1fUvkTERERkaimAhghMlPjGD88i2DI5v31xW7HiWgtBw9w4Mc/pOXAfmJychn4vQc1vYOIiIiI9AkqgBHk8vYpId5bX0QgGHI5TWSq37iBAz99mEB1FfEjRjLw376Pv18/t2OJiIiIiPQIFcAIUjgonbzMBI7Ut7JuZ4XbcSJO9eJFFP/q59gtzSRfNI2Cb34Lb3Ky27FERERERHqMCmAEMQyDORPDo4DvrjnkcprIYYdClD33F8rnPwO2Tcb1N9Lvc1/AExPjdjQRERERkR6lAhhhpo/pR6zfy46DRzhUVu92nF4v1NxM8a9/yZFF74DXS79/+jxZN96s7/uJiIiISJ+kAhhh4mN9XDIm/J21xWs1CngqbdXhJ302bFiPJyGR/t/8FikXX+J2LBERERER16gARqDZ7beBLttymMbmNpfT9E7NB/Zz8OGHjj7p89//gwRzlNuxRERERERcpQIYgQqyEikclE5rW4gPNx12O06vU79xPQcfeZhAdbWe9CkiIiIi0oUKYISaM7EACN8GGrJtl9P0HuEnff4Cu6VFT/oUERERETmOz+0Acm7Gj8giPTmW0uomtu6tYszQTLcjucoOhSh//lmOvPsOABnX30jmDTfpYS8iIiIiIl1oBDBCeT0eZk/oGAUscjmNu0LNzRQ/9otw+dOTPkVERERETkoFMIJdOi4fn9dgw64Kyo80uR3HFYEjR8JP+ty4QU/6FBERERE5DRXACJaS6GfKqBxs4L11fW8UMFBXy6H/+W896VNERERE5AypAEa4OZPCU0K8v6GY1ragy2l6TrCxgaJHf0ZrSTH+gv560qeIiIiIyBlQAYxwQ/NSGNQvmYbmACu2lbodp0eEmpso+vmjtBw8QExuLv2/+a960qeIiIiIyBlQAYxwhmFwefvE8IvXFGFH+ZQQodZWin71C5r37MaXmUn/B76NLzXN7VgiIiIiIhFBBTAKTC3MITHOx/7SOnYV1bgdxzF2IEDx44/RZG3Hm5pG/29+m5iMvj39hYiIiIjI2VABjAL+GC+z2qeEeO2DvS6ncYYdDFLy+/+jcfNGPElJ9H/gW/hzc92OJSIiIiISUVQAo8RVUwcS5/eyZV811oFqt+N0KzsU4vCTf6R+zWo88fH0/8a/Eptf4HYsEREREZGIowIYJZLiY5g7ZQAAL7+/J2q+C2jbNmV/eZq6j5ZhxMZS8PUHiBs02O1YIiIiIiIRSQUwisydMpDEOB87DtWwdV/kjwLatk3Fi89T848lGDExFHz168QPG+52LBERERGRiKUCGEUS4nxcfdFAAF5eGvmjgFWvv0r1wrfA6yXvS18hYVSh25FERERERCKaCmCUuXxSf5ITYthTXMuG3ZVuxzlnVW8voPK1V8AwyPv8F0kaO87tSCIiIiIiEU8FMMrE+X1cN20QAK+8v4dQBI4CHlmymIoXnweg32c+R/LkKS4nEhERERGJDiqAUWjWhALSkvwcKKtnrVXudpyzUrvsQ8r+8mcAcj51DynTL3E5kYiIiIhI9FABjEL+GC+fmD4YgFc+2EsoFBmjgHWrV3H4iT8AkHX7HaTNnuNyIhERERGR6OJzO4A4Y+bYfBYsP0BxRQMrtpVy8QX93I50SvUb11Py+/8D2ybj+hvJuOoatyOJiIiISIQKBgNUV5cTCLR2y/HKyjyEQqFuOda58Pn8pKdn4/Wef31TAYxSMT4PN1wymCcWbOfVD/YytTAHr6d3Dvg2bttKyeOPQTBI+tyrybzhJrcjiYiIiEgEq64uJy4ugcTEfhiGcd7H8/k8BALuFEDbtmloqKW6upysrLzzPl7vbATSLaZf2I+c9HjKqptYtumw23FOqLWkmKLHfokdCJB62Wyybr+jW/6QioiIiEjfFQi0kpiYEhW/VxqGQWJiSreNZqoARjGvx8ONM4YA8NqHe2lz6W8tTibU3ETx449htzSTNHkKOZ/6dFT8IRURERER90XT75Xd+c/i6C2gpmnOA74PxAA/tyzr18dtN4HfAunAYeBOy7KqnczU11xUmMvfP9pPcUUDSzcWM2dif7cjAeGh7MNP/onWkmL8+fn0u++fMHrpLaoiIiIiItHCsd+4TdMsAH4MzADGA/ebpjm6y3YDeA34qWVZ44B1wHedytNXeTwGN7WPAr6xbB+tbUGXE4UdeWch9atX4YmLI/+fv4onLs7tSCIiIiIijigpKWbWrGncd9+8Y35KS8Nf01q1ajlf+9qXeiSLkyOAVwCLLcuqAjBN8yXgNuCh9u0TgQbLst5qX34YSHMwT5810cxmYG4SB0rreW9dEXOnDnQ1T+MOi/KXwhO9537mc/j7nf+XWUVEREREerOsrGyefHL+MetCoRDPPvsMTz/9BEOHDuuRHE4WwHygpMtyCTC1y/Jw4LBpmn8EJgDbgK86mKfP8hgGN80cyi9f2sjfl+/n0vH5xPndeQBs4MgRSn77OIRCpF99LcmTJruSQ0RERET6jp+/uIGNuysdOfbYYZl8/fZx5/TZffv2sn//Xr7zne/x4ovPdXOyE3OyBXiArjOQG0DXp5D4gFnApZZlrTZN84fAo8B9Z3qCzMyk80/pgOzsZLcjfMwVWUm8vfIg1oFqlm8v5/bLR/Z4hlAgwJZHHyFYU0PqhWMovP8+DK+3x3NEut54fUl00TUmTtL1JU7S9SUdyso8+HxHv+3m5PNgDINjznUiXq+HiopyPvOZeZ3rrrrqGu6++16+//0HWbNmNYZhnPI4Ho+nW65xJwvgIWBml+V+QHGX5cPATsuyVrcvPwu8dDYnqKysJxSyT79jD8rOTqa8vM7tGCf0iemDsA5U89fFO5k6MpuEuJ4dBSx7bj61W7fhTUsj877PU1HV2KPnjwa9+fqS6KBrTJyk60ucpOtLugqFQsfM2/e1285thK7D6eYBPN0cgcFgiKysbJ544thbQDs+FwyGsG37lMcJhULHXOMej3FOA2JOPnZxEXC5aZrZpmkmALcCb3XZvgzINk2z47/G9cAaB/P0eaMHpWMOSKOhOcDCVQd69Nx1q1ZyZNFC8HrJ/9JX8KWm9uj5RURERETEwQJoWVYR8D1gCbAemG9Z1krTNN80TXOyZVlNwM3A703T3ALMAR5wKo+E5w+5+dKhACxcdZD6prYeOW9LcRGHn/wjANl33EX8sOE9cl4RERERETmWo/cAWpY1H5h/3Lpru7xfwbEPhhGHjRyQxgVDMtiyt4q3VhzgtlnOPm0o2NRE8eO/wm5pIfmiaaTNvtzR84mIiIiIyMlp5u0+6OaZ4VHARWsOUtPQ6th5bNum9Mk/0nb4MP6C/uTQjrVjAAAgAElEQVTe8xkMJ7+BKyIiIiLSC+Xl5fPSS6+fdPvEiZN57LHf9UgWFcA+aGh+CuOHZ9HaFuLNj/Y7dp7qhW9Rv2Y1nvh48v/5K3hiYx07l4iIiIiInJ4KYB9108whACxZV0R1XUu3H79x+zYq/voiAP0++zn8uf26/RwiIiIiInJ2VAD7qIG5yUwelUMgGOKNZfu69dht1dWU/PY34cner7mOpAmTuvX4IiIiIiJyblQA+7CbZgzBMOD9DcWUH2nqlmPagQAl//drgnW1xI8qJOumW7rluCIiIiIicv5UAPuw/KxEpo3OJRiyef3Dfd1yzPIXn6d59y586Rnk3f8lDK+3W44rIiIiIiLnTwWwj7thxhA8hsGyzYc5XNV4XseqXbGcI+++A14veV/6Mr6UlG5KKSIiIiIi3cHReQCl98tNT2DG2H68v6GEV5bu4Ys3jjmn47QUHaL0qT8BkHPnPOKHOju/oIiIiIhIpCgpKeauu25h8OChx6x/5JFHeeSRH1NVVYnHY/DlL3+dSZOmOJpFBVC4fvoQlm0uZeW2Mi4bV0Xh4Iyz+nywsTE82XtrK8kXTyd11hyHkoqIiIiIRKasrGyefHL+Meseeug/uOSSmdx66yc5cGAfX/nKF3j55TfxOvg1Kt0CKmSmxnH99EEA/Plti7ZA8Kw+X/bsM7SVluLvP4Dcu+/VZO8iIiIiImfg0ktnceWVVwNQUDCA1tYWmpq65+GMJ6MRQAHgmmmDWLGtjOKKBl5ftp9bLh16+g8BDZs3UffRMoyYGPK/+GVN9i4iIiIivU7jgkcJHtzoyLG9A8aScM03T7tfRUU59903r3N57tyrmTfvns7lZ599mhEjTJKSkhzJ2UEFUADweT3ce7XJT55Zy4Ll+7moMIeC7FNffKGWFkqfeQqAzOtvxN9Pk72LiIiIiJzIiW4B7fDCC/N59dWXeeyx3zqeQwVQOo3on8as8fm8t76Yp96y+O7dE/Gc4nbOyldfJlBRQeyAAaTPvboHk4qIiIiInLkzGaE7FZ/PQyAQ6qY0x3r88V+wbNmH/PrXvyMnJ9eRc3Sl7wDKMW6bNYzURD+7imp4f33xSfdr3reP6nfeBsMg957PYPj0dwkiIiIiImfjhRfms3btGn7zmz/2SPmD04wAmqYZC9wP3AKYQBDYDrwIPGlZVqvjCaVHJcTFMO/Kkfzmlc28+N5uxo/IIi3p2O/12cEgpX9+AmybtCvmEjfkzL4vKCIiIiIiYbZt86c//Z7ExES++tUvdK7/2c9+QVZWtmPnPWkBNE3zMuBXwAfAL4G9QBswBLgaWG2a5tcty1rsWDpxxWQzm7HDMtm4u5L5i3byzzcdOzdg9aKFtBzYjy8zk6ybbnEppYiIiIhIZMjLy+ell14/Zp1hGLz11pIez3KqEcDrgRmWZdUet34L8IZpminAg4AKYJQxDIO7547k+39YwertZWzYVcG44VkAtJaXUfnqywDk3n0vnrg4N6OKiIiIiMhZOOl3AC3L+lfLsmpN0yw4fptpmqMty6q1LOsBZ+OJW7JS47l5ZvjWzmcWWjS3BrBtm7KnnwpP+D51GokXjnU5pYiIiIiInI1T3QKa0f72TdM0ZwEdj4OMAf4GjHI2mrjtisn9Wb6llP2ldbyydC/XJpTRuHULnsREsu+cd/oDiIiIiIhIr3Kqp4A+C1QAFwKV7e8rgIPAWuejidu8Hg/3XTMKw4APlu/k8LPheUuyP3knvpQUl9OJiIiIiMjZOukIoGVZVwGYpvkny7I+23ORpDcZ1C+ZKycPwP/Gs9DYQPyoQlKmz3A7loiIiIiInIOTjgCapnkFwKnKn2mac50IJb3LVZkNjKnbQ5vhZee4uRinmBxeRERERER6r1M9BfQK0zS/AvwOWGhZVgDANE0/MBv4KuE5ARc6nlJcE2ppofq5ZwD4MGMs6zbWMnZqE1mp8S4nExERERGRs3WqW0C/a5rmRcAPgOdN0ywhPGKYCywFfmhZ1kc9klJcU/nqywQqKogdMIDghbNo2VHJMwt38LXbxmokUERERETkDJSUFHPXXbcwePDQY9Y/8sijPPzwQxw5Uo3f7+fb3/53RowwHc1yqhFALMtaYZpmE/BZwhPBh4DdlmXVOJpKeoXm/fuofudtMAxy7/kMd2UXsHn/CjburmS1Vc6UUTluRxQRERERiQhZWdk8+eT8Y9Y9/PB/MXv2HG666TaWL1/G//zPI/zf//3J0RynLIDt/gp8BcgjfDvoPicDSe9gB4OUPvUE2DZpV8wlbshQ4oDbZg3j6bct5r+zgwsGp5MQF+N2VBERERGRU3p8w5/YUrndkWNfkDmKfx53bs/M/O53/6PzfUlJMcnJyd0V66ROWwAty/oL8BfTNEcRHglcZZrmMuBXlmWtdDqguKN60UJaDuzHl5lJ1k23dK6/bHw+H20+zK6iGl76xx7uucrZIWoRERERkWhQUVHOffcdnUt77tyrmTfvHgDmzbuVw4dL+MlP/sfxHGcyAohpmh5gBDCy/TNlwOOmaf7dsqwHHcwnLmgtL6Py1ZcByL37XjxxcZ3bPIbBvVeb/OCJVby3roiLL8hlRP80t6KKiIiIiJzWuY7QdfD5PAQCofM6xoluAe0wf/5f2bnT4hvf+Arz579ESkrqeZ3rVE41ETwApmn+iPDk798GngeGW5b1AHAZ4VtDJYrYtk3Z009ht7aSPHUaiReO/dg+BdlJXDNtIABPvWURCJ7fHwYRERERkb5o2bIPaGxsBGDECJN+/fIoLi5y9JynLYBADnCtZVkzLct61rKsNgDLshqAuxxNJz2ubvkyGrduwZOYSPad80663/XTB5ObHk9xRQMLVhzowYQiIiIiItFhwYI3eO21vwGwd+8eqqoqGThwsKPnPG0BtCzrfsuyNpxkm+YAjCKBulrKnn8WgOxP3okvJeWk+8b4vJ3f/3v9w30crmrskYwiIiIiItHia197gJUrl3PvvXfx8MP/xYMP/piEhARHz3lG3wGUvqH8+WcJ1deTUDialOkzTrt/4eAMLhnTjw83H+bPb23nW3dN0NyAIiIiIiLHycvL56WXXv/Y+qysbB599LEezXImt4BKH9CwZTN1yz/CiIkh5+57z7jIfXLOcJLiY9h+4AjvrXP2fmURERERETk/KoBCqK2VsqefAiDz+hvx5+ae8WeTE/zcPXckAM++u5O9JbWOZBQRERERkfOnAijULFlCW0U5/oL+pM+9+qw/P7UwlzkTCwgEbR5/eTP1TW0OpBQRERERkfOlAtjHhZqbqFrwBgBZt9yG4Tu3r4XeMWcEQ/JSqKxt5vevbyVk290ZU0REREREuoEKYB9XvegdgnV1xA0bTuLYced8nBifh3++aQxJ8TFs2lPJG8v2dV9IERERERHpFiqAfViwoYHqtxcAkHXzref9BM/M1Djuv340BvDq0r1s2VvVDSlFRERERKS7qAD2YdVvLyDU1ERC4WgSRhV2yzHHDM3khhlDsIHfvraFqtrmbjmuiIiIiEikKikpZtasadx337xjfkpLDwPQ2NjAJz95I2vXrnY8i+YB7KMCNTVUL1oIQObNt3brsa+/ZDC7i2rYvLeKx1/ZzHc/NRGfV3/XICIiIiJ9V1ZWNk8+Of+E2x599L+pq6vrkRz6rbyPqlrwBnZrK4njxhM/dFi3HttjGHz++tFkpsSyp7iW5xfv6tbji4iIiIhEi3ffXUhCQgLDhg3vkfM5OgJomuY84PtADPBzy7J+fZL9rgMesyxriJN5JKytqpKa95YAkHXTLY6cIznBz5duupCfPLOGd9ccYnhBKheNPvP5BUVEREREukvRLx6lYdNGR46deOFYCr72zdPuV1FRzn33zetcnjv3aubMmcsLLzzLL3/5Gx544F8cyXc8xwqgaZoFwI+BSUALsMw0zSWWZW09br9c4GfA+T2BRM5Y1RuvYQcCJE+ZSuyAgY6dZ2h+CnddMYJnFu7gyQXb6Z+TREFWomPnExERERHprY6/BTQUCvHNb36Fb3zj28TGxvVYDidHAK8AFluWVQVgmuZLwG3AQ8ft9wfgv4CfOphF2rWWllLzwVIwDDJvvNnx882eUMCuohqWbynl8Zc38f17JhMfq6+eioiIiEjPOZMRulPx+TwEAqFuShO2f/8+DhzYz09/+kMAiooO8sgjP+I73/k+EydO7tZzdeXkb+L5QEmX5RJgatcdTNP8F2AtsPxcTpCZmXTO4ZyUnZ3sdoST2vH0nyAUIueKORRcOLJHzvnApybzwC/f58DhOp5bsptv3T3pvKec6Mt68/Ul0UHXmDhJ15c4SdeXdCgr8+Dzde/jTs7neN72ByJ2PcaIEcN57bUFnctf+tLn+dznvsCkSScufx6Pp1uucScLoAewuywbQGdtNk1zDHArcDnQ/1xOUFlZTyhkn37HHpSdnUx5ec88wedstRQdovz9peD1knjltT2a8wvXj+ahp1azdH0RA7ISuGLygB47dzTpzdeXRAddY+IkXV/iJF1f0lUoFOrWEbvzHQEMBsOfPdUxbNsmGDx57lAodMw17vEY5zQg5uRTQA8BeV2W+wHFXZZvb9++GngTyDdNc6mDefq8yldeBtsm7bJZxGRm9ei58zIT+ey14bkGn1+8i11FNT16fhERERERt+Tl5fPSS6+fcp/HHvudo7d+dnCyAC4CLjdNM9s0zQTCo31vdWy0LOtBy7JGWpY1HrgWKLYsa6aDefq05r17qF+3BsPvJ+O6613JMGVUDldOHkAwZPObVzZT29jqSg4RERERkb7KsQJoWVYR8D1gCbAemG9Z1krTNN80TdP5aivHqHjlbwCkzbkCX2qaazlunz2M4QWpVNe18LvXtvS6W3hFRERERKKZo49jtCxrPjD/uHXXnmC/fcBgJ7P0ZY07LBq3bMYTH0/G1R/719+jfF4PX7ppDD94YiVb91Xz6gd7ufnSoa5mEhERERHpK5y8BVR6Adu2qXz5rwCkX3kV3iT3n5yanhzLF2+4AMOA15ftY+PuCrcjiYiIiEiUse3oudOsO/9ZVACjXOOWTTTt3IEnKYm0K69yO06nwsEZ3DwzPPL3+9e3UnGkyeVEIiIiIhItfD4/DQ21UVECbdumoaEWn8/fLcfTjNxRzLZtKl4Of/cv45rr8MbHu5zoWNdePIjdRTVs2F3JY3/bxLfnTSQhTpekiIiIiJyf9PRsqqvLqa8/0i3H83g8hELdOxH82fD5/KSnZ3fPsbrlKNIr1a9dQ8v+fXhT00ibfbnbcT7GYxh87vrR/PDJ1Rwoq+eXL23gG3eMJzbG63Y0EREREYlgXq+PrKy80+94hqJpnkndAhql7FCIylfDo3+Zn7gej797hoy7W2JcDA/cOZ705Fh2HKrh1y9vIhB0729XRERERESimQpglKpbsZzW4mJ8WVmkzrzM7TinlJ0WzwN3jCcpPobNe6r47WtbCLo4xC4iIiIiEq1UAKOQHQhQ+drLAGRefyOGr/ff6ZuflcgDd4wnPtbHGqucJxdsJxQFX9oVEREREelNVACjUM2HS2krLyemXz9Spk13O84ZG9Qvma/fPhZ/jIcPNx3muUU7o+LJTSIiIiIivYUKYJQJtbVS9cZrAGTdeAuGN7IeqDKifxpfvWUsPq/BojWHeHnpXrcjiYiIiIhEDRXAKFOzZAmB6mpiBwwgadJkt+OckwuGZPCFG8bgMQzeWLaPBSv2ux1JRERERCQqqABGkVBzM1UL3gAg86ZbMTyR+593kpnNZ68bBcCLS3bz3roilxOJiIiIiES+yG0I8jHVixYSrKsjbthwEseOczvOeZs+Jo+7544E4Om3LZZvOexyIhERERGRyKYCGCWCDQ1Uv70AgKybb8UwDJcTdY85E/tz62VDsYE/vLGNdTvL3Y4kIiIiIhKxVACjRPU7bxFqaiKhcDQJowrdjtOtrrt4MNddPIiQbfObV7awdV+V25FERERERCKSCmAUCLW2cuS9JQBk3nCzy2mccculQ5kzsYBAMMSv/rqJXUU1bkcSEREREYk4KoBRoG7VCkL19cQOGkzc8OFux3GEYRjMu3IkF1/Qj5a2ID9/YQMHSuvcjiUiIiIiElFUACOcbdsceXcRAGlzroia7/6diMcw+Ox1o5g4MpvGlgCPPr+ew1WNbscSEREREYkYKoARrnn3LloO7MeblEzy1Klux3Gc1+PhCzdcwAVDMqhtbONnz62joqbJ7VgiIiIiIhFBBTDCHVkcHv1LvfQyPDF+l9P0jBifh6/cfCHD+6dSVdvCz55bT3Vdi9uxRERERER6PRXACBY4coS6NavB4yF11my34/SoWL+Xr982joG5SZRVN/Hjp1dzqKze7VgiIiIiIr2aCmAEO/KPJRAMkjRhIjEZmW7H6XEJcT4euGM8wwpSqKpt4Sd/WcOWvZoiQkRERETkZFQAI5QdCFDz/nsApM2+3N0wLkpO8POtOycweVQOTS1Bfv7iBpZuKHY7loiIiIhIr6QCGKHq1qwiWFODv6A/8eYot+O4yh/j5Ys3XsA10wYSDNk8sWA7f3t/N7Ztux1NRERERKRXUQGMUEcWvwtE/9QPZ8pjGNw+azj3XGXiMQzeWLaf37++lbZAyO1oIiIiIiK9hgpgBGret4/m3bvwJCSQMu1it+P0KrMmFPAvt40l1u9l+dZS/ue5ddQ3tbkdS0RERESkV1ABjEBHFr8DQOolM/HExrqcpvcZOyyTf/vURNKTY9lxqIYfP72GsmpNGC8iIiIiogIYYQJ1tdStXAGGQeqcvvvwl9MZmJvM9z49if7ZSZRWNfKjP69hV1GN27FERERERFylAhhhape+jx0IkHjhWPzZOW7H6dUyUuL4t7snMmZoBvVNbfy/Z9exenuZ27FERERERFyjAhhB7GCQI+8tBiDt8itdThMZ4mN9fO22scwan09bIMTjr2xmwYr9ekKoiIiIiPRJKoARpH79WgJVVcTk9iOhcLTbcSKG1+Ph01eZ3D57GAAvLtnN0wt3EAzpCaEiIiIi0reoAEaQI+8uAiBtzuUYHv2nOxuGYXDNRYP40k1j8Hk9vLeuiF/9dRNNLQG3o4mIiIiI9Bi1iAjRcuggTTssjNg4UqbPcDtOxJoyKodv3zWBpPgYNu6u5JG/rKW6rsXtWCIiIiIiPUIFMEJ0TPyeeskleOPjXU4T2Yb3T+V790wiNz2eA2X1PPTkKjbvrXQ7loiIiIiI41QAI0CwoYHa5csASJutqR+6Q256At+7ZzLmgDRqGlp59PkNPLtoJ22BoNvRREREREQcowIYAWo+eB+7tZWE0Rfgz8t3O07USIqP4Vt3TeCWS4fi9Ri8s/ogP3xqNYfK6t2OJiIiIiLiCBXAXs4OhajpmPphzhUup4k+Ho/BJ6YP5t8/Hb4l9FB5Aw89tZqFqw4S0lQRIiIiIhJlVAB7uYZNG2krLycmK5vEsePcjhO1huSl8IPPTOWy8fkEgiGee3cn//v8ej0gRkRERESiigpgL3dkcXjqh9TZczT1g8Ni/V7uvXoUX73lQpLiY9iyr5oH/7SStTvK3Y4mIiIiItIt1Ch6sdbDJTRu2Yzh95N6yUy34/QZE0Zm89A/TWXMkAzqm9p47G+beHLBNppbNWegiIiIiEQ2FcBerGP0L2XaxXiTklxO07ekJcXy9U+O464rRuDzenh/Qwk/eGIVe4pr3Y4mIiIiInLOfE4e3DTNecD3gRjg55Zl/fq47TcC/wUYwF7gM5ZlVTuZKVIEm5qo+fBDANJm6+EvbvAYBldOHkDhoHR+99pWDpXX8/DTa7hxxmCuu3gwHo/hdkQRERERkbPi2AigaZoFwI+BGcB44H7TNEd32Z4C/Aa4zrKsccBG4AdO5Yk0tR99iN3STPxIk9gBA9yO06f1z07iP+6dxNwpAwjZNi8v3ctP56+l/EiT29FERERERM6Kk7eAXgEstiyryrKsBuAl4LYu22OAL1uWVdS+vBEY6GCeiGGHQp23f2rqh94hxuflzstH8MCd40lL8rPrUA0P/mklH20+jK3pIkREREQkQjhZAPOBki7LJUD/jgXLsioty3oZwDTNeOC7wCsO5okYjdu20nb4ML70DJImTHQ7jnRxweAMHvqni5g0Mpvm1iC/f2MrP39xIyWVDW5HExERERE5LSe/A+gBug6NGEDo+J1M00wFXgY2WJb11NmcIDOzdz4YJTs7+bw+X/7BewDkX3c1Of3SuiGRdKds4MH7L2bRygP84bXNbNpTydZ9VXxixlDunGuSFB/j7PnP8/oSOR1dY+IkXV/iJF1f4qRoub6cLICHgK5zF/QDirvuYJpmHvA2sBj4xtmeoLKynlCod91+l52dTHl53Tl/vrW8jOrVazB8PnwTp53XscRZ44dm8PDnp/Hy0j28v76YV9/fzburDnDLZUO5dGy+Iw+JOd/rS+R0dI2Jk3R9iZN0fYmTeuP15fEY5zQg5mQBXAT8wDTNbKABuBW4v2OjaZpe4HXgBcuyfuRgjohSs2Qx2DbJUy7Cl5Lidhw5jZREP/dePYpZ4wt49t2d7Dh4hD+/ZfHe2iLuumIE5sB0tyOKiIiIiHRyrABallVkmub3gCWAH/iDZVkrTdN8E/hPYAAwEfCZptnxcJjVlmV9zqlMvV2opYWaD94HIO1yPfwlkgzql8x35k1gtVXOC4t3cqCsnkfmr2PKqBxunz2MrNR4tyOKiIiIiDg7D6BlWfOB+cetu7b97Wo0Ef0xald8RKixkbihw4gbPMTtOHKWDMNgyqgcxg3L5K2VB3jzo/2s2l7G+l0VXHPRQK65aBCxfq/bMUVERESkD1MB6yVs2+bIux1TP1zucho5H/4YLzdcMoSH75/GtNG5tAVCvPbhPv7998tZvkXTRoiIiIiIe1QAe4nm3btoLTqENyWF5MlT3Y4j3SAjJY77b7iAf7t7IoP6JVNd18LvXt/KT55Zy96SWrfjiYiIiEgfpALYS9R+9CEAKdNnYPgcvTNXetiI/mn8x72T+cy1o0hJ9LOrqIYfPbWaP/19GzX1LW7HExEREZE+RE2jFwi1tVG3aiUAKRdPdzmNOMFjGMwcm89kM4c3lu1j4aqDfLCphJXbSrl0XD5XXzSQjJQ4t2OKiIiISJRTAewFGjauJ9TYSOzAQcQW9Hc7jjgoPtbH7bOHc+m4fF5Ysot1OytYtOYQS9YVMe2CXK6dNoi8zES3Y4qIiIhIlFIB7AVqP1oGaPSvL8nNSOCrt47lYFk9by7fz8ptpXy46TDLNh1m4shsrr14EEPyNA+kiIiIiHQvFUCXBepqadi0ETwekqdOczuO9LABOUl84YYLuHnmEN5acYAPNpWwZkc5a3aUM3pwOtdNG8SoQekYhuF2VBERERGJAiqALqtbuQKCQRIvHIsvNdXtOOKSnPQE7rl6FDfMGMI7qw6yZF0RW/dVs3VfNUPyUrju4kGMH5HldkwRERERiXAqgC47evvnJS4nkd4gLSmW22cP59qLB7F4bRHvrDrI3pJaHvvbJvKzErnjypEU9k/F59UDfEVERETk7KkAuqiluJiWfXvxxMeTOH6C23GkF0mMi+H66YOZO2UASzcU89bKAxRXNPC/z64jMyWWq6YOZOa4fGJjvG5HFREREZEIogLoorrl4dG/pMlT8Pj9LqeR3ig2xssVkwcwa0IBK7aW8vaqgxwqq2f+op28vmwfM8fmM3NsHrkZCW5HFREREZEIoALoEjsU0u2fcsZ8Xg+XXJjHDbNGsHDZXt5cvo+9JXW8uXw/by7fz8gBacwcm8dkM4dYv0YFRUREROTEVABd0mRtJ1BdhS8ri/jhI9yOIxHC4zGYZGYzcWQWOw/VsHRjMau2l7Hj4BF2HDzCX97ZwUWjc5k5Np8hecl6eqiIiIiIHEMF0CWdo3/TpmN49EAPOTuGYTByQBojB6Qx74qRrNpextINxewuruUf64v5x/piCrITmXlhHheP6Udygm4xFhEREREVQFeEWlqoW7Ma0OTvcv7iY31cOi6fS8flU1Rez9KNJXy05TBF5Q08t3gXL763mwkjspgxNp8xQzLweDQqKCIiItJXqQC6oH7dGuyWZuKGDcef28/tOBJFCrKTuPPyEdw2axgbdlWwdGMJm/ZUstoqZ7VVTnpyLJdcmMeMsXnkpMW7HVdEREREepgKoAu63v4p4gSf18MkM4dJZg7/v707jZHkvO87/q2jr7n3mOWeXO2K5LMUKXNJkaJJy7YSHYAEx4ZhJYKkOFYASVEQ50KcA7AARwKcGEggGwaMIIClOIGgOIjt5IVsCY6ow5TJpSSeJld8luK13EO7c+323N1d9eRFVXdX9/Ts7O50T+9M/z7Y2nrqueqZmWd66t9VXTVbXuGJF3/C9164yOUry3ztiTf42hNvYI5M8ICZ5IE7J9kzXuz3kEVERERkCygA3GK1K3MsnX4JgoDRh97d7+HIANg9VuQXHn0bH37kKGfOXuHxFy7wQzuFfesK9q0r/M9vvsLR20Z54K69PHDXJAf3DuvmMSIiIiI7lALALVZ+6hQ4x8h9JwlGRvo9HBkgvudx4uguThzdxSc+UOP5V6d55swUf/PaDG9emufNS/P8n8df57ZdJe6/a5IH7prk+MExfAWDIiIiIjuGAsAt5Jyj/MRfA3r2n/TXUDHkkXv288g9+6lUI06/McczZ6Z47sfTXJpb5htPneUbT51lfDjP/XcmZwZPHN1FGOiOtSIiIiLbmQLALbT61lkq58/hj4ww/M6f6vdwRADI5wJO3rmXk3fuJYpjXnnrKs+8MsWzZ6aYKa/ynecu8J3nLlAqhNz39j3cf9ck7zy+m2JeLx8iIiIi242O4LZQ/eYvow89jBfqWy+3nsD3G5eJfux9d3L20gLPnJnimVemOD+1yKnTlzh1+mMUcF4AABWwSURBVBJh4HPn4XHuPrqLu4/u4m0HRgn0PEsRERGRW56ikC3iooj5p54EdPmnbA+e53F0/yhH94/yyz93nEtzSzxzZopnz0zz6vmr/OjNOX705hwAxXzAXUcmGgHh4X0j+uygiIiIyC1IAeAWWTr9ElG5TG7/forHjvV7OCI37LZdQ3zo4aN86OGjlJcq2LNXGkHgpdklXnh1hhdenQFguBhyIg0G7z66i/27h3RnUREREZFbgALALVJ+snnzFx0Iy3Y3NpTnoRP7eOjEPgBmyyu8fHauERDOlld52k7xtJ0CYHwknwSDtycB4V49hF5ERESkLxQAboFoaYmFZ58BYOynH+nzaES6b/dYkUfvPcCj9x7AOcfUleVGMPjym3NcXahw6qVLnHrpEgB7xgocOzjO8QNjHDuQXGaqm8qIiIiI9J6OuLbAwtM/wFWrlMwJcnv29ns4Ij3leR77dg2xb9cQP3/yEM45LkwvNgJCe/YKM+VVZsqX+eHLl9M2cGjvMMcOjHHs4BjH9o9xaHJYj50QERER6TIFgFugfvdP3fxFBpHneRyaHOHQ5Ajvf/AIcey4MLPI6xfKvH6xzGsXy5y7vMi5qWR5/IWLAORCn6O3jaZB4SjHD4wxOVHSJdQiIiIim6AAsMeq01Msn7F4+Twj73qw38MR6Tvf9zg8OcLhyRF+9r6DAFSqEWcvL7QEhZfnlvnx+av8+PzVRtvhYsixg2Pcvm+UQ5PDHNo7zIE9Q+TCoF9fjoiIiMi2ogCwx8qnkkc/jNz/AEFJN74Q6SSfC7jj0Dh3HBpv5C0sV3kjDQbrgWF5qcqLr83y4muzjXq+53Hb7hKH9g4nZxr3DnNocph9u0p6NqGIiIhIGwWAPeScy9z989E+j0Zkexkp5bj3+B7uPb4HSH6fZsurvH6xzLmpBc5PL3J+apFLc0tcnEmWH6Z3HQUIA5+De4aSM4WZwHDPWFGXkYqIiMjAUgDYQyuvvUr10iWC8XGG7r6n38MR2dY8z2PPeJE940UeTB8/AVCtRVycWeL81CLnphc4P5UEhjPlFc5eXuDs5QXgUqN+IRewb1epsdy2a4h9E0l6YrSgB9iLiIjIjqYAsIfKp9Kbvzz8CF6gzyiJ9EIuDLj9tlFuv220JX95tcaF6UXOTy8mZwynknR5scJblxd46/JCh778RjCYLENJkDhRYvdYEd9XcCgiIiLbmwLAHomrVea//xSgu3+K9EOpEPL2Q+O8PfO5Qkg+Wzh1ZZlLc0tcnlvOLEuUl6rJpaXTi2v6C3yPyYlSchZyrMDu0SK7xgrsGSuye6zI7tEC+Zze6BEREZFbmwLAHln8mxeIFxfJHz5C4ciRfg9HRFIjpRwjpRzHDoytKVterSXB4JUkILyUCQ6vLFT4yewSP5ldumbfu9PgMAkMC80gcbTIxGheN6YRERGRvlIA2CO6+YvI9lMqhBzdP8rR/aNrylYrEVNXlpkurzBXXmF2fpWZ8gqz5VVmyyvMza+ysFxlYbnK2UtrLy8F8IDRoRxjwwXGR/KMD+fTdSFJN7bzlAqhblYjIiIiXacAsAeihQUWX3gePI+xhx/p93BEpAsK+YDD+0Y4vG+kY3nsHOXFCjPlFebSoHCmvMrs/AqzaaB4dbFCealKeanKuamO3TSEgd8SEI4P5xkZyjOansEcGco1zmaOlHIU84ECRhEREdmQAsAemP/BUxBFDN37TsKJiX4PR0S2gO95TIwUmBgpwMHOdWpRzPxSlfJihauLq1xdqHB1sbmUF1Yb6ZVKxEx5hZnyynXtP/A9RoZyzQCxJVDMM1IKGSrkGCqGDBVChoohpUKowFFERGTAKADsgebln7r5yyBxzhG7mBiHc3GSdo6YuFnmYhwuyXcxLrOdrNO6aZtGWdrXxajI3JVFHK61PLPtGn1k85rrjvkd8nAQ4yDTjkZ9WuqT1Gq2zZbjSKpk+6DRb7oHkn/Zvhqptv7I9FGv08xr/O8yZR32Ue8zW6veT2teff+NGq37ds10tqy+v/ZRtLR3DkJw48B4Um8UGAXi2BHFjiiO03RM5BxxDHGczq04nUtx0nYFWAEaJxcjB/Mkyzoxngd4XvKYDc9LAknSPN/30nKvpU57m41kv+cdy117ueuQurm+OxVv2GYTfd/Mfm6y+9623lTz9Rt7ntfh5909m/2e99P2Hfmto9fzq6u2yzilYf/oPv7FyX9MIcj3eyibpgCwy5bOnWfltdfwi0VGTt7f7+HccpxzRC6iFteoxRE1l6yjuEbkYmqulh7oRkRxRM1FxC5K6qR5UXY7kxe7mMglbeM4k86u4+Z2lNaPXdQIyLLpKA3QokYwl1lozc8GPiI95QHB2qzNnsNz6RJfq1BERGRAzVcWiF3U72F0hQLALpv6zncBGHnXQ/iFwg21ddVViGsQ5MAP8Xp0t8BaXKMSVanEFSpRlWpcpRJVqMZVqnGNalyjFlWpxDVq2by4SiWuUotrVKNmXrJOt12tEdxFcY1aI9hLlx3yi7MeDw/f8/E9D8/z8UnSvufjeR4+ab7n41PPz5alefU0aZnn4eFRLOSoVmJ8zwePZnla31uz9httW/LTdPLQ86R98i9Tb01fSTmNtN+SX7+MsNkePM9vfF/qdaF+Fimtlabr/+O1bLX0R0sfzZCnY3nalkxfNPKaOa1joKV+0rytxFu/XvP/5tfZoQeyFbMjoS3Vfmlmth861GmvvSavQ91a5FitRqxUaiyvRuQLOS5PLza2Vyo1llZrrFYiVqtxuk6XNF2pRLgbDEGdu776ge+RDwNyoU8Y+ORCj1xmOx/4hLmAfOARhj650CcXBORDnzD0yAU+Yb1tkKyDwCcXeM106JPzPcIgIAjTMt8jCH0Cb4PX4Q2+jE4/sxvR19abat658d69I0xPd75JUrds7wuat/fo+21y7whTPZ5f3aTL77eXg/t2MTOz9jFR21FPA0BjzMeBzwE54PestX/QVn4S+ENgDPgr4LPW2lovx9RLLo4bAWCnu3+6uIZbmCUuTxHPT+Hmp4nLl4nnp3HzU7iV+UbdGKj4PpUwZDUIqYYhq0HAahBQCQJWAz8p930qHqz4HhXPo+pBxYOq56jgqAJV4jSdrOM+v5Xvez6hH5LzQgI/IPRDQi8g8AOCzHptnk/ghQS+31qWpn3PT7eTgzbfCwh8v5HfUp6pnwRo9e00OMssgefj4Tf6SoK67NIatPXS5OQoU1PzG1cUuUk3M8di56hUI1Yq9aUZPNbzKmnQWKnGjQCy0mE7SceNdORgeRWWG3tzwNb9mQh8jyDwCH2fMPAIAp/AT4LH+nbop+t6UJnZDvz64jf6Cnwf30+DzMx2s7y1jZ8u9XTgtW77XqYss/a8ejmNes016RtVW3sAOpQrUQq37Z95ucUVc0WKYbXfw5Adyt9Bj3HqWQBojDkE/DbwLmAVeMIY821r7elMta8An7LWnjLGfAn4NPBfejWmXlt+5QyrU9N446Ncrr3OyvefZmVpluWlOZZX51mtLrHiQcX3WPF9Vn0vWfIeK5MhFX9Psu15VP1r/VF2QJQuN853jlzsyDtHzkEuduRcusSO0EHOOcJ0yTkI40w6rRum+fW6QVoWOgjq5c4R0MwLHDR+fTwPPB+8AHyf5JRWckasWdaa9rL5/to69TNUa9unF8m19NGa32yXniFqpNvrNvuPgXhNfW9t/fQsU3uZ16mssW6vD+XREpWF1ca2l61bb9tId8rL7GudPG9Nnfp/bdstB45t+21Pr1vedpaqY/tr5LWe3qJFx/ad6nYqXz/fu6627WUbHWR3KO/YpFO9GzmA37h9rVAhXuzwDucG+8kDeR/GikCxvq9cuqyz7w36d84RRTGVyFGtxlSjiGotplqLqdUclSiiWo2pRTGVWkytFlOJ0nUtphZFVGqOWi2mFsfUIket5qhFMVGc9hOl+S3rZr7LxJs1bib03FyA1eu36zwywWEaGNavEvD95Cx2PVCsfy7UT8+0N9L1Om31fY/MOvndyRdComqU/vo36zTbNfcP9e3k5Z56vfrZ9Uy6cTWC12xTb082v6VO/coB0jfwmuVk07Sm1yvLvuxmz/R7LW0y48n+HDLjbX359chWXFOPtfUyr9aZvxXtL68ebVXw6n+d219e27Y7v8yvneftv9LZ71dm4B3yvPbijv219wuwWlnhylzb61f7ODYYaHv5NV/a1/n9vt6X5c2+AXOt5hv2vMG+NzeyDXa96Q56Obr1TUwU+7LfXujlGcD3A9+y1s4CGGP+BPgI8IV0+yhQstaeSuv/EfB5tnEA+OLXvswE8NThiCcvf6dZUASKPtD59vHryft5imGegp+n4OcoBDkKXo68H1L0cxS8gIIXkPcCCvgU8Ml7Pjnnkccjh0fOQd6Rrj1yzhE4l1xqGke4OIL64iJcHKfp5rpRp54XxziXrEnXrl7u4jS/mXYuSj7sXK/v4mTbuWSbqCWWvdYBz6B/DGm13wOQHe9WvbglTJdSr3bgp0tuo4rS0PEDoyL9UwGG+j0I2bFOuxLjH/sdhsfG+z2UTetlAHgQuJjZvgi8e4Pywzeygz17biyg6rXcgT3MXpjhzeND7CdHKSxQypUYKo5QKo4xMjRBqVCiGBYZyhXXrnNFSmGRYligEOaTz3ntUM65ZoAYJwGiywSPrWmXBqLNwLKZdkka19xutHXN/dTzW7ZdM1B1tJQndxFzmUC1Po5mvmsJZJt3y7zW9rXyG3egbKlLc6wteXHmrpT1u022tWnbztZzHfPb+qj/nNrL16nfGH8jP9uGzv1mijvVy/bXfrfNjndQy+6/tfO17dfrY02+65hs7fsab02s19c1sjrXW5u37l47jqe7+74x127f07v2barvTY5r01/W1rzl1Xm6XOu2Vq2/E9czynV+da5V6zp/dO4aWxvt73p77V4D1yHVM+u9vHXIW/98ysbj3NFvzG7TL26bDvuWt+IPc+f+XQwND/d7KJvWywDQp3UOerS+X7hR+YZmZhaI41tnmt//8X/L5D8f5adv5jNa6YmwGrBAlQUG+Rr2+lvxG1QZQNnPZ93ERYUiG9LnTKWXNL+klzS/pJfekc6vxaVbZ475vndTJ8R6eRh9DjiQ2d4PXLiBchEREREREemiXgaA3wTeZ4yZNMYMAb8CfKNeaK19E1gxxtSflv6rwNd7OB4REREREZGB1rMA0Fp7HvhN4NvAc8BXrbXfN8b8hTHmwbTaJ4DfNca8THKHlN/v1XhEREREREQGXU+fA2it/Srw1ba8D2fSz9N6YxgRERERERHpkQG9lYaIiIiIiMjgUQAoIiIiIiIyIBQAioiIiIiIDAgFgCIiIiIiIgNCAaCIiIiIiMiA6OldQHsoAPB9r9/j6OhWHZfsDJpf0muaY9JLml/SS5pf0ku32vzKjCe4kXaec677o+m99wCP93sQIiIiIiIiffazwPeut/J2DQALwEPARSDq81hERERERES2WgAcAH4ArF5vo+0aAIqIiIiIiMgN0k1gREREREREBoQCQBERERERkQGhAFBERERERGRAKAAUEREREREZEAoARUREREREBoQCQBERERERkQGhAFBERERERGRAhP0ewE5ijPk48DkgB/yetfYP+jwk2QGMMWPAE8AvWGvfMMa8H/giUAL+l7X2c30doGxbxpjfAv5euvnn1tp/o/kl3WKM+QLwEcABX7LWflHzS7rNGPOfgb3W2k8aY04CfwiMAX8FfNZaW+vrAGXbMsZ8G9gHVNOsfwS8nR1wrK8zgF1ijDkE/DbwHuAk8BljzDv6OyrZ7owxDwPfA+5Kt0vAl4FfAu4GHjLGfKh/I5TtKj0Q/yBwP8lr1ruMMR9D80u6wBjz88DfBn4KeBD4p8aY+9D8ki4yxrwP+LVM1leAX7fW3gV4wKf7MjDZ9owxHsmx133W2pPW2pPAOXbIsb4CwO55P/Ata+2stXYR+BOSdz5FNuPTwD8BLqTb7wZesda+nr6r+RXg7/ZrcLKtXQT+lbW2Yq2tAj8i+WOn+SWbZq39LvC30nm0j+SKowk0v6RLjDG7SQ7G/0O6fRQoWWtPpVX+CM0vuXkmXf+lMeZ5Y8yvs4OO9RUAds9BkgOquovA4T6NRXYIa+2nrLWPZ7I0z6QrrLUv1Q+UjDF3klwKGqP5JV1ira0aYz4PnAYeQ69f0l3/FfhNYC7d1vySbtpF8rr1y8D7gM8Ct7ND5pgCwO7xST7nUOeRHEyJdJPmmXSVMeYe4P8B/xp4Dc0v6SJr7W8Bk8ARkjPMml+yacaYTwFvWWsfy2Tr76N0jbX2SWvtP7DWXrXWTgNfAr7ADpljCgC75xxwILO9n+ZleyLdonkmXWOM+RmSdzj/nbX2v6P5JV1ijDmR3pADa+0S8GfAe9H8ku74KPBBY8xzJAflvwh8Cs0v6RJjzHvSz5jWecAb7JA5pruAds83gX9vjJkEFoFfAT7T3yHJDvQUYIwxdwCvAx8nuamCyA0xxhwB/i/wUWvtt9JszS/pluPA540x7yF5x/yXSC7Z+0+aX7JZ1toP1NPGmE8C77XW/kNjzIvGmJ+x1v418KvA1/s1Rtn2JoAvGGMeJbnj568Bfx/4yk441tcZwC6x1p4nuRb928BzwFettd/v76hkp7HWrgCfBP6U5HM1L5N8CFnkRv0GUAS+aIx5Ln0n/ZNofkkXWGv/Avhz4FngaeAJa+0fo/klvfUJ4HeNMS8DI8Dv93k8sk1Za79G62vYl9M3FnbEsb7nnNu4loiIiIiIiGx7OgMoIiIiIiIyIBQAioiIiIiIDAgFgCIiIiIiIgNCAaCIiIiIiMiAUAAoIiIiIiIyIPQcQBERkTbGmL8keU7d/wB+w1p7us9DEhER6QoFgCIiImt9AMBa++F+D0RERKSb9BxAERGRDGPMfyN5YPmLwDuAh0keKv0fgbOAARaB3wH+Wbr9p9baf5m2/zvA54A8sERyBvHJrf0qREREOlMAKCIi0sYY44BJ4IfAR0gCwG8CD1lrnzXGfB0YB94LjAEXgLcBw8CfAe+11s4YY+5J291hrV3c6q9DRESknS4BFRERuT6vW2ufTdOvAlettRVg2hhTBnYDPwccAB4zxtTbxcAdwPNbPF4REZE1FACKiIhcn9W27WqHOgHwmLX2o/UMY8wRkjOEIiIifafHQIiIiKwVAbmbaPcY8EFjzAkAY8yHgReAUhfHJiIictN0BlBERGSt/w18l+Szf9fNWnvaGPMZ4I+NMR5QA37RWrvQgzGKiIjcMN0ERkREREREZEDoElAREREREZEBoQBQRERERERkQCgAFBERERERGRAKAEVERERERAaEAkAREREREZEBoQBQRERERERkQCgAFBERERERGRAKAEVERERERAbE/wfChIzsERIR9wAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "# plot results\n", "plt.figure(figsize = [15,5])\n", "plt.plot(t,y1,linewidth=2,label='F1')\n", "plt.plot(t,y2,linewidth=2,label='F2')\n", "plt.plot(t,y3,linewidth=2,label='F3')\n", "plt.plot(t,y4,linewidth=2,label='F4')\n", "plt.xlabel('time')\n", "plt.ylabel('y(t)')\n", "plt.legend()\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": 11, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAA4wAAAFFCAYAAACnsmP3AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJzs3XmUXOV95/93VVev6tbekloSWhDwgNjEYsBsBgNe8B6wHcOAHcdmiMeTZJJMMvMzzow98XFm4sTYx5DJhOTYDiZ2bMfxCrYxigGziFVi0yOQhPaltbZ6X6p+f3R1q2kk1JL69q3qfr/O0emq5966/S3rOZgP3+c+N1MoFJAkSZIkabhs2gVIkiRJkkqTgVGSJEmSdEgGRkmSJEnSIRkYJUmSJEmHZGCUJEmSJB2SgVGSJEmSdEgGRkmSJEnSIRkYJUmSJEmHZGCUJEmSJB2SgVGSJEmSdEgTJTDmgEXFn5IkSZI00RxTJpooAWoh8ApwGbA55VokSZIkaazNBx4CTgLWjvRDEyUwNhV/PpRqFZIkSZKUriYMjK+zDWDv3jby+ULatbzOjBn17N7dmnYZGqecX0qS80tJcn4pSc4vJakU51c2m2HatElQzEYjNVECYx9APl8oycAIlGxdGh+cX0qS80tJcn4pSc4vJamE51ff0Zw8UTa9kSRJkiQdJQOjJEmSJOmQJsqSVEmSJEk6LoVCgdbW/XR0tJLPH35l586dWfL5/BhW9lq5XBXTpjVSUXH8cc/AKEmSJEkjsHdvM5lMhunTZ1NRkSOTyRzyvFwuS29vOoGxUCjQ1tbC3r3NzJzZdOQPHIFLUiVJkiRpBLq7O5k6dQa5XOVhw2LaMpkMkyZNpre3e1SuZ2CUJEmSpBEpkMmUfoQazTBb+t9WkiRJkpQKA2MZyxcK9PQe1WNUJEmSJI0jl156fqLXT3TTmxDCDcBtQCVwe4zxjmHHlwF3AZOBB4FbY4y9IYQFwN3ALCACN8YYW0MI04BvAfOALuCWGOOzSX6HUvZX9zzD7pZOPv+7F1BT5f5FkiRJkkZXYh3GEMI84AvApcAy4JYQwtJhp90NfDrGeAqQAT5ZHL8TuDPGeCrwJPDZ4vgfAc/FGM8G/hfwtaTqL3WFQoG1W/eza38nz76yK+1yJEmSJI1DSbalrgYeiDHuAQghfA+4Hvh88f1CoDbG+Fjx/K8Dnwsh3AVcDrx/yPivgT8DKoCG4vgkoCPB+ktaT2+e3r4CACte3MlFS+ekXJEkSZI0sdz+3ZWsWrs7kWuftWQGf/jBsxO59tFIMjDOBbYNeb8NuOAIx+cDM4GWGGPvsHGALwGPhRC20r+M9ZqjKWjGjPqjOX1MNTY2HPmkIfa0dA6+fn79bmonVVNfVzXaZWmcONr5JR0N55eS5PxSkpxfOlo7d2bJ5Q4u0kzyyRqZDK/5XW/kUOdls9lRmeNJBsYsUBjyPgPkR3B8+DhDPvc14Gsxxq+GEN4MfCeEsDTG2DqSgnbvbiWfH37p9DU2NtDcfOCoPrN1V9vg696+Ar94dD2XnTV3tEvTOHAs80saKeeXkuT8UpKcXzoW+Xye3t6DkeYPrj90BzCXy77mvGM10msc6rx8Pv+aOZ7NZo6pgZbkLqmbgaYh7+cAW0dwfCcwJYRQURxvGvK59wH/CBBjfBTYAZw26pWXgfau3te8X/HSzpQqkSRJkjReJRkY7weuCiE0hhDqgOuA+wYOxhg3AJ0hhEuKQzcB98YYe4CHgA8Xx28G7i2+Xknx3sYQwsn0L2tdk+B3KFntnf2BcXFTAxXZDC+9upeW9u6Uq5IkSZI0niQWGGOMW4DPAMuBZ4F7YowrQgg/CyEMPCzkRuDLIYTVQD3w1eL4p+jfVfVF4DL6H80B8FHg4yGE54FvAx+NMe5P6juUsvauHgAap9Zy+uLp5AsFnorNKVclSZIkaSw9/PCTiV4/0Yf3xRjvAe4ZNnbtkNcree1GOAPjG4ArDjH+MvDWUS+0DA10GOtqKjlryQxWrd3Nihd3cOU581KuTJIkSdJ4keSSVCVoMDBW5zjn5EZyFVnWbNrH3gNdKVcmSZIkabwwMJapgU1v6mpy1FbnOHvJDArAE6vd/EaSJEnS6DAwlqmhHUaAC5bOBmDFSztSq0mSJEka3zIUCsf/uIykFQqj9yhBA2OZGtphBDhryQyqKytYt7WF5n0daZYmSZIkjUtVVTXs27eL3t6eUQ1lo6lQKNDW1kIuVzUq10t00xslp6Ozf5fUgQ5jdWUF55w8k8de3MGKl3bwrjcvSrE6SZIkafyZNq2R1tb97Nmzg3y+77DnZbNZ8vn0OpG5XBXTpjWOzrVG5SoacwMdxtqag3+FF5w2uxgYdxoYJUmSpFGWyWRoaJhKQ8PUNzyvsbGB5uYDY1RVslySWqYG7mGcVFM5OHb64unUVefYtLOVrbva0ipNkiRJ0jhhYCxTg/cwVh/sMFbmspwb+lvPbn4jSZIk6XgZGMtQoVAY7DDWVr92VfGFpw3slrqzZG/ElSRJklQeDIxlqLs3T1++QFUuS2XutX+Fpy6cSkNdJdv3tLNpZ2tKFUqSJEkaDwyMZWiwu1jz+j2LKrJZzg+zAHjcZamSJEmSjoOBsQy1D3ukxnAXnNYfGFe86LJUSZIkScfOwFiGBje8OUSHEeDkE6Yytb6K3S2drNvaMpalSZIkSRpHDIxlaGBJal115SGPZzMZLihufuOyVEmSJEnHysBYho7UYQQGA+MTq3eSz7ssVZIkSdLRMzCWoYMdxsMHxsVNDcycUsP+1m5e3rxvrEqTJEmSNI4YGMvQSDqMmUyGC5cOLEvdOSZ1SZIkSRpfDIxlqKPzyIERDi5LfXL1Tnr78onXJUmSJGl8MTCWofauN36sxoD5jZNomlFHa0cPqzfsHYvSJEmSJI0jBsYyNHgPY82hd0kdkMlkuNDdUiVJkiQdIwNjGRq8h/EIHUaAN502C4Cn1zTT0+uyVEmSJEkjd+TEcRxCCDcAtwGVwO0xxjuGHV8G3AVMBh4Ebo0x9oYQFgB3A7OACNwYY2wNITw5pOZaYAkwL8Y4odpnbSO8hxGgacYkFsyuZ+OOVp5ft5tzTmlMujxJkiRJ40RiHcYQwjzgC8ClwDLglhDC0mGn3Q18OsZ4CpABPlkcvxO4M8Z4KvAk8FmAGOP5McZlMcZlwOPAn0+0sAhDNr0ZQYcRcFmqJEmSpGOS5JLUq4EHYox7YoxtwPeA6wcOhhAWArUxxseKQ18HPhhCqAQuL54/OD70wiGEq4Czgf+dYP0la2BJau0IOowAbzq1f1nqs6/soqu7L7G6JEmSJI0vSQbGucC2Ie+3AfNHcHwm0BJj7D3M5wA+B3wmxjjh0k+hUDi46c0IO4wzp9ayZN5kunvyrFy7K8nyJEmSJI0jSd7DmAUKQ95ngPwIjg8fZ+jnQginAzNjjD852oJmzKg/2o+MmcbGhhGd19HVS75QoLqqgqY5U0Z8/beev4C1W57n2bW7edflJx1rmSpTI51f0rFwfilJzi8lyfmlJI2X+ZVkYNwMXDbk/Rxg67DjTYc4vhOYEkKoKHYQm4Z97v3Ad46loN27W8nnh2fR9DU2NtDcfGBE5+5p6QSgtqpixJ8BOHX+FDLAky/tYMOmvSPaMEfjw9HML+loOb+UJOeXkuT8UpJKcX5ls5ljaqAluST1fuCqEEJjCKEOuA64b+BgjHED0BlCuKQ4dBNwb4yxB3gI+HBx/Gbg3iHXfXPx+IQ0+EiNIzyDcbhpDdWEBVPp7SvwzMvNSZQmSZIkaZxJLDDGGLcAnwGWA88C98QYV4QQfhZCOL942o3Al0MIq4F64KvF8U/Rv6vqi/R3KW8bcukT6e9OTkjtR/FIjeEucLdUSZIkSUch0XWJMcZ7gHuGjV075PVK4IJDfG4DcMVhrjn80RwTymCHcYQb3gx1Xmjk7l+s4cX1eznQ3k1DXdVolydJkiRpHElySaoS0HEcHcaGuiqWLp5GvlDgqeiyVEmSJElvzMBYZo6nwwhwYXFZ6gqXpUqSJEk6AgNjmWnr7AGOrcMIcM7JjeQqMsSN+9h7oGs0S5MkSZI0zhgYy8zgpjfVR7dL6oC6mhxnnjiDAvBk3DmKlUmSJEkabwyMZebgYzWOfb+iC5e6LFWSJEnSkRkYy8zgpjfHeA8jwNlLZlJVmWXtlhZ27esYrdIkSZIkjTMGxjIz0GGsPY4OY3VVBctOmgn4TEZJkiRJh2dgLDPto9BhBLjo9DkAPLhyK/lC4bjrkiRJkjT+GBjLTHvX8e2SOuCsE2cwY3INzfs6eX7d7tEoTZIkSdI4Y2AsMwMdxkk1x7ZL6oBsNsOV584D4IGntxx3XZIkSZLGHwNjGSkUCgfvYayuOO7rXXZWE7mKLM+t3c2Ove3HfT1JkiRJ44uBsYx0dvdRKPRvWlORPf6/uoa6Ki5cOosCsNwuoyRJkqRhDIxlpKNrdDa8Geqt584H4OFV2+jq6Ru160qSJEkqfwbGMtI2sEPqcW54M9TipsmcOHcy7V29PP6ij9iQJEmSdJCBsYy0dxZ3SB3FDiPAVcUu4wNPbabgIzYkSZIkFRkYy0h7AktSAc4/dRYNdZVs3NnKK1v2j+q1JUmSJJUvA2MZaU9gSSpAZS7L5WfPBeBXT20e1WtLkiRJKl8GxjJysMN4fM9gPJQrz5lHJgNPxWb2tXaN+vUlSZIklR8DYxnpKHYYa0e5wwgwfXIN55zcSF++wIPPbh3160uSJEkqPwbGMpLUPYwDrjp3HgDLn91Cb18+kd8hSZIkqXwYGMtIUvcwDjh14TSaZtSxv7WbZ17elcjvkCRJklQ+DIxlZKDDOCmhwJjJZHhr8REbbn4jSZIkKZnkURRCuAG4DagEbo8x3jHs+DLgLmAy8CBwa4yxN4SwALgbmAVE4MYYY2sIYTLwt8DS4iV+N8b4dJLfoZQk9RzGoS4+Yw7f//Va1mzax+adrcyfVZ/Y75IkSZJU2hLrMIYQ5gFfAC4FlgG3hBCWDjvtbuDTMcZTgAzwyeL4ncCdMcZTgSeBzxbH/wbYFGM8B/jv9IfHCWPwHsaa0d8ldUBtdY6Lz5gDwANP22WUJEmSJrIkl6ReDTwQY9wTY2wDvgdcP3AwhLAQqI0xPlYc+jrwwRBCJXB58fyh4xngOuAvAWKM9wEfT7D+ktOe4C6pQw0sS33khe2DXU1JkiRJE0+SyWMusG3I+23ABUc4Ph+YCbTEGHuHjc8CuoBPhRDeA3QA/+VoCpoxo3SXVzY2NhzxnI7uPgAWzJtKQ11VorWcffJMVr68i5Xr9/Ley5ck9rs0NkYyv6Rj5fxSkpxfSpLzS0kaL/MrycCYBQpD3meA/AiODx+nOJ4DZgP7Y4xvDiFcA/wAOHGkBe3e3Uo+P/zS6WtsbKC5+cAbnpMvFGjv6O/2tR3ooLOtK9GaLj2jiZUv7+JHD67lwlMbyWYyif4+JWck80s6Vs4vJcn5pSQ5v5SkUpxf2WzmmBpoSS5J3Qw0DXk/B9g6guM7gSkhhIrieFNxfBfQC9wDEGP8JVAfQpiVSPUlprOrjwJQU1VBRTb5zW2XnTyD6ZOr2bG3gxfX70n890mSJEkqPUkmj/uBq0IIjSGEOvrvP7xv4GCMcQPQGUK4pDh0E3BvjLEHeAj4cHH85uJ4F/BL4LcBQggXAW30B8lxr72ruENqwvcvDqjIZrnynHkAPPD0ljH5nZIkSZJKS2KBMca4BfgMsBx4FrgnxrgihPCzEML5xdNuBL4cQlgN1ANfLY5/iv5dVV8ELqP/0RwAvwu8M4TwPP07pP52jHHoMtdxa2DDmyQfqTHcZWfPJVeRYeUru2je1zFmv1eSJElSaUg0fcQY76G4hHTI2LVDXq/ktRvhDIxvAK44xPg24L2jXmgZ6Oga+8A4ua6KN506m0df2M7yZ7bwoStPGrPfLUmSJCl9yd8Mp1Ex2GFM8BmMh/LW8/qXpT60civdPX1j+rslSZIkpcvAWCbaBzqMY3QP44ATmyazaE4DbZ29PP7SjjH93ZIkSZLSZWAsE2ncwwiQyWS46rz5ADzw1BYKhdJ7LIkkSZKkZBgYy0RaHUaAC06bRX1tJRt2HGDd1pYx//2SJEmS0mFgLBNtncXHaoxxhxGgMlfBZWf3PzLzV09vHvPfL0mSJCkdBsYy0VFcklqbQocR4Mpl88hk4MnVO9nf1p1KDZIkSZLGloGxTAwuSa0e211SB8ycWsvZS2bS21fgwZVbU6lBkiRJ0tgyMJaJg4/VSKfDCAxufvPvz2yhL59PrQ5JkiRJY8PAWCYOdhjTC4ynLZrGnOl17D3QxbMv70qtDkmSJEljw8BYJkqhw5jNZLjy3HkA/OopN7+RJEmSxjsDY5lI87EaQ11yRhPVlRWs3riPLbvaUq1FkiRJUrIMjGUgXyjQWQyMtVXpBsa6mhxvPmMOAA/4iA1JkiRpXDMwloHOrl4KQG11jmw2k3Y5vLW4LPWR57YPPh9SkiRJ0vhjYCwDg/cvprjhzVDzG+tZumgaXT19/HzFprTLkSRJkpQQA2MZaCuBDW+Ge/+lJwLwyyc3caC9O+VqJEmSJCXBwFgGSuGRGsOdNH8KZ544g67uPu59fGPa5UiSJElKgIGxDJTCIzUO5QOXLwbggac2s6+1K+VqJEmSJI02A2MZaO/q31imlDqMAIvmTObcUxrp7s3z00c3pF2OJEmSpFFmYCwDHcUOY22JdRgB3n/pYjLAr5/dwp6WzrTLkSRJkjSKDIxloBTvYRwwf1Y9bzptFr19BX78yKtplyNJkiRpFBkYy8DBexgrU67k0N536WIyGXh41TZ27utIuxxJkiRJo8TAWAZKucMI0DRjEhefMYe+fIEfP7w+7XIkSZIkjZJEE0gI4QbgNqASuD3GeMew48uAu4DJwIPArTHG3hDCAuBuYBYQgRtjjK0hhLcA/woMPC3+mRjj7yT5HUpBqe6SOtR7L1nMYy/s4JEXtnPtmxfSNGNS2iVJkiRJOk6JdRhDCPOALwCXAsuAW0IIS4eddjfw6RjjKUAG+GRx/E7gzhjjqcCTwGeL4+cDX4oxLiv+GfdhEQ52GCeVcGBsnFrLZWc1USjAD+0ySpIkSeNCkktSrwYeiDHuiTG2Ad8Drh84GEJYCNTGGB8rDn0d+GAIoRK4vHj+4Hjx9ZuAt4UQVoUQfhRCOCHB+kvGQIextkSXpA5498WLyFVkWfHSTjbtbE27HEmSJEnHKckEMhfYNuT9NuCCIxyfD8wEWmKMvcPGAfYB/xJj/NcQwq3At4FLRlrQjBn1R/UFxlJjY8Nhj3X29AFwwtypNE6vG6uSjlpjYwPvvHgRP35oHfeu2MhnfufCtEtS0RvNL+l4Ob+UJOeXkuT8UpLGy/xKMjBmgcKQ9xkgP4Ljw8cZ+FyM8daBgRjj/w0h/GUIYUqMcf9ICtq9u5V8fvil09fY2EBz84HDHm9t7wago62L5r6+sSrrmLz17CZ+/uirPPb8dlas2sLipslplzThHWl+ScfD+aUkOb+UJOeXklSK8yubzRxTAy3JJambgaYh7+cAW0dwfCcwJYRQURxvAraGELIhhM8MGR/QyzjWl8/T2d1HBqipHv7VS8+U+mquOq+/IfyDh9alXI0kSZKk45FkYLwfuCqE0BhCqAOuA+4bOBhj3AB0hhAGlpTeBNwbY+wBHgI+XBy/uTieBz5QvA4hhJuBx4v3R45bHV39HcXa6hzZTCblakbmnRctpKaqgufX7eHlzfvSLkeSJEnSMUosMMYYtwCfAZYDzwL3xBhXhBB+FkI4v3jajcCXQwirgXrgq8XxT9G/q+qLwGX0P5oD4KPAH4YQXgB+B/hEUvWXisFnMJbwDqnD1ddWcs35/fsR/eBBu4ySJElSuUo0hcQY7wHuGTZ27ZDXK3ntRjgD4xuAKw4x/gJw8agXWsI6Bp7BWOI7pA739gtO4FdPbWb1xn289OoeTls0Pe2SJEmSJB2lJJekahS0d/YA5dVhBKirqeQdFy4A4F8fWkehUHqbDUmSJEl6YwbGEjewJLXUn8F4KFefP5/62krWbmnhuXW70y5HkiRJ0lEyMJa49s7yu4dxQE1VjmsvWgjADx5cb5dRkiRJKjMGxhI30GGcVFOZciXH5spz5zGlvooNOw7w9JpdaZcjSZIk6SgYGEtce5luejOgurKCd795EQD/9tA68nm7jJIkSVK5MDCWuIHAWFuGS1IHXH72XGZMrmbLrjZWrN6RdjmSJEmSRsjAWOLau4q7pJZphxGgMpflPZcsBuCHD79KXz6fckWSJEmSRsLAWOLKedOboS4+Yw6zptayY087jz5vl1GSJEkqBwbGEjew6U05dxgBchVZ3ndpf5fxR79ZT2+fXUZJkiSp1BkYS9xgYCzTXVKHunDpbJpm1LFrfycPr9qWdjmSJEmSjsDAWOLKfZfUobLZDO+/7EQAfvzIq/T09qVckSRJkqQ3YmAscQc7jOUfGAHOC42cMKuevQe6WP70lrTLkSRJkvQGDIwlrC+fp6u7j0wGqqsq0i5nVGQzGT5weX+X8Ye/Wc/eA10pVyRJkiTpcAyMJayjq3/JZl11jmwmk3I1o+fsJTNYdtJMOrr6uOf+NWmXI0mSJOkwDIwlrL2z/xmMtePg/sWhMpkM/+Ftp1BdVcFTsZlnXm5OuyRJkiRJh2BgLGED9y9OGgc7pA43fXINv1XcAOfuX6yho/hdJUmSJJUOA2MJa+scXxveDHfVefNZNKeBvQe6+MFD69IuR5IkSdIwBsYS1jGOHqlxKNlsho++41SymQy/emoz67e1pF2SJEmSpCHeMImEEKqBW4DfAgLQB6wGvgt8PcbYnXiFE9jAktTacdphBFg4p4G3vekE7luxkW/cu5rPfux8KrL+dwxJkiSpFBz238xDCG8BngBOA74KXAu8A/gKcBbwZAjhrWNR5ETVPs47jAPed+liZkyuYePOVn75xOa0y5EkSZJU9EZJ5D3ApTHG4esEXwB+EkKYDPwP4IGkipvo2rv6d0kdr/cwDqiuquCmtwdu/+5K/u3hdZwfGpk5tTbtsiRJkqQJ77BJJMb4JwAhhHkxxi1Dj4UQlsYYXwT++I0uHkK4AbgNqARujzHeMez4MuAuYDLwIHBrjLE3hLAAuBuYBUTgxhhj65DPzQdWAefGGF8d4XctOxOlwwhw1pIZXHDaLFa8tJN/+sUa/vCDZ5EZR8+elCRJksrRGy1JnR5CmA78LIQwbeB9CGE28K9HunAIYR7wBeBSYBlwSwhh6bDT7gY+HWM8BcgAnyyO3wncGWM8FXgS+OyQ62bpD5lVI/2S5WrgHsbx3mEc8JGrTqauOsdz63bzxOqdaZcjSZIkTXhvtLvIPwO7gDOB3cXXu4BNwNMjuPbVwAMxxj0xxjbge8D1AwdDCAuB2hjjY8WhrwMfDCFUApcXzx8cH3LdPwXuL9Yyrh3sMI6/5zAeypT6aj545RIA7vnlGto6e1KuSJIkSZrY3mhJ6tsBQgj/GGP8+DFcey6wbcj7bcAFRzg+H5gJtMQYe4eNE0I4D3gr/ZvvfPoYaiorE63DCHDZ2XN55PntvLx5P99dvpaPvfPUtEuSJEmSJqzDJpEQwtUxxvvfKCyGEN4WY/zFYQ5ngcKQ9xkgP4Ljw8cB8iGEOvqXqn4wxpgPIRyurMOaMaP+qD8zVhobG1431t3b/z/XvKYphzw+Xv2XG87j9/96OQ+u3Mq1l57I6SfOSLuksjeR5o/GnvNLSXJ+KUnOLyVpvMyvN2pdXR1C+DTw/4BfDHT8QghVwJXAf6b/mYyHC4ybgcuGvJ8DbB12vOkQx3cCU0IIFTHGvuI5W4vXmg38qBgW59J/f+UHYoxxBN+V3btbyeeHZ9H0NTY20Nx84HXjB9r6H3PZ1d51yOPjVU0W3nnhQn78yKt85dtP8z9/5wIqcz6b8Vgdbn5Jo8H5pSQ5v5Qk55eSVIrzK5vNHFMD7bD/Fh5j/G/AF+kPhntDCGtCCK/Qfz/jHwBfGNhJ9TDuB64KITQWu4PXAfcNuf4GoDOEcElx6Cbg3hhjD/AQ8OHi+M3F8Z/HGBfFGJfFGJfRHyKvHWlYLEcD9/BNpCWpA9598UJmT69j2+527n18Q9rlSJIkSRPSG7ZtYoyPAx3Ax4EbgA8B82OM18YYHz3CZ7cAnwGWA88C98QYV4QQfhZCOL942o3Al0MIq4F64KvF8U/Rv6vqi/R3Fm87pm9Xxnr78nT35MlmMlRXVqRdzpirzFXw0bf3Lzv+ySOvsm13W8oVSZIkSRPPSFpX36d/g5km+penvjrSi8cY7wHuGTZ27ZDXK3ntRjgD4xuAK45w7UUjraMcDd3wZqI+j/DUhdO49MwmHn5uG//088h//cg5E/Z/C0mSJCkNR7wxLMb4rRjjW4D3ArOAJ0II/xRCeF3Q0+jpGHykxsRbjjrUh956EvW1lazeuI+Hn9t25A9IkiRJGjUj2kkkhJAFTgZOob8ruRO4M4TwuQRrm9AGOoy1E/D+xaHqayv5yFUnA/AvD7xCS3t3yhVJkiRJE8cRA2MI4S+ATcCfAt8BToox/jHwFibAsxDT0m6HcdBFp8/m9EXTaOvs5Tu/ejntciRJkqQJYyQdxln070Z6WYzxn4u7mBJjbAM+kmh1E9jQexgnukwmw01vD1Tmsjz6wg5eWL8n7ZIkSZKkCWEk9zDeUtyc5lDHDvcMRh2n9oFHathhBGDWtDree8kiAL7589V09fSlW5AkSZI0Afg09BJlh/H13n7BAuY3TqJ5Xyc//s2raZcjSZIkjXsGxhLlPYyvl6vI8tF3nEoG+PmKjazdsj/tkiRJkqRxzcBYog52GCtTrqS0LJk3hWvedAJ9+QJ3/tvz7poqSZIkJcjAWKIGO4wuSX2d669YwpJ5k9l7oIu//9EL5POFtEuSJEmSxiUDY4lySerh5Sqy/N77zqC+tpIXXt3Lj36zPu2SJEmSpHHJwFii2ruKu6TaYTyk6ZNr+I/vO50M8OPfvMqqtbvTLkmSJEkadwyMJcoO45Gdvmg6779sMQXg73/8Arv2d6RdkiRJkjRN4Md9AAAgAElEQVSuGBhLlJvejMy7Ll7EWUtm0NbZy50/eJ6e3nzaJUmSJEnjhoGxRHXYYRyRbCbDJ969lJlTanh1+wG+/auX0y5JkiRJGjcMjCWopzdPd2+eimyGqkr/io6kvraS33v/GeQqMix/ZguPPr897ZIkSZKkccE0UoI6istRa6tzZDKZlKspD4ubJnPDNacA8I37VrO5uTXliiRJkqTyZ2AsQQfvX3Q56tF4y9lzufiMOXT35rnjB88PBm9JkiRJx8bAWILcIfXYZDIZbnp7YH7jJHbsaecff/YShUIh7bIkSZKksmVgLEEDz2CcZIfxqFVXVvCfPnAmtdUVPBWb+eUTm9IuSZIkSSpbBsYSNNBhrPWRGsdk9vQ6Pn7taQB899/XsmbTvpQrkiRJksqTgbEEuST1+J0XZvH2C06gL1/gb3/4PPvbutMuSZIkSSo7BsYS5KY3o+O6tyzh5PlT2N/azd/98Hn68vm0S5IkSZLKSqKJJIRwA3AbUAncHmO8Y9jxZcBdwGTgQeDWGGNvCGEBcDcwC4jAjTHG1hDC0uL5k4A9wMdijBuS/A5psMM4OnIVWW593xl87utPsHrjPv7tofVc95YlaZclSZIklY3EOowhhHnAF4BLgWXALcXAN9TdwKdjjKcAGeCTxfE7gTtjjKcCTwKfLY7fAXw+xng28B3gi0nVnyY7jKNnWkM1t773dDIZ+OmjG3jm5ea0S5IkSZLKRpJLUq8GHogx7okxtgHfA64fOBhCWAjUxhgfKw59HfhgCKESuLx4/uB48fU1Mcb7QghZYCGwN8H6U9Pe2b9Lqh3G0XHqwmmDncW7fvISO/d1pFyRJEmSVB6STCRzgW1D3m8DLjjC8fnATKAlxtg7bJzictWpwItAHXDF0RQ0Y0b90Zw+phobGwZf9xYfHdg0e/JrxnXsbn736WxqbuPxF7bz/378Iv/nP19GdWVF2mWNGeeRkuT8UpKcX0qS80tJGi/zK8nAmAWGPjU9A+RHcHz4OEM/F2PcB8wNIbwD+FEIYXGMsW8kBe3e3Uo+X3oPcm9sbKC5+cDg+/0tnQD0dPa8ZlzH56ZrTmb9lv2s27KfL33zCT7xnqVkM5m0y0rc8PkljSbnl5Lk/FKSnF9KUinOr2w2c0wNtCSXpG4Gmoa8nwNsHcHxncCUEMJA+6dp4HMhhA+FEDIAMcb7gFpgWiLVp2jgHsZa72EcVXU1lXzqA2dQXVnBYy/u4Fu/XEOhUHr/AUGSJEkqFUkGxvuBq0IIjSGEOuA64L6Bg8XdTTtDCJcUh24C7o0x9gAPAR8ujt8M3Ft8/SfABwBCCFcCu2KMuxL8Dqlwl9TkLJjdwO9fdya5iizLn97Cvz64Lu2SJEmSpJKVWGCMMW4BPgMsB54F7okxrggh/CyEcH7xtBuBL4cQVgP1wFeL45+if1fVF4HL6H80B8DHgD8KITwL/E+GbKIzngx0GCfZYUzEaYum86n3n0E2k+Gnj27gp4++mnZJkiRJUklKNJHEGO8B7hk2du2Q1yt57UY4A+MbOMSGNjHGF+l/TMe41dPbR09vnlxFhspckg3giW3ZyTP5xHtO4+9/9CLf//U6aqtzvPXc+WmXJUmSJJUUE0mJGbocNTMBNmRJ00VL53DTOwIAd/9iDY88v+0In5AkSZImFgNjiTm44U1lypVMDFcsm8eHrjwJgH/86Wqeis0pVyRJkiSVDgNjiXHDm7H3jgsX8J6LF5EvFPi7Hz3P8+t3p12SJEmSVBIMjCVmoMNY54Y3Y+r9ly3m6vPm09tX4Gvff46XN+9LuyRJkiQpdQbGEmOHMR2ZTIbfvvpkLjlzDt29eW7/7ko2bC+th61KkiRJY83AWGLsMKYnm8nwsXeeyvmhkY6uPv76O8+ydVdb2mVJkiRJqTEwlpj2zh7ADmNaKrJZbnnv6Zxx4nRaO3r46+88S/O+jrTLkiRJklJhYCwxdhjTl6vI8p8+cCanzJ/C3gNdfOnbz7D3QFfaZUmSJEljzsBYYjq8h7EkVFdW8AcfPJuFcxpo3tfJX3/nWQ60d6ddliRJkjSmDIwl5uBzGA2MaautzvFHHzqbuTMnsXVXG3/zLyvpKP79SJIkSROBgbHEtBU7jJNqKlOuRAANdVX88YeX0Ti1hg3bD/CV762iq6cv7bIkSZKkMWFgLDE+VqP0TGuo5k9++xym1lexZtM+7vjBc/T0GholSZI0/hkYS4yb3pSmxqm1/Mlvn0N9bSXPr9vDX/2z9zRKkiRp/DMwlpgOH6tRsubOnMSffuQcpjVU88qW/Xzhn55ix572tMuSJEmSEmNgLCGFQsEOY4mbP6ue224+nwWz69m5t4Mv/NNTrNm0L+2yJEmSpEQYGEtIT2+e3r4CuYoslbmKtMvRYUxrqOa/3XguZy2ZQWtHD1/69jOseGlH2mVJkiRJo87AWELsLpaPmqoc//m6M7nynHn09hX4vz98gZ8++iqFQiHt0iRJkqRRY2AsIe6QWl4qsln+w9tO4UNXnkQG+P6v1/GN+yK9ffm0S5MkSZJGhYGxhNhhLD+ZTIZ3XLiA33v/GVTmsjy4citf/d4qOop/l5IkSVI5MzCWEDuM5ev8U2fxpx85h4a6Sp5fv4cv3v0Ue1o60y5LkiRJOi4GxhLS3lV8pIYdxrK0ZN4UPnPTecyeXsfm5jb+4ptPsnHHgbTLkiRJko5ZoskkhHADcBtQCdweY7xj2PFlwF3AZOBB4NYYY28IYQFwNzALiMCNMcbWEMJpwN8Vz+8Afi/G+GyS32EsDXYYaypTrkTHata0Oj5z03l87furWLN5P1/81tP83vtO56wlM9MuTZIkSTpqiXUYQwjzgC8AlwLLgFtCCEuHnXY38OkY4ylABvhkcfxO4M4Y46nAk8Bni+N/D/zvGOMy4DPAN5KqPw0uSR0f6msr+ePfPoeLls6mq7uPr3xvFcuf2ZJ2WZIkSdJRS3JJ6tXAAzHGPTHGNuB7wPUDB0MIC4HaGONjxaGvAx8MIVQClxfPHxwvvr4LuK/4ehWwIMH6x5yb3owflbksn3zPUt598SIKBfinn0f+Zfkr5H3shiRJkspIkslkLrBtyPttwAVHOD4fmAm0xBh7h40TY/z6kPM/D/zb6JacLjuM40smk+G3Lj+Rxik1fPPnkfse38iufR184t1LqaqsSLs8SZIk6YiSTCZZYGg7JQPkR3B8+DhDPxdCyAB/BVwEXHk0Bc2YUX80p4+pxsYG+oqv5zQ20NjYkGo9Gj2/dXXgxAXT+OI3nuDJ2Mye1mf5kxvP44TZY/d37HxSkpxfSpLzS0lyfilJ42V+JRkYNwOXDXk/B9g67HjTIY7vBKaEECpijH3Fc7YChBBywDeBecCVMcb9R1PQ7t2t5POltySwsbGB5uYD7NvfAUBvdw/Nze6uOZ7Mm1bLf7vxXL76vVWs27KfP/ybf+dDbz2JK8+ZRyaTSfR3D8wvKQnOLyXJ+aUkOb+UpFKcX9ls5pgaaEnew3g/cFUIoTGEUAdcx8H7D4kxbgA6QwiXFIduAu6NMfYADwEfLo7fDNxbfP0l+ndIfdvRhsVyMHAPY633MI5L8xvr+dzHL+DiM+bQ3Zvn7l+s4SvfW0VLW3fapUmSJEmHlFhgjDFuoX8n0+XAs8A9McYVIYSfhRDOL552I/DlEMJqoB74anH8U/Tvqvoi/V3K20IIjcCngQA8HkJ4NoQwbh6pAd7DOBHUVuf4xLuXcuv7TqeuOseqtbv58394nFVrd6VdmiRJkvQ6iSaTGOM9wD3Dxq4d8nolr90IZ2B8A3DFIS45rpPUwV1SfQ7jeHfBabM5ad4U7vrJi6zeuI/bv7uKK8+dx4euPIlqN8SRJElSiUhySaqOQqFQGNJhNDBMBNMn1/AnHzmHD165hIpshuVPb+HzX3+CDdtLa727JEmSJi4DY4no7s3Tly9QmctSmTMwThTZTIZ3XriQ224+n6YZdWzb3c5ffPNJ7n18g89slCRJUuoMjCVisLvohjcT0sI5Dfz5x97ElefOoy9f4LvL1/Klf36GPS2daZcmSZKkCczAWCLaO3sAN7yZyKorK7jpbYE/uP4sJtdVsnrjPv78H1aw4qUdaZcmSZKkCcrAWCIObnhjYJzozj5pJp/73Qs5a8kM2rt6+b8/fIF/+MmLdBTniCRJkjRWDIwl4uCGN+6QKpgyqYo/uP4sbnrbKVTlsvzm+e38j39cwSubx93jRyVJklTCDIwlwg6jhstkMlx57nz+/GNvYsHsenbt7+SL33qKb/48cqC9O+3yJEmSNAEYGEvEwQ6jgVGvNXfmJG67+XzeedECMmT492e28N//7jF++eQmevvyaZcnSZKkcczAWCLsMOqN5CqyfPCKk/jc717A6Yum0d7Vyz/f/zL/4x9X8Py63WmXJ0mSpHHKwFgiOuwwagTmzZzEH314Gb9/3VnMmlbLtt3t/M2/rOQr313Jjj3taZcnSZKkccZ0UiLau/ofq1Frh1FHkMlkWHbyTE5fPJ37n9rEj3/zKivX7ub59Xu45vwTePfFi+xUS5IkaVTYYSwR3sOoo1WZy/LOCxfyxVsu4tKzmsjnC9y3YiP/3/97lAdXbiWfL6RdoiRJksqcgbFEeA+jjtWU+mo+fu1p3PbR8zlp3hRa2nv4+r2r+V/feJI1m/alXZ4kSZLKmIGxRLQVO4yTanwOo47N4qbJ/Pf/cC63vHcp0xqq2bDjAH/5raf5q7ufZE9LZ9rlSZIkqQzZzioRbnqj0ZDJZLho6RzOOamRex/fwL2Pb+TBZ7bw2HPbeOdFC3nHhQuorqxIu0xJkiSVCdNJiRhYkuqmNxoN1VUVvP+yE7n0rCZ+9MgGHl65lR8+vJ7lz2zhmvPnc+U5813+LEmSpCNySWoJKBQKbnqjRMycUsuf3fwm/uyGc1g0p4GWtm6+/+t1/Ne/fYTv/3otLW3daZcoSZKkEmY6KQGd3X3kCwWqKrPkKszwGn1hwTQ++9HzeeHVPfzs0Q2s3riPnz66gV88sYnLzmriHRcsYObU2rTLlCRJUokxMJaAto7+ZzDaXVSSMpkMZyyewRmLZ7B2y35++ugGnn1lFw88vYV/f2YrFy6dzbUXLWBeY33apUqSJKlEmFBKwGBgdIdUjZEl86bw+9efxebmVu59bAOPv7iTR1/YzqMvbOeck2dy7ZsXsmTulLTLlCRJUsoMjCWg1Q6jUjK/sZ5Pvud03n/Zidy3YiMPrdzGMy/v4pmXd3Hqgqm8682LWLpoGplMJu1SJUmSlAITSglo6xzoMPrXoXQ0Tq3lprcF3nvJYn75xCaWP7OZ1Rv3sXrjsyyc08C7LlrIuaGRrMFRkiRpQkk0oYQQbgBuAyqB22OMdww7vgy4C5gMPAjcGmPsDSEsAO4GZgERuDHG2Drkc78LXBZj/FiS9Y8V72FUqZgyqYrrr1jCtRctZPkzm/nlE5vYsP0Ad/7b88yeXscVy+by5jPmMLmuKu1SJUmSNAYS25IzhDAP+AJwKbAMuCWEsHTYaXcDn44xngJkgE8Wx+8E7owxngo8CXy2eM2aEMJfArcnVXcaBgKjz2BUqairyfGuNy/i//zexdx4zSnMmFzDjj3tfOeBV/jjr/2GO37wHKvW7iKfL6RdqiRJkhKUZEK5GnggxrgHIITwPeB64PPF9wuB2hjjY8Xzvw58LoRwF3A58P4h478G/qw4ngX+FLgwwdrH1MA9jJMMjCoxVZUVXHXefN6ybC4rX9nNw6u2smrdbp6KzTwVm5nWUM0lZ87h0rPmMsvHckiSJI07SSaUucC2Ie+3ARcc4fh8YCbQEmPsHTZOjPEXwC9CCB87loJmzCjNxwUMdBgbp9fT2NiQcjUaj0ZjXjXNmcI7Lj2R3fs7eODJTfxyxUa27WrjJ49s4CePbODMJTO55sIFXHzWXKorK0ahapUL/7mlJDm/lCTnl5I0XuZXkoExCwxdr5YB8iM4PnycYZ87Zrt3t5bkErqBwFjo66O5+UDK1Wi8aWxsGPV5dcVZTbzlzDms2bSPh1Zt48nVO3lu7S6eW7uLv/3+Ki5aOpvLzm5i4ewGd1gd55KYX9IA55eS5PxSkkpxfmWzmWNqoCUZGDcDlw15PwfYOux40yGO7wSmhBAqYox9xXOGfm7c8bEaKkeZTIawYBphwTRuuPoUVry0g4dWbWX9tgMsf2YLy5/Zwgmz6rnsrCYuOn0O9bU+Z1SSJKncJLbpDXA/cFUIoTGEUAdcB9w3cDDGuAHoDCFcUhy6Cbg3xtgDPAR8uDh+M3BvgnWmzk1vVO7qanJccc48PvvRN/H5j1/ANeefQH1tJZt2tnLP/S/zR197mDv+9TkefX774GNkJEmSVPoSSygxxi0hhM8Ay4Eq4K4Y44oQws+AP48xPgncCPx9CGEy8DTw1eLHPwV8I4RwG7AR+EhSdZaCwecw2mHUODB/Vj0fufpkrr9iCc++souHVm3lhXV7eGpNM0+taaYimyEsmMq5pzRyzsmNTGuoTrtkSZIkHUamUCi9e/oSsAhYX6r3MP5/f/8Y23e388X/eBGzp9WlXY7GmVJYQ7+npZNnXt7F02uaiRv3kR/yz53FTZM595SZnHtKI00zJqVYpY5FKcwvjV/OLyXJ+aUkleL8GnIP42Lg1ZF+zpZWCWjzHkaNc9Mn13DVefO56rz5tHb0sPKV/vD4wvo9rN/WwvptLXz/1+tomlHHuac0cu4pjSya44Y5kiRJaTOhpKxQKBy8h9HAqAmgvraSS85s4pIzm+jq6eOF9Xt4ek0zK1/Zxbbd7fz00Q389NENTGuo5tyTGznnlJmccsJUchVJ3nItSZKkQzGhpKyzu498AaorK/wXYk041ZUVgx3F3r48azbt45k1u3j65Wb2HujiV09v5ldPb6auOsdpC6dx6sJpnLZwGk0z6uw+SpIkjQEDY8o6unqB/l0mpYksV5Fl6aLpLF00nRuuOZlXtx/g6TXNPL2mmW272wc3zQGYMqnqNQGycWptytVLkiSNT6aUlLV1FgOjy1GlQZlMhsVNk1ncNJnr3rKEnfs6WL1hL6s37OWlDXvZ39bNYy/u4LEXdwAwc0rNYHg8dcE0d16VJEkaJaaUlLUPPFLDDqN0WLOm1jJrai2Xnz2XQqHAtt3tvFQMkKs37mXX/k4eXrWNh1dtA2DO9DpOKwbIsGAqDXVVKX8DSZKk8mRKSVl7lx1G6WhkMhnmzpzE3JmTuOq8+eQLBTbtaO0PkBv3EjftY/uedrbvaWf5M1sAmDdzUn/Hcu5kTmyazLzGSd4zLEmSNAKmlJS1d3oPo3Q8spkMC+c0sHBOA++4cAG9fXk2bD/AS8Xlq69s2c+WXW1s2dXGw8/1dyBzFVkWzq5/TYicNa3WjXQkSZKGMaWk7GCHsTLlSqTxIVeRZcm8KSyZN4V3X7yInt4+Nu5sZf3WFtZvO8D6bS1s39PO2q0trN3aAk/1f25STY5FxfsmFzc1cGLTZKbUey+kJEma2AyMKesodhhr7TBKiajMVbBk7hSWzJ0yONbe2cP67QeKIbKFdVtb2N/WzQvr9/DC+j2D502fXM3iOZOZP6ue+Y2TmNdYz6yptWSzdiIlSdLEYEpJmfcwSmOvrqaS0xdN5/RF0wEoFArsPdDVHx63tfBqsRO5p6WLPS0HH+cBUJnL0jSjjnkzB0LkJObNrGf65GqXtEqSpHHHlJIy72GU0pfJZJg+uYbpk2s4L8wCIF8osH13O69ub2FLc/89kFuaW9nd0sXGHa1s3NH6mmvUVlcwd2Z/eJzXOIn5jf0/J7tDqyRJKmOmlJTZYZRKU3bIbqxDtXf2snV3f3gcGiRb2ntYu6WFtVtaXnN+fW0ls6bV9v+ZWsvsaXWD7+trK+1KSpKkkmZKSZnPYZTKS11NjpPmTeGkeVNeM97S1j0YHvt/trFlVyutHT20dvSwbmvL665VW13BrKkHA+SsaQcD5ZRJVYZJSZKUOlNKygY7jAZGqaxNnlTF5ElVnLZw2uBYoVCgpa2bHXs72Lm3g5372tm5t6P4vp2Orj427DjAhh0HXne9qsoss6bWMqO4VHb65GqmT67pf99QzdSGap8lKUmSEmdKSdngPYwuSZXGnUwmw5T6aqbUV3PKCVNfc6xQKNDa0dMfJPd2sGNvOzv3dQy+b+3oYXNzG5ub2w59bWBKfdXgvZfTGwYCZfXgWENdJVm7lJIk6TiYUlJ2cNMbn8MoTSSZTIaGuioa6qpYMmx5K0BbZw/N+zrYvb+LPQc62dPSWdy1tZM9B7rY19rFvtZu9rV2H3K5K0BFNsPkSVVMGfhTX1V8X82UYkd0Sn3/sZoq/+9AkiS9nv+GkKJ8oUBHcUlqbXVFytVIKiWTaiqZNKeSRXMOfby3L8++1q7+EHlgSJgs/tzd0klbZy97D3Sx90DXEX9fdWVFf4isPxgwG+qqqK+t7P9TV0l9TSUNdf3vqyr9Z5YkSROBgTFFnV19FIDa6hwVWe9FkjRyuYosM6fUMnNK7WHP6entY39bN/vbumlp7R58vb+tm/2tXbQUX+9r7aarp69/Sey+jhH9/qpclvq6SqY21FBTmaW+tpKG2iom1eZoqKtiUk2OupocddWV1NbkqKvuf1+Vy7qZjyRJZcTAmKL2rv4dUifVuhxV0uirzFUcMVRC//2UHV197G87GCL3t3ZzoLjDa2tHD63t3bR29NLa0U1rRw/dvfliN/PI3cuhKrKZYpDMDf6sral8zfu6mhy1VTlqqir6/1QPvO7/WV1V4b2ZkiSNEQNjigbuX6w3MEpKUSZTDHE1OZpmTDri+YVCga6ePlrbe6isqWLT1n20dvRwoKOHtiE/27t66ejspb2rl/biz57ePAfaezjQ3nNcNVcPhMliiKwdeF1dQU1lBVWVFVRXVlBVmaV68HVF8XX24PGq4s9c/3nZrEFUkqShDIwpGrh/0Q6jpHKSyWSKQS1HY2MDU2pGfj9jT28f7V19tHf2DIbIgz97Bt93dvfR2VX82d1HZ/fB1109fXR19//ZT/eofrfKXJaqXJbK4p+qXMWwsf4QOvh68Lz+9wOfy1VkyFVkqazIkstlh7zODI5VVvSP5yqyVOb6z3e5riSp1CQaGEMINwC3AZXA7THGO4YdXwbcBUwGHgRujTH2hhAWAHcDs4AI3BhjbA0hTAW+BZwINAMfijFuT/I7JGmgwzjJHVIlTRCVuQqm5Po32DlW+Xx/h7Ozu4+OwVB58GdXdx9dPXm6e4rhsqeP7p78kNevPT70fU9vnp7e/Ch+46NTkc1QUZEhl+0PnRUVB8NnxeDY8OPF19kMFdksFRUZstn+97mB95n+zw2ekx32vvg6myn+zL7+Z/ZwxzOvPd7/s/8/LAyel8mQKY5JkspLYoExhDAP+AJwHtAFPBJCWB5jfHHIaXcDn4gxPhZC+Afgk8DfAncCd8YYvx1C+CzwWeDPgL8AHooxviuEcBPwFeDDSX2HpLUPdhht9ErSSGWzGWqrc9RW55jWUD1q1y0UCvT05ukuhsbu3oMBcmiY7B5+rDdPT29/KO3py9Pbm6e3L09PX4He3uFjeXr7CvT09hV/9o/3Fsf78v1/ukkvtCYpk+FgqBwImJn+v9PM0NdkyGYhl6ugkC/0f25I8Oz/2X9O/+f+//buNEaSsgzg+L/3QBa5YQjLIajIwxEFI4eRQ5QjERU0IEQQWCMgETQhojGBBFnjkWjAkPDBBBANQYxA/CBgkIVwCIgajsDKE6OgApvIFZAlHDtTfqjqneremp3d7e6tmdn/70vXe9bbu0966qmjezJJ7dRe53V6k9Xua09f+sd2t7vtvX0n9w9U23T7UBtfbdNfX5unw2Q7a+xrirbV5WoN3bmrTs3t1bjumsqt2n6q8T37au5X1dS2V29OrrdnfO+Jgk59jtr+Vk9Wa+s/v7DG2Ia23rrafnsaypc3xgteeXllzwT1ftPNuca8DRVrtk89fvq51u19T2WQEzbTDZ125mkmWGvrgOeZBhk+yL/ZVm+vGmDPM8soM5VjgLsy82WAiLgJOBlYWpX3ABZl5kNV/+uAyyLiauBI4HO1+nsoE8ZPV20AvwKuioiFmTnYwzAtWX2F0VtSJal1nU6HzapnHdtQFGWyuDp5HJ/oKa/qK49PTPart49PFExM1MZ3yxMTjNeS0vFan+6Yie520VdX1Nr6+vXUF91X+soFRQFFAePV3JI0ly2YP48LTzmAfffYru2lDGyUCeMuwIpaeQVwyDTtuwE7Aq9l5qq++p4x1a2rrwFjwPPrsqAddthyPd/CaB38wcXc+/gKDtlvZ8bGtmp7OZrDjC+NkvGldVHUksf+xLMo6EkuJ9ZW1y0XBcUEk9u18uqxRUFRJbBFrb5oSGSLotavWu/qsRSr+6zRtz5moju2b5563+6/RVXoluvb9Tno2+6umVpfuu+BNeehtv76/8XkHNW4arvbXg2t5u72mJyv3n+qefr3Va+rz1++9K2hPnGt3DNX39g12vv2N+U8U7yftbX3rHGKiqKvYq3j+xr75256X82VDeuauus0o9Zt7ABTV80bvu9pdz3Q+MF2vsXmC9l1563nxN/IUSaM8+j9l+5Azz02U7X311Mb139duH/OtXrppdeZmEFnNbfdfAHf+8ohjI1txQsv/K/t5WiOMr40SsaXhq0DdK/xTsZXB+b7/KOGy88vjVI3vmZSjM2b19mgC2ij/LX4Z4HFtfLO9F4JnKr9v8A2EdH9e7G4Nu65qh8RsQDYCnhp6CuXJEmSJI00YbwTODoixiJiC+Ak4Pfdxsz8F/BmRBxWVZ0B3F49j3gfk19mcyZwe7V9W1Wmar9vtj6/KEmSJEkz3cgSxsx8DrgYuBt4FLghMx+OiNsi4qCq2+nAFRHxFLAlcGVV/zXg3IhYDhxB+dMcUDyxZjIAAAXmSURBVH5b6kcj4smqz/mjWr8kSZIkbeo6xaBPk84OewJPz7RnGLu8h16jZHxplIwvjZLxpVEyvjRKMzG+as8wvhd4Zp3HjWpBkiRJkqTZzYRRkiRJktTIhFGSJEmS1MiEUZIkSZLUyIRRkiRJktRoQdsL2EjmQ/nNQDPVTF6bZj/jS6NkfGmUjC+NkvGlUZpp8VVbz/z1Gbep/KzG4cB9bS9CkiRJklp2BHD/unbeVBLGdwEHAyuA8ZbXIkmSJEkb23xgMfBn4K11HbSpJIySJEmSpPXkl95IkiRJkhqZMEqSJEmSGpkwSpIkSZIamTBKkiRJkhqZMEqSJEmSGpkwSpIkSZIamTBKkiRJkhotaHsBm7KIOA24BFgI/DQzr2p5SZoDImJr4AHgM5n5TEQcA1wOLAJ+nZmXtLpAzVoRcSlwSlW8NTO/bXxpWCJiKXAyUADXZOblxpeGLSJ+AuyYmUsi4kDgamBr4F7gvMxc1eoCNWtFxN3ATsA7VdVXgfczB471vcLYkojYFfg+cDhwIHBuROzX7qo020XEocD9wN5VeRFwLXAisC9wcER8qr0VaraqDtyPAz5M+Zn1kYj4IsaXhiAiPg58EvgQcBDw9Yg4AONLQxQRRwNn1aquBy7IzL2BDnBOKwvTrBcRHcpjrwMy88DMPBB4ljlyrG/C2J5jgLsy8+XMXAncRHlmVRrEOcD5wPNV+RDg75n5dHXW9HrgC20tTrPaCuCbmfl2Zr4D/I3yj6PxpYFl5j3AJ6o42onyDqhtMb40JBGxPeXB+w+q8h7Aosx8qOpyHcaXNlxUr3dExGMRcQFz6FjfhLE9u1AegHWtAHZraS2aIzLz7My8r1ZlnGkoMvPJ7oFVRHyA8tbUCYwvDUlmvhMRlwHLgWX4+aXh+hlwMfBKVTa+NEzbUX5ufR44GjgPeA9zJMZMGNszj/I5ja4O5cGXNEzGmYYqIvYH/gB8C/gnxpeGKDMvBcaA3SmvYBtfGlhEnA38JzOX1ar9+6ihycwHM/PMzHw1M18ErgGWMkdizISxPc8Ci2vlnZm8jVAaFuNMQxMRh1GeQf1OZv4C40tDEhH7VF9AQma+AdwCHIXxpeE4FTguIh6lPIg/ATgb40tDEhGHV8/IdnWAZ5gjMea3pLbnTuC7ETEGrAROAs5td0mag/4ERETsBTwNnEb5JRLSeomI3YHfAqdm5l1VtfGlYXkfcFlEHE55Rv5EylsIf2x8aVCZeWx3OyKWAEdl5pcj4omIOCwz/wicAdze1ho1620LLI2Ij1F+I+pZwJeA6+fCsb5XGFuSmc9R3kt/N/AocENmPtzuqjTXZOabwBLgZsrngp6ifOhaWl8XAZsDl0fEo9WZ+iUYXxqCzLwNuBV4BPgr8EBm3ojxpdE6HbgiIp4CtgSubHk9mqUy83f0foZdW52ImBPH+p2iKKbvJUmSJEna5HiFUZIkSZLUyIRRkiRJktTIhFGSJEmS1MiEUZIkSZLUyIRRkiRJktTI32GUJGlAEXEH5e8E/hK4KDOXt7wkSZKGwoRRkqTBHQuQmce3vRBJkobJ32GUJGkAEfFzyh+YfwLYDziU8kfAfwj8GwhgJfAj4BtV+ebMvLAa/1ngEmAz4A3KK5QPbtx3IUlSMxNGSZIGFBEFMAb8BTiZMmG8Ezg4Mx+JiNuBbYCjgK2B54E9gXcDtwBHZeZLEbF/NW6vzFy5sd+HJEn9vCVVkqTReDozH6m2/wG8mplvAy9GxGvA9sCRwGJgWUR0x00AewGPbeT1SpK0BhNGSZJG462+8jsNfeYDyzLz1G5FROxOeQVSkqTW+bMakiQNbhxYuAHjlgHHRcQ+ABFxPPA4sGiIa5MkaYN5hVGSpMH9BriH8tnFdZaZyyPiXODGiOgAq4ATMvP1EaxRkqT15pfeSJIkSZIaeUuqJEmSJKmRCaMkSZIkqZEJoyRJkiSpkQmjJEmSJKmRCaMkSZIkqZEJoyRJkiSpkQmjJEmSJKmRCaMkSZIkqdH/AUHWe7lj4ukGAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "# plot results\n", "plt.figure(figsize = [15,5])\n", "plt.plot(t,I,linewidth=2,label='I')\n", "plt.xlabel('time')\n", "plt.ylabel('y(t)')\n", "plt.legend()\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Simulation" ] }, { "cell_type": "code", "execution_count": 26, "metadata": {}, "outputs": [], "source": [ "PTC = [{\"Ni\":0*1e-3,\"Ci\": 2*1e-6},\\\n", " {\"Ni\":100*1e-3,\"Ci\": 2*1e-6},\\\n", " {\"Ni\":0*1e-3,\"Ci\": 0*1e-6},\\\n", " {\"Ni\":100*1e-3,\"Ci\": 0*1e-6},\\\n", " {\"Ni\":0*1e-3,\"Ci\": 0*1e-6},\\\n", " {\"Ni\":100*1e-3,\"Ci\": 2*1e-6},\\\n", " {\"Ni\":100*1e-3,\"Ci\": 0*1e-6},\\\n", " {\"Ni\":100*1e-3,\"Ci\": 2*1e-6},\\\n", " {\"Ni\":100*1e-3,\"Ci\": 0*1e-6},\\\n", " {\"Ni\":0*1e-3,\"Ci\": 2*1e-6},\\\n", " {\"Ni\":100*1e-3,\"Ci\": 0*1e-6},\\\n", " {\"Ni\":0*1e-3,\"Ci\": 0*1e-6},\\\n", " {\"Ni\":100*1e-3,\"Ci\": 0*1e-6},\\\n", " {\"Ni\":0*1e-3,\"Ci\": 0*1e-6},\\\n", " {\"Ni\":100*1e-3,\"Ci\": 2*1e-6},\\\n", " {\"Ni\":0*1e-3,\"Ci\": 0*1e-6}]" ] }, { "cell_type": "code", "execution_count": 37, "metadata": {}, "outputs": [], "source": [ "PTCt = [{\"Ni\":0*1e-3,\"Ci\": 2*1e-6},\\\n", " {\"Ni\":100*1e-3,\"Ci\": 2*1e-6},\\\n", " {\"Ni\":0*1e-3,\"Ci\": 0*1e-6}]" ] }, { "cell_type": "code", "execution_count": 38, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "16" ] }, "execution_count": 38, "metadata": {}, "output_type": "execute_result" } ], "source": [ "len(PTC)" ] }, { "cell_type": "code", "execution_count": 45, "metadata": {}, "outputs": [], "source": [ "Tc = 50\n", "TT = [i*Tc for i in range(len(PTC))]" ] }, { "cell_type": "code", "execution_count": 46, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAnUAAAHoCAYAAAAv5xlUAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJzsvXl8VdW58P/dJydECIZqEgYHrNq6FFEJpIgaGcTmqq3XtpZKkEHwIEGR6/Veh7fawbbW29v+eKuvolWQWNNQp9veeh0uKjKEYgXa4oBsxSqDoBkFkpiQc87+/XGG7HMSBjH7rMN+1vfzySdnr7P3Wc/aa6+1nv08az3LchwHg8FgMBgMBsORTUC3AAaDwWAwGAyGL45R6gwGg8FgMBh8gFHqDAaDwWAwGHyAUeoMBoPBYDAYfIBR6gwGg8FgMBh8gFHqDAaDwWAwGHyAUeoMBoPBYDAYfIBR6gwGg8FgMBh8gFHqDAaDwWAwGHyAUeoMBoPBYDAYfEBQtwAZIg/4GrALiGiWxWAwGAwGg+FA5ABDgHVAx6FeJEWp+xqwWrcQBoPBYDAYDJ+DC4HaQz1ZilK3C6C5uZVo1PE0o8LC/jQ2tniaR7Yiuewgu/ym7DLLDrLLL7nsILv8Xpc9ELA45ph8iOsvh4oUpS4CEI06nit1iXykIrnsILv8puxykVx+yWUH2eXPUNk/15Qxs1DCYDAYDAaDwQcYpc5gMBgMBoPBBxilzmAwGAwGg8EHSJlTZzAYDAaDIQuJRMI0N9cTDu/TLcohU1cXIBqN9spvBYN9OOaYYnJyvrhKZpQ6g8FgMBgM2mhurueoo/qRnz8Yy7J0i3NIBIMBwuEvrtQ5jkNr6x6am+spKhryhX/PuF8NBoPBYDBoIxzeR35+wRGj0PUmlmWRn1/Qa1ZKo9QZDAaDwWDQikSFLkFvlt0odQaDwWAwGAw+wCh1BoPBYDAYDHHKykp7THcch/vv/zVTplzJ5MlX8sYbf8+wZAfHLJQwGAwGg8FgOAgrVrzC1q0fUF39FB9//BE33zyf6uqnCAazR5XKHkkMBoPBYDCI59dPbeSN9xs9+e2zTy3kpknnHNa1a9euYeLEcgKBAEOHnsSgQYN56603GDFiZC9LefgY96vBYDAYDAbDQWhoqKewsCh5XFhYRF1dnUaJumMsdQaDwWAwGLKGw7WkeU00Gk1Zqeo4DoFAdq3aNUqdwWAwGAwGQxqLFj1Ebe0qAEKhOQwcOIjGxobk901NjRQVFesSr0eMUmcwGAwGg8GQRihUSShUmTzet6+T5577Exdf/E/s2rWD7du3ccYZwzRK2B2j1BkMBoPBYDAchAkTJrJp01vMmFGBZcHtt/+AvLyjdIuVglko4SM6313DZy/djxMJ6xYl43Q2N/PRvQto2/yOblG08Oaa5/jbiv/SkrfTGWHnHzex78NPteRfX1/HK6+8wN69e7TkX7dtDQ0f/gHHcbTkr5NPmtr48SNr+WCXnnuvm6q/PslzH7ykWwwt7PiwmaWLXqe1pUO3KL1Obe36HtMty2LevJuorn6SpUufZvToMRmW7OAYpc5HtK94hPAH6wlvWatblIxTX1NN65tvsONXv9AtihY2/eV/efevK+js+CzjeXe800DrlibaVm7NeN4AL7zwRz76aDurVr2iJf/tm/9IW/ObdLR8qCV/nTz032+zYXMdP32s50HQz7R1tvH8e6/yvFCl7tnfv8F779RR+9IW3aIYXHjqflVKTQHuBHKBX9u2/UDa9yOARUABsAqotG07rJQaClQDAwEbuNq27Ral1DHA74DjgQ7gOtu2sy+ks2acznbdImScSGuLbhGygmg0mvE8nc7M59kTn33WpjV/J9qpNX8d7GnrnU3Ij0SiyLPM9kR7m7znPpvxzFKnlDoeuBsoA0YA1yml0mcUVgPzbNs+DbCA2fH0hcBC27ZPB9YDP4in3wy8adv2OcBPgfu9kv+IxsmOQTajCN4MOhUNA4259QA4Agf5LIvmkFEC5sEHQOCsg6zGS/frxcBy27abbNtuBZ4Gvpv4Uil1EtDXtu3X4klVwCSlVC4wNn5+Mj3+OQc4Ov45H8i8r+lIQGIjM0odgJ55XVly67XPaRPY7izB7U5y2d1IfJnJZrx0vx4H7HId7wJGH+T7E4AiYI9t2+G0dIBfAa8ppXYSc9l+/fMIVFjY//OcftgUFx998JM8YG/8f35+H76kSQZdZf8kLzep4euSQXfeAIXH5tPv6MzK0JjfRGKqtM7yBwKWlvy3xf8XFORxjOb6zzQ5OV12Ad3PfqZp6+waPouK+otV8oLBnC9c93V1AYLBI2+Kf2/KHAgEeqUNeanUBUh9d7WA6CF8n56O67r7gftt275PKXUe8IRSapht24c0oaqxsYVo1Nu3iuLio6mv33vwEz2kteUzOjXIoLPsneGuR0uXDNlQ9w0Ne+nbntnOsd01r0pn+SORqNb89+z5jHBAb/1nGsfVn+p+9jNNe7hr7nJd/R4C1pGnlPQGnfvCX7juo9Eo4fCRNW0oGAz0qszRaGr/FQhYh2WI8vIp3AEMcR0PBnYewvd1wAClVE48fYjruiuARwFs214LfAKc0euSH+Fod0PpQOhbcjp63K/Zce+1P/e689eAJXlSnQvtz55GBBc9K/FSqXsZmKiUKlZK9QOuBF5MfGnb9lagXSl1QTxpGvCCbdudwGrgqnj6dOCF+OeNwLcAlFJfJebCfdfDMhyZiGxlZnCJIbHuY+geWCXOLZLc6pyUz/LqPoEfh5uystIe08PhMOXl47jmmilMmzaZa66ZQiQSwXEc7r//10yZciVTp07ijTf0BeXwzP1q2/ZHSqk7gFeBPsAi27ZfV0o9D/zQtu31wNXAI0qpAuCvwH3xy68HHlNK3UlsykpFPH0G8Bul1O3EQprMsG17t1dlOGIRuPrVWAxiSF4ooV2h9ePodhCyxEirnajAuk+g+2Uqk2zZ8h7Dh5/FggX3p7hfX331ZbZu/YDq6qfYsWM7t956E9XVTxEMZn7TLk9ztG27BqhJS7vM9XkjqYsnEulbgfE9pL8HXNTrgvoOOY0siRldANlKnf7BRXf+mScgut05rk/y6j6JB0Vve2EBke1v9P4PAzknnk2/S28+rGs3b36bTz9t5tprpxEM5lBZeSMlJaNYu3YNEyeWEwgEGDr0JAYNGsxbb73BiBEje1n6g2P2fvUjEvsX0YOLC8Fz6vQjr+GZqo/hCPSOJJCl0FpceOF4pk2bydat/+Bf//VGfvvbJ2hoqKewsCh5VmFhEXV1dVokNEqdHxHYwVhCV56lo99apQ/dZdedvw6khvFIR7b7tfd/83AtaV7zrW9dmfys1OkMG3Ymb775d6LRaEpbcByHgKYpQUap8yMClTrRoe1TkGup069U6c4/82RJ1WvB/bjJslalor/deceiRQ9RW7sKgFBoDi0tLZx11jkcf3wsdK7jOASDQQYOHERjY0PyuqamRoqKirXIbMwbBn8geXRxIXtOnXYJdAuQcYylLoafFZuD4uOih0KVVFXVUFVVQ1nZOLZseY+lS6sB2Lr1Q957713OOaeEMWMuYNmyF4lEIuzYsZ3t27dxxhnpu6JmBmOp8yMiOxgzuICeuT3Zc+fN6tdMkz11r4Ou+o4i0DsSR5JCO3NmiHvu+QnTpn2PQCDAnXfeRb9++UyYMJFNm95ixoxYoI7bb/8BeXlHaZHRKHV+RKD71YQ0iSHbUqd7cNGdf+YxlroY+p89ffix6LW163tMz8/vz89+9p9A6o4SlmUxb95NzJt3U8Zk3B/G/epDRHYwZnCJIXj1q+7nXuK8KsnvUib4cAzd7c6QilHqfInARpYlioVuBNZ89iDw5htLXQyj2BiyBaPU+RGBHYwJaRJHh+vdDOxxJLY73RJkB1GBU14SOFF5z302Y0ZCPyKxgxE8uritBFo6WM23PnusRfIGt+y595nHMTtKABKf+uzGKHUGfyB5co8LyQsldCPRBWeaXQwTfNiQLRilzo8YS50wXBYDLXWv995nj7XIjG6icAcfltjnJjBaXVZhlDo/IrCRmTl1MXRYi3TrVFmj1Ilsd1ly7zUj2v3qw6KXlZUe8Pv3399CRcV3U9KWLq1mypQrmTz5O6xcuTyZvmzZi0ydOonJk7/NM8886Ym8bkycOj/ix1Z2MMzgEkOg+zV7FAt57U7X/pbZhmz3q6yyv/DC//Cb3zxAMNilPr3zztssW/Y8S5bU0NbWypw5MykpGUVHRwePPLKQxYsfJze3D5WVsxg5spSTTz7FM/mMUudHhDUyQLRS52h3A8m996nIa3eSa94slIjhxXCzcOOjvN24ufd/GDiz8HSuP2fWYV3b0tJCbe1Kfvzju/nZz36UTF+7dg3jxl1EXl4eeXl5lJSMYs2a1QCMHFlKQcEAACZMmMiKFa94qtQZn5UvETi/Q7BS50biQolssdRJs1hA9tx73RhLnQz69+/P3Xf/kkGDBqekNzTUU1hYlDwuLCyivr6ux/S6ujpPZTSWOj8iqJElkD24uBdKyFPqskCAOBLbnW4JsgNH4ot0HC+6nMO1pOnCcZyUthA7DhCNRlLGJsdxPJ+yYJQ6HyJQpzOxFRIIDD6cPQq9vIYXyJp7rxdJ1qru+LfsixY9RG3tKgBCoTmUlY3r8bzi4oE0NDQkj5uaGhk69CQANm78W0p6UVGxhxIb96tPEfjWKHn1a8qcOv92sPsja/QKgfdeMu55dLLdr7ol8I5QqJKqqhqqqmr2q9ABjBlzPitXLqe9vZ3m5mY2bFhHaeloSktHs2HDOpqbm2lvb2fFiuWce+55nspsLHV+xM+tbD9kj7VGLxLdr9lS9xIny2fJrdeOxLpPIPFFMp1hw4ZTXn4ZodB0IpEwoVAlxcUDAZg9+3rmz59DZ2eYyy+/gmHDhnsqi1Hq/IjERmZGF0D/6tfYXJJM10W21L28dmfcrzFEKzY+LHpt7foDfj9kyHH88Y/PEQ539bcVFVOpqJja7dzy8ksoL7+k12XcH4J9Vj5GYgcjeHBxUv2vWiTo6WOmyJqqF9jusubea8Bd3VHBO0qIVmizEKPU+RGJHYzk0cWF9g5WS/6m7nWRLa5v3ch2v+qWwODGKHW+RF4rEz24uHpViXtQZk/dS2x3uiXQieZQQlmC5LJnI0ap8yMSG5kJaQJkQQerxf2aHXWv/d5rwMypiyHZUmfILoxS50cEDi6iQ5q40VH37iw15J8tSp1ES53xfMeQPadOtwQGN2Yk9CMCO5jsGdj1ot1aJLqDl1d4S7BWZ/Z+NWQjRqkz+AOj1AH659TpUCqzRqHXrVBrwMx6iKH9ZUoD2dLsDKkYpc6H6B7YtSC4h3GcLJqwLXlOnUBrTbbce91IdL/6ue7LykoP+P3772+houK7KWmTJl3BNddMSf598snHACxdWs2UKVcyefJ3WLlyuWcyJzDBh/2I7oFdA1bAvJ9ANih1gi11AjG3PoZEhV6q5/2FF/6H3/zmAYLBLvVp9+5Pyc0NUlVVk3LuO++8zbJlz7NkSQ1tba3MmTOTkpJRFBQM8Ew+o9T5EoEdjAEQaqXNFnQr1BqQrFBnlYVcA17W/Ef3LqD1zTc8+e38s87m+H+5+bCubWlpobZ2JT/+8d387Gc/Sqa/884mHMdh7txZdHR0MHXqTC666GLWrl3DuHEXkZeXR15eHiUlo1izZjWXXvrN3ipON4xS50cEdjBm9Wsc3XUv2P0q8WXKzKmLERVY97F2J6vc/fv35+67f8muXTtT0vft28e5557H3Lk30tTUxA03zOaUU06loaGeM844M3leYWER9fV1nspolDo/ontg14FrdNGz/2h2oN1iIHlHCd33XgNS21k6Ii3kHlb94VrSdDF27HjGjh0PxPaFHTt2Aq+//lp8LOo6L3bsrQHCmDf8iMQOxo24wVWzG8idp2BLnZlXJRftL1MayJZ25yWLFj2UXPhQW7tyv+etWbOazZs3uVIcgsEgxcUDaWhoSKY2NTVSVFTkocRGqTP4EYEdbALdFgM9IU0ynuV+kPfcBYxWB0h1v+qWwHtCoUqqqmqoqqqhrGzcfs/7+OOdPProI0SjUZqaGqmtXcX555cxZsz5rFy5nPb2dpqbm9mwYR2lpaM9ldm4X/2IREude9JyNIqVk6NRmMySokfpVmgFW+q033sNZMut14Gj20Kumaxpd1nAFVdcyZYtW5g+/Sqi0Shz597I4MFDGDx4COXllxEKTScSCRMKVVJcPNBTWYxS50cEdjApCC6/9sFFS/bZMrjIe+7MwB7DWOr8RW3t+gN+P2TIcfzxj88RDscMKMFgkNtuu6PHcysqplJRMbXXZdwfxv3qR3QP7DrIJmtVxnFbDDRYaZ39HmSE7FEspD13/h7YPw/aX6a0YCo/GzFKnQ8ROWHbhRMV6H6Oo31w0XDrs0Wp037vNRDIknuvg5T1QQKnvAiu+qzGKHV+RGAHk4Lk8mtXLHTnrxN5ZTcDewyR7lcTpDAr8XROnVJqCnAnkAv82rbtB9K+HwEsAgqAVUClbdthpdRQoBoYCNjA1bZttyil1rtk7gucChxv2/YnXpbjiEP7wK4DV5mjwsqfTZHtzUIJWXSLwZUldZERsqjdaUBSTR9JeGapU0odD9wNlAEjgOuUUsPSTqsG5tm2fRqxZ2R2PH0hsNC27dOB9cAPAGzbLrVte4Rt2yOAvwA/NAqdoRsCO9gEut1AekKaZMvwIve5A9HNTuaUl6xpdwY3XrpfLwaW27bdZNt2K/A08N3El0qpk4C+tm2/Fk+qAiYppXKBsfHzk+nuH1ZKTQTOAX7hofxHLhLdj2khTSSRukZEnqUuWxA5sKcYyGWV313aqMA+1+h02YmX7tfjgF2u413A6IN8fwJQBOyxbTuclu7mLuAO27Yjn0egwsL+n+f0w6a4+OiM5JPO3vj/YE5Amwy68u3IzyMRt7vw2Hz6HCun/B3tXTH5+vXNzbgMe+o/oy3++Zgv9eOoDOefl5eb/Kzj/m9LyNEnqO3510W/fn2Sn4uK+pMbFBQfsqUj+Tk/v4+4ug8Gu2xCX7TsdXWBlN/TzZgxI3nttb92S29ra+Puu+/iww8/AOCaa67l61//JwB+97vH+dOf/kA0GuX6629kwoSJAPzv/77AkiWLCIfDTJ48he9+96oe8wwEemfc9lKpC5D6MmORujZuf9+np+O+Til1JlBk2/b/fF6BGhtbiHo836q4+Gjq6/ce/EQPCXeGtcigs+ytrV0dbEPDXnIjmQ/BqKv8+9rbkp9bW9szLsO+Pe3Jz83NrQQzfOs7O7ve7XS2vY6OTu1tP9O0tu1Lfv6kbi95uXKUusbPWpOf97Zkvt3pJhLpGku/aNmj0Wgy5lu20JM8VVWPUlw8iLvuuoc9ez5l2rTJnHPOSD755GNefPE5Hn30d7S1tTJnzkzOOWckHR0dPPTQAyxe/Di5uX2orJzFOeeM4uSTT+n229FoNOU+BgLWYRmivOx+dwAXuo4HAzvTvh/Sw/d1wAClVE7cEjck7bpvAU94IrFfEOYG6YYw96sbme7XbPEDyW532p+9jNNVXuN+7V2ee+pNtr3f5MlvDz31WL4x6azDurakZCQnnvhlAI499lgKCgpobGxk7do1jBt3EXl5eeTl5VFSMoo1a1YDMHJkKQUFAwCYMGEiK1a80qNS11t4ae98GZiolCpWSvUDrgReTHxp2/ZWoF0pdUE8aRrwgm3bncBqIGGjnA684Prd8+LfG/aLvA4mBYEdbALtA6vkhRK6771mJBdf4nxKiSFNvva1MQwePBiAl19eRmdnJyeffAoNDfUUFhYlzyssLKK+vq7H9Lq6Ok9l9MxSZ9v2R0qpO4BXgT7AItu2X1dKPU9s1ep64GrgEaVUAfBX4L745dcDjyml7iQ2ZaXC9dOnELPyGfaHwN7VyaawHjrRUnanx4+ZImt0OoEDe0rVC2t3juCyg7f28cO1pGWK5ctf5t57f8WvfnUfwWAwHs6n6/vYcYBoNJLy0uk4DgGPlWFPZ7/Ytl0D1KSlXeb6vJHUxROJ9K3A+P38ZnpYFEM6AjuYFKTFqXOh2w2k59HLEq1OolLnQnCzE6nQZ42F3EMWLXqI2tpVAIRCcygrG8fTT/+emprHuffeB/jyl08FoLh4IA0NDcnrmpoaGTr0JAA2bvxbSnpRUbGnMmfPchNDryHxrTH1tVmW+9VJNRnoE0RT/lkztui+9xpwUuaVySt/At0vU1rIlnbnIaFQJVVVNVRV1VBWNo5Vq1bwxBNLefDBxXzlK19NnjdmzPmsXLmc9vZ2mpub2bBhHaWloyktHc2GDetobm6mvb2dFSuWc+6553kqc+aXCBq8R2IH48IRbDLQrtBL3lFCoLXGje5HL9M40neUyJp2lzkWL/4N+/a1c9ttN2NZsWf+9tvvZNiw4ZSXX0YoNJ1IJEwoVElx8UAAZs++nvnz59DZGebyy69g2LDhnspolDqD/xDYwSbQH3hZ7kIJiQO7G6/DRWUzMt2vuiXwjtra9T2mP/bY0uTnYDCQEvakomIqFRVTu11TXn4J5eWX9L6Q+8G4X/2IcEud6JAmOgYXZz+fM0a2jC7yBnbZiwVku56z5WXKkIpR6vyIwA7Gje79T7Wiu+41z6nTq1hIb3e6JdCHsdQZsgWj1PkSeR0M2bRYIOPondvj7PcgQ/lnjbVI2nOXijRLXVbtuawFo9VlI0ap8yMiOxgXkt2vmsuue3DTmr/wdifRBZlA4upXY6nLToxS50cEdjBudCsWmSYl8LLIOXVu9Akg0QUn2kDuQmLdG6UuOzFKnR+R3LuC7PKLLLvb/axRDOGIs9QJ38XGLJTIToxSZ/AfkkMr6LbSGverWCQXPyrQUmem1GUnJk6dD9E+sOsg5a1ZYPnj6Am8rHfvVzdmoURmkbyjROpCCXl9jp8tdWVlpT3Gqmtra+Oee37Ctm0fAjBt2kwuvvifAJg06Qry8/OT5/7iFwsYNGgwS5dW8+yzfyAadZg7dx7jxl3kqexGqfMjwjrXbkheKKFbsTCWOrFILr5ES51bpYttYO9fJS9BdXUVgwYN5qc//Q/27PmUadMmM3JkKTk5OeTmBqmqStnqnnfeeZtly55nyZIa2tpamTNnJiUloygoGOCZjEap8yMCe9es2v8002TTvrfadUqzUCKjZE04Gb1ILrsXrPrDQ+z6YJMnvz3k5GGM/XblYV1bUjKSE0/8MgDHHnssBQUFNDY20tjYgOM4zJ07i46ODqZOnclFF13M2rVrGDfuIvLy8sjLy6OkZBRr1qzm0ku/2YslSsUodX5E98CuGd1hPXSie99bLXHyskaxkD2wS3O/pi7QkdvnQGyLuJwc/1vqvva1McnPL7+8jM7OTk4++RR27drJueeex9y5N9LU1MQNN8zmlFNOpaGhnjPOODN5TWFhEfX1dZ7KaJQ6PyKuc01DcAerfXARbKmT3u4Ev0sJVGhT6e3iH64lLVMsX/4y9977K371q/sIBoOMHTuesWPHAzBkyHGMHTuB119/Le6W7rouduzt+lSj1PkSwb0r+q1VmUb7hG3tceqyxWIi67mD1BJLU2xSyi6wz5Wwo8aiRQ9RW7sKgFBoDmVl43j66d9TU/M49977AF/+8qkArFmzmsLCQk4/fVj8SodgMEhx8UAaGhqSv9fU1MjQoSd5KrMJaeJHfNrADhnBJgPtnav2hRI685bd7iSXX3LZwb8v0qFQJVVVNVRV1VBWNo5Vq1bwxBNLefDBxXzlK19Nnvfxxzt59NFHiEajNDU1Ulu7ivPPL2PMmPNZuXI57e3tNDc3s2HDOkpLR3sqs7HU+RGJHUw2LRbINO5wLroVWu1z6oylLqO4ihz16cC+P9yKnMRtwlIXyegTI5MsXvwb9u1r57bbbsayYuW+/fY7ueKKK9myZQvTp19FNBpl7twbGTx4CIMHD6G8/DJCoelEImFCoUqKiwd6KqNR6vyIxA7GhXbFRiO6LQa6O3e9c+rkPncgOua3yJAmbnT3O71NTzHqAB57bGnyczAYIBzuavO33XZHj9dUVEylomJq7wp4AIz71eAT3CYDf3Uwnwf9CyV0K5UmpEkmkRx82I32dqeBlLoX3OdmG0ap8xMCgj8eChI72ATarZTaF0qY1a+68Ou8qkNBskILiJx5kK0Ypc5XmOoEZC+U0N27CrbUSR/ZBOt0OAJXv7oRr9RmEUYL8BMBwZa61NgK2sTQgaN7oYT2kCau7I1Sl1FkhzQxrucEkq202YZR6nyFqU7IAhekRrRPWBZsqZPs9gfZA7vI1a8utPc7hiRGC/ATAcHV6e5UxCl1uoPvuvPXkL0LE9Ikw6SsTxJY/jgilRqBIU2OBARrAT7ELJQAhHawcbRbKbXHqTMLJXQh2FBnLHU+e/bLykp7TG9ra+OOO25hxozJTJ9ewbp1f0l+t3RpNVOmXMnkyd9h5crlyfRly15k6tRJTJ78bZ555knPZTdx6nxFl1IX22NOqJKnW7HRiJbO1cypi+Wtu/AakLBV1KEgvu6FaPS//301J5wwlLvv/iU7dmzlhhvm8N///SLvvPM2y5Y9z5IlNbS1tTJnzkxKSkbR0dHBI48sZPHix8nN7UNl5SxGjizl5JNP8UxGo9T5FccRa7mTPLdJ+8CqJX8T0iQbkBarTPyOEi56u+pbXv4H4Y/29u6PxgkefzT9Lz48pWrWrOsIh8MA7Ny5k6OPLgBg7do1jBt3EXl5eeTl5VFSMoo1a1YDMHJkKQUFAwCYMGEiK1a8YpQ6w2HgRBHrXRdtqdNcdu3eX52KhdznDmTPqZNcdpBjqQMIBoPcfPM8NmxYxy23fB+AhoZ6zjjjzOQ5hYVF1NfXYVkWhYVFKembNr3trXye/rohw8jd/9QRvFBCe0gTN8JWvzqpE/oymndW4LZWyWp2qeFcJCr0jncW8sO1pGWKBQvup77+Y2bPnsnw4WfHpzt1fR87DhCNRlKmQTmOQ8Dj0GNGqfMrwpQ6N9pdkBrRbanTfef1zqmT2+ZAtrVKd7vTjV+rftGih6itXQVAKDSH/Pz+nHgCMKnPAAAgAElEQVTiSRQVFTFkyHEMH342H3zwPsXFA2loaEhe19TUyNChJwGwcePfUtKLioo9lVmof04Afm1l+0OwpS4F3fWu2Q1j5tRlFrNQIobEskuo+1CokqqqGqqqaigrG8fatbVUV1cBMZfr5s2bOOOMMxkz5nxWrlxOe3s7zc3NbNiwjtLS0ZSWjmbDhnU0NzfT3t7OihXLOffc8zyV2Vjq/ETqciRtYuhGuwtSI9o7Vx2Lb7MlpIl2O6VeBE2riuNyPQuveylz6q65JsQ99/yU6dOvIhgMMn/+vzF48BAGDx5CefllhELTiUTChEKVFBcPBGD27OuZP38OnZ1hLr/8CoYNG+6pjEap8ymOE0Xm2ldEW+q0K7Ta59Rpdj8LDiUkbfWrG7P61V91X1u7vsf0fv3y+elP/wOAYDBAONxV7xUVU6momNrtmvLySygvv8QbQXvAuF99hXFBQhYoNpkmZcKyxL1fdYY0Sc9P1rOXMuvBZwP7wRC/96vcdXlZjVHq/IrkViaxg42j21Kl+97rdz8LfvYEW+q0tzvNaG93hiRGqfMr0hqZ4IUS+icsZ9Per5qVSnFzq9zWKo1iaEb8nDrdDd+QxCh1fkXwm6M496sL7Z2r5L1fQXa7033vM03Kcye33kH/y5yhC6PU+RVxnYzgCR5ZFXxYc/a659QJG90kz6lzI7HsKS9Tks20WYZR6vyEI1ixcSO4g9FuMdDu/tXtfpXb7qStfk1dKCG33kF/uzN04WlIE6XUFOBOIBf4tW3bD6R9PwJYBBQAq4BK27bDSqmhQDUwELCBq23bblFKFQAPAsPiP3Gtbdt/9bIMRyraB3eNmLLrFEBv9lHtlkq5g5swnS4FkduEuRD82GcdnlnqlFLHA3cDZcAI4Dql1LC006qBebZtnwZYwOx4+kJgoW3bpwPrgR/E0xcA223bLgH+DzEFz5AkxReiTwwdpGzEKK2DdbtfdUT/dX3UbSnLdP7dspP17KU2O2F9jgvdz70e3Pv++qv8ZWWlPaa3tbVxxx23MGPGZKZPr2Ddur8kv5s06QquuWZK8u+TTz4GYOnSaqZMuZLJk7/DypXLPZfdS0vdxcBy27abAJRSTwPfBX4SPz4J6Gvb9mvx86uAu5RSi4CxwLdc6SuVUrcDVwInA9i2/aJSaruH8h/hyBpc3Eie36F9cNF+64UplVmEvJW/XUiud5BT/t//vpoTThjK3Xf/kh07tnLDDXP47/9+kd27PyU3N0hVVU3K+e+88zbLlj3PkiU1tLW1MmfOTEpKRlFQMMAzGb1U6o4DdrmOdwGjD/L9CUARsMe27XBa+kCgA7heKXU58Bnwr59HoMLC/p/n9MOmuPjojOSTTotlJbvVYwb0JU+DHNrK3jeX5vjno/JytMmhI99cq8115FBU1D+juxo099/DZ/HPfXIzf+9zc3OSn/Pz8zKafzQaxv1mWXhsP/r01fPs6eCovNyuz0f10dbudLA70C/52bEcUWUHCAa72t3R/Y/6QuWvqwsQDHY5Dpcte54dO7Z9Ifn2xwknDKW8/LKDnueWJ8F111USDocJBgPs3LmTgoICgsEA7777DgBz515LR0cHM2bMZOLEr/OXv/yZCRMmkp/fl/z8vowaVcprr63hssu+2e23A4FArzxDXip1AVJfmy1SzUf7+z49nXh6EBgE7LZt+zyl1NeBPwCnHKpAjY0tnpuJi4uPpr5+r6d57A+3haq5qYWcnMzKobPsn7Xtc33u0CKHrvLvbW5NOa6r200gkLOfs3uf9paO5Od9HeGM34N9+yLJz3v2fJbR/J1oOOW4sXEvwTw5uy+2t3cmP7e06ml3umje09XuItGIqLIDdHZ2tbvdX7DdRaPRlC23vLT8OY6Tktf+2P85AebPv54NG9Zxyy3fJxyO8tlnHYwePYa5c2+kqamJG26YzUknnUJdXR1nnHFm8reOOaaQjz/+uMffjkajKfcwELAOyxDlZe+zA7jQdTwY2Jn2/ZAevq8DBiilcmzbjsTP2Qk0AGGgBsC27ZeUUv2VUgNt267zrhhHKLonzGeYlC5AnPs1tbxONAoZVOrSoh9nLt8e0O0GkuaCdHw8r+rzIDGkiZvenvIyceKlvfp7vc2CBfdTX/8xs2fPZPjwsxk7djxjx44HYMiQ4xg7dgKvv/5afC/orutix94GHfHy118GJiqlipVS/YjNh3sx8aVt21uBdqXUBfGkacALtm13AquBq+Lp0+PpHcBLwGQApdQYoJWYsmcARMdqcyNuoUQqWhUbcXu/piG43UlTbFKDXsutd9D+LucZixY9lFz4UFu7kr/9bQMNDTGVY8iQ4xg+/Gw++OB91qxZzebNm1xXOgSDQYqLBybPB2hqaqSoqMhTmT1T6mzb/gi4A3gV+DtQY9v260qp55VSiaUlVwP/Vym1GegP3BdPv57YatlNxKx9d8bTrwUuVUq9RWzl62TbtmW3pv2gfXDTiPYAvJrRWfe6nzvd+eteqKET7bdeI9K3CfOrQh8KVVJVVUNVVQ1lZeNYu7aW6uoqABoa6tm8eRNnnHEmH3+8k0cffYRoNEpTUyO1tas4//wyxow5n5Url9Pe3k5zczMbNqyjtHT0gTP9gng6+cO27Rri7lJX2mWuzxtJXTyRSN8KjO8hfRfwz70uqB+R9uZoAi8n0Wo1EKbUpbtb9SuVGcbsKAEIrHdIcwzJKP8114S4556fMn36VQSDQebP/zcGDx7CFVdcyZYtW5g+/Sqi0Shz597I4MFDGDx4COXllxEKTScSCRMKVVJcPNBTGeXM6BWBUWxATgezP7RaKnV7P3U/97rz14i8dudy++PE50tlbtV5NuE3nba2dn2P6f365fPTn/4HEFsdm1jwEAwGue22O3q8pqJiKhUVU70RtAfMNmF+RZwLUq5Cm24lyLxi4/T4UQf6LSa6888sKcGHtd97vUjbKszJprmshiRGqfMTKasQZXUwbsycOp0LJTKft7u8ugcX7ZZCjUhrdt3jbslVbORZabMXo9T5FclvTtJGl3SEKXWp2et+7nXnrw/9914vkssvuezZhlHqfIVcF6TsvV9Tybi1KCW0Q2azTifze786Bz72Oe7iynO/ppZXmvtV8lbj2YxR6nyKZDeQ9LdG2QsldAcfltvupA/skute2stMNmOUOr8iTalLMRkILjsmpIlmAfTmn3HMjhIJpFkqU50jssqezRilzk+YhRKAWSghbaFEava66153/voQp9SkFVec+9WF36q+rKz0gN+Hw2FCoRk8//yzybRly15k6tRJTJ78bZ555slk+rp1f2HGjMlMnvxtHn54oWcyJzBKnV/xWys7KHItdVkVAFd3mLiMFz2L7r0GsmjbX+1I2/fXjbTVr1VVi9i2bVvyuL6+jkceWcjChYtYsqSGP/3pD3zwwT/o6Gjnnnt+wj33/H9UVz/F5s2bWLt2jaeymeDDvkLwQgk3wkcXnZZKHUpNakgT3c+93GdPmgsuXYmTbanr3bqve7+G9j1bevU3ExxV8BUGnjrlsK9/882NbNnyLmVlY5Np69e/zsiRpRQUDABgwoSJrFjxCiNGjOTEE4dy3HHHA1Befimvvvoy5513QY+/3RsYS51fEWatcmPcr3LLr91SJvjeS3O/pqP92dOIlKK3trZw330LuPXW1N0jGhrqKSwsSh4XFhZRV1fXY3p9fZ2nMhpLna+Qa6lzJM8n7BZVQ19IE91LIHUPrOJccKJDmqQizVLnZdDvL2JJ85IFC37B9OkzOfbYwpT0aDSaskWc4zgEAlZ86zj3mQ6W5a0tzSh1PkW0tUaYGygdreXXbSjTrVjozl8j0tudZKXWr673RYseorZ2FQDXXnsd69ev4/3332fx4oepq/uY9etfJxgMMnDgIDZu/FvyuqamRoqKiikuHkhDQ2MyvbGxkaKiom759CZGqfMr4joYs1Aieax171dplrp0M6m0Z68Ln47rh0xU9yohjWh/mfKIUKiSUKgyeXzhheOTn3/+87sYMWIk5eWXUl9fx6OPPkxzczN9+/ZlxYrl3Hrr9zn11K+yfftWduzYzpAhx/HSS//LN77xz57KbJQ6PyHZBelGctnxbwe7P7Jq71fdpkqN6L73mSa9vNLK70Zw0QEoLh7I7NnXM3/+HDo7w1x++RUMGzYcgO9//0fcccet7NvXwXnnXcCECRM9lcUodX5FsGIj3Q2k0/Wu+95rH1h1568Rv7rgDhVpc+rc6G73vU1t7fqDnvPDH95FONxV5+Xll1Befkm380pLR/PYY0t7Vb4DYVa/+orscYNlHMl7v3bzAEreJkzjIhFAe6C+DONWon02rh8CaZY63Q9/phE83GQzRqnzK9IUGxcmpInkHSVk568TyQsFQLilTnjdZxNGqfMr0joYyXu/6l4o4eznswa0L5QQZqlzY9yvsure7P2anRilzqeIfnOSXHaMpU6zAHrz14i0oqcXV5z71Y3gomcbRqnzK8LeGt29ijT3a7fBRfScuswKIH1gTzGQS9Pq0hBX/pT5lMLKnsUYpc4ndBvMxCl1LiSXHf2BpzOuWGVRSBPJz568gd3s/ZpAe7szJDFKnV+R3MEIs9Slo72D1Zi9/rk9uvPXh/57rxfdL1M68VtIkyMZo9T5BuGWOtEhTdIWSmS6/N18kPo6eN0LJaQN7JJ3lEh/1KRZKlPWR/ms6GVlpQf8PhwOEwrN4Pnnn02m3XjjHKZO/R7XXDOFa66ZwttvvwXAsmUvMnXqJCZP/jbPPPOkp3KDCT7sX6QpNi6MpU5z+bWu09Bddp+Nbp8D6dYayduESbPSVlUtYtu2bcljx3HYvn0bTz/9LMFgl1pVX1/HI48sZPHix8nN7UNl5SxGjizl5JNP8Uw2o9T5Fd2DW8YxIU2SR7oVC1GWunSEPXtmsnwSceVPCT7cu2V/7N2PsHe39epvJlAD+jHjtOMP+/o339zIli3vUlY2Npm2bdtWAG6+eR67d+/mn//5W1x55VWsX/86I0eWUlAwAIAJEyayYsUrnip1xv3qF9K9r+KUui7EW+qiEc3563OB6lbqJLe7iDBrTfpK56ijt93pRIqVtrW1hfvuW8Ctt96Rkr537x5GjfoaP//5r7j33gf54x+fYd2612hoqKewsCh5XmFhEXV1dZ7KaCx1fkWYYpMymIsre/qx5gC8ktyv3e61jMEtQeqcOlllT0dy+Xu77F/EkuYlCxb8gunTZ3LssYUp6cOHn83w4Wcnj7/5zStYu3YNAwZ8CcuykumO4xAIWHiJUep8Q/rAKvmtUZZSl05Us6Uu0+7XVH3eWOp0ofve60Z0SBOf1v2iRQ9RW7sKgGuvvY7169fx/vvvs3jxw9TVfcz69a8TDAYZNGgInZ37KC0dDcSUt2AwyMCBg9i48W/J32tqaqSoqNhTmY1S51cEvzVKs9Slo9sFqXNamfayCx7YjVInue51S+ANoVAloVBl8vjCC8cnP//853cxYsRIyssvZc2a1Sxe/BAPPvgokUiYF154jltu+T8MGjSYRx99mObmZvr27cuKFcu59dbveyqzUer8il9b2f5IW1/vRKNYASlTRk1Ik66sdT/3whSblB0l9ImRDeh/9jJLVgX91swFF1zIpk1vMWvW1UQiUb7znUlJd+zs2dczf/4cOjvDXH75FQwbNtxTWYxS5xuM+zWFaBTEKHWp6HY/6+zgdQ8u0gZ2N9IsdekLJSKC695v7tfa2vUHPeeHP7yLcLirzmfPnsvs2XO7nVdefgnl5Zf0qnwHQuaoJwB5g4vgILCaVz53685F7SghO+i3WSjRRVSaldaF9LrPJoxS5xfS25Q092s6Ebnl122py7wfzu0G0myllDyw+8xac1C67Sgh1zviN0vdkYxR6vyKMItBOtoVG43oVmxEBx/Wfe81Ik6pS0OytUrwY591GKXON8h2A0m2VHYLgqp9Tp2+/IxSl1lSJssjS7FJb3faX6YyTEooIUH1nu0Ypc6vCFJqesJY6rQKoC1r7QqtYPcryLbWmYUShmzAKHV+RVoHk65I6A7Am0m6GWk1hzTR2MHr301DWLtLQ7ulVCPaX6Y0IlmZzzZMSBO/kNaZ6t7/UzfGUqdTAH0Z6i677vx1E4k65OoWIkOkP+aiLXU+U+bLykp7DGsSDoe57LKJHHfc8VhWbNhdvPhxAoEADzxwL3/+82oCgQC33noHZ589AoClS6t59tk/EI06zJ07j3HjLvJUdk+VOqXUFOBOIBf4tW3bD6R9PwJYBBQAq4BK27bDSqmhQDUwELCBq23bblFKjQP+C9ge/4m/2bY908syHLGI62DSLXWSyq85+HD3mCoZzt+ddYa3KDuEFD/TzUAuqdml4TfF5qC451MKsdRt2fIew4efxYIF9xMMBpJx6l599WW2bv2A6uqn2LFjO7feehPV1U/x3ns2y5Y9z5IlNbS1tTJnzkxKSkZRUDDAMxk9U+qUUscDdwOjgA7gz0qpV23b3uQ6rRoI2bb9mlJqMTAbeBBYCCy0bfv3SqkfAD8AbgNKgV/Ztn2PV3L7BnFKXSqO4JAmurcrynT2qXu/6rZSyn3uQNqEeRN8OEFv1/uvn9rIG+839upvJjj71EJumnTOYV27efPbfPppM9deO41gMIfKyhspKRnF2rVrmDixnEAgwNChJzFo0GDeeusN/vrX9YwbdxF5eXnk5eVRUjKKNWtWc+ml3+zlUnXh5Zy6i4Hltm032bbdCjwNfDfxpVLqJKCvbduvxZOqgElKqVxgbPz8ZHr889eAcqXUG0qpPymlTvRQ/iMb3YObbgR3sNpdz4Isdd3zl/vcgey5VVGdmx5rRs5jb3HhheP5zW+WcOut3+dHP/o+n376KQ0N9RQWFiXPKiwsoq6ursf0+vo6TyX00v16HLDLdbwLGH2Q708AioA9tm2H09IBPgWetG37v5RSlcDvgQs8kP3IR04ri9FtlzQ55U9XZLQvlJCk1HXLT85z1xOyLHWpSFPoU3YT6WVl/nAtaV7zrW9dmfys1OkMG3Ymb775d6LRKJZlJb9zHIdAwMJxHFzJ8WNv16d6qdQFSK13i9Qeb3/fp6eTuM627cpEgm3bDyml/kMpNcC27d2HIlBhYf9Dl/4LUFx8dEbycRPt7KDFdRzMsbTIoSNPgOajUqdnH/Olo8gXUv7IZ/kpx0cdFcyoHA35TXS4jgcc3Tej9z4Y7OokHcehqKh/SgfrJfvao+x0HffJzdHWBnTQp0/qEPKlL+VTfExfTdJkloLO1HLm9c1su9NNTk6qcvJFyl5XF0hpx9lAMBjg4YcfZPXqlUBsb9fW1r2cddY5nHBCl5MwL68PgwYNorm5MVmG5uZGBg0ayKBBg2hqakqmf/ppEyeffHKPZQ0EAr3y/Hip1O0ALnQdD4aU/m8HMKSH7+uAAUqpHNu2I/FzdiqlAsD/Af4jnp4gzCHS2NjiuXuguPho6uv3eppHTzjhjpTjzn2dGZdDV9kB2tv3pRw3NeylLV9G+ZubW1OO21rbMyrHZ62pz97uT9toq8/cGkj3ptoAdXV7CAQyM0CE97WkHHd07NPWBnTQ0dGZctzQsBfCh9wlH9Hs3t2Wctza2iGq7iMub0gkEv1CZY9Go93asW7C4SizZs1h1qw5ybT77/81Gze+wb//++189NE23n3XZvjwc2hpaeO55/7ERReVs2vXTrZt28Zpp51BTk4uv/zlz/ne96bw2WefsW7d68yaNafHskajqfcwELAOyxDlpVL3MvBjpVQx0ApcCVyX+NK27a1KqXal1AW2ba8BpgEv2LbdqZRaDVwF1ADT4+lRpdS3gfeAJ5VS04G/xOfrGbq5wGSHNJE8p1C3G0h3rDi98+rkPncAEcHu14jgPlf3XNZMMXNmiHvu+QnTpn2PQCDAnXfeRb9++UyYMJFNm95ixowKAG6//Qfk5R3FsGHDKS+/jFBoOpFImFCokuLigZ7K6JlSZ9v2R0qpO4BXgT7AItu2X1dKPQ/80Lbt9cDVwCNKqQLgr8B98cuvBx5TSt0JbAMq4ukz4uf/iJhFb7pX8h/xCJvf0T1On6TyZ5NSg/aoHpktf5bd+wzTfXc+WeV3I67u3dvz+azee4pRB5Cf35+f/ew/AVJCmliWxbx5NzFv3k3drqmomEpFxVTvhE3D0zh1tm3XELO2udMuc33eSOriiUT6VmB8D+lvA+f3uqC+ID1glL8a2edGlFKXSlR34OkMP3vpY6lWS6W0l6k0JHU73fZcFmyllazMZxvZNTPR0HsIdgWANEtdKtotBpqz1znAOIIHdpA9uEuOU+c4WdDvGACj1PkW3fOqMk53P5AWMXTQzVKVaUtdN0uZ7lhxGvOX1u7SkKzUietzu0090CSGIQWj1PkFs19PCqItdboHVu0LJXS6X2WPbJLj1Em21IH+lzlDDKPU+RVxHUzawB6R5H5OV2oyXfb0F4oMZ5+G1oUSugufYbq/S8od2MUpNabusxKj1PkV3ZPldWMsdRoF0Ju/1v1fpQ3saUR0P3sZJF2JkxzSBLKg3zEARqnzL+IsdalIdr9GdQ8ugufUyZtXlYo4a5ULyWUHU/5swdOQJgaNSGtgghdKdFuooPmNWbdek9HBpVtWgp47egjrIdhaE9X94GcYP3e5ZWWlPcaqC4fDXHbZRI477ngsKzbMLl78OI7jJNMTLF78OIFAgAceuJc//3k1gUCAW2+9g7PPHuGp7Eap8w3pwXdluwJ81cN8TrRbizL8QtF9Xpdxv+pC0o4S3ZQa3e1OMxIsdVu2vMfw4WexYMH9KcGHN29+J5nu5tVXX2br1g+orn6KHTu2c+utN1Fd/RTBoHeql1Hq/Iq4DkauUptuLcm467nbFnWS3K9ZtPJWB2bRfRJxwYc9XCjR9sICItvf6LXfc5Nz4tn0u/Tmw7p28+a3+fTTZq69dhrBYA6VlTdSUjIqJT0nJ4e5c2Ppa9euYeLEcgKBAEOHnsSgQYN56603GDFiZC+Xqguj1PmFxEBm5cQCD0vtXXNyIBKRFdo+jhUI4ESj+hQLi1hHn/FbH8vQsiwcx9FU/gAx16vMdpcTsIhEHWHu11hZA1aAqBMVG9IkELCIRh3t0z4yg8WFF45n2rSZbN36D/71X2/kt799IiX9H//Ywr//+7/w298+QUNDPYWFRcmrCwuLqKur81RCo9T5jYAFEeRZ6uJKrRUI4EQioix1ibIHrBwiRPW5HwMWRBxtCrVlBXCciBY3UCxvjQq1JhJ3OpBQ6gS44NJJKHUS3I9uEh4CK2BB1OnV8h+uJc1rvvWtK5OflTqdYcPO5M03/56SftppXenRaBTLspLfOY5DIGDhJUap8xtWDhCWp9QlCOQAnSJXv1qBAEQ0ugCtmKlO19gWsxhomttjBeIajqyBPUFioJJkqUsqtFYsiITUOXWBgEUEf9b9okUPUVu7CoBQaA4tLS2cddY5HH/8CUCsrwkGg7z44nM9pg8cOIjGxobk7zU1NVJUVOypzCakid8IxKtU4JsjxN8aASLyOlgrXvfaFNrEC6im585KDK46yh9/G5dmqUsQiJdfoqUuR7hSlzBE+dH9GgpVUlVVQ1VVDWVl49iy5T2WLq0GYOvWD3nvvXc555ySlPRt27rSx4y5gGXLXiQSibBjx3a2b9/GGWcM81RmY6nzHRbJtdaO09XifE5yLMnJiR/L62ATSk3Gy56Yzhmw4sYqPZ17wlqUyZeZpAuKQLzswp67+L3OEWipSyDWUhev6qSVVkDVz5wZ4p57fsK0ad8jEAhw55130a9ffkq6ZVnJ9AkTJrJp01vMmFEBwO23/4C8vKM8ldEodX7BPZBZgdhiCSeCNGNswlqFqG3CYgS0W+riLxCZ1ikT8ymTSq0m9yvIU+riJAZ2STtKdM1lFarUxUl4R/xkqespRh1Afn5/fvaz/wRICWniTndjWRbz5t3EvHk3eSdsGrJGfAlYltABJj6wJyx1gubUOa5FIqCx7BosZW4SE5J1uF+TCqWwOXVJA3lArvs1QEKpk1V29yIZkBGn7kjAKHV+JBBTbESGNUlY6nz01nioWAHNrmfNc+qSlkpjqcs4yTl1gtpd94US8rwDILPusxmj1PmRrpmreuXQQJe1Sl4Hq9NS5c5fV6g2y9JnMbCSXWnvhnY4UhD8LtW1UEKYlTaBH92vRzJGqfMhliVwsUBiIBVspQxos9Ql7n3iZUK3pS6D5U8pa9JUmbn8NZNsdoHEymM5ZU+QLLuk/hZSFkiBPPdztmKUOp+QMpcn+dos0FqVo3lemUa0z6nT5H5NbqZiaVQsLEu0hVzi6lcnuaNE7GVKnFIXJ2AsdVmFUep8hoV7oYTARibSUpe++lPT3q/JCdOZzT6B7gnbCRestMUSYOLUgUSlLqHU6m33hlSMUucX0kOagExLXcJaZUKaZJyuOXWZ7t01K7UJBC+WkGipS2BCmviv7svKSntMdxyHqqpFzJw5he9979u8+OJzye+WLq1mypQrmTz5O6xcuTyZvmzZi0ydOonJk7/NM8886bnsJk6d37CslF0lxJAI6xEPaSKp7En3o6Y5dcmuPDmnLqPZJ9GzUML9MiXP/Zq41wHB86rEhjRJCz4swf26bNkLrFv3Fx5++DHa2vYydepkLrhgLDt2bGPZsudZsqSGtrZW5syZSUnJKDo6OnjkkYUsXvw4ubl9qKycxciRpZx88imeyWiUOj9iSXRBxkla6uSV3UpOVpc1py5BQGv5reSuEqIWKMXJERh82OkWfFiedwC82fd34cZHebtxc6/9npszC0/n+nNmHda1r7zyEhUVU8nNzaWwsIiFCxeRl5fH2rVrGDfuIvLy8sjLy6OkZBRr1qwGYOTIUgoKBgAwYcJEVqx4xSh1hs+HFZA7uCQtdQIV2oBu92NAl/s1nr3OOHXQ5fx6aA4AACAASURBVH6VOKfOhy64QyUnINNSl8ASFHz4o4+28+GHH/Dkk0tpaWnh6qtncOKJQ2loqOeMM85MnldYWER9fR2WZVFYWJSSvmnT257KaJQ6P5Kw1El6c+wKbx47FKXUJXaU6NpNw3GcrjluGcLSPGG6K39dcwo1K9UaEb2jhFn9CvSuQn+4ljSviUQivP/+FhYsuJ/du5uYPXsmp52m4v1t13mx4wDRaCSlH3YcJ3m/vMIodX7B3ZmKDmkSV2wiYc2SZB7LsrAsC8dxcJxoMl5hxtBkqUvf+1Wf+zkxl1Veu0u6XyOSlLpYWROrXyMC6x38baVdtOghamtXARAKzeHYYwsZP34iwWCQQYMGc+aZZ/HuuzbFxQNpaGhIXtfU1MjQoScBsHHj31LSi4qKPZXZrH71G5YlMqyHk7b3K5Lm1LkUei2x6tJCmuifU6dnoUQypIkgi02iqiXOqUvQpdTJqXfovlDCj0pdKFRJVVUNVVU1lJWN44ILLmT58pdwHIfduz9l06a3+OpXT2PMmPNZuXI57e3tNDc3s2HDOkpLR1NaOpoNG9bR3NxMe3s7K1Ys59xzz/NUZmOp8w1uS11CqZP35mjlxB5piduEYcV2lYhGIkSjUTJsp9O+XZCWHSWSWF3THnTtk6aRnBx5O0okSpoj3f2aDPju/7q/6qqrWbjwXqZNuwrHiTJzZihpkSsvv4xQaDqRSJhQqJLi4oEAzJ59PfPnz6GzM8zll1/BsGHDPZXRKHU+ROTcni6TQey/KCtlF13z6jKp1CZiqiQOBVnqUiKayGt3XVNZ45Y6iXPqAlKVurRwNj5S6mpr1/eYHgwGmT//3+KfA4TDXXVeUTGVioqp3a4pL7+E8vJLvBG0B4z71Y+IttQl5tQJLDuWXsUiOacu81mD5oUSFqKDDwd9OLAfKoFYy8PBEajYyYpTdyRglDo/InpHCblKHWiO1ZYlIU10LZSQaKlLEJC4UMJllQwInVcHXdMuJM6nzEaMUuc7XAslJK7GkhinrseFEhms++SOFnrjVemYU5e6z6u8HSVI21FC6sAud/9XY6nLNoxS5xdSBvaueGXSkOx+ha6612Kt0rT3a9dWVfom66e4vgUulAgKjFPnLmlXrDp5/U4gR67rPRsxSp3fsCyZ7tfk3q/yFkq4CWhZKJHIPP5fW/Bhze5PS+C+w/H/iV0VJFrqLMsSGdYkGdLEMpa6bMIodX4kIG9wSZIIaSLRUmdZeuLUJbLXPqdOr8VAu1KpkRyB1hq36z0g2f0qsO6zmQOGNFFK5QHXAd8BFBABNgNPAVW2be/zXELD50fyjhIJpUaiUodrcNHoftU9p07/jhICB/bkQgl5ZQfXQgmBfa7ulylDKvtV6pRS44D/B9QC9wEfAJ3AycAlwHql1E22bS/PhKCGz4HEkCZJ96u8hRJODwslolrcr7otdTosZcJ3lIj/zxEcpw4gR2ysOr1zWb2irKy0x1h1juPw2GOLWblyOR0dHUyfPotLLvkGAJMmXUF+fn7y3F/8YgGDBg1m6dJqnn32D0SjDnPnzmPcuIs8lf1AlrrLgTLbtvekpb8N/I9SqgD4EWCUumzAPbDH5/ZIGlySCF4oYWHpXSST3CYss9lmx96vrrmsAhdK5Ai31kgOaSJp9euyZS+wbt1fePjhx2hr28vUqZO54IKxRKMRcnODVFXVpJz/zjtvs2zZ8yxZUkNbWytz5sykpGQUBQUDPJNxv0qdbdv/DqCUOt627Y/c3ymlhtm2vQn4N88kMxweliXb/Zojd+UvuK1VGkKaaFr9mqCr7GZOXcZI2/tVolJnYYkMaeKkhbPpzbr/6N4FtL75Rq/9npv8s87m+H+5+bCufeWVl6iomEpubi6FhUUsXLiIvLw8/vrX9TiOw9y5s+jo6GDq1JlcdNHFrF27hnHjLiIvL4+8vDxKSkaxZs1qLr30m71cqi4O5H49Nv7xeaXUeLo2AcoF/gs4/WA/rpSaAtwZv+bXtm0/kPb9CGARUACsAipt2w4rpYYC1cBAwAautm27xXXdCcAbwEjbtj88eDEl4N6vSKD7NU7S/SrQUpfY+xX0Bh/W9cau11KH6Dl1ib1fw4KUOqfH4MPy+h1JCyU++mg7H374AU8+uZSWlhauvnoGJ544lH379nHuuecxd+6NNDU1ccMNsznllFNpaKjnjDPOTF5fWFhEfX2dpzIeyP26FPh6/HOjKz0MPH2wH1ZKHQ/cDYwCOoA/K6VejVv4ElQDIdu2X1NKLQZmAw8CC4GFtm3/Xin1A+AHwG3x3w0QUwT7HEL5BOK21AkaXNLm1GkJ6aGNnoIP63S/6rbUZbDsKdMe5FnqEitAE2EtJAzs3bFkr371wFJ3uJY0r4lEIrz//hYWLLif3bubmD17Jqedphg7djxjx44HYMiQ4xg7dgKvv/4ajuMkw3cC8WNvg44cyP36TwBKqUdt2551GL99MbDctu2m+O88DXwX+En8+CSgr23br8XPrwLuUkotAsYC33KlrySu1AG3Ai9zCJZCsZgdJUDkKjxLz44SidyTU8p0r37Vkb9FV4Qoec9eIqSJxDh1ADmWWSjhxzl1ixY9RG3tKgBCoTkce2wh48dPJBgMMmjQYM488yzefddm27atFBYWcvrpw+JXOgSDQYqLB9LQ0JD8vaamRoYOPclTmferMiqlLgY4kEKnlCo/wG8fB+xyHe8CTjiE74uAPbZth9OvU0qNAi4CFhwgX/GI3lFCZ/DdLCAZ2V7njhKOnnltVjKkiu69X+U9e5Ln1AEigw8nsHxc96FQJVVVNVRV1VBWNo4LLriQ5ctfwnEcdu/+lE2b3uKrXz2Njz/eyaOPPkI0GqWpqZHa2lWcf34ZY8acz8qVy2lvb6e5uZkNG9ZRWjraU5kP5H69WCk1D3gYWJZQspRSfYAJwI3EYtYt28/1Aei2MWL0EL5PTweIKqX6EXPLTrJtO6qUOlC5eqSwsP/nvuZwKC4+OiP5uAn3aacVCOQE6Jffl31Aft8gx2RYFh1lB6jvk0MrcPSX+lEHWI6jRRYdebY09AWgT59gsoM9+ug+GZPlk6Ny2Qf0P/ooPrMAB4qLju4KRuwxCffPgAH9AAgGAxkre9uePXwcz7NvvzxaGqF/fubuvW6CwdhLxJcGxJ5BK2CJKfvRLUcB0PeoXPIiuQAUDMgTU/4EAwbE7kOfPjmHXfa6ugDBYHbthdCTPFdfPZX777+X6dOvIhqNcu2113HKKSczdOiJ/OMfW5g+/Socx+GGG+ZzwgnHc8IJx3PJJd9g9uzpRCIR5syZy5Ahg3vMLxDonX7rQO7X25VS5wI/Bp5QSu0ipnANAlYDP7Vte+0BfnsHcKHreDCwM+37IT18XwcMUErl2LYdiZ+zM/5bg4A/xRW644gt4vi2bdv2IZSVxsYWz98miouPpr5+r6d59ER0b2wdSTQSpa09ZuRsbWkjnEFZdJUdoKMjZh1paesEINoZzrgsusq/e3cbAJ2dYQKBWJP+tLk1Y7K0t8dikLe0dMSsdY5Dfd2eri3bPCYR8La1dV9cns6MlX1fW+zehyNR2ttjcuzd24alqR1kmnBnrN21xe99x76Itj4g0+zZ+xkA7R2dRMKxcaWxqYV6ZJQ/YXpJtLu2tn2HXffRaJRwOHusnLW16/cjTw7z5t3MvHkxpS8cTsgd4JZb7kg5M3H9VVddzVVXXd0tPZ1oNJpy/wIB67AMUQfcUcK27b8opT4DZhELPhwF3rdte/ch/PbLwI+VUsVAK3Alsd0pEr+9VSnVrpS6wLbtNcA04AXbtjuVUquBq4AaYHo8/X+BLyeuV0p9CFxmVr+mkbL3a/Y0Eu+RvFCiC0vLYgHX54AVm1MXdSAncyKA/uDDXe5nOe0uPfiwtpXHGomFNJE3py5R94nVr36cU3ckciiv0s8A84gpWBdxiF11PLbdHcCrwN+BGtu2X1dKPa+UKo2fdjXwf5VSm4H+xHauALgeuE4ptYmYhe7OQyyPYFwNKm6tkajYWMm9X+V0rl1YeneUsEgGPtKxAFZ38GGJO0okkLn3axeJFwqRIU18PKfuSOSAljoA27Z/B/xOKXU6MYvdOqXUn4H/Z9v26we5toaYMuhOu8z1eSPQbdagbdtbgfEH+e0vH0x2qSQGdkkWgySCAy9DV5w6XYtkrIAVG+w0dPC6gw9L3lEiEdJE7upXE9IkKnSLuGzjkCa9xGPDfRU4jZgiWAcsVErd5aFshsPF7P0qa+VvD7HStFjqYgLE/mvo4JPbFekaWCXGqUvsKCHQUucmIND9muxzBW0TdiRwUKVOKfUzYDux+HBPAF+xbfvfgHHE3LKGbCBlHbFga1XCUheJ6LPY6MKCgG6lNrHiNaMdfDbs/dq157JEC3lOwv0oaWA3O0oALktdRFDdZzEHdb8S26rrsrirNIlt261KqQpvxDIcNpYlOviwFbCSKzCJRruCEfuYFH0+YS3StlUW8fz1uV8zay2SvaNEovyJhRKilLo47r1fJcWpSy6SyUm4nuXVfTZyKHPqrjvAd/uLUWfIOO6togS6IJNYWDk5OOEwTjTStResACz3QglNCn1yTl0G+/fEWKJn9WsC144Sggb2BBKDD7tLKnH1awI/LpQoKyultnZ9t/Rf/vLnvP32W0DMdvD++1u4666fM2HCxSxdWs2zz/6BaNRh7tx5jBt3EQDLlr3Ib3+7mHA4zKRJFVx55fc8lf1QLHWGIwpLpvvV/ZaYkwPhcGyrsFx9ImUMtxtIh0KfHtIEtC6UyKT7NaWUCUudoIUSXQq1XEsdFgQQuFAiXtWS5tTdcsv3k59feOFPvPTSS4wfP5F33nmbZcueZ8mSGtraWpkzZyYlJaPo6OjgkUcWsnjx4+Tm9qGychYjR5Zy8smneCajUer8SEDu3B6smLXKQWJIl2wIaaJvoUSX+1ND3q78JbY7s1BCnvs1gReWuueeepNt7zf12u+5GXrqsXxj0llf6Dd27/6Uhx9+iAcfXIxlWaxdu4Zx4y4iLy+PvLw8SkpGsWbNagBGjiyloGAAABMmTGTFilc8Veqya18OQ+8gPKxHYh6dE5FX/q7NtXWFNCGev77Vr9oC4IqcUxcjYFnEd4gTpNh1lTPHLJQQVO/wxBM1fP3r/8TgwbFNsRoa6iksLEp+X1hYRH19XY/pdXV1nspmLHV+xJIY0qTrYzJOn8AAxHrmU/awq4KGW69lTl0P4WQkWepS5pXlWIQjDpGokxzoJWBhJS11otyvcbxQ6r6oJc1LotEozz33J5YseTyZ5jhOsuvrOg4QjUawXF84jvdtwyh1fsHpYaGEwA4GXLtKSFJqITa3R8O8shQCOtyv6SFNtGxnAYJ3lID44B5xxFhsnBRLneCFEjkJ74D/6n3RooeorV0FQCg0h7Kycbz11huceOJQBg4clNzHtbh4IA0NDcnrmpoaGTr0JAA2bvxbSnpRUbGnMhv3q9+wLKHuV/dCCVmWOqcHa1FGFdrUmCqx/1p3lNAVp07gjhIpzU7mYgmLrjl1kva+TfQ7fna/hkKVVFXVUFVVQ1nZOADeeutNzj57RMp5Y8acz8qVy2lvb6e5uZkNG9ZRWjqa0tLRbNiwjubmZtrb21mxYjnnnnuepzIbS51vcLvApC+USLggJSm1iZAmcYuBlrJbXSvhMmipSw9pon9OnaznDuLvkpaw7aLcCq1ZKOFLS11P7Nz5EV/5yldS0oYNG055+WWEQtOJRMKEQpUUFw8EYPbs65k/fw6dnWEuv/wKhg0b7ql8RqnzHTItdSnWqoTFRsxCCXdIk8wvlEgN6xH/LySkSU/BhyW9TKW4IBOWOiEW8i4ss6ME/rLQ9hSj7v9v783DJDnrO8/PG5FXnX1UVV+6D/RKQghhBMggCQFCGGzANmAezpHHAgPGs8/s4GMHWA88Y8947QXGzyLjMQZ5rBVmjfEsNogFCSwJnQgkgdTSq7Nb6ru77ivv3D8iIjMqKzKrqjsjIut9388/3RWZVfG++WbE+43fGfCJT/xh5PH3vOf9vOc97191/Lrrfonrrvulno1tLaz7VUccQ2PKwAtKDQoOm7a5iFY8ZSqxPZ4fyvt/KmFtLSthKi3iDM5+FYhmZwGdNvf14qZ53aWM45plqet3rKjThajsT4PiO8II10z3K6Rf0iTd3q9ihbBLGmFwRwkIuV8N2dzDVkqb/WrOuvc7VtTphsDMkiZhjHO/tki1+DAha1kqN/iwqEvDUmmgpc4mSvjxhOa5X1fHspq17v2KFXXaELqgHLNFXbPfq3GiTrTahKUwdwGpWOrCRrn0kiUEIniYwrTvnUfQVaJqyOYe1fvVJFEX0LTUmRbu0qdYUacdIiTqqukOJUnCiRKGdZRoRPR+TdRSF9X7NaUMyORFXUQ5GYMsdSuEjcGJEk1RZ+CDdDOmrpFO2INlJVbUaYipJT0Ar6RJUHzYEFEXIMKJEin1fk2zTRikXNbE4FJCQoDrmJwoYW5JE7Bxdf2EFXWa0LDu1xZNS51BlkofJ+0kkVQSJVqIFIvACoPr1EHL/WqOqFvdUcJE9yuEXbCmrH3/YuvUaYYX12SgqAu7XzOmxtSF3K9pzT3FjhIijTZpK7ppBO35zPneNaISJQxzvwow2v0Kvgu2qo+l7sorL4+sVfdnf/YnPPbYo4B3v3nmmaf5zGf+hNe97lre9a63MzQ01Hzvn/7p59i5cxdf+9rN/PM//xP1eoOPfvTjvPa1r4917FbU6cKK4JagTp15liowL6auhWha6lLLfm1Wl0/yrFExhQkLCxFKlDDUBZfx69QZkyixQtAGljrz1l4I0SpArLmg/73f+4/N/99667f4/ve/zzXXvIHZ2Rmy2Qw33XTLivc//vhjfO973+GrX72FpaVFfvu3f5OXvezljI5uiW2MVtTpRqhNllGWuhDNmLqqIaJ2RTeNNETdin5J+ANI8PwBok8SJUy67qI6Spgh6pqIUKJEw4x7TntChOP2tqzJnf/0JQ4/t7cnf6ud3edczNW/9pFT+huzszP89//+Jf7yL/8GIQSPP76XRqPBRz/6bymVSrz//b/J619/Lffeezevfe3ryefz5PN5Xvayl3P33Xfx5jf/So9msxor6rTBxtQBfqKEoZY6IVJ2v4qQ+zW5s/ZHSRMMT5QQoTp1psw/HFMXJImYMveVmBZT9/Wv38Ib3/gmdu3aDUC5XOZVr/pFPvrR32Vqaorf+Z0Pce6553HixHEuuujFzd8bGxvn+PFjsY7NijrtEGaKuhW+ELNE3YqKImlkPq/oZpJe8eFwTF1SZUVWhNQZaKlbEfURtAkzZGMPEIiQ+9WctQ/Tyn7tzXV3qpa0OKnX63z729/iq1/9u+axq6++hquvvgaA3bv3cPXVr+OBB+6j0Wg0n3MB/+d481Nt9quOhCwGxtUNClnqTEuU8JqJpBRTFgwg5ezXtCx1AhFKlDDPWiMwt6MEGJz96l/uvXa/9gtf/vKXuP7693L99e/lRz+6A4BHH/0ZZ5xxJjt27Gy+7+677+KJJ8Lu4gaZTIaJiR2cOHGieXRqapLx8fFYx2xFnS6ErSWitcEYZa3zMdb9SljUmFTSJJwokX5JE0zb2H1aJU3MELWRHSUMmXs7urpfb7jhI9x00y3cdNMtXHnlawF49NGfc+mll61435Ejh/jKV/6aer3O1NQkP/rRnbz61VdyxRWv5o47fkCxWGR6epqf/OTHXH75K2Mds3W/6kZg6nVdqNY8UeeatczGiboVHSW8ta6nVaMvcL+mYiFOO1HCvJImYWWT8T/7qmYb+1oIIGOY+7X98jap+PChQwc5//zzVxx7+9vfwdNPP80HP/hu6vU6H/3o77Jr12527drNdde9hRtu+CC1WpUbbvgIExM7Yh2fWbu91rRfZeFWYfnER5MqhrpfwyVNGilZDISh7ldgVdiDCAfT6I4wr/hwuOC7Y6j7NfiGB63CUrnuYiCqRl3AJz7xh6uOZTIZ/uAPPhn5/ve85/285z3v79nY1sK6X7XDu7iEE9SqM+Mm04goPmyKpa4RWact5d6vRnaUELS2OVOETQtTiw97vV9tmzDQz/26GbGiThs6WerMEDYBItz71ZQ6dQECRGCpSrmjRJLZr/1R0sQXcwZmwEKQKGFY79eo4sOmFXz3r/fWdWfI2vcxVtTpygr3qwGELXXGxdS1/ptKSZMQwjXM/doWXGRaV4mwhdw092uYVvareesOobU/BSutcZUaQvRy7lbUaUczx9z71zBLnR/cAxgk6nxE2m3CUi5pkkadutbJgzGYaakjXHzYMPerEOHiw2atezOm7hQTJTKZHIuLc0YKu0ajweLiHJlMrid/zyZK6ELbtSAclwbmxNSFCdyvpJUBmiLptAkLDyCN4sP9UdIEMLqrhGl16horWqSZmSgR0KxTd5Ixddu2TTA9fZyFhZleDitWHMfp2X0mk8mxbdtEb/5WT/6KpX9oPjr5S2ukqDPNUrc6UaJRTzIDM9x7Nk33a6ukSS2ltTeyrIlPxjWzpAkIMoa5X9s5VUud62YYH9/dyyHFzsTECMePz6c9jFVY96s22EQJwEsWME7U+Qg/USTNAsTBHSXlmLrkuzrYRAnTig+HMbWkSSvax5w6df2OFXXaYWhMne392iTxZImokiYpxcakWXwYTEyUaP2/mf1qmKXOy/xNOewhaToVHzYsnrIfsaJOU0TKWZCpIUSzTp2JxYchVKsu6fkL0XS/JhtTFx5CujF1piZKCGFeTF2YcJ06k4L9RfOeY0ua9AtW1OnCqr4tgbAxMFkgsNSZVqfOJ9VkicBSl6C1ptGwiRL9gGnu15UdJRwEggYN6gaufdP9apiVth+xok43Wn1bvH+NsdSF69T5xYdNsdS1CfpUW4Wl6H4VIo3er+1jMMtSF91RwqyNPbBWtTJg9Rd17SvsNhMl9J97v2NFna4YJ+pCGBZTFyBEu/s1eUtly/2a+KkB21EiTTL+Z181xQXXXoC36YI1aO1tokTfEWtJEynle4FPAVngC0qpL7a9fhnwZWAUuBP4iFKqKqU8E7gZ2AEo4H1KqQUp5cX++4eAKeB6pdT+OOew+Wjv/WqICzJ0Lwncr+bF1Hmk0v+1efK0O0oEc7eJEokQ1VHCtGB5/yufERlKlM0Qde3eAefU6tRZekdsljop5WnAHwNXApcBH/ZFWZibgY8rpS7AuzQ+5B+/EbhRKXUh8CDwaf/4F4HPKqVeCnwd+C9xjX/zYUuaAF6wvqGWugCRaleJ5EVdf/R+9TC1Tp0Qwrjer6uzzoP5myNq2ztKmLL2/Uyc7tdrgR8opaaUUovAN4B3Bi9KKc8CBpRS9/mHbgLeJaXMAlf7728e9///RqXUd6WUDnAWMB3j+DcpZpY0aTTMjalrz7ZL3FoVtpKmmP2aSpuw9mn6LjhTLHUrYuoM7f3ajKnzBX3dAEG/StC6NqauX4jT/boHOBz6+TDwyjVePx0YB+aUUtW24/iu2a3AXmAQuCaWkeuAsR0lzO39ilhZXiDx+afc+xXSj6lrWerM2twEkDG092uAa2IB4sCG0KxTZ5ag70fiFHUOKwW9AOrreL39OOHfU0rNAHuklL8EfEtKeY5Sal1X0djY8PpHfwpMTIwkcp4wpeogS0Am43jtS4YKzAPDQ1lGExxPGnMHOJR1KQHbtg0iMi4vABnRSHw8acz/+EgBgEIhy8TECLm81xh6dDSfyHgq+QwVYHR0gMGJEebwLuakPovA4zs+PsLk5CAAuZybyPnnnAGOhc43dyjP8iyMjOTYntK1kCSu3xps+/Yhar6gdtxkPvu0GZrMAzAwkPOuu2wGirBla4GJUb3nXym3ttyJiRFGRwcAKOSzRqx9QD/ONU5RdwC4KvTzLuBQ2+u7I14/BmyRUrq+WNsd/J6U8jeAf1BKNXw37ACwDTixngFNTi7Enp2TVj+42vQSANVanePH5ymWPB08P7tAKaHxpNkLr1LxbjLT00s4vqiplMqJjiet+c/PLwNQLFY4fnyewFAyNTVPZjD+8ZRLnlF9bm6Z5WnvllL3v4dJELifT5xYYGGhDMDSUimR8y/PeddduVzj+PF5ymXvw5+bXaTm9l9fyF5Tq3rznZpaZGHR++yXS5W+7InZaxYWiwAsL/v3mbr3dHF8cp5saSjNocVOWNQdPz7P0rK39gsLyVx3/UDc93vHESdliIozpu424A1Sygkp5SDwDuC7wYt+1mpRSvka/9AHgFuVUhXgLuDd/vEPArf6//8E8GsAUsrXASeUUusSdPrTJlZdM92vQoRi6qqGzd3/13ECYWVG8eEwgeUo6SSRwFpoaqKEF/Vgtvs14699tWFIxQFaZZRcG1PXN8Qm6pRSB4FPAj8EHgZuUUo9IKX8jpTycv9t7wM+L6V8AhgG/sI//jG8bNm9eNa+T/nHrwf+Vynlw8B/IpR4YQkISpoEBWgNucGEEyUyhnWUWJX4nHCduhXpp62YujTaJaWd/Yphoi7cVSFjbO/XoPiw9zBVM+xBGkKZv4atfT8Sa506pdQtwC1tx94S+v8jrEyeCI7vJyIJQim1F69EimUtAkudgckCIhNkvxoi6pr4iRKplTTxe78KPKHZoGU+TIgg8zd5a9HKRAkMEXUBAsj41pqKoZa6IFGiaqKoa7YJM3Pt+wnbUUIXVvV+DdyvpgkbEJksYJClrg3XDdyvKc0/4QzY6N6v6WysTfergRt7JmOWtabdEp1xDMx+9QnCHkxZ+37Gijrd8GMcWrXaDBE2ke5XM26ujTb/a9Ixdatu4ynVqhNC4Lopd5RwDHO/hqbfahNmmLUm6CjhX3dVAx6k2wWtay11fYMVdbpirKVONF3PjWol5bEkTLPudEq9X4NkAf+pnRRu8Glb6loxdYZdd0I03a/VqpkbeytRwgxB2bgvmgAAIABJREFUD60EIcda6voGK+p0pRlTZ9jmwsrer2kE6yfOqkQJ31KXlqBPsQBxejF1HqbG1EHLBVc1bGNvJUr43z3jHqStpa6fsKJOG9rdQIZa6gQIx2l2lTApUUQ0EyXSKWkS5EQ0W4WlsLm7KSaJgLklTQSQdc1yv7aHPbTcryatfVDSxKy+v/2MFXXa4W+tBsfUAa0MWCNcsB16vyYl6trv426alrqES5q0f+8MS5QIT79Vp86sjT14mDKpTt2qvLxA1Bnqeu8nrKjTFWMtdW2JIoYkSwCt3q/+3GspZ78mnSgB4XpZKa27Y6b7VQBu0Pu13qBuQthDGybWqQti6lrFh81b937DijpdaL+JNt2Phok6H5ENRJ1583fTSpTwEQl3lVhZ0iTp7NeVCOF/7wwTdeBlH2eaAfMGWGzavt4mJkoEOCate59jRZ1uBBmIjl+rzRhLXZsbzKAM2Hbp1EqUSGlzabpfk73BryxpknKdOmM29ra4siAD1jAXLJhV0qSdpqXOul9Tx4o6XTE4+xXCMXWmbK4tUi9p4rtA03S/Jm6pE2YnSjRrtTUzYPXf3NsTJVrZr4atPSFLnXW/po4VddrRLBzk/WvIDabd+2xUq7C2yTuJd5Rod/0n634NE46pS6aczaqIce9fU667tp9NtNQFTe3Ncr+2Fx/2H6YMEPP9jhV1mtLqKKG/+zFMcINtirqKAaLOJ5h70h0lQiPw/kk5USL4HNKoUSgMLT7czAA1yFLXTmCpM8n9GlxrpmY+9yNW1OmC6b1fV5U08TdXEyx17SVNko4ra7eSplh8GFJ0wWKe+3V1fpY5om6VldKg7NeOJU0MWPd+x4o63RAr69SZVHzXo72kiQmibiXJu199gurDCWe/Nk/ftFSmkSzR7JEGmCPqAkQzps5ci41Z7teVtDpKNMzo4tPHWFGnCe2XUdBRwpzs15WIjJ/9a6Coc1OuUyfcwP0a/1N71AbSrG6fwgNNq02YmRaLwP1aMcJi0ylRwrx7jhACx7G16voBK+p0xfjsV99iYoCoW+0KSbijRDu++zONRAnv9MnVqmvPgGx1lND/exeFiZa6VkcJ755bMcD9GoUTstZZ0sOKOm0wPKYuIHABZgx0v65yP6bbUcLImDrHTPdrQMYxKKauvfhwYKkzKEkmcLtDOPNc/7XvZ6yo0wzRHlNmiqWuvayHSSVNViVK+DF1KQnalvs1HVEXFCBOx/1qVkeJdvd3JmOOqAsIdI1JbcKiwubcjLXU9QNW1OmKY6j7tRmxbWBJE//f1DtKpJQo0Tx901KX3PxFW6KEcb1fg1ptjnl16gKCRImKQZa6MK7NgO0LrKjThU69X01zv/oYZ6kMkVpHCR+RYJ26qESJNPu/mlbSpB2z6tTZjhJhgkQJK+rSxYo63RBhi4GARj2RLMT0aQtYz/qWShNi6lYlSiRcfLhdV6XQ+1WEgnsSje1ZVaMvSJQwY2NfVastcL8acc/xESsTJapGrH1E1nkm6CphnpW2n7CiTlOEEIZa60yuUxdUd0+ppEnT/5tu9qubdPFlaM3d0I4SAU33a1X/jb098zlIlKgatPbhhynXJkr0BVbU6YxJGbCder8aIOraN5fki++2WUlTbBMGScfURZc0oZFU79mU6dRRwiBLXStRwhz3a9RX2zGwnE0/YkWdxgjXL8BrUlxZ4H02SNQ1CUqapNVRIsBNu6RJkP2apLDwLcTCobXNGyRs/ClnA1FXNWfuAS33q0n3nNZ/A0Fft5a6VLGiTheic8y9f00Qdat6vxok6jrG1CU093YraVDSxCT3a4hWNxcDLDZtP2cy5mS/dur9amKbMGiVNLHu13Sxok43VlSDDERdJZ2xpIJvMcn6VsqKOXNvuoF8QVtLXND6I2iWNEnn5p5mmzAIZ8Aa8EDh06z5bVSbsABv9lnHPEtdaLdpXXcGWmn7CSvqtGH1k3Gz/6lRos4jEHX1qnlzb1nq0pm78G/uSVjqokuaBIkiaYk6c/suZ/0MyIoJG3t74WUnaBNm3j0HWtmvNqYuXayo0xk/ps4E9+uqZIFA0BpRfHh1ooRwHBqNRjouyJRLmiTbUSIqYty3kJvghmsTNlmDO0oYZamLjPaxlrp+wIo6nXHNs9Q1y/Rlg44S5sydFcImBRdsKw0QSD+mLllLXeuzN9JSF9Rqc82x1HWKqasYuO4Qynw2SND3I1bUaUdoc2la6gwQNu3B+kFMnQHu13YrJYDjr30aLtggUcKoOnUhmokSBsTUta+wUe7XAF/YuMJFIKg36tqXNYm651hLXX9gRZ3ONEVdOd1xJIl/g23GE5pkqQvRSpZIYf5NS11aiRKBpS6lNmkGJ0pkDbLUtSOEMC4DNpwokbGJEn2BFXW6EBEwbladuvY2YQbF1EUWAk2wrEmHkibJWOoiLAZugnXqoq47xxz3a/v0zYqpW732xsTVdYupM2Lt+xcr6nQjXNLEtSVNTHC/BqxIFvAtlanE1AUdJVIo/gvh4sNJzt30mDrvHxMtdSusVaZlwK4oPuzXqTNo7fsRK+p0JnC/GiRsAgL3a91U92uz/2saMXX+bSWBjhKRNbdTSZRo0YqpM8MFF6YZU2eAtSbq222MpS4Ca6nrD6yo05im+9WEG0x7VwWjig9HJUqk2CoshUSJFQbqtBMlDIypC8gY2CZMhMxVGce77+gu6qKu7FaihK1TlyZW1OlM4H41yVLXXtLEpLmHSMX96iP6JlEi7TZhem/sYZqJEgZZ6qLI+K5/U8qaRHWUMEnQ9yNW1OlC10QJA4RNe+9XgxIlol2QgahLYu3bzaQtS11Ux4e4SVLURZV2aMbUGeB+bV9fo0qaRHy3s76lTntRZxMl+hYr6nRjhR/KoDp1TQwuaRJa+3Tcr/5n74iWsEsgrq6dVjxhcqIq7ILDREudgcWHW7TWvhVTZ8h9J3TPydg6dX2BFXXaELF5ZkwUdR4mZr+GadapS6ucTUJxdVGWwNTdr35MHQbG1JlU0iTqm93KftXfStuOa+vU9QVW1GlHqLSCY5D7dVWdOpPahEUkSmSSc79GeVhFUNYkIUtduJxLq6SJjalLGqPcrz5h50jGGEtdl44SBgj6fiYT5x+XUr4X+BSQBb6glPpi2+uXAV8GRoE7gY8opapSyjOBm4EdgALep5RakFJeBPyV//5l4KNKqYfjnMOmxkRLXXtHiWqVRqOxYtPXlbAL0G26XxNc+xVR0w5QgxRu8Olb6sxpE9aOSe7XqHjKrGH9X6MSJaylLl1is9RJKU8D/hi4ErgM+LCU8uK2t90MfFwpdQHe9+ND/vEbgRuVUhcCDwKf9o//NfCnSqnLgE8CfxvX+HXAqI4S7V0NhEBkDMmA7dJRIrVWWW5QgDi9RIlkSpp0SZQwYGM3u6NEQDimzpCSJl0SJcxa+/4jTvfrtcAPlFJTSqlF4BvAO4MXpZRnAQNKqfv8QzcB75JSZoGr/fc3j/v//zLwXf//PwPOjHH8mx8TO0qEHh1NyoBtJyhpUk+wpMkKW2hQgNgYS1047MGc7NeAwBAe7iiRRuZz2phW0iSqo0Td1qlLlTjdr3uAw6GfDwOvXOP104FxYE4pVW07jlLqptD7Pwv8z94OeRMT+ehkjvs1srREJgssax9XFzX3ZgZoElbKqJi6xCx1KSdKRM09KD5swMbePn3HEbiOoFZvUKs3yLj6hz2EMaX4cBQZa6nrC+IUdQ4rr3kB1Nfxevtxwr8npRTAnwFXAK/byIDGxoY38vaTZmJiJJHzhFmaG2AZyOUzzfMvzY5yBMg6jcTGlMbcAfY7DlW8NS74Y9hfyFGbh22jueaxuElj/i8M5QEYHMo3zz+yZQiAfN6JfUzFnEsN2Lp1kAH/XMVClhrLbB0pNI/FQblcBjx3ezDPYtG3UtZrsc/drRaYBPKFbPNcTmWYqRcgnxOpXQ9JEVSuGR8fYXjA+9xzWYflUo0tWwcZLGRTHF28DB7OATAcuu62HBoEIDcQ/3WXJln/wSl83QVWShrJ7Tdp04/zjFPUHQCuCv28CzjU9vruiNePAVuklK5Squa/5xCAlDID/A/gNOB1SqnZjQxocnKBeszZeBMTIxw/Ph/rOaKozi4BUC7XmuevLnhPi+Xl5UTGlNbcAWp1T/dPTS6Sdbwx1H032OTRGXJiIPYxpDX/hcUSAMtL5eb5i0Xv85ifW4x9TJWyZxGbmVlioeA/rfvrMT25yEI+viiPQNQ1GjTnGcQRVqvV2Oe+OFcEoFSsNM+15F93ywldd2kS3E4nT8yz7As41/GSZA4fmWN0KJfe4GJm0b/uFhdLzXWu+NfdzHz8112azM8Wm/8P5rm44H0e5VJN67kHxH2/dxxxUoaoOGPqbgPeIKWckFIOAu+gFQ+HUmo/UJRSvsY/9AHgVqVUBbgLeLd//IPArf7//xwv8/W6jQo6ExEZ74baqJVTHkkCRLifnZxvsSlrPv+oWm0JljSJJOGYunByc1DSpF6vU6/Hff4It79jUvbr6vnnst7al6tmxBSKiESJiuYhL1HRPpmMd91VDVn3fiU2UaeUOoiXofpD4GHgFqXUA1LK70gpL/ff9j7g81LKJ4Bh4C/84x/Dy5bdi2ft+5SUcgL4OCCB+6WUD0spbTmTbgQxdbpnf4ZZkSjhi1rdRV0ELVGX4NzDn30QU5dCeQMhRPLJEiKiPqT2tcrChISNv7mbUNaknawbtAkzae09mjF1FfPWvZ+ItU6dUuoW4Ja2Y28J/f8RViZPBMf3A9dE/MlYx7upiQrYznixVo0kN/a0iJp/zhd1midKRJHxBW0qvV8B/Bt8GiVNADKZDLVajVqtSjabbFyXSaIuymKT99e+bODmnvPXvqy5pS6y4LkrEALq9Qb1egPHMStJpl+wHSV0Y0V58yD71QBRFxDuLOCLOu3drwGhuaftfhWB+zUla02Q/ZuGK8ikOnUB4dtO1jD3a9hC3XS/GiDoYeW6CyFaXSUMtNL2C1bUaUNUSQ+DLHURNOvUmSLqQrh+PGWtks7cRSYoaZJ8TBt4ljpIp/iySZa6KHK++1V3S11UKaGc6113ZUPXPoirq1QMEfR9iBV12hFhqTNC1EUlSniiVndLXWSdOn/tqynVqWsmSiRWiHSlq6dlqYtb1EUlSviizoBEiajVzWXMstStTJTw24Rp736NJpO1lrq0saJOG7oXH240TLnIQgHrObMsdSv6MCYaUxcMIPTZN2Pq4v3edWpakLylbnWiRN1Qa00ua0aiRDdLnSnu1/aHqVayhBmCvh+xok5jhHDM6SoRlSgRZL+m5IJMjMjyAn1S0iShjV20xWQH2a/xW+oixhKUNDFhY4/47gWWupKBG3vWkESJzg9TwXWnt6DvZ6yo0x0/tsqYsibhemU2UcKAmLpomm3SEu3/6hEkStCoGWMhX5koYYalLoqcYSVN2h+mAverLWuSHlbU6UKHJ6dmAeJqKcHBpEFEbFPOlDp1UTF1KbhfQwTZr43YY+q6J0rEbamLalovhDAmAzbSBWlKSZOIr17LUqf7PSeaIPvVFiBODyvqdGOVH8owSx0GlzRhtaWumlaSTCbZjhKrEyWC4sPJx9SBeRmw4WQB4zpKhO65plnq2sn4VlprqUsPK+q0IdpiYVSrsDaMiamLwPUtVfVaNRUXYMtSl3aiRDrCwqQM2HZyhnSUiLRSBpY6zS20nWgmSmi+9v2MFXW6Y0hZk0g3WDP7Ve+n5ihhI4SD6wbJEjFvMFHCyo+pI6WOEsmVNInGmGQJmyixglbvV/Put9ASdbakSXpYUac5xhUgXtFRIqhTp3s8occqz3s24f6vK3q/JmOp60TSJU3aP3vT3K/htTelpEkUGSeDQFBt1KgbkCSzOlEicL+aJ+j7BSvqdKGTQcQ1w1IXXdLE31i17/0avfiptgpLqE5dJ5IradIh7MGYRInVZJuJEuZt7EKIZgFi3cuaRGHdr+ljRZ12tAVsGxZTF35ydPJmdJRo0ZYs4K99NYWYwqR7v66yGKRYfBjMs9SFZx9Y6sqGbOyibe2zRiVLtBUfztriw2ljRZ0mRAXtAqE6dboLm4iYumbvV83drx3iW1qWupjXPiqeMahTF/vG3mHuqcfU+V0lGnpv7FFfPVNKmnS65+Ycv/+rxpa6TglKrrXUpY4VdZohRAdLnfaiLiDcJswvaVIyZe4ryWQDS11Sm0vouxdY6hJLlGhvV5S0pW4ljmmWuoiYOlPcr21G4lBZE/3vO51j6qyoSwsr6rShw+aZLXj/VjS3VkXg5L25a2+p60Am67mfaymsvcimXdIkqNOXlqXOf5gyRNSFKeS8jV377NcO3728Xxu0ZEjIS5isv/YV3de+j7GiTnOa2a+VYsojiZmI3b0ZU1fUW9R1soW1LHUpbC6O8EwY9UYiyRKdYuoqsVspOyRKBO5X7Tf2iFptWUNEXUDbly9nsqgLMp/Lhqx9H2JFneaIbFDSRG9h0yScKFEwraRJe9CyN/9qUpa6cEkTIVoZsCnE12SzSVvqVn72jhtY6kzZ2FvzL5gm6trIu951Z2KrsKxNlEgdK+p0oaO5xrvBaO9+jSpp4tepa5RKNOoax3h0KgSaS1jUtSH8zT3eDNhORVADUZdS71vHDFEX9dXLB+5Xza01HRMlmpY6je+5He45Tfer5mvfz1hRpxvtiRJ+TJ327tcmoUQJx2kmS+hfq241bpruV0AElrpEgqajEyXid79G03S/GhJTF77t5H0xXzRkY28vaWJSTN1q70AQU6fxQ3SfY0WdNphe0iSaVlydKaK2RTZwv6bkfhZGuV9X4jhm1YcMk3EFriOo1RtUUyo+nQSdnCMmibp2bKJE+lhRpx2mWuo6BKznzYqrC5OYpa6j6z89URdY6mJ3v3ZwQ5mc/SqEMC9ZIkQQU6ezqOt0yWdtm7DUsaJOd4xLlGgLWA/KmmieAQusmnuaiRLQKmuSVFeJMOHiw52aj/eU9rAHN3C/6ruxd6NgSFxdFK3iwwasfds1b2Pq0seKOl3oZDEwOFECwMn7BYg1ttR1Ei2pljQhHFMX3w2+k15zHCfB/q9R5zc3UQJMKWsSPfl8xoREiejDrZg6nde9v7GiTnNMLmkCLUudCTF17ZXtE7fUtZFkTF17nTpIN67ONPdr++dvUlmTVYkSBlnq2i871xUIAfVag5rG8ZT9jBV1utPsKKG3qOlUXiCIqTOxq0RgqavFbqnr9NiehPu1s2s1mbImaxQf1t1S18la5bvedXa/drJSmpAo0emqE0I0XbC2VVg6WFGnG6t6vwaWOn1vMGHan5pN6SoBrI6p8+vUVVKrU5dkosRqU10rWSIJS12H4sMab+xhVlmrct5nb4Klrv2rZ1RHiQgTue0qkS5W1GlDJ2tJK6au0TDvyckxOPs1OUtdNCLj3dyTqVO3mjQLEJvmfm0nsNTpXavO4OLDXbBxdeliRZ3mCMcBNwc09K5V1ylRxIiYuk6JEr6lrhzv3DsmlyZafDji9L6oi7MAcaepB4kS9brmG3snF2TOnALEq4sPByVNNF77LhnlNgM2XayoMwCRM6VWHatLmhRMEHUBK+eezQ8ACRYfbnf95/zbS0pP7LlcIOrif5hZ5YQSLggHGnUa9XQKICdK2wcw4LtfTRB17Qz43pGizqLOJyI/iVw+EHUGfO/7ECvqdKFbKa6ct7lTXk5kKKnQYf7uwCAA9WV9597poTnri/m4LXWdCHq/xmmp61aDLuu7n+NtFdbBQiwEjus/UGi8uXf69At5X9SV9N3YO84946170ZSKA23kfEFfLpkn6PsBK+p0IyJwVWQ9UdfQWdR1wBkILHX6z7196R03g+O41GtVagnElbV/81qiLv6be3sPSgiLuiTCDlaf33H8eE7dXbCsnv2Ab61ZNsBa0z73ghuIOgO8AxGmusBSVzbQStsPWFGnDZ0tFsK31Ontfo2evzPgzV1nS103a1EmRWtd4H5NK6YuqFMXq6Wui4Vc+LFVDY0tdZ0Y8C11yxpbazqWc3FzCATleoVaXc/5d2vSks0Hljr9BX0/YkWddkRZ6vyYOhMsdati6kwQdQER1io/UaSahqhL0FIXRWCpK5cTsNRFWCyCsiY6u18D2i2lrZg6Ezb2tlhSIShkDEiWoENMnZ8oYUVdOlhRZwJBTF3FBGGzkqalzgD3axRpxtU169Sl5IZJMlEiCmGQ+7WdpvtVY0tdNwIX7LKBcXU531Jns1/TwYo6XehiDzcipq7D/FuWOo1dz13LC/iirhTj/DucPplEic6vpZkoAeBo7n7tlqRS8C11WsfUdZt/MwNW4/tOB1qWOivq0sCKOgNoxtTpLOoC2pMFDEqUiPKFJBpT137+bKtNWKPeLT07HlJPlPBFnQnu13YGDMh+DYjqO9xKltB87aMSlIKYOp0FfR9jRZ0JmCTq2nBMKGnS5bUgpi4V96sQK4RdPHQraeK5XxOJqYugmf2qqajr9r0byOnvfu06/6CsiaaWum5W2mb2q8Zr389YUacbUSVNjIqp69T7tUijrnebtPbK9pCQ+7ULTRdszPE13UuapNOqq5n9qnlMXaSlqpn9aoK1ZvUHELhflzUvaxKdKGGzX9PEijpt6BZTF2S/6n2DiUI4jmFdJVYSiLo0sl8hlCyRQgZsq6RJSpY6g92vhVyrTl03q87mpktMnUm16tpodpSwlrpUyMT5x6WU7wU+BWSBLyilvtj2+mXAl4FR4E7gI0qpqpTyTOBmYAeggPcppRZCv/dbwFVKqevjHP/mJMpS57kgG+WlpAeTHF02DmdgkHqxSH15CXdwMMFBJUS3RAm/VVi5FKOVtluSTi4ZS10UuVwSJU26JUoEHSU03di7aLWM65DLOpQrdYrlWjPGTkeirFWB+1V3S10UQfZryVrqUiE2S52U8jTgj4ErgcuAD0spL257283Ax5VSF+BdGx/yj98I3KiUuhB4EPi0/zcLUsr/CnwhrnFrSX4IgEZpMeWBJECEL8jxhVxt0YD5t5Hzs3/LxSQEfdQDRXqiLpwoUY/d9R6VKKG5qPOJEjUAQwXPUrpUNG9zH/QrDixVNQ95iVj8fMEXdcV0wh5MJ07367XAD5RSU0qpReAbwDuDF6WUZwEDSqn7/EM3Ae+SUmaBq/33N4/7/7/aH/PvxzjuzUk3a0net05pLOq6eXjcIU/U1pc0tlR2IOevfaWUztxFPl5R18215zhOwhmwbecPRJ2mG3unjgoBg/7mvqjp5t7tnjOY8UVdxcR7Tqv3az2FrHfTiVPU7QEOh34+DJy+jtfHgTmlVLX995RS31NK/T6g512yF0Q8OYnAUqez+9UnKmjbGEtdxORzBW/u5SRKukR99wJLXczxNVGJEgB5P1EmjQxYx/VrJGpuqYu86IAhf3PX3VIXlaDUFHWaCvqAqOvOcURI2Om99v1InIEODiujLgRQX8fr7cdp+72TZmxsuBd/Zk0mJkYSOU+Y+SMFikAhn111/ka1wHN47tfx8eGOG2AvSGPuAE8L70szNjZCZnBgxWszY9tYBIbcWuzjS2P+AwOem2t4KL/q/PXiuPdvrRTb2BZdhzqwbdsQ+fGVMYuT2waZYpKBjMtYDOfPZj2x6Dgicn6DgwMsLMwzOOjGMv/GUp5pYGAgt+rvV0oNDj8BNOL77NOkWmvdlqPmt23LAByYxc1ltJx/4Tlv+xweKaya3+6qd91VRFnLuVeKrYe0yOtuKEe5VGVoMM/28aEkh5Yo/bi2cYq6A8BVoZ93AYfaXt8d8foxYIuU0lVK1fz3hH/vpJmcXIjdHDwxMcLx4/OxniOKyrxnDSiWqtHnd3NQK3P88IlmNmyvSWvu0HLDTU4u4CyufDqsOJ7omTkyiRPj+NKa/9KSZ4VaXCyvOv+ibyhYWpiPbWw1f3Ofnl7Ebay0yBWr3s8L08vUYzj//Lxnfa3XG5Hzc/y1P3p0Gtft/eaysOBdd8vFyqrzN/xm7tXyEseOzcX6MJUGgagTEPnZZ3w/0OFj8X330qTou5UXFoqr5ldZ8u5Hs0t6zn16uuX1iVx7P+v98KEZag09S0nFfb93HHFShqg43a+3AW+QUk5IKQeBdwDfDV5USu0HilLK1/iHPgDcqpSqAHcB7/aPfxC4NcZxGkEQV2dEskQbjskxdU33azpzd1JMlADI572YulIKdfqEk0GIDFCnUdczrqwbRidKNGPq9Ha/dqKVLGHe2qdNbKJOKXUQ+CTwQ+Bh4Bal1ANSyu9IKS/33/Y+4PNSyieAYeAv/OMfw8uW3Ytn7ftUXOPUhjVqQYm8p/i1FXVd5h8kSmgfUxdBUNKkUirSiOuJuctXr5Uokc7NPZeLN6ZuLbu/kwni6szb3Id0T5To8tpg1nuYWtQ2pq77N7/gh4RYUZc8sRYPUkrdAtzSduwtof8/Arwy4vf2A9d0+bs34WXFWtaJMZa6yJImgaVO97mvPuQ4Dtn8AJXSMpVSsWm5S+r8SSVKdKIl6tIpAOy4A9Qq836yxJZUxhA3nbzKg76lblHzjb1bokSxWqTeqOMIPev8d1p7a6lLDz2/aZZVCJNq1bXhDumd/bpWaYlcUIC4mPz8hd8yKI2SJtDKfi2V4hJ13c+vc1mTtRpFNC11y7pa6rqU0xEOBbdAg4aWBYjXWvtA1BU1Xft+xoo63ejw6CQKvvu1uBD5us64w16GUm1B97lHr31uwFv70nIKoq7gW+pifmLvXNLEE1Vxx9RFWWsAnIz3QFGv6hzPGT334UHPUreg/cbeYf6+C3ahoufDZDcC96sVdcljRZ02rBFTV/CETWN5LonB9BXuiC/q5vXLQlsPhaaoS17UCr9eVaNUo5FCIdKC3/e3mESdvghcX9TVtBZ10YwMeEkq80vp9N6NnTW+zsM577pbKBso6nxBv7xkRV3SWFGnHR0sdQMH0CIHAAAgAElEQVS+qCtqKmy6JUoElrr5OT2bi6/lghz0NpfiUlyirks3E0e0kiVSKERa8NukFYvpuMBaljodN/bu37sRf2Of13xj71SoZjjrhbwsVHT3EKxmwF/7ouZr349YUWcIojAKQGNZU1EXEJUokc8j8nka1Sr1lDb3JOgYtBxY6pbiXvtOrn/fWldMPlmiJerSstT5mdcaW+o6fe9GBj1L3cJyhbqOD1NrMGKCpa7D4g/4a28tdcljRZ0urFXSxGBLHZjtgs0PeMIitpi6NfbrwAVbj8VSt1aiRNwxdWskSviiTkdL3Vo6LZtxGMi71OoNLWvVrZWgFFjq5jW01K219gPW/ZoaVtTpRoen5qalrqh5TF2H+YddsPrSIVlgMKGYug6fvdO01MW5sUefPJvN4jgO1WqVajV5YWFCTF23Phnax9VBR2vVcM53v2psqeu09oWm+7WsZ8hLH2NFnSE0LXW6u187kNHYUrdmWY+m+zUdi4GIUdSttV8IIeJ1wa5xfp2zX9ezVZsSVxfFSNa77nS01K21+tmsSybrUKs1qKTUTcZUrKjTjk5xTS33q95PTtHzd0c8S2VtTmNLXYfH5sKgt/bFlAR9IOrqMVrqurVVTSYDtsP3LhvE1Olrren24QdxdTpb6jomShgRU9f5JRtXlw5W1OnCWjF1mRzkBqFeo1HS8cmxO5mtWwGozs6kPJLkKQx5gra4kI6gdQZ8S11KN/eBAc9atrycvLXMi6kT1KuLNBrmWSy2Dnsb+8yCfqJurZi6LTnvYWq2rPGDZBcG/bVfXEinm4upWFGnHZ0fnZwhT9g0FvUTNmtZH92mqJtNYjh9RVPULc1Rr8fQ/3UtF6TvgqsvpxMsP+AniiwtxSHq1niYEk4rA1a3IrTrMPhvHfY6eswu6ruxdyo8vSXvXXdzJf1CPtaz9kO+qFvSUND3M1bUGYQY9EXdkn6iLqBTZ4HMFq/vZnVmOsnhJEqnzcV1M+QHhmk0GukUIPary6dlqRscTMBS18UN5WT9eM6Khps73V3fWwJL3bx5G/tQdhBHOCxWl6jU9HRBdrrnAAwOBZY689Y+Tayo04S1XAFghqjrRGbrNgCqMxrOfR0xkgPDnqhdXkjeUukM+jF1MbQMWk98aOB+jcdStzauHzCvm6hbzz0nsNTNaOiCW+ur5wiHLTnPWjdb1m3t12ZoxFv7JQ3Xvp+xok4zuj01O76oqy/qa63qRBBTV9M6pq7z4g8MeaKuuBijqOtUTiew1C1XY2sV1slCC2FLXZzuz87nz+huqevyms6ibj0ELtjZkqZxdV0WvxVTZy11SWJFnUGIIc9apaWlbq2YulHf/To7SyOOuLIUWY9MKgx7m8tSCpY64YhWWZMUGnynGVMH4Goq6taTRL/Vt9ZMz+so6tb+AAJRN1PSLJZ3AzF1i1quff9iRZ12dH50EkPbAagvTCY1mOTpYLFxslmvrEm9rqcLFtZ4avYslcvz6VhpnSE/WWIxeVE3NOSJusXFdLK+Xd8FV9U0C7Kbd2BkMEvGFSwWq5Q0rVfWzUq8Ne89TE6X9LzndLPSDvuCfmHOiroksaLOIJyRMQAa8xqLui5kx8cBqE6aN/+hLd7aL85OpXJ+x39qry8m74opFAZwXZdyuUSlkvz5MzlvY69V9NzYu+EIwfZRr07g5Jy+fZc7MVbwvCNTRfPWfniLt+4LcyXNa6P2F1bU6cI6Lhpn2BM19YUT+l1k65hPZswTNpXJ43GPJlnWMffBEW9zWZzvvahbz3fJ8TPh6j2Pr1lHgpAQDA35hWAXkrfWuTm/nI5uLrh1Mu5v7idm9RJ160kUaYm6dB6m4mI9c89mXQqDWer1ho2rSxAr6nSjmz08PwTZAlSKUNKsZlZAF1dIdswTtZUTJ5IaTaJ0Ky+QiKWuy2fvDPvu14W43K/dvvgw7Pf+XVyMK66tW6KE382kMkejoVc8p0f3zz4QdZOzcXb06E+2D3iibnJZ0+S0br53YNRf+3nNBH0/Y0WdQQghcEZa1jrTCNyvlUnd5r4OS93wVoQQLC/MUqslXwS46X7t8RP7eg3OgahbWOitqFuPlVI4GT9ZokFNo7i69X72Y1u83ru6WerWw1jBi2OeKk7r5R1Z51RGrKhLHCvqtKP7k5MzMgFAfe5YEoPpK7Lj3twrxzVzvwZ0K2fjugyMbAMaLM3FY63r9s1rirqYMuHWMBg0Rd38fDyiao3TkwlcsGW93HCw9mc/4W/sx2b0tNR1LcCbGaDgFijWSizo1lGEtdc+EHVz03qufT9iRZ0urPMpUGzZCUB95kico+lLsjt3AVA5at7cAUa2eaJ2bupo4ud2/Ey4+kKZRi15F+SoX9Jmbi6duLZM3o/nLJqXpLNzu1cn8MhUOsWf42JdVloh2DHoeQiOLenmIVibrds9K+2MFXWJYUWddqxhqdu6G4D6rGbCZh032Oz4OCKToTo9Tb2ojztgvU6dkW2eoJ+fTt5KKzKOZ61rxJEssTajo56lbDal4tPZgi/qSjqJuvV983b5ou7o1DL1mIpPp8laVtqdg97D1NElTT0EXdjqr/3MpBV1SWFFnWE4WzxrlbaWui7+AOE4ZHfsAKB8LHlrVfysEbS8zZv7/FQ6rndn1LfWzSZft2pkZBQhBIuLC9RqMdRLW8MPlcn75XQ0tNStJWoG8hm2DOeo1upMGVjWpGWpM1DUjfmWuqklvWIK+xgr6rRhfRdM01I3c9jIiyy3y5t/+dDBlEfSQ9a5jiPbPUtdz92v6/waOVs8UVfradD0+k7uui7DwyM0Go0eu2DXd/6mpa6ojwtuI3eP3b7F5tCkfnFla7Fz0HuYOrKkz4PkereOwkCWfCFDpVyzZU0Swoo63VjjsdkZGEUMjEJlmca8PhtMQLfq7gC5004HoHTgQBLDSZS1gpa3jnuCdubEwXhKa6xxfnebFzRdm47DWrOWvQi2bfMyEaenk7eWZfLbESJDrTJLvaqZtWqtLx5w+oRXJ/CFY+l09YiVNeZ/2rDnHTm4oJ93ZK2lF0IwNuF1dJnUce37ECvqDMQZOxOA2tTzKY+kN2zE4pg/w5t76QU95r4RCkOjFAZHqJZLqXSWcP2g6dpU7+JrNmJs3rbNs5ZNTSUv6oRwyA74rv9lPTb3jXz2Z+zQT9StpwAvwI7BCbJOlqniNEsVXZJF1r/4Y/7aTx4zz0qbBlbUacfaT82uL+rqJwwUNoGoe/55Dd3Pa6/91onTAJg5nrz72d1aAAH1uRKNam8theswFjVFXTyWurUHkB3w3N/lZX3ccLCemcOZO72SMs8f1UfUBaw1f0c47PGtdQcWDsc/oERZe/XHd3qi7oSGa9+PWFFnIM742QDUjj+b7kBSIDM+jjM0RG1+jqp2RYjXZttOT9ROHt6X+LmF6+BsLUADapPJWyzG/I4iJ04cS0XQ5wb9eM4ljeI518me8SEyruDI1BKLxbi6ivQvZ454YR/75sx7kJ7Y5Ym6Y4f0Kbzdz1hRpwsb2KTcnecDUDv6tB5tizYwdyEEA+d5819+5um4RpQoGxEo43vOAeDE4ed6OYB1vzWzw4uvqR7rlahb/7mHhoYZHByiUqkwM9Ortk0bcP0PnQFAafGFHp1785DNOJy9y2uX9sxBPTb3jTwWnLvlLACend0fz2ASZkNhD+ND5PIu83MlFmIqPm5pYUWdbqzDD+UMb0cMbYfyMvVpjawG6/HBAQPnvwiA5aeeinM0ybOO+Y/tPhuAqSPPU6v22mKy9vkzE4Go63V8zfrWfscOzwV67FjycW3ZwgTCyVErz1LVqF3YOi87zj/NKwD91IF0agXGx9ofQEvU7aOuw4O0z3rW3nEEO/d4gv7IgXSKf5uEFXXasDF3krvnQgBqBx6LYzB9zYD05r70uHlzzw8MsWV8D/ValeMHk3e/Z3b5ou7oQk86S2zUi7pz5x4ADh9OPvtZCIf8sLe5F+c3f+jDRl3Y8kyvAPTefZo0t9/A/McK29mW38piZYkDC4diHFR/ssdf+wO6rH0fY0WdoWROvwSA6gs/T3kkyVM4+xycwUEqR49SPq5PD9x1GkzYffZFABzZ93jiA3CGcl69ukqd2vHk4+r27PFimw4fPkS93kuLyTqtxKPnAVCce6aH506b9c39wjO3kXEF+w7PMb+kT82ytcooBe+5aPsFAOydfDLuIfUdZ5yzDYAXnpvWMEGtv7CiThc2eJ24p18CQlA7/ASNklmp5sJ1GXqxJ2oXfvqTlEeTPLvPuRiAA08/ksoNNnu654opP5+8K2ZkZJQtW7ZSqZQ5ciR5i0lh1I/nnHuKRr2a+PnTJJ9zkWdspQE8/NTmT1La6JXz4nHPQ/DI8Ud7P5g+Z3znMINDORbmShw/YrNg48SKOu1Y31OzMzCKu+ciqNeoPPdgzGOKmZMQJsOXvwKA+Qfu7/VokmeD8x8/7TwKQ6Mszk4ydTT5bLzsWZ4rprJvhkaPeoGuN64L4KyzzgVg374eWMs2+Nln89vJDuymUS+zPLe5YzpPZuUuv9Cr1ffA43qVdVkPF2+X5N0cz88f4NjS5he1G0EIwbnSyz5/+nF9vCP9iBV1BpM9/xcBqDx+R8oj6REb2NmHXvJSnMFBSvv3Udy/L74xJco6Bb3jcKb8BQCe/fm9cQ4oEnd8AGc0R2O5SvVA8gkD55zjWcv27XuGcrlHbsANiMqhbb6V+MRPe3PulNmIoH653EHGddi7b5pjM2Y1ec+5WS6beAkA9xx6IOXR9Ib1uJ4DXvRiL0lJ/fwotR7XqbS0sKJOGzb+3Jw575WQG6R+/FmqRzZxnMdJWOqcXI7RV18JwPT3vtvrESXKeivbhznv0tcAsP/xH7O8eIrCaoOnF0KQu8B7ai8+dvwUXcAb/90tW7aya9ceqtUqTz659xTOfXIMjb0UITIU55+hvLyJrRYnsWzDA1leedEOGsD3H9jspV02/gFcddoVANx96H6WN3G7uJO5ZnfuGWF8xzDF5QrqUfMstUlhRZ1ubOCpWWTy5C65FoDyA9/Y/AGsGzEZANveeB24LvMP3E/xeQ3qR21g+qPbd3LaeS+hVq2w974eidoNnD//ou2InEvt2CLVg/PJnhx48YtfCsBjjz1CqdSLzXX953czgwyNvQyA2UM/6MG502WDlx2/9KozEcAdjxzk2PTmb5slNrD2Z4+eyflbz2Gpusxt+/81vkH1IUIILrvCq9X4k3v2UynXUh6RnlhRZzi5l7wJURihduRJKns3/wazEbJj42y95vXQaHD0b79KvbJJK92fpBa/5DW/jBAOTz9yN8cOJFuIWeRc8i/x4quW7jtA/SRv8Cf7HLJnz+ns3LmbUqnEAw/cc9IPNCdjJQXYsutKhJNjee5JFqc3Z+D8yT4Cnj4xzBUv3kW11uCmW5+g3qO4yqQ5mVELIXj7eW8G4PvP38EL82aVNznvwgnGdw6zMFfivn/d/GV9+hEr6rRjY4/NIj9E/jUfAKB07y3GlTgZ+9VfJzM2Rmn/Po7e9BUatc379LgRiwHA1vE9XPiKa4EG9/zLV5mbStYlkr94AndsgMZihaV/3XdK/WA3ai0SQnDFFVeRyWR47rmnefTRh0/63CeDmx1h62melXzq+X+mtLB5XZEb/d4BvPv15zM6mOWJ52f4v297kvom9hJsdPbnbjmbq077RWqNGn/1s5uYKm7e2m0bve4cR3DNmy/AcQSP/vQQj/7ULFGbBLGKOinle6WUe6WUT0kpfyfi9cuklA9KKZ+UUn5ZSpnxj58ppbxTSvmElPL/lVIO+8e3Sim/LaV83H99V5zjN4Xsea8k+5I3Qb3G8v/3Bco/+y6N+uYVNxvBHRhgz0d/F5EvMH//vRz8b5+jMhlHw/f+5JJXv5ldZ11IaWme27/2eV548qHE3PDCEQy+9ixEIUP18AIL332a2kxycUZbtmzl1a9+LQAPPfRj7rnnTsrl5NoYDY+9nKHtL6VRr3Ds6b9j/sRPNn8IxDoZHcrx0V+9hIwr+OFPD/LFb/6c2QVzWki94/xf4ZzRM5kuzfBnD/5fPDb5RNpDSoyJXSNcdZ2XrHTX957i7tuftq7YHhKbqJNSngb8MXAlcBnwYSnlxW1vuxn4uFLqArwHng/5x28EblRKXQg8CHzaP/6fgbuUUhcBfw38t7jGv+k4xc0gf8W7m8KudN/fs/j1P6T0029RO/4cjVp/19M61Y2wcPbZnP7v/wPu8AhLex9j3yf/gCN/89csPPwQtYXNUFPp5OfvOC6vedtvsefcSyiXlrjnX77K92/5c5586A7mJo/0uEDvatyRPMNvOg9nOEdtcpn5bykWf7iP8r4Z6svrcYef2tqfffZ5XHnl63Ach6effoJvfvNrPPDAPRw+fJBKJd4CuUIItp/5VobGXkajUWX6hW9z+PEbmTt6N+WlI/3fl/kUrzt55jb+3TsvZSDv8tBTJ/iDL93LTbc+zs+emWSpuAlCIU5h/lk3y8de+ltcsPU85srz3PjIV/g/f/JF7jp4H8eWTvR9K7FTffa4+LI9XHXd+QgBP/vxQW7+0v3c+8NnOPT8DJWKFXinQibGv30t8AOl1BSAlPIbwDuBz/o/nwUMKKXu899/E/AZKeWXgauBXw0dvwP4A+CX/dcAvgZ8UUqZVUptgjtAQmzUHt78NYfCL76HzJ6LKN73NRqzRyk/+E3KD34ThIsYGccZ2orID0NuAJHJgZtFOC44LgiHud17aOx5BUKk49XfSHp9OwPnv4iz/tNnOf7/fJ35H9/P3L13M3fv3QC4I6NkxsbIjI7iDAziFPKIbA7hugjXBcdBZLOM/vIbgUKPZnMSnOT0M9k8V779Bp752d08es+tTB99gemjnjvQdbMMbdlOYWgL+YFBMtkCbjaH62ZwHAfhuOyp7MbFPelhu1sLDL/1Aoo/PUz5yUkqz89S8QsTi7yLM5RDDGQQOReRdcB1EI4AR1Ba9vqINmonv8uce+6L2LZtjAceuJujRw/zxBOP8sQTXpzb4OAQQ0PD5PMFcrkcmUwG13VxHBfHcRjkEANA/SRbngnhMHbmWymMnMvMwduoliaZOXQ7cDsIl0xuK252GCcziOPkEY53zQnhAg4IQTY/xtD2l5z0/E+Zk7/suOScMT7zm6/ka7c/xUNPneDORw5z5yOHAdgynGNstMDIQJbBQoZ81iWbccm4AscRuI4gl3V5zSW72DKc79FkToKTvO8MZgf4+GU38MMDP+K7+37As7P7eXbWS9jKuTnGC9sZzY0wmB2g4ObJujkyjosrXFzhIBCct/UcLtz+ol7OZoOc/OJf8gunMbFrhB/d9jTHDs3z8P0HePh+r33f8Gie4dE8hYEsuXyGbNbFzQhc/9oXQlAYyHDRpbvJ5k7+3qMjcYq6PcDh0M+HgVeu8frpwDgwp5Sqth1f8TtKqaqUcg6YAPrCMf+Ff3iEnz2TjuvuivxzvGcIfvz0DFe9+uT/Tuasyxg641JqB35O9bmfUD2iaMweozF3lNpc55irr+7egprLg/rHkz/5SSLqDT7mQOMUr+3M1m3s/vBHGPvVX2f+/ntZ2vsYxf37qM3PUZufo5tz6OHd1zL583QKGQ+6Jyi4njg7WYRwOP+lV3H2xa/i4NM/4+Czj3Li4LMsL8wwN3W0a7zd2NbfYMAZOumkAQAn5zJ4xekULt1J+dlpqofmqZ5YolGqUSt1rmdWZwkEOEunZtnYtm07b3rTWzlx4jj79j3DkSOHmJmZYmlpkaWlzh1Xzhhf4LzdMD07x/bTO75tTYa2vZjBrReyPPsUSzNPUF58gWp5mmppkmqp8z3lO7Wreb5RgGeSL2Rcr3ifeekUrbnjWwf43XdcyqETi9y39yiP75vi+WMLzC6UmV3obi398PDtOD87SC/ypzeK2DkKIwXyTu6k/4bruFx75mu5cs+reOjYz/n55OM8N7ufufI8hxaPcGjxSMfffdu/zuAcKpNGMaqZwg44/ZcpPvsU0O6AWz8794zy6x94GUcOzvHME8c5tH+GqROLLMyVWJjrfMcdztxNzjnC43ee9KlPmbPkpVzxyzekN4AOxCnqHFb6RgRQX8fr7ccJ/V77Y0H73+zK2Njwet96UuRycX6c3XmodDY5qtTHL+fXJ0ZO/Q/uvBJe7tVxq1dKVGeOUVucobY8T720TKNSolGreK7Zep3s5MNQmjr1854EDUfwnSu3cM7E2VzTi7lPjMCLvR6djXqd8uQUpclJKrNz1BYXqZWK1MtlGtUajVqNRr2OejIPyXe9AqBYu4DtO3dz6SuvIJs7dYvF7j2v4/KrXwdAqbjE/PQJlubnKC4vUikXqZbL1GpV6rUa9XqN6lwNd7TArjPHTvncAJy1HfDc6rXFCtX5EtWlCvVSlXq5RqNa9yxz9Trb6g1efSzPGRefzUQP1n5iYoSLLvI6TtRqNebn51lYWGB5eZlSqUS1WqVarVKr1Wg0GtRry9Trh7jksjczPNKD796OVwBet5NatUR5eZpKeY5qZZl6tUi9VqFer9Jo1Gg06uQPDENKXf6crMOWS8Z40Y6Rnn32L73IC5Ou1RtMzixzYnaZucUyi8sViuUalWqNSrVOrd6gXm9w2tP3w8wpn/qkuHZqkfOHd/HqCy4j45yqtWiEM3a/nrfxegAWSoucWJpipjjHQnmJYrVIuVahUqtSa9So1Wtsf/D7wIFTnsfJsKV4HHnsHs46d3tP1n7HjlEufZn3VFSr1pmZXmJ+rsTyYplSsUqlXKNarVGrNajX6zz78MMspLPdrKAXc+81Iq7AXCnlvwGuUkrd4P/8aUAopcLu19uVUuf7P18FfAZ4EzAJbFNK1aSUZwB3KKXOlVI+5//NA35SxRQwtg7369nAc5OTC7Gnz09MjHD8eBrPjelj8tzB7PnbuZs5dzB7/ibPHcyef9xzdxwRGKLOAfat+/fiGhBwG/AGKeWElHIQeAfQrHKqlNoPFKWUr/EPfQC41RdodwHv9o9/ELjV//93/J/xX7/LxtNZLBaLxWKxxCjqlFIHgU8CPwQeBm5RSj0gpfyOlPJy/23vAz4vpXwCGAb+wj/+Mbxs2b3AVcCn/OOfBq6QUj7mv2dVmRSLxWKxWCwWE4k1CEwpdQtwS9uxt4T+/wgrkyeC4/uBayKOTwFv6/lALRaLxWKxWDY5tqOExWKxWCwWiwZYUWexWCwWi8WiAVbUWSwWi8VisWiAFXUWi8VisVgsGmBFncVisVgsFosGWFFnsVgsFovFogFW1FksFovFYrFogBV1FovFYrFYLBpgRZ3FYrFYLBaLBlhRZ7FYLBaLxaIBVtRZLBaLxWKxaIAVdRaLxWKxWCwaYEWdxWKxWCwWiwZk0h5AQrgAjiMSOVlS5+lHTJ47mD1/O3dzMXn+Js8dzJ5/nHMP/W13I78nGo1G70fTf1wJ3JX2ICwWi8VisVg2wFXAj9b7ZlNEXR54BXAYqKU8FovFYrFYLJZuuMBu4MdAab2/ZIqos1gsFovFYtEamyhhsVgsFovFogFW1FksFovFYrFogBV1FovFYrFYLBpgRZ3FYrFYLBaLBlhRZ7FYLBaLxaIBVtRZLBaLxWKxaIAVdRaLxWKxWCwaYEWdxWKxWCwWiwaY0vs1dqSU7wU+BWSBLyilvpjykGJBSjkK3AP8ilJqn5TyWuBzwADwdaXUp/z3XQZ8GRgF7gQ+opSqpjTsniCl/CPgN/wfv62U+n1T5i+l/CzwTqAB/I1S6nOmzD1ASvnnwLhS6vpOc5RSngncDOwAFPA+pdRCaoPuAVLKH+LNp+If+m3gPCLud52+E5sZKeVbgT8ChoDvKaX+FxO++1LKG4CPhw6dA/wd8D/RfO4AUsr3A/+b/+OtSqlPbIbr3lrqeoCU8jTgj/F6zF4GfFhKeXG6o+o9UspX4fWgu8D/eQD4CvB24CLgFVLKN/tvvxn4uFLqAkAAH0p+xL3Dv4lfB7wMb41fLqV8DwbMX0r5WuD1wKXA5cDvSilfigFzD5BSvgH4N6FDneZ4I3CjUupC4EHg04kOtMdIKQXe9f5SpdRlSqnLgANE3O/WuB9sSqSU5wJfAn4V7/v/C/6ctP/uK6W+HFrz9wHHgD/FgLlLKQeBvwBeC7wUuMrfA/r+ureirjdcC/xAKTWllFoEvoFn1dCNDwG/Axzyf34l8JRS6jn/iexm4F1SyrOAAaXUff77bgLelfRge8xh4D8opcpKqQrwON5mp/38lVJ3AK/z57gDz8K/FQPmDiCl3I4nYv7E/zlyjlLKLHA13vXfPJ7oYHuP9P/9npTyESnlx+l8v4u8H6Qy6t7xa3jWqAP+df9uYAlDvvsh/hL4j8C5mDF3F08fDeFZo7N4luq+v+6tqOsNe/A2/YDDwOkpjSU2lFI3KKXuCh3qNG/tPg+l1GPBxSylfBGeG7aOOfOvSCk/A+wFbsegtQf+CvgkMO3/3GmO48BcyOWkw9y34a33rwFvAD4CnIk5a38+4EopvyWlfBj4GGZ99wMvxYBS6h8wZO5KqXk8a9sTeJbpfUCZTXDdW1HXGxy8WKMAgbfh606neWv7eUgpXwx8H/g94FkMmr9S6o+ACeAMPCul9nP344peUErdHjq83u89bOK5Ayil7lVKfVApNauUOgH8DfBZDFh7nwyeZfK3gF8EXoVnrTJl/uDFUH7O/78R93wp5aXAvwXOwhOsNbzwm76/7q2o6w0HgN2hn3fRclHqTKd5a/l5SClfg2e1+EOl1N9iyPyllBf6AcIopZaAbwLXYMDc8dxt1/lWms8CbwNuIHqOx4AtUkrXP76bzT13pJRX+vGEAQLPamHC2gMcAW5TSh1XSi0D/4Qn8oyYv5QyhxdX9i3/kBH3POBNwO1KqWNKqRKeS/UaNsF1b0Vdb7gNeIOUcsIPsHwH8N2Ux5QE9wNSSgzMEgYAAAMlSURBVHm+/4V+L16W0H6g6IsggA8At6Y1yF4gpTwDL+vrvUqpv/cPmzL/c4G/llLm/Zv82/FcktrPXSn1RqXUJX6w+P8OfEsp9ZtEzNGPuboLTwgCfJBNPHefrcCfSSkLUsoRvGSR9xN9v4u8HtIaeI/4F+BNUsqt/pzejBc7pf133+dS4Ek/dhLMuec9AlwrpRzyk4XeCtzBJrjurajrAUqpg3gxNz8EHgZuUUo9kO6o4kcpVQSuB/4RL9bqCVrBou8DPi+lfAIYxssk2sx8AigAn5NSPuxbbq7HgPkrpb4DfBt4CPgJcI8vbK9H87l3odMcP4aXDboXuAqv7MemRSn1L6xc+68ope4m4n63xv1gU6KUuh/4P/Cy/vcC+/GSBq7HjO/+uXhWOMCce75S6nvA1/C+8z/DS5T4r2yC6140Gu2uYIvFYrFYLBbLZsNa6iwWi8VisVg0wIo6i8VisVgsFg2wos5isVgsFotFA6yos1gsFovFYtEAK+osFovFYrFYNCCT9gAsFoulH5FSfg+vDtf/AD6hlNqb8pAsFoulK1bUWSwWSzRvBFBKvSXtgVgsFst6sHXqLBaLpQ0p5Vfxiqw+ClyM1/NzGPgvwPOABBbxCpL+O//nf1RK/Xv/99+KV4A0ByzhWfruTXYWFovFNKyos1gslgiklA1gAngQeCeeqLsNeIVS6iEp5a3AFryekKN4/R7PBobw++MqpSallC/2f+/8ULsli8Vi6TnW/WqxWCzr5zml1EP+/58BZpVSZeCElHIO2A5cjdfU+3YpZfB7deB8vJ6SFovFEgtW1FksFsv6KbX9XIl4jwvcrpQKGnwjpTwDz5JnsVgssWFLmlgsFks0NbxG3hvlduA6KeWFAFLKt+A1BR/o4dgsFotlFdZSZ7FYLNH8A3AHXizdulFK7ZVSfhj4eymlAKrA25RSCzGM0WKxWJrYRAmLxWKxWCwWDbDuV4vFYrFYLBYNsKLOYrFYLBaLRQOsqLNYLBaLxWLRACvqLBaLxWKxWDTAijqLxWKxWCwWDbCizmKxWCwWi0UDrKizWCwWi8Vi0YD/H35ser3lYI5+AAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plt.figure(figsize = [10,8])\n", "for i,T in enumerate(TT):\n", " \n", " ptc = PTC[i]\n", " Ni = ptc['Ni']\n", " Ci = ptc['Ci']\n", " \n", " F_3n = Ni**2.5/(Ni**2.5 + (17*1e-3)**2.5)\n", " K_con1 = 0.1*1e-6\n", " K_con2 = 20\n", " K_coff1 = 0.05*1e-6\n", " K_coff2 = 0.3\n", " K_12 = 0.2\n", " K_43 = 0.1\n", " K_inact1 = 0.2\n", " K_inact2 = 25*K_inact1\n", " \n", " k = {\"Ni\": Ni,\"Ci\":Ci, \"F_3n\": F_3n,\\\n", " \"K_con1\": K_con1,\"K_con2\":K_con2,\\\n", " \"K_coff1\":K_coff1,\"K_coff2\": K_coff2,\\\n", " \"K_12\": K_12,\"K_43\":K_43,\"K_inact1\": K_inact1,\\\n", " \"K_inact2\": K_inact2}\n", "\n", " y0 = [1,0,0]\n", " t = np.linspace(T,T+Tc)\n", " '''solve ODEs'''\n", " y = odeint(model,y0,t,args=(k,))\n", "\n", " y1 = np.empty_like(t)\n", " y2 = np.empty_like(t)\n", " y3 = np.empty_like(t)\n", " y4 = np.empty_like(t)\n", " I = np.empty_like(t)\n", " for i in range(len(t)):\n", " y1[i] = y[i][0]\n", " y2[i] = y[i][1]\n", " y3[i] = y[i][2]\n", " y4[i] = 1 - y1[i] - y2[i] - y3[i] \n", " I[i] = y2[i]*F_3n\n", " \n", " y0 = [y1[Tc-1],y2[Tc-1],y3[Tc-1]]\n", " \n", " '''Plot the simulation'''\n", " plt.plot(t,I,linewidth=2,label='I-'+str(T))\n", " plt.xlabel('time')\n", " plt.ylabel('y(t)')\n", " plt.legend()\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "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.6" } }, "nbformat": 4, "nbformat_minor": 2 }