{ "cells": [ { "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": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "t = 0;dt = 0.001; s4 = 0;\n", "ci= 0 \n", "s5=0.7; factive=0.2; s7=0.01; s6=1-factive-s7-s4;\n", "tnow =[]\n", "ixnow = []\n", "F1=[]\n", "F2=[]\n", "F3=[]\n", "F4=[]\n", "for i in range(1,250000+1):\n", " t = t + dt\n", " ni = 0\n", " if t>30 and t<40:\n", " ni = 100\n", " ci =0.1\n", " if t>50 and t<60:\n", " ni = 100\n", " ci = 0.3\n", " if t>70 and t<80:\n", " ni = 100\n", " ci = 5\n", " if t>90 and t<100:\n", " ni = 100\n", " ci = 20\n", " if t>110 and t<120:\n", " ni = 100\n", " ci = 30\n", " if t>130 and t<140:\n", " ni = 100\n", " ci = 40\n", " if t>150 and t<160:\n", " ni = 100\n", " ci = 50\n", " if t>170 and t<180:\n", " ni = 100\n", " ci = 60\n", " if t>190 and t<200:\n", " ni = 100\n", " ci = 70\n", " #if t>310 and t<340:\n", " #ni = 100\n", " #ci = 70\n", " #if t>340 and t<370:\n", " #ni = 100\n", " #ci = 80\n", " #if t>370 and t<400:\n", " #ni = 100\n", " #ci =90\n", " #if t>400 and t<430:\n", " #ni = 100\n", " #ci = 95\n", " #if t>430 and t<460:\n", " #ni = 100\n", " #ci = 100\n", " \n", " \n", " f3n=ni**2.5/(ni**2.5+17**2.5)\n", " kcon1=0.1\n", " kcoff1=0.05\n", " kcon2=20 \n", " kcoff2=0.3\n", " kinact=0.05\n", " #kinact=1\n", " #kinact=0.001 \n", " #kinact=1\n", " #kinact=0.2\n", " \n", " s5=s5+ (s6*ci*kcon1-s5*kcoff1+factive*f3n*kinact-s5*0.3) *dt\n", " factive=factive+(s7*ci*kcon2-factive*kcoff2+s5*0.15-factive*f3n*kinact)*dt\n", " s7=s7+(factive*kcoff2+s6*0.1-s7*ci*kcon2-s7*f3n*kinact*25)*dt\n", " s6=1-s5-factive-s7\n", " \n", " F1.append(s5)\n", " F2.append(factive)\n", " F3.append(s7)\n", " F4.append(s6)\n", " \n", " incx=factive*f3n\n", " tnow.append(t)\n", " ixnow.append(incx)\n" ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "scrolled": true }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAl8AAAE8CAYAAAD+N9R0AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzt3XecXHd97//XlN3ZXrRa9bIq9lduarZly7iCwRhc6MRACO0mlAcEyA3lFxIg4YYbCDjlUu4NoRPABGwHgwvEtizLloskq1jWV7233ZW2a9vM/P6YWWmRd6WV9pw53zPn/Xw8/PDszOrsZ+acmfOebzuxbDaLiIiIiBRGPOgCRERERKJE4UtERESkgBS+RERERApI4UtERESkgBS+RERERApI4UtERESkgJJBFzBWzc2dvq+JUV9fwfHjPX7/GTlH2i/u0T5xk/aLe7RP3FSI/dLYWB0b7TG1fA2TTCaCLkFGoP3iHu0TN2m/uEf7xE1B7xeFLxEREZECUvgSERERKSCFLxEREZECUvgSERERKSCFLxEREZECUvgSERERKSCFLxEREZECUvgSERERKSCFLxEREZECUvgSETmLTDZL30A66DLOaO+RTtq7+oIuY1TdvQOs3nyYbNb3K8Wdt407W9my53jQZUgEKHyJSKDSmQz/+fgOmttOBF3KqH7w4BY+9LUVdHT3B13KiE70DfKF7z3H+/7ukaBLGdWPHrb8v//azO+e3x90KaO6+571fOWn64IuQyJA4UtEArVuawu/Xb2Hz/7f1UGXMqqVGw4BuZYRF6Uzudak/sGMsy1Lre29ADxvjwZcydkd73S3BVGKg8KXiASqoiwJ5Lr2XFWazH1U7jzUEXAlZ+dqcJgztQaAXQfdfw13HmwPugQpcgpfIhKoZOLUx5Cr46rmTssFh52hCA5u1lhVXgLkWulcbZ0b4uprKMVD4UukyN37xE6+dd+moMsYkz2HO4MuYUQ1laVArj4Fh/E7ctzN8X0lQy2cIXgNJdwUvkSK3K+f2s1zW47S7uhg8eHCcNJrzo9dclUYusxcrbG8NAHA7sOdZDJuh2wJN4UvkYjYccDNE95wYRhT5WpwGLL7SCfpTCboMs7I9ZDdN5DmYEt30GVIEVP4EomI7SEIX7scDzbgfnDoH8hwoNnt4OD6awjh+CIg4aXwJVLkJtWVA+EIX60dfbQ5vFAowI4D7p+UXd/X+4520dfv5uSKIdv2twVdghQxhS+RIpcaGsdyqJOBQXe7oxLxGABb97l70ptYW8aew52c6BsMupQRNeWXc7B73X0N50ytJp3JOhsQJ0+ooLIs6fRrKOGn8CUSEYPpDHuPuDmbEMDMqgPcDg4LZtWTybobHKZOrKSuqhS797izszLNrHoAtux18zI+8Viuxpb2XlocvuqChJvCl0iEbNvvZmgAmD2lmlRpwtmTMpwKiC7XuGBWPR09Axxq7Qm6lBFdML2WeCzmdMg++UXA4VZYCTeFL5EIcTk0JOJxLphRy6HWHmeXxZg/o5ZEPCTBwdF9XVaaoGlqNbsOdTg77muB461zEn4KXyIRUJ5KMmVCBXZfG4Npd8d9DZ30nA0OJbngsPuQu+O+TgUHtwOiy+O+pjdWatyX+ErhSyQiLm6qp68/zS6Hp9Cf7Nbb42b4glPjvlztwp1UX05dVSlb9h539nqZQwFx8+5jAVcysngsdnLc15HjbnbfSrgpfIlExMVNEwDYvNvdYNM0pZryVIIXdx9zdsD4xbNzwWHjztaAKxlZLBbjkjkT6OwZcPZyTRfOrKMkGWeDo68hwGVzc++XDTvcrVHCS+FLJCIWzKojFnO3tQFy474uaZpAc1svh4+52eJwwcw6ykoTbNjR4mxAXDRvIgDrt7cEXMnIUiUJFsyq50BzN62OXq5pYf413KjwJT5Q+BKJiIqyEpqm1LDzYIez45UALpvXALjb4pBMxLlkjtsB8eKmCSTiMWdb5wAWDu1nR2usr04xc1IVW/a2OTsxQMJL4UskQi6bO4F0JsumXe62fi2c63b4glM1utoqUlGW5IIZtew61OnszNGT4cvR1jnI1TiYzvCSw2MQJZwUvkQiZMkFjQCs29YccCWjq61K0TSlmq372pxtoRsKDusdDV/gfrdZY105UxsqeGnPcfoH3GxZGtrPLzgcECWcFL5EImTW5Com1KTYsL3V6SUnFs5rIJ3JOtttVluVYs7UXEDsOjEQdDkjWnxBLnytsUcDrmR0Sy5opH8w42wr57xptdRUlrJ2azPpjLvvFwkfhS+RCInFYiyeP5GevkG2Obx69+VmEgDPbXE3OFyxYBLpTNbZcDNlQgUzJ1WxadcxunvdDIjLLnJ7P8fjMa4wjXSdGFDXo3hK4UskYoa6Htduc7crZUZjJVMbKtiwo9XZrscrF+SCw7MvuRkcIBdu0pksa7e62c08c1IVkydUsH5Hi7OD2pddNBlwez9L+Ch8iUSMmVVHZVmS57YcdbYrJRaLceWCSQwMZpwdbzOxtpy502rYsve4s4Par8wHh+ccDQ6xWIxlCybRP5Bh/Q439/P8GbXUV6dYa5ud7qqXcFH4EomYZCLOsosm09Hdz0sOL7h6ssVh85GAKxndsgWTyGbheUe7zSbVldM0pZrNu90NiENdj6tfdHM/x/NfBHr6Bp0dmybho/AlEkHLL5kCwNMvHg64ktFNm1h5csySs8Hh4snEYzFWbjgYdCmjesVlU8lkszy18VDQpYxoemMVTVOqWb+jheOdfUGXM6JrL5sKwBPr3d3PEi4KXyIRNG96DY11ZazZ2kxvv5tjqgCuXzSNdCbLKkeDQ11VikXzG9h7pIvdh928ZubVl0ymJBlnxfqDzq7If8PiaWSzOBtiZ0yqYu60GjbuaHV2RX4JF4UvkQiKxWIsv2QK/QMZZ8cDASzPB4cn1h909iLRNyyeDsCKF9wMDpVlJVxhJnH0+Am27HVzhuuyiyaTKkmwcv1BMhlH9/OiaWRxNyBKuCh8iYzDsY5evvzjNWzf3x50Kefs+kXTiMXgv9fud7ZFpKKshGUL8sHB0an+l86ZQENNitWbjzg7M/OGxdMAeGzdgYArGVl5KslVF0+mtaPP2YH3yy6aTFlpghXrDzIwqIH3Mj4KXyLjsGXvcbbtb+fvf7wm6FLO2YSaMpZc0MjeI13sOOhmlxmcaln6/fP7A65kZPF4jBuXTKevP+1s69cFM2qZNbmKNfYoR4+7eT3Km6+YAcCDz+wNuJKRpUoT3LB4Gu1d/U6PlZRwUPgSGYey0uTJ22H8Nvyqpblg8+gaN4MN5ManzZtewwvbWzjQ3BV0OSO6ccl0UqUJHnlur5PHQSwW49arZpPNwkPP7gu6nBHNaKxi4bwGtu9vZ9t+N7tHX3PlLBLxGA89s9fZbnAJB4UvkXGIDbv90h53L1Y9mgWz65k+sZLnthylue1E0OWMKBaL8bqrZgPwkKOtIpVlJdywaBptXf2sdrRV5IoFjUysLePJDYdo73JzVuGtV80C4MHVbu7n+uoUyy+ZwuFjPaxzdOFaCQeFLxGPPG/D92Eci8V43fLZpDNZfrt6T9DljGrRBROZ2lDB6s1HaHE0JL7mypkk4jEeeHq3k4txJuJxbr1qFoPpDA885ea+vnBmHfNn1PLC9hZnx1HeevUsYjH41RM7nV2kWNznW/gyxsSNMd82xjxtjHncGDP/tMf/pzFmjTHmOWPMG/2qQ6RQ1m1tdrLL6WyWXTSJyfXlPLnhkLPT6OOxGLdf00Q6k+VXK3cGXc6IJtSUcdOS6TS39To7sP26RdOYVFfO4y8c4IiDY79isRhvvXEeAPc8vt3JiSBTGyq5buE0DrX28OQGN5dAEff52fL1BqDMWrsc+AzwtaEHjDF1wMeA5cBrgH/ysQ4R38Vi0N07yLpt4Wv9SsTjvH55Ltjcv2pX0OWMatnFk5k1uYrVLx5hz+HOoMsZ0W2vaKKsNMGvV+2mp9e9mY/JRJw33TCXdCbLL1e4GWIvmFHHkgsmsn1/O2u3ujnz8c5r51BaEue+J3c5O8NV3OZn+LoWeAjAWrsauGLYY93AHqAy/1/4mgtEhrluYW4q/8qQfhNefulkpjdWsmrDIWeDTTwW42035RrQf/bf25xsFampKOV1V8+m68QA9znaQnflgknMnVbD81uOsmmnm5fLecuN80jEY/zH77c6GW7qq1PcetVs2rv6+ZWjIVbc5mf4qgGGd9qnjTHJYT/vAzYDa4F/8bEOEd9Na6hg/oxaNu865uyYpDNJxOPc9aoLyAI//f1WJ4MNwMVNE1g8fyJ2X5uzXT63LJvJlAkV/Pea/ew46N64pVgsxrtvMcRjMX74sKWvPx10SS8ztaGS1y+fzfHOPn65YkfQ5YzodVfPZmpDBY+u3e/s+DRxV/Lsv3LeOoDqYT/HrbVDX2FuBaYCc/I/P2yMWWWtfXa0jdXXV5BMJvypdJjGxuqz/5IUnKv7pfpwbumDyqoybr9uLnf/dB1Pbj7C/7jzsoArOyWZjBOPnf01vKGxmic3HWb1psO8sOs4r8nPMByNV/vkaGfuuo2VFaVj2uaf37WUD3/lUe55fAc3LZtNfU2ZJ3WcSSpVAkBDQ9WY/t6f/9ESPvvNVfzoka3c/fEbKC3x97MrNezal2N5DRsbq3njjfP45WPbeeCZvXzwTQv9LA+AisoUALV1FWOq8T13XMrabS08tu4AN105i8UXTvK7ROLxOIlEfMzH9sfvWsqn/8+TfO+hLfzTJ26goqxkxN9z9fMr6oLcL36Gr1XA7cA9xpirgY3DHjsOnAD6rLVZY0wbUHemjR0vwODQxsZqmpvd7HKJMpf3S0dHrpWru7uPq8xE6qtTPPz0Hm5eMp2q8pE/iAttcDBDJsuYXsO3XD+X9dua+bf7NjKzoZyJteUj/p6X+6StLffe7u7pH/M233T9XH7yu6185YfP8fG3LSIei539H41DX98AAK2tXQzmb5/J5JoUNy2ZzmPrDvCvP1/Hu28xvtbXdeJUTWN9DV+9dDpPbzzEb1btYnZjJVcs8Dfc9HTnlrdob+sZc43vee0CvvzjNXz1R8/zhfcto64q5WeJZDIZ0unMmOtrrCrl1qtm8eAze/nHHz3PB++8hNhpx6LLn19RVoj9cqZw52e3471ArzHmKeBu4BPGmE8aY+6w1q4EngNWG2OeBrYCv/OxFhHfJRNxbrlyJn0DaacXLT2TCTVl3PWqC+ntT/PvD7zk7FT6m5ZOZ+G8BjbtOsZvnnZz2YS3v3I+MydV8fi6A6ze7N7aX6UlCT70hkspLYnzvQdf4lBrd9AlvczcaTW87ab5dPQM8K37Njk5m/iN189l/oxanttylIcdXcBW3ONb+LLWZqy1H7TWXmOtXW6t3WKt/bq19r/yj3/eWntV/rG/tNa6OchE5Bxcv3gaVeUlPPzcXjp7+s/+Dxz0isumsPTCRuy+Nn7xmJvjbeKxGB+47WIm1KS4b+VOXtjm3qy4oXBTVprgu7/ZwtZ97q3aPn1iJX/y2gWc6Etz9z3rnVx89eYrZnDlgkls29/Odx7Y7NzK8slEnA/deSl1VaXc89h2ntl8JOiSJAS0yKqIh8pKk9z+iiZO9KX59ardQZdzXmKxGO9//UVMmVDBI8/tY9VGNwe2V5WX8JE3XkZJMs6379/EjgPuDXqeMqGCD7/xUrLZLP/6yw3sd/DySMsvmcIbrptDS3svd9+z/g+6MF0Qi8X4wG0XcWG+delHD1vnAlh9dYpPvG0x5akE33lgs1a/l7NS+BLx2E1LpjOpvpzH1h3gYIt7XTljUZ5K8tE3X0Z5Ksn3frvF2fXL5kyt4UN3XspgOsvd96x3MoBdOqeB99y6gO7eQb7yH+vYdci9i5jffk0TNy6Zzt6jXfzDT9ZyvNOtFrCSZIKPvmUhMydVseKFg3zngc3OXUVg5qQqPvbmhSQTcb5x7yae2uTmlxZxg8KXiMeSiThvv2k+6UyW7z+4hUzGrW/pYzW1oZJPvHURyWSMb933Ii9sd69rD2DR/Im8/7aL6O1P848/e4EXd7t3jc1XXDaV9966gO7eAb7603Vs2OHW+lqxWIx3veZCbr5iBgdauvn7H61h92G3QmJlWQmfescS5k2vYfWLR/j6z1+gvdutrn0zq56/+KPFlJUm+M4DL/GLx7eTdiwkihsUvkR8sOTCRq5cMIntB9r575AOvgeYP6OWj755IfEY/OsvN7DiBTcvm7P8kil8+I2Xks5kuPvn63nomb3OrVV23aJpJ1vp/vkX67n/yV1O1RiPxbjrVRfwpuvncqyjl7//0VrnJgpUlpXwP9++hKUXNrJlbxtf/N6zHD7m1mWS5k+v5bPvWsrk+nIeXL2Xv/vuM07tZ3GDwpeIT9756gupKi/hF4/vcHbV+LG4pGkCf/mOJVSWlfCDh2xuYc4B9xbmXHphI3951xKqK0q457Ht/OoJ91Yev2LBJP7qjy+nuqKE+5/cxcadbrXSxWIxbrumiT9/6yIG0xl+8sjWoEt6mVRpgo+88VJetXQGbV39rNxwMOiSXmZ6YxV//SdX0lBTxpotR529ZqoER+FLZBzO9IW2prKU97/+IgbTGb5x70a6e90ayHwu5k2r5a/efTkzGnNLJ3zqX1Y6Oe3/ghl1fP69VwKwybFgM2T2lGruuDa3vrSLswsBFs5roGlKNf0O7mPIhcRrLpsC4Gy3fkVZkgWzc8tXulmhBEnhS8QDoy3xuWj+RG6/pomW9l6+8auNDAy612I0VpPrK/jcuy/n0rkT2Hmwnb1H3GzNq6tKkSr1/2oY45HyecV7L5y+WKiIeEfhS8Rnd1475+QYlW/f/6KzC5eORWlJgpmTqoIuQ0Qk1BS+RHwWj8f4szsuZsGsOtZta+Fb970Y6hYwEREZH4UvkQIoSSb46JsXsmBWHWu3Nju5mKWIiBSGwpdIgZSnknzibYu4fNg0+Z0H3VpLSURE/KfwJVJAJcnc9f7ecO0cjnX08eUfr+HeJ3aqG9IH2RDMMXO9wjAsTxWKGoMuQJyj8CVSYPF4jDuuncMn/2gxtVWl/Pqp3fzNvz/LGns0VIsxulyp6/P0wjCR0PUaXa8PIOb8kShBUfgSCcglTRP4u/dfxc1XzOBo2wm+ce8m/vb7z7PGHnV27SIRERk/hS+RcRlfSCpPJXnHzRfypQ9cxVUXT2bvkU6+ce8m7l3p3ursQ/RtXkRkfBS+RLwwzjwytaGSP7vjEj79zqUAHOvQ5UhERIqVwpeIQybWlgVdgoiI+EzhS0RERKSAFL5EpDhpzoIH9CJ6IkSzmKUwFL5E5Pw4fD5xfRmCMExacL3CMLyGYShRgqHwJSIiIlJACl/itGw2S2v7iVAtPuoFl5+t661KIiKuU/gSp63b1sJ7/vYRHl17IOhSRhSxTCgiIh5Q+BKnvbCtBYD7n9wVcCVnpsYgEREZK4UvcVrfQO6C010nBsiomUlEQkifXHI6hS9xWm9/+uTtI8d6AqxEwiYMJzzXxzI6Xh4QjhpFTqfwJU7r7R88eXv34c4AK5HTZZ2ON453BDteHhCOGh2nl1BGo/AlThve8rX7UITCl8u5RkRExkXhS5zW15+mqryEeCzG7sMdQZcjIiIybsmgCxA5k3QmS0VZkrqqFHuOdJLJZInHi7cxP6ZFtEREip5aviQUmqZW0z+Q4VBrd9CliIiIjIvCl4TCnCnVAOyK0rgvESkOGsMpp1H4klBomloD4Oy4L3UXuicMSxCEoETnuT3rVmRkCl8SCjMaq0jEY1puwiEuhxvXo7Dr9QHEHK8yDN93wlCjBEPhSxyXO8OXJOPMaKxi75EuBtOZgGvyn8O5RkRExknhS9yX//rYNLWawXSGgy0adB8kfZsXERkfhS8JjaaTg+7dHPclIiIyFgpfEhpz8oPudxxQ+BIRkfBS+JLQmDGpisqyJFv2Hg+6FBFxRQgGSIagRCkwhS8JjXgsxoUz62hp76W1vTfocsR5ITjlOV6iyzNaRcJM4Uucdvpnv5lVD4Dd50brl85NbnJ+UoDr9UE4anSeXkQZmcKXhMqCWXUAbNnbFnAl/sqqyUFEpGgpfInzhn93nNFYRUUqiS3ScV/Ot9gA+jYvIjI+Cl8SKvF4DDOrjua2Xo4e7wm6HBERkXOm8CWhc9m8BgA27GgNuJJwUA+miIhbkn5t2BgTB74JLAL6gA9Ya7cPe/xW4PP5H9cCH7HW6jQhZ7Vw7qnwdfMVMwOuJhzUUSjFKgwnDY3hlNP52fL1BqDMWrsc+AzwtaEHjDHVwFeB26y1VwO7gYk+1iJFZEJNGTMnVbFlbxt9/emgyxFHheF053qNWecrdFs4xnBKEPwMX9cCDwFYa1cDVwx77BpgI/A1Y8xK4Ii1ttnHWiSkRvvCuHBeA4PpDC/tKc6B92Ggb/PnLxaCtkjXK4wp2UiI+dbtCNQA7cN+ThtjktbaQXKtXDcBi4EuYKUx5mlr7dbRNlZfX0EymfCx3JzGxmrf/4aMXSIeIxZ7+X65/vKZ/ObpPby0r41XXzMnoOqgOn+po+qqlCfHTiJVAkAqVeLZsZhMxonFY55tr7KyFIC6ugpPtnm0sz+33YpSz2qMx2MkE3HPtpfK75eGhirqa8rGvb2amtxHo1fHTaq7/+Rtr55zSUmCmIfbq6hMAVDr0XHTNZABoLzcu/dKPB4n4eFxU1aWO24mTKjUucVBQe4TP8NXBzD8mcXzwQugFXjOWnsYwBjzBLkgNmr4Ol6AmW2NjdU0N3f6/ndk7NKZLEl42X5pqCyhrqqUpzYc5K03zCWZCGbuSEfHCQA6u/o8OXbau/oA6Osb8OxYHBzMkM1kPdteT/5E39bW48k229py7+3unn7PasxksgymM55tr69vAIDW1i4G87fHo6Mjd4UGr46brhOnavLqOQ8MpMl6uL2e7tyx3e7RcTN0Tjhxwrv3SiaTIe3hcdPbm9svx451k1JDnVMKcb4/U7jz84y1CngdgDHmanLdjEPWAJcaYyYaY5LA1cBmH2uRIhOPxbhiwSS6ewfZvPtY0OVobIeIiIyZn+HrXqDXGPMUcDfwCWPMJ40xd+THd30WeBh4BviVtXaTj7VIEVp20WQAnn3paMCViIiIjJ1v3Y7W2gzwwdPu3jLs8Z8BP/Pr70vxmzethoaaMtZubaZ/IE1pif9jAkUiJQxzKsJQo8hptMiqhFYsFuOqiyfT259mzVZNlpXThOGkrBmjRU2jEWQ0Cl/ivDNNy79+0VQAVrxwsFDlFITOyePj+jIEjpcHuB8cXK9P5EwUviTUJtVXcHFTPVv3tXGwpTvocsYvBGflEJQoIuI0hS8JvRsWTwfg8RcOBFyJiIjI2Sl8SegtuWAitVWlrNxwiO7e8a/BJCIi4ieFLwm9ZCLOa66cSV9/msfWBtP6pZ44EREZK4UvcdpYrx944+LplKeS/O75ffQP6GLbIl4Iw7yPMFz8WxNo5HQKX+K+MTQrlaeSvHLpdDp7Bng0oNYvL+mzevzC8BqGoUYZB81OkVEofEnRuGXZLCpSSX7z9G6N/SoAfZsvcq4HB8fLEzkThS8pGlXlJbz+mtl09w7ym6f3BF3OedH5RESk+Cl8SVG5+fIZTKhJ8fvn9xXHul8iIlJ0FL6kqJQkE7zj5gsZTGf5wUNbyKhvTEREHKPwJU47n+i09MJGLjeNbNvfHtjSEyIiIqNR+BLnnc84qHe++kIqy5L8/NHt7D3S6XlNQ9SwJuOh42f8wvAShqFGKawxhS9jzL+OcN8PvC9HxBt1VSnef9vFDKYzfOu+TZzoG/T3D3o9M0xn5XEb6xpxEk5hmJwShholGMkzPWiM+Q4wF7jCGHPJsIdKgFo/CxMZr8XzJ/Laq2bx0DN7+dZ9m/jYWxaSTKix1ysuRxvnV0lwvUAUHET8dMbwBXwJaAL+GfjisPsHgZd8qknEM2++YS4HW7rZsKOVHzy0hfe97iK3T3wOlzbE6ddPRCQEzhi+rLW7gd3AImNMDbnWrqFP3irgmJ/FiYxXIh7nQ3deyld+upZVGw+TiMd59y2GeFwBQkREgjHWMV+fBfYDTwAr8v897l9ZInke9G2lShN8/K2LmDW5iifWH+T7D20Z/0ZFRETO09m6HYd8AJhnrW32sxiRkXjRy1VdUcqn7lrCl364hic3HOKuV11AeWqsh7+IiIh3xjr6eC/qYpSQqygrYcqECsD9yYSOlycR4vysUcfLA9z/wJGCG+tX/23Ak8aYx4DeoTuttX/rS1Ui4j6dUCRIYZj4EYISJRhjbfk6ADwE9JE7nIb+kyLQ0zvAum3NZDI6mZ6rrMdfu/Wm8obrr6Pr9UE4so1IWI2p5cta+0VjTCUwD9gElFtrddXiIvHLFTt5bN0BbrumiTddPzfockIpSuepKD1XERE/jHW24yuB9cD9wCRgjzHmNX4WJoXz4q7ccL419mjAlbycf21xauUTEZFgjLXb8cvAtUCbtfYwcD3wVd+qkoLqPDEAwKHWHrp7BwKuZiTetbWoK0VERII21vAVz4cuAKy1m32qRwrsRN/gH1z3cNehjgCrERERKX5jne243xhzG5A1xtQBHyG3/ISEXFe+1auyLEl37yA7D3Zw6ZyGgKsS9YpGg/PLOJA7FF1uMHb/FQxHjVJYY235+jPgncBMYAewGPhTv4qSwunrTwOwYHY9ALsPdQZZTsHow3D8XH8NQ5BrZBxcDoRDwlCjBGOsLV8fs9be5WslEojegVz4mlRfTn11il2H1e0YJF202iN6HUXEYWNt+brdGKNPsyI01PJVVpKgaUo17V39HO/sC7iqYVxvvnC8PF/ok0BEZFzG2vLVCmwxxqwFTgzdaa19ny9VScH05sNXqjTJnKk1rNvWwu5DHdRXNwZc2SlqxBARkWIy1vD1A1+rkMD057sdy0oTTJuYu+7hrsMdLLnQnfAVCgqIIiIyRmMNX++w1t7iayUSiKExX6mSBE1TaoBoDLp3vTdTRESK11jHfJUbY2b6WokEYmjMV6okQVV5CY11Zew61BGKKfDnIywD2ovz1S+sMLyGYajR+SLD8FkVghKlsMba8jUJ2G3kJMF9AAAaxElEQVSMOUpuzFcMyFprdSHAkBtIZwBIJnOhpGlKDc9tOUpLey+NdeVBliaO0/mkuLn+RcXx8gCIaTyCjGKs4UtdjkUqPRS+4rlG0DlTc+Fr16EOJ8KXTvByPlw/5YUhOIiIf8Yavm4Y5f4felWIBCOdycWbZCIXvpqmVAOw+3Anyy6aHFhd4i7lBhGR8Rlr+Lpp2O0S4DrgCRS+Qm8w3/KVSOROqbOnVBMDdjt0jUe1EoiISDEZU/iy1r53+M/GmAnAz32pSAoqnc61fCXiuYRTnkoypaGC3Yc7yWSzxIss+Xj9bNQtKiIi52qssx1P1wU0eViHBGQw3+2YSJw6FJqm1NDbn+Zwa09QZYWO56EuDDO4RETkvIyp5csY8xinvuTHgLnAb/0qSgrn5ID7xKn4MH96DU+/eJit+9qYNrEyqNJExicMATYUJWZxeaRfCF7CUNQohXXW8GWMqQe+CRzN33UD8DfW2pV+FiaFMXhatyPAgtn1AGzZe5wbl0wPpK4hYTh/eilUvbwR2zci5yxM72cpqDN2OxpjlgCbgU5r7Qpr7QogBfzMGLOwEAWKv9KZoZavU4fClAkV1FaWsmVvW9F2fxXr85I85096zhcoIj4625ivfwTustY+NHSHtfavgPcBX/ezMCmM0wfcQ25xxQWz6+no7uegE+O+PDxR6Zw3fqFqnhMRcc/Zwle9tfbx0++01j4MTDzTPzTGxI0x3zbGPG2MedwYM3+U33nQGPPBc6paPHP6Ol9DFsyqA2DLnuMFr0lERKSYnS18lRhjXvY7+ftKz/Jv3wCUWWuXA58BvjbC73wJmDCWQsUfJ9f5iv9ha8ZF+XFfm3cfK3hNIiIixexs4WsF8PkR7v8c8PxZ/u21wEMA1trVwBXDHzTGvAXIAA+OqVLxxWgtX5PqK5hcX87m3ccZGEwHUZqvNOJLRESCcrbZjp8FfmuM+RPgBaAXWEpu5uMdZ/m3NUD7sJ/TxpiktXbQGHMp8A7gLcDfjKXQ+voKksnEWH51XBobq33/Gy6JxWPEYjB5cs3LHlu+cBr3rdjBofY+Ll8QzKWGhi7u69V+SaVyh/zEhipqq1Lj3l51dVv+/2We1Fh+YgCA0tKkZ885mYwTi8c8215lZa7Ru7a23JNtHu3sz223otSzGhPxGPFE3MPjpgSAhoYq6mvKxr292sOdAFRWpTypMdXdf/K2V8+5tDT3eTtxYvXLvpydj4rK3Puttq7Ckxr78t+gyspKPHvO8XichIfHTXl57ripr/fmOYu3gtwnZwxf1tpOY8z15C4vtIRcS9U3xrjMRAcw/JnFrbWD+dvvBqYDj5JbrLXfGLN7+MD+0x0/7v/A78bGapqbO33/Oy7p7RskEY+P+LzN9FwgW/H8PmY1VBS6NODUrESv9kt/X+4QbGntov9E/1l+++w6O3oB6Orq86TGnt5c+OrvH/TsOQ8OZshmsp5trzt/om9r7/Fkm21tufd2d0+/ZzWmM1nS6Yxn2+vry+2X1tYuBvO3x6O93dvjpuvEqZq8es4DA+mT2/MifPV09wHQ3ubNcTN0TjhxYsCz55zJZDw9bnrz++XYsW6qSsb/Gop3CnG+P1O4O+s6X9baLLmQ9Og5/t1VwO3APcaYq4GNw7b5qaHbxpgvAIfPFLzEP4PpzMnrOp5u/oxaKsuSrN/RwruyF55shSo0Ta6Tc+X6IaNjWiTa/Izi9wK9xpingLuBTxhjPmmMOVt3pRRQOpMlGR/5TJCIx7lsXgPHOvrYdShaLYIyOuUGEZHxGdPlhc6HtTYDnL6ExJYRfu8LftUgZ5dOZ18203G4ZRdNZvWLR1j94mHmTnv5uLDQ0oh7EREJiDqhIy6TzRI7Q/i6dM4EqspLePalIydXwxf/aQF+EZHipfAVcZlMlvgZBqAkE3GuXDCJjp4BXtpd+AVXPb8MkAbbiIhIwBS+Ii6b5YzhC+DqS3LLTDy58VAhSoo4hUPPhKH1UE2cHtBrKOGj8BVxmWz2rI1B86fXMn1iJWtsM21dfYUpbBiX40g2yh/8EX7q4+XyMT3E9Rpdrw8ISZESBIWviMtks8TPMOYLcgudvvLyGaQzWVa8cLBAlflLueH8haHnNqhlUURExkLhK+LG0u0IsPySyZSnEjy+7gADg+EdeK9TsoiIBE3hK+IymbN3OwKUlSa5buE02rv7WaWxXyIiIudN4SvismPodhxy61WzKEnGeeDp3Qymw9v6JSIiEiSFr4jLjLHbEaC2KsVNS6ZzrKOPleuLY+yXiIhIoSl8RdxYZjsOd+tVs0iVJLh35S66e8d/geGx8GXwtKNT/DVO3DthmInqfoXOvlVOcr0+CEeNUlgKXxGXPcsiq6errUpx+yua6DoxwH0rd/lYmT8Ubryj88k4hOE4dP3N4nh5ALEwFCmBUPiKuEz23FuWXn3FTCbXl/Po2v3sPNjhU2UiIiLFSeEr4nID7s/t35Qk47z7tQvIZuH//fpFevsH/SkuBNSdICIi50rhK+Iy2XPrdhxy0ex6blk2k6PHT/CT3231/hqMPgtXtSIiUkwUviIsm82SPY9uxyFvun4esydXs2rjYX733D6Pq8sJS6bzemRH2MKsiIiMncJXhA2d38e4zNfLlCTjfPTNl1FbVcrPH93O81uOeleciIhIkVL4irBMPn2NdZHVkUyoKeNjb15IaWmCb9//ogKYOCMMjYdhqNH1Tnq3qxMZmcJXhA11bY13Ha05U2v4i7ctprQkzrfvf5Gt+9q8KE8cF47g4KYwLEHgeoWu1wfur9YhwVH4irBM/gpB5zPg/nTzZ9Ty1hvnkclmQxG+ohUcvH2yvix667EQlCgiEabwFWEnux09OlE11pd7syEREZEipvAVYV51O75su45uC8LRagN+PG+PNygiIudN4SvCMkOzHT1q+vJrHIuCg4iIFBOFrwjzutsx0jx6DRU0RUSKn8JXhGUz/nQ7Rmw0e4RpP0eB829n1+sDsmEoUgpK4SvCvO52DMPc7xCUKFEQhgPR9RrVTCwhpvAVYVl1O8p5CMPhEoYaRSS6FL4iLONTt6O3DexqrhcRkeKi8BVhpwbcezXb0R9hWA1cRERkrBS+ImxoIG0Uh05kHR1FrKApIlL8FL4izIsLa4/IzVyTo4mdIiISMIWvCDs529Hxbkdxk+vB0/X6IBxLELheYSheQ/dLlAJT+IqwU+t8BVxIEVB3oWvc3h9uV5fj+jHtdnU5+myV0Sh8RZjXA+6HeHptR31jdI9OKCIi46LwFWFZzxdZ1VnZKwqdIiLFS+ErwoZavpzPTB7W5/pTFRGR4qfwFWFDi6x63u2oZpvzp3QoIlL0FL4izOtuR+UGERGRs1P4irBMhFuoIvzUPeP+S+h+hWEo0fkaXa9PZAQKXxGWPTnmy902K9c/V8OwxlAUOXxIA+7XB+7X6Hh5gPvLdUhwFL7Es48H1z+sc0JRpNN0QhERGR+FrwjLhOTajq7XB97XqBY1EZHipfAVZT51O2o8lYiIyOgUviJsKCOFoWXJa662LEVwV4iIRI7CV4QNtVDphC8i4h/1Bsjpkn5t2BgTB74JLAL6gA9Ya7cPe/wTwB/lf/yttfaLftUiIzvZ+uPwtR29bqCKYiufbxw/oTheHhCWGt2u0u3qREbmZ8vXG4Aya+1y4DPA14YeMMbMBd4JXAMsB15jjFnoYy0ykqFFVj27tKOSjcjY6L0ybmF4CcNQowTCz/B1LfAQgLV2NXDFsMf2Aa+11qattRmgBOj1sRYZQSYkXxkj+fkVkn0jIiLnzrduR6AGaB/2c9oYk7TWDlprB4AWY0wM+Cqwzlq71cdaZER+LbIaguTgVYkheKoiIuIWP8NXB1A97Oe4tXZw6AdjTBnwXaAT+PDZNlZfX0EymfC8yNM1Nlaf/ZeKRE1zNwDVVSlPnvfRzn4AKiq82R6cGqPl1fbKykoAmNBQReOEinFvr7q6DICamjJPauwfSANQWpr07DknkwlisZhn26uqSgFQW1vu6XFTWVHqWY2JRBwyWc+2l0rljpuGhirqa8rGvb3alh4AKiu9ea+kuvtP3vbqOZeW5k4PEydWU54a/6miojJ/3NRVeFJjNpE7H5SVlXj2nOPxOIlE3LPtVZSXArnzV5TOLWER5D7xM3ytAm4H7jHGXA1sHHog3+J1P/CotfYfxrKx48d7fClyuMbGapqbO33/O65oazsBQHd3vyfPu60tt496evo8ex2HZgl5tb2+3gEAjrV2EU+nx729zs5cb3lHR68nNQ4M5mrq7x/07DkPDqbJZrOeba+rqw+A9vYTnh433T3eHIcA6XSGjIfPua8vd9y0tnYxmL89Hu3tQ+89b94rXSdO1eTVc+7vz31XbmnppKx0/KeKnu78cdPW40mNrfnXsLd3wLPnnMlkSKcznm2v50QuFB8/3kNzuZ+nWzlXhTjfnync+Xk03Au82hjzFLlhO+81xnwS2A4kgBuAlDHm1vzvf9Za+7SP9chp/JrF5OW0atdnWokUO9eXSXC9PtDnmLycb+ErP5D+g6fdvWXY7fG33cv4eD7b0ZvtFII+CsfP9RNKKE7KYShSRDynRVYjLBOWJe5dr0+c4/oh43p94P4s4zBc4N39CiUoCl+Rlp/tGHAVBRWSJ+tye0gYgoOIiMsUviIs6/EC92H4JioiIhI0ha8IO9Xr6PHlhVxutnGeAqyISLFT+IqwocG+Tncj+RTktMaqiIgEReErwk52O3q1QZdDnM+cDrAiRS0EX4FCUKIUlsJXhJ1cKsDrbkePP2m8rE7j0jykE4qIyHlR+Ioyf7KXSOBcD9luV5fj/TVfveV4eTlhqFECofAVYV53O+pzxkMOz1rQfhYRGR+Frwgb6h70/Buuu7nhFEfDTSi+zYuIyLgofEWY5wPufeBmRBIRETl/Cl8RdrLxJ0r9jmGoUUREiprCV4Sd7Hb0OJF43Vqlrjg5H2G4aHUISnS+RsfLA8JRoxSWwleEeX15IYkWnVBERM6PwpdE8tqOnrfOOd566KkwpHXXS3S9PvFEmD4TpbAUviLs1OWFojPbUR+FIiISNIWvCAvDbEcREZFio/AVYScbqLzqdvQhxbk+2FdERORcKXxFWT7ZxJ2/tqMfqc77TYqIiIyFwleEZRRApIiF4fAOQ43OV+l4eUA4apSCUvgS5y+g66UIPVXfqUtYROT8KHxF2MnZjp5v1+MNOsyvhTxdfg3DkF9drzEMSxC4/kUlDF8aQ1CiBEThK8K8XmRVHzQiIiJnp/AVYacaV1xOTT61LHm9QQVYEREZI4WvKDs52zHgOs7G9fpERETOgcJXhJ2c7RipywuFoUYRESlmCl8SktAk7nF4VgA4Xx7g9syKPNcrdL0+8H7tQwk/ha8IO3VtR6+36+32REREionCV4Sd7HWMYMOXX0tEREIYjhfXD2rHyxMRfyl8Rdipli93zwRRy0jqAhYRKX4KXxF2cp0vr7fr8PgGz7tYvd2ciIhEgMJXhJ0MDo6vURWGtqAw1CgiIm5Q+Ioyv7od1RwkIiIyKoWvCPOr21GiwfXxeI6XB4SkRseLDMPkmRCUKAWm8BVhHq+xKh4KwwlFRETOj8JXhPk129Hl2KCgOX5heA1dr9H1+sDtWdDg/moiEI4aJRgKXxF2sttRHxDu0L4QESl6Cl8R5nULlV/flF3/Bi4iInIuFL4izLdFVl3ud8zzbEhVCJ6riIi4ReFLotXtGKknKyIiLlL4irBTS014E0giHWsi/eQdFYYZoyEoUUS8p/AVYae6HT3erodnlDCcP0VERM6FwleEeX15oTBxNdOFYVeEYQKE6yU6Xh7gfo2u15cTjiql8BS+oiyfQOJenan0OSMiInJWCl8RlvGpT8/VViVQPhQRkeAl/dqwMSYOfBNYBPQBH7DWbh/2+P8A/gwYBL5krX3Ar1rkzFzvohERESkmfrZ8vQEos9YuBz4DfG3oAWPMFOBjwCuAW4AvG2NSPtYiI9BsRxERkcLzreULuBZ4CMBau9oYc8Wwx5YBq6y1fUCfMWY7sBB4zsd6zqind5D7n9jBseM9QZVQcDsPtQPet3ztPtTJA0/t9mRbfl1gesW6A9RUlo57O7sOdXhQzcu1tPd69hp29Ax4sp3TPW+PcvjY+N8vLe29HlTzcr0Dac9ew4Mt3Z5s53R2X5snNfYPpsdfzCgefnYvqZLEuLdj97V5UM3LHWjp9mw/9w1kqCz3ZFN/4KlNh9jq0/OXc5eIx3jddfMCrcHP8FUDtA/7OW2MSVprB0d4rBOoPdPG6usrSCbH/wEwmifXH+A792/ybfsumzm9jsb6inFvJ16aJBbLBRIvQ0lNZSmNjdWebGvSxEoAfr9mvyfbGzJ9Sq0nNWazWarKS2hp7+VXT+z0oLKcqRMrPXsNp07KbefZl47y7EtHPdkmwJTGKs9qrKsu41Brj6evYTIRZ9rUWirKSsa9ra6BDAAv7TnOS3uOj3t7Q7x8r0yckPtM+M3TezzZ3pCZ0+s8qbGqppxkIs6B5m5+1ezdfq6rLvPsNZzSWAXAE+sPebI98U6qvIS332wC+/sxv1oWjDFfB1Zba+/J/7zfWjsjf/sO4LXW2g/nf74X+F/W2udH215zc6ev47gzmSwt3QMcben08884p64yxYxJVZ5tb39zF21dfZ5tD+DKy6Zzosub1pGBwQzb97eR9vC4L08lmTu1xrMlGFraTnDY4xbY6ROrqK/2pmc/k8nS2jPAkWbv3isliTjzpteSTHgzEqKjp5+9R7x9LzfUlDG1odKTbWWzWfYe6aLzRL8n24PcEiBXL5xOZ8cJT7bX159m+8F2T1ufq8tLmTW5yrP3yqHWblo7vG05nTW5mpqK8beKAwymM7R0D9DS2uXJ9sQbiViMqxfPoL3N356uxsbqUQ90P1u+VgG3A/cYY64GNg577FngfxljyoAUcBEQaLNTPB7jkrkNTKr25k0XVTMaq5jR6F2YA6gqL/EsfJUk41zUNMGTbfllYl05E+t86PvwSDwe4+I5DTRWufteqako5dI5DUGXMapYLMbsKd60rgxXlkriVeRMlSa4xPH3ytSGSs8CsR+SiTiXzZtIc42GNLum1IOu9PHwM3zdC7zaGPMUubHY7zXGfBLYbq39L2PMvwAryQ36/ytrrT8DP0REREQc4lv4stZmgA+edveWYY//G/Bvfv19ERERERdpkVURERGRAlL4EhERESkghS8RERGRAlL4EhERESkghS8RERGRAlL4EhERESkghS8RERGRAlL4EhERESkg367tKCIiIiIvp5YvERERkQJS+BIREREpIIUvERERkQJS+BIREREpIIUvERERkQJS+BIREREpoGTQBbjAGBMHvgksAvqAD1hrtwdbVXQZY9YB7fkfdwH/F/hnYBB4xFr7xaBqixpjzFXAP1hrbzTGzAe+D2SBTcBHrLUZY8zngdeT2z8ft9Y+G1jBEXHaflkK/BrYln/4W9ban2u/FIYxpgT4LtAEpIAvAZvReyVQo+yX/TjyXlH4ynkDUGatXW6MuRr4GnBnwDVFkjGmDMBae+Ow+14A3gzsBH5jjFlqrV0bTIXRYYz5FPDHQHf+rq8Dn7PWPm6M+TZwpzFmD3ADcBUwE/glcGUQ9UbFCPtlKfB1a+3Xhv3OUrRfCuVdQKu19o+NMQ3AOuAF9F4J2kj75W9x5L2ibseca4GHAKy1q4Ergi0n0hYBFcaYR4wxjxpjrgdS1tod1tos8DDwqmBLjIwdwJuG/Xw5sCJ/+0HgZnLvnUestVlr7V4gaYxpLGyZkTPSfnm9MeYJY8y/G2Oq0X4ppF8Afz3s50H0XnHBaPvFifeKwldODae6uQDSxhi1CgajB/hH4Bbgg8D38vcN6QRqA6grcqy1vwQGht0VywdgOLUfTn/vaP/4bIT98izwl9ba68m1Dn8e7ZeCsdZ2WWs78yfy/wQ+h94rgRtlvzjzXlH4yukAqof9HLfWDgZVTMRtBX6c/xayldybYsKwx6uBtkAqk8yw20P74fT3jvZP4d1rrV0zdBtYgvZLQRljZgKPAT+y1v4Heq84YYT94sx7ReErZxXwOoD8mK+NwZYTae8jN+YOY8w0oALoNsbMM8bEyLWIrQywvihbZ4y5MX/7VnL7YRVwizEmboyZRe6LS0tQBUbUw8aYZfnbrwLWoP1SMMaYycAjwKettd/N3633SsBG2S/OvFfUtZZzL/BqY8xTQAx4b8D1RNm/A983xjxJbqbQ+8h9i/wJkCDXN/9MgPVF2V8A/2aMKQVeAv7TWps2xqwEnib3Ze4jQRYYUR8C/o8xph84DPyptbZD+6Vg/j+gHvhrY8zQGKM/B/5F75VAjbRfPgn8kwvvlVg2mz37b4mIiIiIJ9TtKCIiIlJACl8iIiIiBaTwJSIiIlJACl8iIiIiBaTwJSIiIlJAWmpCRELLGPMN4BVAKTCf3AWNIXcx9qy19tse/I2lwNuttZ82xnwP+IK1do8x5ofAZ621B8b7N0QkWhS+RCS0rLUfATDGNAGPW2sX+/Bn7gbemL99E/DF/O3/nX/sbT78TREpYgpfIlJ0jDFfALDWfsEYcxi4D7iK3MKK3wU+BswA3mOtXWGMmQ98C2ggdy3Rj1pr1xljXgkcstYeM8Z8BpgG/NYYc521drMxpskYM89au6PgT1JEQktjvkSk2E0GHrTWLgHKgDdaa68DvgB8PP87PwA+Za1dCvwp8LP8/XcATwBYa/83cBB4nbW2Nf/4k8BthXgSIlI81PIlIlHwYP7/e8gFpqHb9caYKuBK4HvGmKHfrzLGNAAXAI+eYbt78r8jIjJmCl8iUvSstf3Dfhw87eEE0Dt8vJgxZgZwjNz1RQfOsOlBctceFREZM3U7ikikWWvbgW3GmHcBGGNeTb6rEdgONA379UH+8EtrU/53RETGTOFLRATeCXzAGLMB+DK5pSWywK/JzXAc8gC5Afdz8j/fkP8dEZExi2Wz2aBrEBFxljFmFXCntbbltPsXAZ+z1r41mMpEJKzU8iUicmYfBz49wv2fAv6iwLWISBFQy5eIiIhIAanlS0RERKSAFL5ERERECkjhS0RERKSAFL5ERERECkjhS0RERKSAFL5ERERECuj/B4MND7Iu8j05AAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plt.figure(figsize = [10,5])\n", "plt.plot(tnow,ixnow)\n", "plt.xlabel(\"Time(t)\")\n", "plt.ylabel(\"Current\")\n", "plt.savefig(\"plot/current.png\")\n", "plt.savefig(\"plot/current.pdf\")\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAmMAAAFACAYAAADwGHt4AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzs3Xd8XHeV8P/Pna4y6s2SVS35yr3FcYmd4vSQhBASCIQOWVjKspvfwy7P/vbZwi4suywLy5OlJJQlCRASQhISQqpT3ONuWWVUrGp1adQ1mnafP0Yjy7baSNMknzfkZWnmzr1f3ytrzpzvueeraJqGEEIIIYSIDF2kByCEEEIIcSWTYEwIIYQQIoIkGBNCCCGEiCAJxoQQQgghIkiCMSGEEEKICJJgTAghhBAiggyRHsB8dXUNhrwnR3JyLHb7SKgPIwIk1yX6yDWJTnJdoo9ck+gUjuuSnm5VpntOMmMzMBj0kR6CmIJcl+gj1yQ6yXWJPnJNolOkr4sEY0IIIYQQESTBmBBCCCFEBEkwJoQQQggRQRKMCSGEEEJEkARjQgghhBARJMGYEEIIIUQELdo+Y5Hy8ssv8tJLL5CVtQyAm2++leTkFB555Ps88sijER6dEEIIIRabsAZjqqpuA/7NZrNdf8njdwF/D7iBn9tstsfme4yn99ZytKpzQePcWprBh/YUT/v8vffez0033QrAuXO17Nv3Dnq9xLVCCCGECFzYIghVVf8a+DgwfMnjRuB7wNbx5w6oqvqizWZrD9fYAvX8889y5MghAP7mb/6OT36ymNOnT0V4VEIIIYRYjMKZzqkD7gWeuOTxVUCtzWazA6iquh/YDTwzn4N8aE/xjFmtYLjnng9OZMaEEEIIMTuv5sXtdePyunF5Xbi9bpwe18RjHs2NR/PinfSfR/OiTfypjT/vwatpvm3wPa6hoWm+VRI1NHz/9z/uf1Qb33ZiK9A0FEXH+2KvR8EUqVMTvmDMZrM9q6pqwRRPJQD9k74fBBJn219ycmxYli9IT7de9L3VasFsNl/2uMmkv+wxETpyrqOPXJPoJNcl+izma6JpGiOuUfodA/Q5BuhzDNLvGKB/bIChsRFGXKOT/nMwOv712HjQFa2SrHF8YPVtETt+NBQ6DQCTfzKtQN9sLwrHQqvp6Va6ugYvemz37psBLnv829/+/mWPidCY6rqIyJJrEp3kukSfxXJNHG4H54faaR1uo2u0h57RXrrH/3N4HLO+XkHBYjBj0VtINCVi0psw6gwYdUbfn3ojhsnf64zoFR06Re/7U6dDh4JO0aNTdOgUZeI5RVHG/9ShU3QoKCiAovi/uvC1goLv/+P/U/yjG/9TUdArerYUlYb8uswUhEdDMFYJlKiqmgIMAdcC/xHZIQkhhBBXBo/XQ8tQKzV956jvb6RlsJVuR+9l25l0RtJiUkmxJJFgspJgsmI1WbGa4kkwxRNnjCPGYMFisGDWm9Api6d7ll4X2YXCIxaMqar6USDeZrM9qqrqw8Cr+Pqe/dxms52P1LiEEEKIpa5vrJ+z3ZWUdVdQ3XcOp8c58Vy8MY7S5BKy47PIiV9GRmw6aTEpWI3xKL7UkgiysAZjNputAdg+/vWvJz3+IvBiOMcihBBCXElG3aMc7zjN4bbj1A80TjyeFZtBcVIhxUlFFCcVkmROlKArzKJhmlIIIYQQIdI+3MkbTe9wrOMULq8LBYWVSStYn76GdWmrSYtJifQQr3gSjAkhhBBL0PmhNv5Y/zpnusrR0EiLSWXHsq1sy9pMsiUp0sMTk0gwFqBLl0Pavv0aDh3aj06nsGHDZu6++wMRHuH8aJrGSPlZXL09JO6+TlLUQgixSPWPDfLSuVc51HYUDY38hFxuzb+BdWmrF1VR/ZVkyQVjv699iZOdZQvax6aMddxbfOe0z09eDumNN17l7//+n/F6vfz1X//VogvGvGNjDBw+SN8br+NsawUgtmQlpmXZER6ZEEKIQGiaxuH24/yu+g84PA6y4jK5t/hOVqeslA/YUW7JBWPhcOlySB6Ph0ce+T4PPPDRCI9s7rwuJ/3vvE3vH1/CMzgAev2k51wRHJkQQohADTqH+FXVM5R1V2LRm/nwynu4JntbxFs2iLlZcsHYvcV3zpjVCobJyyF1d3fzwx/+Fx/5yMcoKVFDetxg0LxeBvbvo+fF53Hb7ShmCyl33EnSnhvpfeVl+t54PdJDFEIIEYCmgRYeLXsc+1gfanIxD5beT2pMcqSHJQKw5IKxcPvGN/4Og8HIb37zJCkpqXz5y38Z6SFNy9FQT8eTjzPWUI9iMpF82x2k3Ho7equ/K/B4Gnt8fS8hhLiS+dc6jOYpvmMdp3ii8mk8Xg93Fd3GLfnXS13YIiTBWIDuuOOui77/wQ9+HKGRzJ13bIzu3/+Ovr1vgKZh3bad9Ps/jCHpkk9O/l84EouJCBt0DhFnjI3aN5XKnmrqBxq5veCmqHyjHnAO8qvK33Htiq2siV8b6eFM6a3m/bzdcoCvbvozUizRl8XpG+vnh6d/TnpMGg+t+3ikhzOld1sO8XT185j1Zh5a/3HWpq2K9JDEPEkwtsQ5Gupp/+mjONvbMGUtI+NjnyC2dOp/sBfeUiQaE5HzbstBflv9PB8suYs9ubsjPZzLnOws46dnnwBgQ/pacuKXRXhEF+sfG+T7J39E50g3zcPn+dbO6AvGXm3Yyx/OvQJAeU8Vu3N2RHhEF7M7+vj+iR/T7ejl/FAbLo8Lo94Y6WFd5I2md3iu9o/EG+P48saHyLXKTVeLWXR+7BQLpmka9tdfpelf/wVnextJN91C3j/807SBGDApMybBmIiMI23H+W318wCc6jwb4dFcrqLHxs/LfzXxfX1/4wxbh9+Ia5T/Pv1TOke6ARgcG2Js0jI30eDdlkMTgRjAuSg7h0OuYR459dOL1mZsGoyuFfoOtr7Hc7V/JMmcyMNbviiB2BIgwdgS5HU6af/Zo3T99jfo4+PJefhrZDzwUXRG08wvlFhMRFB5j40nq54hxhADQOtwO17NG+FRXdAw0MRjZ59Ap+j48EpfC5toCiScHhc/Kfsfzg+1sTtnB3tyd+PVvDQNNEd6aBOOd5zm6ernsRrj+fvtXyPGEBNVAe2Yx8mPTv+C9pFO9uTu5jNrfHfIn+tviOzAJjnddZZfVz1LnDGWr2z8HJmx6ZEekggCCcaWGPfgAM3//q8MHj6EpbCIvL/7R+JWr5njq/0TlRKNLWXlPVUcaTse6WFcpGmwhZ+WPY5e0fGF9Z/i6qzNjLpH6RjpivTQAOh12PnR6V/g8rj4zJoH2ZWzDYvectH6fpGkaRq/rnqW2r56NmWs50Mr309RYgEQPQFj00ALj1f+FrPezJc2fpbM2HQKE/LoGu1h0DkU6eGhaRq/qnyGhoEmrs7azAeK30dhYj4A9QNNER6dT+tQO/9T8RRGvZEvbvgMWXGZkR6SCBKpGQvQpR34CwoKqampJi4ujhtuuIlt2yJX++Dq6aHlP7+Dq6OdhB3XkPGJT86eDZssCguRRXANOAf54emfA756J4vBHOER+aaFHit7ApfXzUPrPk5xUiFtw+28136C+v5GlkX4DcflcfHTsicZcg3z4ZX3sCHd9+GmMDGPyt5qhpzDxJviIjrGd1oOcrTjBAUJeXxy9QPoFB2FiXkAUREwDjmHebTscTxeDw+t/zi51hzAdw4rem00DDSxLm11RMe4t3kfxztPU5RYwIOl96FTdCSbk0g0JVDf34imaRG9WWPENcqjZb/E6XHyubUfpyAhL2JjEcG35IKxrmeeYvDY0QXtw3rVVtLvf2Da5yd34D927D3uueeDeDxeHnnkexELxpzt7bR8999x23tJvu0O0j54//x/ccg85ZLk8rp5rOyJie8bBpooTSmJ4IjAq3n5ZflT9Drs3FF4MxvSfcXmk7M6O7OvjuAI4ZmaP9A42My2rC0XFZoXJviCsfqBxogGEuf6G3i29kWsxng+t/ZjGHW+X+tJ5kTSY1Oo72+KaCDh1bz8T8VvsI/1cWfhLRfd8efPPJ3rj+w5rO2r5/m6l0k0Wfnc2o9hGD+HiqJQlJjPya4yehz2iC2orWkaT1Y+TddoD7fk38CmjHURGYcIHZmmnIfnn3+Wb37zH/nmN/+RjRs309XVxde//nDEAjFXd9dEIJb2wQ+Rft+HovJ2exE5mqbxtO05zvU3YNL7sqW1ffURHpXvjrCKXhurU1VuL7hx4vFlcZlY9OaIT7Gd7jrLgdYjLI/P5gH13ov+XfkDxvr+yE1hOdxj/LL8KTRN4zNrH7xs8eeStCKGXMN0jXZHaITw7vlDVPZWsya1lFsL9lz0XEFCHgpKROvGHO4xnqj47fg5/BiJ5oSLnr8QMDZEYHQ+h9uPc7q7nJKkIu4qujVi4xChs+QyY+n3PzBjVisYJnfgP336FKtWreYnP/kFX/nK57nllttDeuxLufv6aPnudyYCsZTb75j/zuRuyiXrUNtRDrYdJdeawxfWf4q/O/At6iIcjLUOtfPHc6+RaLJOTK356RQdBQl5VNlrGHINE28M/zTgoHOIX1c9i0Fn4NNrPoLpktYGBYm5KCjU9p0L+9j8nqv7I92OXm7Ou56VySsue740bQUHm45R03eOjAgUeneMdPF87cvEGWN5sPT+y/rGxRgs5MQvo36gCafHddk5Dofnal+i29HLLfk3UJxUeNnz/sdq7Oe4OmtzuIdHz6id31W/gEVv5uOrPhy1vffEwshVXaChoUG++c1/4N/+7V/Yvn1nWI/tdTg4/1/fxdXVScqddy0sEGNSl2mJxZaUjpEunql+gRhDDA+t/QRJ5kSWxWVSP9CE2+uOyJg8Xg+PVzyFW/PwkdIPThlsFScVAVBtrwv38AD4re05hlzD3F1025SF0jGGGPKsy6kfaIpI+4iq3hr2nz9MdlwW7yu6Zcpt1mb6lmiz9daGc2iAb3ryycqncXldPKDeS6LZOuV2anIxbq87IpknW28t+1uPkBO/jDsKb55ym1xrDjEGC9X28J9DTdP4je1ZHJ4x7lv5flniaAlbcpmxULu0A/811+zmmmvC35hS83pp+9mjjDU3k3jt9aS+/95g7j2I+xKR5Pa6+Z/yX+P0uvjMqg9N/DIvTiqkdbid5sHzE9Mw4fRm87s0D7WyPeuqaWuF1JRiXqp/FZu9ls0Z68M6vvIeGye7yliRWMANubum3U5NKaZxsJnavnrWpIZvbVq3183T1S+goPDxVR+aqBO7VI41i0STlWp7Xdjrxo60n+BcfyMb09fNeP1WJq/gzeZ3sdlrw1rD6DuHz6Og8LHS+6c9hzpFR0nSCs50l9Mz2ktqGOvGTneXU9lbTWlyCduztoTtuCL8JDO2SPW88BzDJ08Qo5aS8dGPBeeX7Pg+NJmmXDJeb3ybpsHzbMvawpbMDROPrxifeolE3VjfWD9/aniTeGMcHyy5c9rt8q3LsejNVIc5q+Pyuvld9Qu+fmLqB2acFlKTiwGw2WvCNTwA3m45QMdIJ7tztpOXsHza7RRFYWVyCYOuIdqGO8I2vhHXKC/UvoxJZ+S+krtm3LY4qQidosMW5szT2y0HaB/p5JqcbTOeQ5h8ncM3RqfHye+q/4Be0fOhle+XOuAlToKxRWi47Ay9f3wRY3o62X/+ZRRDkBKc8m99Sekc6eKVxr0kmqzcv/Lui57z1xdV9YY3iAB4vvZlnB4ndxfdRqwxdtrt9Do9xUlFdI52Y3f0hW18bzXvo3O0m905O2Zd6qgoMR+Dog9rwNg/NsDL9a8TZ4zlzjkUc6sp4Q8kXm54nUHXELcV3HjZTQWXshjMFCTk0TTQwqh7NCzjG3AOTpzDu4tum3X7SJzD15vewT7Wx57c3WTGZYTtuCIyJBhbZNz9fbT//DHQ61n2hS+hj48P4t6lgH+p0DSNp2zP4fa6uW/l+ye62vslmKwsj8+mtr8eZxjrnc71N3K04yR51hx2ZG+ddXv/m2BVmN4Eh1zDvNqwl3hjHHdOU0M0mUlvojAxn5ahNoZcw2EYIbzS8CZjHid3Fd1K3AzBrJ86HniHK3vX67Czr+UQqZYU9uRdO6fXqMnFaGhU28NzM8SrDXsZ8zh5X+EtczqHWbEZJJqs2Oy1YVkVYsg5zJtN72A1xnPbpLuMxdIlwdgionm9tP/sMTyDg6Tf9yEs+QXBPYDcTblknOg8jc1ey5rUUjalT92TaFXKStxeNzVhvBvwxXOvAvDBkrvndFeYf3qoqrc6pOPye6PxHRyeMW4t2DNj1m6y0pQSNDSqekI/xp7RXg60vkdaTCo7l82t/1qKJZnM2HRs9jpcHleIRwh/qn8Tt+bhfYU3T1uHdalVKSsB3+oQodbrsLP//GFSLclcM8cedoqiUJqykkHnEM1hWKfytaa3GPM4ubVgT1Q0ZhahJ8HYItK/7x1GKsqJW7eepJumvntKCJfXzQt1r6BX9NxfMn2tyepU3xtgRY8tLOOqttdRba9lVcrKKVsITCU7LoskcyIVPTY8Xk9Ix9c/NsjbLQdIMieyO3v7nF/nvwGhrKcyVEOb8HLDG3jGAx29Tj/n161NW4XT46S6L7R3pnaOdHG4/RhZsRlszdo059cVJuYRb4zjbHdFyDNP/mDx9sKbJ5q7zsXEde6uCNXQAF9N5bstB0kyJ7Ire1tIjyWihwRjAXr55Rf54hc/xze+8X/4xjf+D4cO7cftdvPlL/8ZZ8+Whey4Lrud7t89jS4mhsxPfjokxZyKZMaWhH3nD9Hj6OXa5TtIj02ddruixAJMehOVYcg6aZrGS+NZsTunacMwFUVRWJ+2mhH3aMhbH7ze+BYur4vbCvZgDKDfVXZcFimWZMp7qkIaMHaN9HCk7TjL4jK5KnNjQK9dn+ZbwulMiAOJVxr24tW8vK/oloD6YekUHWtSS+l3DoY089TrsHO4/RiZselcnTn3YBFgVUoJBkUf8nP4ZtO7uLxu7ii4KaCfQ7G4LbnWFgf31nGuqnNB+ygqzWDnnssbKPpNXg4J4Mc/foTMzNCtn6dpGp2/ehzv6CgZn/gUhqQQ9ZqRu3UWvRHXKK/Uv0mMwTJrrYlBZ0BNXkFZdyWdI91kxKaFbFw1fXXU9TewLm1VwGvqrUtbzbvnD3Gmu4KSKRqbBsOwa4QDrUdINiexY9nstWyTKYrCurRVvNNykNq++ok6t2Db27wPDY3bCm4MuPFnYUIeccZYyroqeGDlB0LyYa5vrJ9jHafIis1g4/iyVoFYn7aaI+3HOdNdQX5CbtDHB/BW8368mpdb8m8IKLMIYDFYWJlcTEWvjZ5Re0h6fo24RjnQeoREUwLblkkriyuJZMbmYfJySL/73VNs2rSF7OyZb41eiOGy0wyfOknMSpXEXXMriF0IaW2xeL3Vsp9h9wi35N0wp67169N8b5qnu86GdFxvNr0LwK35e2bZ8nIlySuw6M2c6a4I2c/mvvOHcXpd3JC7K6CpK79QT2ENuYY51HaUFEvytDWAM9Hr9OOZp4GQZZ7ebj6AR/OwJ2/3vLrEl6asxKDoQ3YOJwc6gWYW/UJ9nQ+0HmHM4+T63Gvm9XMoFq8ld7V37lkxY1YrGCYvh/TVr36R5cuXU1VVSV1dLd/61neCeizN7abr6adAUch48BMouhDGz5IZW9QcbgdvN+8nzhDLtcvnthrE+vTV/Mam42RXGTfnXx+ScbUPd3K2p4qixPx5NZg16gysSlU52XmGtuEOsuOzgjo+l9fNOy0HsOgt816UvCSpCIvewqmus9xbcmfQl6zZ13IY13iwGGhGx2992hreaz/Bya6yWftqBcrhdrC/9TBWYzxXZ85vySCLwUxpSglne6roGO4MejsHf6BzW/6N8w501qev5unq5znReYbrc68J6vjcXjdvtxzArDexK4CaRbE0SGZsgf7rv37I1772t+zcuYuPfvQTQd9/37tv42pvJ/Ha6zHn5AR9/1OSzNiitL/1CCPuUa7PvWbOd2DFG+NYmbSCxoHmkPXy2tu8D4A9ufPP6vqnvY53nArKmCY71nGKAecg12RfTYzBMq99GHQGNqSvwT7WF/TFzV1eN++cP0CMwcLOAKdQJ1uTWopFb+Zo+8mgF8kfbjvOqNvBdcuvWVCd05bxjNXRIF9nr+blnZaDmPQmduXMvyg+yZxIcVIhdf319IzagzhCONVZRt9YPzuXXU2sMWb2F4glRYKxAN1xx10X1Yv5ffazn2ft2sCnD2biGRmm5w/Po4uJIfWeDwR132JpcXlcvNn0Lma9ieuXB/aJfWOGL9A5FYKpymHXCO+1HyfVksKG9DXz3s/6tNWY9SaOdpwK+lTlvpZDKCgLznT47x482nEyGMOacKarnEHnEDuWbcUyz2ARwKQ3sjF9XdADRk3T2N96GL2iX1CgA77snUln5GjHyaBe54oeG/axPq7O3DTnliXT8V/nYH8w2N96BIBrl+8I6n7F4iDBWBSzv/4a3qEhUu64E4M1IfQHlLspF62j49md3Tk7An6z2ZC+FgUl6EEEwHvtJ3B53Vy7fMeCpu5MehMb0tfS4+ilfiB4gUTLYCuNg82sSS0lxbKwgmw1uZgEk5WTHWeCugD7wdb3AObcE2smV2X5M0/Bu9b1A020DXewPn0NVtPCmlBbDGbWp6+he7SHhoHmII0QDkycw4W3itiUvg6Dog/qOewY6aKm7xwlSUVkxKYHbb9i8ZBgLEp5Robpe+M19FYrSXtuCssxL9xhJcHYYqJpGu+0HECn6ALOioGvG/+q1JU0DjQHdf1CTdM42PoeekXPtiAscrx1vBXB0fbgvQkeCGKgo1N0bMncwLB7JGi927pHe6iy17AisYCsuIXfsT05YAxWA9gD4xmdYJxDuHCdj7QfD8r++sb6OdtTSW58dlBq5WKNsaxJW0XrcDtNgy1BGOHkgFv6il2pJBiLUn1vvI53dJTkW25HZw5zB2aJxRaV+oFGWoZaWZ+2etZ1AKfjb+dwqO1o0MbVMNBM63A769NWLzhjAhcCiaMdp4KyhJPT4+RoxwkSTQmsSS1d8P4AtmVdBVwIUBbqYKvvegTrTVqn6NiWtYVh9wgnuxbeF3HU7eBEx2lSLSkTqyUs1KqUlSSZEznafgKH27Hg/R1uO4ZX83LNAqdQJ/PX7u1rObzgfbm9bg63HSPOEDuvliBiaZBgLAp5Rkawv/Ea+ngrSTcE3gpg3sYzY9LaYnF5p+UgANfN8Q7KqaxLW02cIZb32k4ErXFpsD/t63V6di7byqh7lGMdpxe8v1NdZxl1O9ix7Kp536F4qVxrNoUJ+ZT32Oge7V3QvryalyPtx4kxWNiUEbx61F0521FQeHf852YhTnaewel1sTN7a9DuINXr9OzK3obDM7bgqUBN0zjSfhyjzjjvdhZTWZ2qkmpJ4WjHSUZcC1vcvLK3miHXMFuzNkmT1yuYBGNRaGD/u3hHRki6+RZ0lvkX7AZOpikXmwHnICc7y8iKy6Qkaf4tXYw6A1dlbWLQNRSUDuNOj5PjnadINicFtQnqRCBx/uCCPzT4pzu3LbsqGEObcO3yHWho7D+/sKxJXV89fWP9bEpfj0lvCtLoIC0mhTWpKvUDTQueZjs+HhRvDbCb/Wx2Zl+NTtHxbsuhBV3n5qHzdI50sy5tFTGG4N2hqFN07M7ZjsvrWvB06rHxGwECWT5KLD0SjEUZzePB/ubrKCYTSdfdEN6DSyy2YEOuYR4re5zzQ21hOd577SfwaB52Z29fcFf1a3N8vY3eGm9FsRBne6oY8zjZmrUpqD23ki1JrE9bTfPgeRoGmua9nyHnMFX2GvKsy4O+8sCm9HXEG+M42PYeYwuYTj3W6Qt0gpnR8dud47tj7+3mA/PeR//YIDZ7LYUJeaTGpARraAAkmhPYmL6W1uF2bPbaee/HHyxuCcE53JG9FaPOwFvN++edTXZ6nJzpriDNkkK+NTSrDojFQYKxKDN06iTunh4SduxEH7/wOpuAyN2UC1beXcWprrN8673vhfxYmqZxpO04ekU/cZfcQmTFZbI6VaWuv4HGBd7J5n8TDEUg4W9o+0bTO/Pex8muMryaly2ZG4I1rAlGvZHdOdsnlliaD4/Xw8nOM1hN8ZQkFwV5hL5ptqzYDI52nKRnntOpJ7vOoKGFJNABuCnvOgBebdg7r9d7NS/HO05j0VtYk6IGc2iAr0ffjmVX0+PonchuBaqsuxKnx8mWzI0hWaJKLB4SjEWZvjdfByDpxpsjPBIxH5Onk0ZcIyE9VstQK63D7axNWzWnpY/mYk/ubuBCo9b5cLgdlPdUkhWbQXZccLvlg6+QPz8hl1NdZ2kdap/XPvw9orZkBD8YA7g+dxcmvYk3Gt+Z112LVfZahl0jbM7YEPRu/uCbZru1YA9ezctrjW/Nax/HO06joLA5Y32QR+eTn5DL6hSV6r466voaAn59fX8T9rE+NqSvCVkt1s3516FX9LzauHdejXSPhzD7KRYXCcaiyFhzM6PVNmLXrMWcHaZu+5MokhlbMG3SHO/xzjMhPZa/ViUYbSP8SpNLyIlfxvGO07QPd85rH2e6K3B53WzJ3BCST/uKonD7+CLorzS8GfDr+8b6qe2rZ0Vi4bzvPp1NvDGOa3N20O8c4FDbsYBf7w8WQ/kmvSVjA2kxqRxuO0avI7Bu8nZHH+f6GyhJKiLRHLoeiP7F7l+ufz3g2rETnaGbovRLsSSzLWsLHSNdAWfHHO4xynuqyIrLDPoSX2LxkWAsivQf8C2mnHR9mGvFJkjR2EJN/nT8XpD6JE3F4/VwtP0kccZY1qQGbwpGURTuKLwZDY2X61+f1z4m3gRDlHUCWJu6iuXx2ZzoPEPLYGtArz3TVY6GxubM0GR0/PbkXotJZ+TlhtcDatHg8Xo421NJkjmRwoS8kI1Pr9Nze8GNuDUPf6h7JaDX+hfK3hjEuzynsiKpgFUpK6my11DeUzXn12maRll3BRa9BTU5tGsV31awB4POwB/qXgmo5UqVvQa31y3tLAQgwVjU8LpcDBw+hN6aQNy60L2JzUgSYws2+dP7uf5GOke6QnIcm72WIdcwWzI2znvR4+lsSFsqf8rQAAAgAElEQVRDnjWH452nAw50xjxOKntryI7LCvpCz5MpisLdK25HQ+PZmhcDypqUdVcCviWWQinRbOWm/OsZdA7xagBTgfUDTQy7RlibWhryOqKrszaTa83haMdJ6vvnfkOE/xyuS1sVqqFNuLf4ThQUfl/70pwL5duGO+hx2FmTqgb938elUmNS2JO7G/tYH282vTvn15V1+QLacJxDEf0kGIsSw6dP4R0aImHHThRDaH95TEsKSBfMnxkrSfIVXb/XfiIkxzk13rAzFPU6iqJwV9FtAAEHOlW9vk/7a8PwBrMmVWVNainVfXVzbsfhcI9Rba8lJ37Zgpc/moub864jyZzI3uZ9dI/2zOk1/qzTuhAHi+CrHbuv5G4Anql5YU51Tw63I6znMDs+i1052+kY6ZpzLaP/5yEcP4cAt+TfgNUUz2uNb9E1Mvt19mpezvZUkmCykmdd+KoAYvELWzCmqqpOVdUfq6p6SFXVt1VVLb7k+f+lqupxVVWPqqp6xa2K3b/f94kqYdfuCI8ESY0tgHd8indTxnrMehNH2k/Mq7B3xmNoXk53lWM1xrMiqSCo+/ZblbKStamrqO6rC2hB5LNhDCTAlzXRKTqeqX6B0TlMBVb1VuPWPGEbn0lv4gMr7sDtdfNk5TNz+lko667EpDMGraP9bIqTCrkqcyONA81zyuxU9dbg1jwhzyxOdmfhLViN8bxU/9qcahnPdlegU3RBW1lhNjEGC/cV34XT6+LJqqdnvc4NA80MuYZZm7oqJDdoiMUnnD8F9wAWm822A/g68F3/E6qqJgF/AewAbgG+H8ZxRZy7r4+R8rNYiooiUrh/gdSMLZQ/i2TRm9mUsZ5eh51qe11Qj1HbV8+Qa5gN6WtC9otcURTuX/l+jDoDz9a+NKcu417NS1lPJfHGOAoSwtMzKSsug1vyb8A+1sdztX+cdftwTVFOtiVzI+vSVlPTd459szSC7RzpomOkk9KUlWHtxn5/yfsnBTszr096JswBN0C8KY4H1A/g9rp5ovLpGacrB5yDNAw0syKxgDhjbNjGuCVzIxvS11LbV8/bzftn3PZC9lOmKIVPOIOxXcArADab7TAwue31MNAIxI3/F9xUQpQbPH4UNA3r9vkvZxMUEostmDb+iVhRFHZl+5qoBlIvNBf+KcpQF0+nxaRwW8FNDDgHecr2+1mnK5sGWxh0DoX90/7tBTeSE7+MA61HODsebE1l8tRQrjV8H3oUReEj6r3EGmJ4rvaPMzYEPhvGWqzJ4k1xPFB6L26vm5+d/RUO99iU23k1L+U9VWE/h+D7eb8qcyMNA008X/fytNuV99jQ0MI2RemnKAoPqB8g3hjHc3Uvc66/Ydptz3ZXYtQZKE0pCd8ARVQLZzCWAPRP+t6jqurk4qhmoAI4AfwgjOOKuKFjR0FRsG4O7rIsAVMkGluotvGsQozBQmFiHqtSVlJtr6W2rz4o+/dqXk51niXOEMvKBSx/NFc3511HUWI+xztPz7rsS6QCCYPOwMdXfRiDzsD/VDw1bc1O8+D58amh0rBPDSWaE/jYqg/h8rp49MwvGZ6mB11lbw3ga8oabhvT13Ld8p20Drfzq6pnpgy+zw+1M+QaZnWKGpHptQfUe8mMzWBv8z6OtU+9bqXNfw5D0Oh1NgkmK59Z8yCapvFY2RP0jfVftk2fY4DW4XZWJBYGdZkrsbiFs1J8ALBO+l5ns9nc41/fDiwDCse/f1VV1QM2m+296XaWnByLwRCcxX1nkp5unX2jBRjr6WG0ppqEtWtYVhLZ5TDc8Ra6gASrhbQQ/70XKtTXZT4GHIMcbHuP1Jhkrl25BYPewEc23s3f7/0PXmt5k+3Ff7ngu+PO9TbS7xzg2oJtZGWGpkfWpR7e9Tm+9to3ebr6edYsX0FRytTtFmoH6tApOq5ZuYlYY/DWAZyL9HSVh/gIPzr6BD+vfJJv3Pj/XTaG/V2+uwWvLlgfkZ+fm9K30+3p4PcVr/C47Td8/dovYZo0FenyuKjrryc3YRkly4NX1B3I3/XzKR+h4+0OTnSeIS91GR9df89Fzx/q8Z3DrfnrIvRv0MrXr/tz/vb1f+PJqmfIzchgbeaFujBN06juryPRksD6guKIdLVPT9+EXbuXJ04/y4/P/oJ/uuFh4s0XmjLvbzwKwJbctVH5e+xKFsnrEc5g7ABwF/C0qqrbgbJJz9mBUWDMZrNpqqr2ATO+09jtoe1uDr4L09U1GNJj2F/3TWFZNmwJ+bFmMzTsm5oY6B9Fi/BYZhKO6zIfL9T9iTGPk7tzb8fe66uxSiWDNamllHdW8UbFoQVPLe6v92UDiuOKw3YOFMx8vPTDPFb2OP/6zn/z11u/QpI58aJtYhP11PY2UpiYx3Cfm2HCf33WWtdxbc4O3j1/iG/ufYQvbvjsRcHO8ZZyALL0ORH7+bkh83rqOps53VnOd97+CZ9d+zH0Ot+Hyhp7HWMeJ8UJK4I2vvn8W/lU6YP85/Ef8nzlq2hOHbfkX+h7eLz5LADLDMsjdg5NxPG5tR/nR6d/zrf3/YivbHyIosR8AFqH2ul3DLA1cxPd3UMRGR/AtpSraVrezjstB/invf/Flzd8buLDQVmHr1/acnNuVP4eu1KF431lpmAvnHnm5wCHqqoHge8Bf6Wq6sOqqt5ts9n2AUeBw6qqHgKqgfl1nFxkBsenKOM3B6+L+vz5PkVqMk0ZsBHXCO+2HMRqimfnsqsveu6DxXeiV/Q8W/sSznksjTNZRW8VCgqrwlxrsiF9DfcU30G/c4BHTv2UQefFb3QVXTVoaGG7A3A695Xczcb0ddT0neOxsscnmnA6PS7O9TewPD4bqynMa75OolN0fHrNR1mZXMzp7nJ+Vv6riZ+JqvEFsdWUyJ5Dqymer2x6iGRzEi/U/Yk/1L2Cpmm4vG5q++pZFpdJojmyGZ3SlBI+vfZB3F43//fkoxNT5FV23xSlGuFaLEVRuK/kLrZlbaFxoJnvnfgRdkcfmqZxpqOSOEMsy+OzIzpGEV3Clhmz2Wxe4AuXPFw16fl/AP4hXOOJBu7+Phx1tcSopRgSE2d/QahJydi8vdNyCIdnjNsLb7ooGwOQGZfBDbm7eKPpHV489wofLLlrXscYdo1Q399EYWI+sWG8S8zvxtxr6Rvr563m/Xz/5E/4i40PTSyF4/+0ryZH9k1Qr9PzqTUf4dEyJxU9Nn5w8lG+sP7TtAy14va6Ix7ogG8h8c+v+wQ/OfNLTned5f+eeow/W/cJbL016BQdxUnBXxg8UCmWZL666fM8cvqnvNq4l15HH1dlbsDldVEa4WvstzF9LQ+t/Tg/L/81Pyn7JR8oft9EvVhphD8UgC/w/tiq+7EYLLzTcoD/OP7f3Fl4Cz0jdjZlrJeWFuIi8tMQQcNlvpna+A1RskisFPDPi9vr5t3zB7HoLezK3jblNncU3kxGTBpvNe+nejwDEqjK3mo0tKAufxQIRVH4YPFd3JC7i/bhDr5z7BGaxzv0n+2wYdIZKUwM3fI9c2XUGfj8uk+yNXMz9QNN/PuxH/Bqw14g8sGin8Vg4YsbP8uWjA2c62/g7w99m/qBJgoScokxWCI9PADSY1P5X1u+RH5CLkc7TvCjM78AIp+5m2x9+hq+uunPiDPE8mzNi5ztqSItJjVka44GSqfouL/kbj5Q/D76xwZ4suoZgIhnkEX0kWAsgobP+Jppxq2PkmBM1kOalxOdZxhwDrIzeyuWad5IzXoTn1zzAIqi8MuK3zLgDLw2oaLHBhC2RpZT8QdkdxXdhn2sj/88/t+81vAWLQNtrEgqDPnSM3Nl0Bn45OoPc0fBTfQ6+qju8/V6K04qnOWV4WPUGfjUmo9wz4o7JqZTo+1N2mqK5682/znXLb9m4rFoyNxNVpiYz/+++i8pTPDVja1Lja7eXYqicFPedfzl5i8QY7Bg0BlYnbIy0sMSUSY6fnNegbwuF8Pl5RgzMzFlZUV6OD6SGAuYpmm83XwABYXrls/cJ64gIY+7Cm/lhXN/4rGyx/mLTZ/HOMfgRdM0KnptJJisEa81URSF2wr2kBWXwRMVv+WFc38CiLqeSYqi8L6iW1iRVMgTlU+zIrEAc5S1EtApOm7Ov54VSQUcOP8eu3K2R3pIlzHqDHxo5ftZn7Yaj+aNmszdZInmBB7e8ufYemspHC/mjzbFSYX8886/xRivYXCE925jEf0kGIuQ0Wob2piDuHXXRnooExSJxgJWP9BE42Az69PWkBaTOuv2N+dfz/nhNo51nOLJyqf55OoH5lQ70jHSyaBziKsyN0bkdv2pbExfy/Krs/llxVM0D52PuoyEX2lKCf+y82+j+saUosQCihILIj2MGUVbsH0pnaJjVWp0Z5xiDBbSrVa6HHIXpbiYBGMRMnzmNBBF9WIgmbF5eKflAAA35F4zy5Y+iqLwYOn99DrsHOs4hUEx8OCq+2YNyKrt54ALC5BHi7SYFB7e/OdYk00M9S3sTtFQUhRl0ocNIYSILlIzFiHDZ8tQzBZiSqLok1yUZFwWi2HXCKe6zpIZm0FJAN3wTXojX9zwGfITcjncfoxflP961pYXtX3jwVhy6LvuB0pRFGKM0Td1JYQQi4UEYxHg6u3F1dFOrKqiGKIvOTnbGoTC52j7SdxeNzuztwY8dRhjiOHLGz5HUWIBJzrP8IOTj9I/NjDltpqmUd1XR6LJSkZMWjCGLoQQIopIMBYBI5UVAMSWro7wSC4l85RzpWkaB9veQ6fo2JY1v4a9scYY/mLjQ1yVuZH6gUa+9d73ONNVftl2nSNdDDqHKE4qipp6MSGEEMEjwVgEjFb5ukXHroqygmeJxeasefA854faWJe2ekEd3Y16I59a/RHuL3k/Ds8YPyn7Jb8o//VFCwxXR/EUpRBCiIWLvjmyJU7TNEaqKtDHWzHlBG8x4OCQaGyuDrb5FvvduWzrgvelKArX517DyuQVPFH5NMc6TnGmu4I9y3dxfe6uC/ViUVa8L4QQIjgkGAszV0c7brud+Ku2ouiiKzE5MQUmNWMzcnvdHO84RaLJyqogNm/Mjs/ia1d9mcNtx/lD3Z94pXEvbzbvw+V1YTXGkxmbHrRjCSGEiB4SjIXZSKV/ijLa6sWQxNgcVfXWMOIeZU/ubvQ6fVD3rVN07MzeylWZGzjYepTXGt+i3+miNGWl1IsJIcQSJcFYmI3YfAsqx6pRVi8G+KOxaG6OGQ2OdfiWsdqSuSFkxzDpTVyfew27c7ZT119PdvyykB1LCCFEZEkwFkaapuGoq0FvTcCYmRnp4VxOMmOzcnqcnOkuJ9WSQr41N+TH0+v0rIyy9QqFEEIEV3QVLS1x7t5e3HY7luLi6JxyUiQam83ZnirGPE62ZG6IzmsohBBi0ZFgLIxGa2sAiCmO1jXepIB/Nsc7fMtYXZUZRctYCSGEWNQkGAsjR914MLYiSqedJDE2ozGPk/KeSjJjM8iOy4r0cIQQQiwREoyF0WhtLYrBgDm/INJDmZIi0diMKnurcXndbExfK1OUQgghgkaCsTDxOhyMNTdhzi9AZzRGejhTWwR9xloGW/n/3/h3zvU3hP3YZV2+ZazWp0dhWxIhhBCLlgRj03C6PBytaMfrDU5g4qg/B5oWxfViXJimjGJH2o9T01PP87Uvh/W4Xs3L2Z5KEk1W8qzRtnKCEEKIxUyCsWkcrujgGz87QkVjb1D2N1pXC4AlWuvFJonevBj0jw0AUNffgBbGDN65/kaGXMOsTVuNTpF/NkIIIYJH3lWm4c+IDY24grI/R2MDAJaCwqDsLzSif5py2DUy8XWPwx62457pLgdgfZpMUQohhAguCcamYTT4To3T7Q3K/sYaG9BbEzAkJwdlfyGxCOr3XV73xNf1/Y1hOaamaZR1VWDSm1ClAasQQoggk2BsGmajb81Bp8uz4H25BwZw9/ZiKSiI8rvwoj8ac3kvZCrrB8ITjHWNdtM52s3qlJUY9VF684UQQohFS4KxafgzY64gZMbGxqcoo7WlxYRFcDely+vCbDBjUPRhy4xV9FYDsDpFDcvxhBBCXFlkbcppmMYzY2NByIwtjnqxxbEaksvjIs4YQ2JsAo2DLTg9Tkx6U0iPWTUejJWmRPGdsEIIIRYtyYxNwxTEzJhjsWTGFsU0pRuj3khhYj5ezUvjQEtIj+f2uqm215ERm0ZqTEpIjyWEEOLKJMHYNEwTNWPBmabUJyRgSEpa8L5Cajw1Fs6WEYFyeV2YxoMxCH3dWH1/E2MeJ6tkilIIIUSISDA2DdPE3ZQLm6a8ULxfGOXF+9HPq3lxel2YdEYKE/IAX7AUSv4pylUyRSmEECJEJBibhj8zttBpyrEmX+bGnJe/4DGFXJQX8Nf3N+H2uslNyibZkkSSOZH6gcaQZvIqe2vQK3pKklaE7BhCCCGubBKMTcN/N+VCC/jHWpoBMOfmLnhMIRflJWMnOk8DsDN3CwCFifkMOodC1vx1yDVM02ALRYn5WAzmkBxDCCGEkGBsGmZjcAr4x877CszNOYsgGIviaMyreTnZeYY4QyxrM0sBKJqYqgxN3VhtXz0amjR6FUIIEVISjE3DoNehKAtv+upsaUExmTBmZARpZKGjRPE0ZV1fA/3OQTakr8Wg800hh7qIv9Z+DoDipKKQ7F8IIYQACcampSgKJqN+QcshaR4PzrZWTMuyUXSL4FRH8f0F/inKzZnrJx5bbs0JafPX2r5zGHQGChIWQ1ZTCCHEYrUIIoTIMRkWFow5OzrQ3G7MOcuDOKorj2+Ksox4YxwrJxXSG3UGcq3LaRlqw+lxBvWYI65RWobaKEjIlSWQhBBChJQEYzMwm/QLmqZ0+ov3ly+WYCw6+4zV2M8x6BpiY/pa9ONTlH6FiXkhaf5a1++rFyuRKUohhBAhJsHYDMxG3YIyY/7ifdPyRTLNFaU1Y8fHpyi3ZG647LniJN8SUzZ7bVCPWdtXP75/CcaEEEKElgRjMzAbDbgW0PT1wp2UiyUzFn08Xg+nuspIMFmnDIxWJhejU3RU9NiCetyavnPoFN3ETQJCCCFEqEgwNgOTUbeg5ZCcLS3orVYMiYlBHFUIRWFmrNpex7BrhE0Z69Apl/+4xhgsrEgsoHGwmUHnUFCO6XA7aB48T751OeYQL0IuhBBCSDA2A7NJj8er4fYEHpB5HQ5c3V2YFlFWLBqXa/JPUW7OuHyK0m9Nqq/vWLCyYw0DzXg1r0xRCiGECAsJxmawkCWRnO3tAJizc4I6prCIksyYx+vhTFc5iaYEimaYLpwIxnqDE4z517uUKUohhBDhIMHYDMzjwdh8ivid7a0AmLKygjqmkFKi627K+oEmht0jrE9fM+UUpd+yuEySzUmU99hwe90LPm7DgC8YKxjv8C+EEEKEkiFcB1JVVQf8ENgAjAGfs9lstZOevx34h/FvTwBfstlsEY0K/Jmx+bS3cLa3+faxLDuoY7qSnOkuB2Bd2qoZt1MUhY0Za3mreT9VvTWsnWX7mWiaRsNAEymWZBLN1nnvRwghhJircGbG7gEsNpttB/B14Lv+J1RVtQLfAe602WzbgQYgLYxjm5LZtJDMmG+a0pi1LKhjCqkoK+Av667ApDdd1Oh1Ov6ashOdZxZ0zO7RXoZcwxRKVkwIIUSYhDMY2wW8AmCz2Q4DV016bidQBnxXVdV9QIfNZusK49imZF5IZqytDcVswZCUFOxhhY4SPQuFdwx30jnSzeqUlXPqgF+QkEuyOYkz3eW4FjBV6V/nsiBRgjEhhBDhEbZpSiAB6J/0vUdVVYPNZnPjy4LdAGwEhoB9qqoestls1dPtLDk5FoNBP93TQeEPxuLiLaSnz33KSvN4qOnsIDYvj4yMhFANL+j6O2NpAWJjTAH9fUPhUM9hAHYWbp5yLFM9dk3+Fl6qfpNWdzNX5ay/7Pm56Gj2ZTQ35ZWSnibTlIGI9M+MmJpcl+gj1yQ6RfK6hDMYGwAm/01144EYQA9w1GaztQOoqvouvsBs2mDMbh8J1Tgn+GvGOrsHybDOvd+Uq6sLzeVCl5ZBV9dgqIYXdKP9owCMDI9FfNyHGk6ioJBnKrhsLOnp1inHtyZhDS/xJq/Z9pFvKpzXcSvb69AreuLdSRE/B4vJdNdERJZcl+gj1yQ6heO6zBTshXOa8gBwB4CqqtvxTUv6HQfWqqqapqqqAdgOVIRxbFOaqBkLsPHrRPH+YrqTEvCvTRlpI64RzvU3UJiYh9UUP+fX5VpzyIlfxpnuCgacgf+jcnpcNA+dZ7k1WxYHF0IIETbhDMaeAxyqqh4Evgf8laqqD6uqevd4fdj/Bl4FjgC/t9lsZ8M4tilN3E0Z4JJIzjb/nZSLqHg/itjsdWhorE5RA3qdoihck70Nr+blSNvxgI/bMtSKV/NKSwshhBBhFbZpSpvN5gW+cMnDVZOefwp4KlzjmYsLBfwBZsY6/JmxRRaMTdxMGdkC/spe3+x0acrKgF+7NXMjz9W+xMHW97gx79oZ+5NdqnnwPAD51sWzaoIQQojFT5q+zsA8zw78zrY2UBSMmZmhGFboBLm1haZp9HYP4/XO/fxpmkZVbzUxhhjyEwIPimKNsWzO2EDnaDflPVWzv2ASfzCWa12EqyYIIYRYtCQYm8GFmrEApynb2zCmpqEzXrmLTGuaxr7Xa/jtT49SfrJ1zq/rGu2hx2FHTS4OKKs12Y151wLweuM7Ab2uZfA8Rp2BzNj0eR1XCCGEmI95B2OqqhpVVb1qvGHrkuTPjI0FEIx5HaN4BgYWX1YMgpoZO36gkfITviCstal/lq0vqJqYoiyZ97Fz4pexJrWUuv56zvU3zuk1bq+b1uEOsuOWodeFtmWKEEIIMdmcgzFVVfNVVX1DVdWrVVW1AIeA94AGVVW3hGyEEeTPjAUSjLm6fL1qjekZIRlTaAWn6WvVmTaO7m8gPsEMwEDf6Nxf21sDwKp51ItNdnPedQD8qeGNOW3fNtyBR/OQa5Xlq4QQQoRXIJmx7wMmoB34KFAMbAN+C/xH8IcWeTFm3/0NYwEU8Ds7OwEwpi++qa5gNOBvbe7j7T/ZMFsM3PnhDaRlxmPvGcHrnX2nHq8Hm72OtJhU0mJS5j8IoDipiJVJK6josVFtr5t1++ZBXxZvudSLCSGECLNAgrEbgC/abLYm4E7gZZvNdhT4T2BrKAYXaRbTeDDmnPvyOq5uX2bMlLEIM2MLnKYcHXHyxh987eFuu3ctyamxpKTF4XF755QdaxxsweFxLGiK0k9RFO4pvgOA5+tenvUO0ZYhf/G+ZMaEEEKEVyDBmAIMq6qqB/YAr48/HgOMBXtg0cBi9k9Tzj0z5prIjC3CYGx8mnI+oZimaex9qYrhQSdXX1tIdp5vTc6U9DgA7N3Ds+6jtu8cACuTiuYxgsvlJ+SyKWM9jQPNHO04OeO2zYPn0Sk6suMWWTsSIYQQi14gwdhh4OvAN4BY4CVVVXOAbwEHQzC2iJtXZqxrEQdjE9OUgYdjtrJ2ms71kluYzKbtF5qmJqfFAtDbPfvyVTXjwVhxkIIxgHtW3IFRZ+TZmhcZck0dEHo1Ly1DbWTFZmCSzvtCCCHCLJBg7CvAFuBLwMPjXfP/BigFHg7B2CLOaNCh1yk4Airg70SfmIjObA7hyEJknkVjI8NODu6tw2DUcd1tKopyYVmllDRfZqx3lsyYV/Nyrq+BjJg0Es3BW1w9LSaFO4tuYcg1zO9rXppym66RbpweJ8tlilIIIUQEzLkDv81mqwauuuThfwS+arPZItuyPYTMRj1jzrlNU2puN66eHiyFwcvsRESAV/PQ3jrGHG6uuakYa6LlouesiRYMRh32rpmDsZbBVhyeMTYnbQh0tLO6YfkujrWf5Ej7cdamrWJzxvqLnm8e8hXvS7NXIYQQkRBQnzFVVZNVVf0bVVV/oapqBr7asYX1IIhyZpOeMdfcpildvb3g9WJcjMX7cCGjFcA0ZVf7INXlHaRlxrN28+XBjKIoJKfGYe8dmbETf+3EFGVhYIOeA71Oz6fWfBSTzsivKn9H92jPRc9PdN6Pl8yYEEKI8Aukz9hKfGtJfgZ4EIgHPggcU1V1Z2iGF3lmo37OBfz+ejHTYqwXAy4Ujc2NpmkcesvXNmLHDSvQ6aZ+fUpaLF6PRr99+jsqa/vqAShJDk1WMSsugwfUe3F4HPzkzC8ZdV8YS8t4W4scCcaEEEJEQCCZse8Bv7PZbCoX7p58EHgG+HawBxYtzCY9Y8651Ywt7jspJ5tbZqylwc75xj5yC5NZXpA87Xb+Oyp7u6Yu4vdqXmr76kmxJJNimX4/C7Vt2RauzdlB63A7Py17Eo/Xg6ZpNA+dJ82SQqwxJmTHFkIIIaYTSDC2Hfi/kx+w2WxefIHYpmAOKpr4MmMevHOYupu4k3KRTlMS4CzlsQO+pYa2XTdzNis5beb2Fu3DnQy7R0IyRXmp+0ruZl3aKqrsNfy8/Nd0j/Yy7BqRZq9CCCEiZs4F/PjSJVOlDjJYon3GACzjSyK5XN6J5ZGm41zMbS0goLsp25r7aG/pJ29FCulZMy9POtsdlXX9vinK4sTQB2N6nZ5Pr3mQH53+Oae6yjjVVQZI8b4QQojICSQz9gfgX1RVjR//XlNVtQjfMkl/DPrIooRpfLHwubS3cHV1oZgt6K2LdO10fzA2h6WLTh5uBmDzpJ5i04lPMGM06acNxur7mwAoTMyf40AXxqw38cUNn7lo/cvChNn/HkIIIUQoBBKMPQykAL1AHL5FwmsAJ/C14A8tOliMc18s3N3TjTEt7aI+W4uJMsfMmL17mMa6HrKWJ7AsN2lO+01Oi6W/dxSP5/KbIRoGmrDozWTFhS+jaNKb+Pz6T3FN9tWUJBVRlFQQtmMLIYQQkwXSZ6wP2KGq6k3ARnxBWLnNZnszVIOLBv6pydmK+D0jw3hHRzGmpoZjWFzTkzkAACAASURBVKGhjMfmsxSNlZ/03X24/qrlc951Slocna2D9NtHJ6YtAUZcI3SMdKEmF6NTAuq0smBGnYGPlt4X1mMKIYQQl5pzMKaq6l7gXpvN9gbwxqTH04FXbDbblhCML+LMxrkFY+4eX+8qw6IOxnx/aDNMU7pcHmxn24mNM1FQkjbnXU/UjXUNXxSMNQz4pjtlmlAIIcSVasZgbLx/WPH4t9cBH1NVdeCSzVYDJSEYW1SYyIzNMk3pGg/GjClzD1CiznhmatA5yO9O/ZT7S+4m85Kpw7rKTpxjHtZtWY5eP/dM1oX2FsOw6sLj9QO+erGCRAnGhBBCXJlmy4x5gZ9yoRvo9y55XgMGgX8O8riihj8z5pglM+bq6QZY1NOU/pqx5oEWKnsV/nDuFR5a94mLtik/1YqiwKoNywLad2qG776P7s6hix5vGC/eL5DMmBBCiCvUjMGYzWY7DJgAVFWtB7babLbucAwsWvhbWzhnyYwtiWnK8Q76Y64xwILNXotX807UcvX1jtDZOkhuUcpla1DOJjbORGyciZ5JwZhX89Iw0ERaTCpWU/wMrxZCCCGWrjnPM9lstsLpAjFVVedeyb3ImOfY2uJCZmwxT1P6gjGXxwnAqNtB+3DnxNM1Fb6vV66e312PqRlxDA2MMeZwAdA10s2Ie1TqxYQQQlzRAingLwL+A1gH+LufKoAZX+PXQBrILhpzLuDv7UUxGNAnJIRjWCFyeQv++oFGsuOz0DSNmooODAZdQIX7k6VmxNNcb6enc5jsvCSpFxNCCCEIrM/Yj/AV6/8KWA48CRwGMoHPB39o0WHuBfzdGJJTUHThbc8QTMr4NKUC5FtzgQsNWbs7hujvHaWgJBWTeX5x96V1Y02DLTDpWEIIIcSVKJDIYQfwkM1m+0egDHjZZrN9BPgG8P4QjC0qzCUz5nU58QwMLO56MbjQgV+DwsQ8THoTDePZq5qKDgCKV2fOe/dp48GYv26sefA8OkVHTnxgNwMIIYQQS0kgwZgRaBj/2gZsGP/6V8DVQRxTVJlLZszd0wss8noxmGhtoQCxxljyrDm0D3cy6hrlnK0bo0lPXmHKvHeflBqDXq/Q0zmEV/PSMthKVmwGJr0xSH8BIYQQYvEJJBirBbaNf10FXDX+dSy+5ZGWJLPRd4pmam2xFNpawIXWFooGMXozBQl5aGhUNjUw2O8grygFvWH+07A6nY6U9Dh6u4ZpH+rE6XXJAt1CCCGueIEU//w38EtVVfXAs8AJVVWHgV341qlckiwm3ylyON3TbrMk2lrAxDSlooHFYCHZkgxATVU7oFBQvPC/X2pGPF3tQ9Sc93Xel2BMCCHElS6Q1hY/Bj4BtNpstnLgs8BOoAP4XGiGF3n+PmMzZsZ6l0BbC7gQjKFh1psoSPAV1vc2O1EUyFsRnGAMoOW875xJMCaEEOJKN+dgTFXVTwAv2Wy2fQA2m+1Jm812NfAhlnABv0Gvw2TUMTo2fWbMtcQyY2hg0BlJtiSRTCr0mVmWm4glZuG1Xf4i/t6uERQUlkvxvhBCiCtcIAVAvwCmaqJVCnw7OMOJTjEmA6MzZMbcPT2gKBiT51/cHg0mWltoYND5MoLLHb6lSTMLglMWmJrh28+YHTJi07AYAuvkL4QQQiw1sy0U/lXgP8e/VYB2VVWn2vTdII8rqljMhhkzY257L/qERBTDIu97O+luSr3iC8bi+lIYAjxpw0E5hNliJNZqxDUSR6pMUQohhBCzFvA/AnThy6A9DnwF6J/0/P9j773DG7vPO9/PQSdAAOy9Nww5fTR9JI3GkmzZKpZsJXFir1vixGsndze7uYlzr5M4N5tkN1nbmzhxnOamuNuSJVl9LM1Imt5YhuSAw947CZIAUc+5fxyA4sywACQAcka/z/P48Yg4OOcHHADne97yfSODwt9IyOo2CSkGLVMz3iUfUxSF4NQUhqI7wLhUWhwZ06EoCr5RLQG9hxHJDeyOy2EM6aDvNZGvK4jL/gQCgUAguJ1ZbVB4CPg+gMPhGAAKgZNOp7Pf4XD8AfBJVBf+1xO8zg0lxajDH5QJhmR02hszu/LcHEowiC49fYNWF0ekG9OUE6NuAl4Zd+Y4PbOzcTtMyOoB9Fjnb/MaO4FAIBAI4kAsNWNHgW8ARQ6H427UOrGTqB2Vf5OAtW0aUowRe4tb68YCU2HD1ztAjEmLuil1ko7+7ikAdDlBemf6CcrLp2pjwWVSOykll6gXEwgEAoEgFjH2CeCjTqfzLPBrwBmn0/l5VIuLDyVicZuFlLC9xVJ1Y8HpsGBJu/3F2I3dlFr6e9TXVlSajl8O0D3TF5fDDOvU/UyPLp36FQgEAoHg3UQsYiwfuBT+9/uAl8P/HmLpLss7BlM4MrakGJuaBkB3m3dSAot8xkCDlpEBF/b0FGoLKwFom2pf9yE8AQ+TyjiKKcDY8CyKoqx7nwKBQCAQ3M7EIsY6gbscDsduoAp4Kfz3R4GOeC9sM5FiXN74NRhOU95pNWOe6QB+X4jcQhvVaRVISDjjIMYG3erAcVOGhMftxz3nX/c+BQKBQCC4nYlFjP0N8CPgFGoR/2WHw/FFVOuLd0XN2IppyjtAjEmShIIqxiaH1RRiXqENi95MsbWALlcv/tD6xNPg3BAAmblmAMaG4tcYIBAIBALB7Ugs45C+DewDfgN4f/jPZ4H7nU7nf8R/aZuHlPB8yvkl5lMGp+6gmjEASU1TTg57AMgtsAPgSK8mpITomO5e1+4jkbGiInV01NiwEGMCgUAgeHcTk0up0+lsApoW/ffxaJ/rcDg0wNeBnYAP+C2n09m+xDYvAM+GZ2FuCkzGSAH/UmnKKTRmCxqjMdnLSgiKJIGiMDU2j0YrkZGtRrAcGVW81nuC5olr1GbWrHn/g3NDSEhUlxZxnkFGhRgTCAQCwbucWNKU6+VxwOR0Og8BXwC+vMQ2/wPYdJXwkciYd8kC/sk7IkV5AwpMT3hIzzCj0agfkeq0CkxaE43jzWsuulcUhUH3MDnmLGypZqx2kyjiFwgEAsG7nmSKsbsJd2CG7TH2Ln7Q4XA8Cci80xiwaVioGbspTSl7vcjz83eUGFMkCGpSCQZk0rPemUep0+jYmulgwjvFoHt4Tfue9rmYD3opsOQBkJ1nxesJMDfji8vaBQKBQCC4HUmmGLNx4yilkMPh0AE4HI5tqLVof5rE9UTNOwX8N6Yp7yiPsTCKJBHQqHVi6VnmGx7bkVUHQONYy5r2PRAu3s9PVcVYTr4VEHVjAoFAIHh3k8zJ1jOAddF/a5xOZyTU9HHUUUuvA2WA3+FwdDudzpdZhvR0MzqdNlFrXSA724qsVY+jSBLZ2e+8hOmhbgDsRXk3/P22RoKARo2IFRWn3/C67rXv5butP6J5uoWPZz8e865nxlVPttqCcrKzrVQ5cjh7opPZaW/M798d837fQYhzsjkR52XzIc7J5mQjz0syxdgpVE+yHzscjoPc2Ajwh5F/OxyOLwHDKwkxgKkpT4KW+Q7Z2VbGxmbxzAfUY7rmGRt7J4oz0zUAgM9gvuHvtzMKEJDUiJiMcsvrcmRU0zLhpKm7gzxLTkz7bhvpASA1ZGdsbBajWYckQdf18Zjev8h5EWwexDnZnIjzsvkQ52RzkozzspLYS2aa8hnA63A4TgNfBX7f4XD8N4fD8VgS17AmTIalTV8jaUr9neC+H0aRIKBJAcBivbVD9EDuHgDOD1+Oed+D7iH0Gj1ZKeqAcINRR0aWhbHhWUIheR2rFggEAoHg9iVpkTGn0ykDn73pz9eW2O5LSVlQDOi0Ggw6zS2mr5Eh4XdazVhQo0bGLKmGWx7fkb0Vk9bI+eHLPFLxXjRSdHo+JIcYcY9SkJp/w3NyC21MjLmZGJ0jJ/+OnqolEAgEAsGSJDMydluTYtThuUmMLRi+3kHdlAABjRmjSYdOf2tNnkFrYFfOdqZ801yf6ox6n2Pz4wSV0EInZYS8QrVZYGRgZn2LFggEAoHgNkWIsSixpOjxeG8SY9PTSHo9GotlmWfdfqjWFgZMZv2y2xzO3w/AmwOno95vxHk/PzX3hr/nFqrRsGEhxgQCgUDwLkWIsSgxm3S4vQHkRQalwalJdGnpSOEB23cCigQhyYDRuHwGu8JeSrG1kIaxZibmJ6Pa71BEjN0UGbOnp2BK0TMy4FrqaQKBQCAQ3PEIMRYlFqMORQFv2GtMkWVCMzPo0tI2eGXxJYQGRdJiNC0vxiRJ4r6iIygonIwyOjY0pxrFFlhujIxJkkRugY3ZGR/uOWH+KhAIBIJ3H0KMRYklRU3bebyqzUVodhYUBa3tzio6D2jVon3DCpExgLtyd2E1pHJq4BxzAfeq+x1yj2DSmkgz2m95LK9IfQ9F3ZhAIBAI3o0IMRYl5nCkyB2uGwvNqGk1nf1WcXE7E9SoYmylyBiAXqPjwZL78IZ8/LL3zZX3KQcZnR8n35K7ZEo3tyBSN5b8VKVrap6hvumkH1cgEAgEgghCjEVJqkmNjLnDkbHgjBrF0drenWIM4J7CQ9gNVk70vc2sf27Z7UY948iKTP5NKcoIOfk2JCn5Rfx9XZN8/5/P8fPv1TPr8ib12AKBQCAQRBBiLEpuiYy57uzI2GppSgCDVs/7yu7HLwd4rmP5+e6RweI3d1JG0Bu0ZOVaGRuaJRAILblNvLneMsKLP1kYAkH39fGkHFcgEAgEgpsRYixKLDdHxsJiTHuHibGQpIowvSG6uZ93FxygwJLH6aELdLp6ltzmnU7KpcUYQEFJGrKsJKVurK15hOPPtaLTa3jgsVpAHckkEAgEAsFGIMRYlFhSVJES8RoLRmrGbHdYN2VYjEU7hF2r0fJrjicA+KHzaQJy8JZtImLsZsPXxRSUqKJ2sDex9VtdbeO8/otWDEYtj/36LqrrcsnJtzLYO40vLLQFAoFAIEgmQoxFiTkSGQsPDQ/doZExWaOKMK0u+o9GVVo5Rwr2MzA3tGS6cnBuCLMuBZth+SGp+UWJF2PD/S5efbYZrU7Dw7+6g+w8dT1l1VkoCvR0ROeZJhAIBAJBPBFiLEosN9WMLUTGrMsLjNsRWVLFmC4GMQbw4erHyDVn83rfWzSNtyz8fT44z9j8BMXWwhXNcY0mPVm5qYwMzRBMQN3Y3IyXl5+5iiIrPPShbQtjmADKq7MAUTcmEAgEgo1BiLEoublmLORyoUlNRdIlbdZ6UlhIU+pj+2gYtQY+tfU30Gt0fKv5+/TPDgLQNzsAQIm1aNV9FJSkIYcURgbjWzcWDIR4+emrzLsDHL6/iuLyjBseT88yY0sz0ds5SSgox/XYAoFAIBCshhBjURLpplyoGXO50N1hthayIi+KjEVXM7aYYmshH6/7CL6Qn39q/Bbj8xP0hsVYsbVg1ecXFKv1d/FOVZ55o5Ox4Tm2bM9j+12FtzwuSRLlNVkE/CH6e6biemyBQCAQCFZDiLEo0Wk1GA1a3PMB5EAA2eO+42wtFEVZEGOx1IwtZk/ODp6oephpn4uvXv4GpwbPAVBuL131ufnF4bqxvviZv/Z0THD18gDpWWbueW/1sqnSMpGqFAgEAsEGIcRYDFhMOtzeIKHZO9PwVY2MRbop1/7ReKDk6IIgG/WMk5OSRYYpfdXnmVL0ZOZYGBlwEQyuv27M4/bzxgvX0GglHni0Dp1++WhfXqEdU4qe7usTKIuGwQsEAoFAkGiEGIsBi0mPxxe4Yw1fZdYfGYvwQMlRft3xIQxaAw+UHI36eYUl6YRCCsP9668bO3W8nXlPgINHK8jKTV1xW41GoqwqE4/bz+jQ7LqPLRAIBAJBtAgxFgMWk455XwjflFrTdMfZWtxQM7b+j8bdhQf58r3/H0cKD0T9nKJyNYLW370+m4nezknaW0fJKbCyfe/qzQPwTqpSGMAKBAKBIJkIMRYDqWZ1VJBnQhUKOpttI5cTdxRFRomIsRVSerGgkWL7iBUUp6HRSvR1rr2QPhAI8eYrbUgSHH2fA41meUuNxRSVp6PTaehuE2JMIBAIBMlDiLEYsJpVe4v5CVUoBE1WRodmCPiTM08x0chxKOBfL3qDlvwiO+Ojc3jc/jXt48qZXmZdXnbuL141PXnDsfVaisrTmZrwMD3pWdOxBQKBQCCIFSHGYsAWjoz5p6YZtFbx0+NT/Ow7l3nq62e43jKywatbPzIyMvFLU66ViA9Yf1fsqUqP20/DhT7MqQb2HimL+fnlIlUpEAgEgiQjxFgMRCJj464QrTmH0ek01O7MR5YVjj/XivPq8AavcH3IiowiaUCRV3TLTzTF4bqxvq7YU5WXT/cQDMjcdbg06mHniymtykSSEKlKgUAgECQNIcZiwBqOjHWGCkDS8P4nt3Pf+x088bHdGIw6TrzkZGz49u3EUxQFBQ2aDbZ2yMxJJcWip697MiabiVmXl+b6Qax2E7U789d07BSzgfwiO8MDM7hnfWvah0AgEAgEsSDEWAxYU/RYAbfWRp6vn/ywY3xmTioPfrAWOaTwxgvXCIVuz5E6amRMQkLeUK8tSZIoLstg3h1gYtQd9fMunupGDinsu6cMrXbtH+0KRzYAXSI6JhAIBIIkIMRYDFgtBrJQ03dlmtEbHiupyKR2Zz4TY24aL/ZvxPLWTUiRUdAgKTLIGysoI6nK3s6JqLafnvTgbBomPctMdV3uuo5dHhZjHc6xde1HIBAIBIJoEGIsBixGLekKmAJzZC3RpHfoWAVGk47Lp3uY96ytE3AjUcI1Y5Iio8gb2yFaXJGh1m61RyfGLrzVjaLA/nvKo7ayWI5Uq5HcQhtDfdO35XkUCAQCwe2FEGMx4JnyopUkMj396NNuNXw1mvTsPVKG3xfi4ts9G7DC9SETrhlDhg1OtaaYDeQV2RkZmFnV4mJ8ZI721lGy81Ipr8mKy/ErHdkoikhVCgQCgSDxCDEWAwM9andfhmdw2VFIW/cUYEsz0dIwyNyMN5nLWzdqAb+EpCgoG5ymBCirUoVVzyrRsfNvdgGw/96KuHWBRkRdp0hVCgQCgSDBCDEWA4N9LkAhfX542SHhWq2Guw6XIocU6s/1JXeB60RBAWlz1IwBlFVnAtDdvnx0anjARU/HBPnF9oU6s3hgS0shOy+VgZ5pfN5A3PYrEAgEAsHNCDEWJYqiMDY8i5YAetmPxmZddtvqrblY7SZa6gdxz90+9ghy2NpCQkYJbfxUgbQMM2mZZvq7pggGbl2PoiicO6lGxQ7cWx53b7QKRzayrNB1Pbq6NYFAIBAI1oIQY1Himpon4A9hYh4Av9Gy7LZarYY9h0oI3XbRsXfSlCgbHxkDKK/OJBiU6e+51QB2oGeKwd5pSioyFmxG4knE4qLzmkhVCgQCgSBxCDEWJaNDqpmrTVb/36M3r7i9Y3seFquRlvpBvPO3R5pL9RbToFFklE3ilRapG7u5kH5xVGz/veUJOXZahpnMbAt93ZMiVSkQCASChCHEWJSMj6giLC2oRmhmMa64vVarYee+IoIBmebLAwlfXzxQIpExNt7aIkJOgQ1zqoGutvEbzHS7r08wOjRLhSOb7LzlU8brpbI2Bzmk0OkUXZUCgUAgSAxCjEXJ1IQHgDTvOAFJy1Rg9fqk2p35GIxami4NEAxuDnGzEgpqZGyzFPADaDQSlY5sfN7gQjerLCucf6sLSYL995Ql9PjVdTkAtDXf/oPgBQKBQLA5EWIsSqYnPJjMelLmp3FrU5h2r562Mhh1bN1dwLwnQNvVzX8xl+VwNyXKpklTghqdAmhvUaceNNcPMDnmpmZbHulZy9fuxQNbWgp5RXYGe6eZE7MqBQKBQJAAhBiLglBQZtblJS0jBckzh1trYjrKLsntdxWh0UrUn+/b0HmP0SDL6vo0igybJE0JkFdow2I10nV9nGAwxNvHr6PRSOw9UpqU49dsjYjBzS+oBQKBQHD7IcRYFMxMz6MoYLcZIBTCrUthOsooicVqpGZrLq7Jebqvb+66o4jRqzoOafNExiRJoqo2G78vxNkTnYyNzFFVl4MtLSUpx6/ckoNGI3G9eXT1jQUCgUAgiBEhxqJgelKtF7Oa1MiRR5fC9Fz0Mwt37S8G2PQ2F8GQ+vokNk/NWISqcKqy6aLaDOHYtr5h4LFgStFTXJHB+Ogck2PupB1XIBAIBO8OhBhbhS5XDz+6/AIA3oCapgqZrVGnKQHSsyyUVmYyPDDDUL8rIeuMB5E0qqRsDtPXxWTnWbGlmRb+O7fAltTj12xVxV+bSFUKBAKBIM4IMbYCwVCQbzf/gKBH7Zy8Ovg2AFKqFdecHzmGGrBdByLRsd74LzROyKFImlIBeXPVt0mSRHXdO9EwvUGX1OOXVmViMGppuzqMvMmihgKBQCC4vRFibAUuD11l3DtJnlQAgM43DYDWZkNWFGY90RuB5hfbycm30n19gqmJzZnqWijg30Q+Y4txbE9eavJm9HotVXW5uGf99HXeOg1AIFgPgWCIaz1TMd3gJRNZVnjjygAnrmxez8SrXRP8z/+4xNWuzTm+bMbt52s/a+R3/vo4HmEiLbgJIcZW4HTfJQBMQQt6g5ZcORyNSVMNX6Mt4gc1srPrQAmweWvH5NA7acrNVjMGYE83c/ShGj72Owc35Ph1O/MBaGkY3JDjb1Y83gDt/dMbvYxlGZpw8/c/beT5U10bvZQl6Rud4w+/cYa/+cEVTjUObfRybsE15+N///AKT73i5LuvOJn1RF8vmwxCsszPTnbwlR810Nbv4pXzm+/31dk7xZ996zxXro8zOO6muVvc0AluJGm5HofDoQG+DuwEfMBvOZ3O9kWP/z7wkfB/vuh0Ov88WWtbCkVRaB5xkma045uTSbUZKZ/PBMaZsI8CpUzP+Sglevf3CkcWaZlm2q6OsPdIGVa7afUnJRFZfidNuZm6KRdTt6uA7GwrY2OzST92dp6VrNxUetoncM/5sKSuPIXh3cCFa6P8y3PNhGSFL358LxVJruVbjVNNQzz1ihN/UKa+fZz3HyxFp90896BvNw7x1KtOAkH1+9bQMcE9Ows2eFXv0D7g4h+ebmLG/Y4Aa+2ZYn/txkWpFzPvC/L1Z5po7p4iO83E2LSX6/3TBIIyet3mOM9vNQ7ynZecAOx1ZHPROUZz1yT7tuRs8MoEm4lkflofB0xOp/MQ8AXgy5EHHA5HBfBR4DBwCHivw+HYkcS13cKEdwqXb5Yycxl+XxCrzYTdr75dXYZuQGZyxhvTPiVJYvfBEmRZoWET3r1FUiQSm6+Af7NQuzMfRQFn0/BGL2VDmfcF+dfnm/mnn18lFE5vN7RvHuuWkCzz3Zev8e8vtKLVStgtBgA6B2c2eGUqiqLws5MdfPPFVvRaDb/34e1k2U209kwR2iQ3Qg3t4/zvH1xhzhPgI++p4k8+sReAq12TG7wylRm3n7/5/hWau6fYVZXFn31yPw/cVYQ/INMxsDkapV4828O3XrxGilHLH/7Gbj77wW1YUvS0dE9uet9JQXJJphi7G3gZwOl0ngX2LnqsD3jI6XSGnE6nDOiB2JROnOmZUQvti7RFAKTajMizcwSNOuY082jsE4zHKMZAHa+TajPS2jCEx725wv1K+BqgFvBvjgvCZqO6LhedTkNrw9C79sd0dHqev3zqEmeaRyjPt/Inn9iLRiPR0r05LtL+QIh/fPoqJ+oHKc5J5c8+uY+PP+QANoeQkBWF777i5IUzPeSkp/Cnn9zL7upstpZnMO8L0j2U/KjvzVxuG+NrP2sC4Pc+vJ337i+hNNeKxaTbFELC7Q3wtz+8Qs/ILPfuzOfzH9qG2aSjrjwDgOZN8Fl89UIfPz3RQYbNyB9/7C5qitPQaCR2VGUx7vIyOj2/0UsUbCKS2ZJmAxbfroQcDofO6XQGnU5nABh3OBwS8LfAFafT2bbSztLTzeh02oQt1qEro3qojK2ZW+imldx8G8rsDIb0dAC0GcPMeUNkZ8c+pPqe+6t56ZmrdLSM8p4P1MZ76WvGMmgI/0vBlmokcw2vLVms5X2PF1t3FdBwsZ/ZKS+VjsSmGubmA7R0TVBTnE6adePTom29U/zVU5eYcft57N4KPvXIVnRaDY6STpw9k6SkmkhN0W/Y+vyBEH/+b2dpbB9nZ3UW/88n92M26SkvCfD1Z65yvd+1oZ8dRVH45vPNnKwfpKLQzp9/5tDCeT20o5CT9YN0j7k5uKsobseM9fU2to/xz881Y9Br+NJnDrG1InPhsV2OHE41DOJHomiD3kevP8jf/rCegTE3jxwp57ef2I4kqR3vR6wm/vHpJto2+DyfuNTHD395nQybkf/1u/eQl/nO2LZdNdmcaRqid9zDtprNke4VqGzkZyaZYmwGbiiw0jidzmDkPxwOhwn4JjALfG61nU1NeeK+wMWkksZfPvhHnHzNGf6LQmBmBlNuNXaDzHT6KP3DrjXVLhVVpJNi1nP+7S5qtudhNCXXpmE5ZmfVSJ+Egmt6DnkD6rKiYaNqxiJUbc2h4WI/b/3yOraM+E8B8AVCXGgd5VzLMC3dUyjA/tocPvvBbXE/Vix0Dc3wv39Yj9cf5D+9z8Gx3YVMTaqdwbtrsmntnuTtS33c5cjekPXJssI//fwqje3j7KrK4j9/cBvuWS/u8Oe6osBGW98U3X2TWEwbIxhfOtfDz092kJ9p5r8+uYOA18+YV42QF6SbkCQ4f3WI+3fFp24s1u/K0ISbv/jORWRZ4f/68A5yrIYbnl+Vb+VUA7x1qY8H9hbHZY2xoCgK//xcM63dkxysy+Xxu8sYH5+7YZvKAhvX+6bp6p3ckBuDrqEZ/u5H9ZiNOv7rkzvRyvIN7+GuGvX7cbZxkH3VWUlfn2BpknFdWUnsJTNNeQr4AIDD4TgINEUeCEfEngUanE7n7zidzk1TsOQOd0ymaEOgKOhsaezO2Y6kCzAW6l/TPnV6LTv2FeH3hWjeRK3i4EwCjQAAIABJREFUCzVjyuYaFL7ZyMm3kVdoo7djcmE6QzyYmvXx/dfa+G//cIpvvthKc/cUJbnql7e5a3LBemQjGBh38+WwEPvMI3Uc2114w+O7atQI4UbaCvz4jXYutY2xpSSN//z41lsKuLeVZ6AocLVzY1JYzV2T/PSNDtKtRv77r+26RSikpuipKLDRPuBibj751gdef5B/fOYqXn+IT3+glq3hlN9itoejZBtVH3j8Uj/nW0epKrTz6Ydr0YQjYovZXpmJAjR2JH+Nsx4/X3+miVBI5nc+uJWinNRbtsnPtJCbnkJL9xSB4Ka51Ak2mGSKsWcAr8PhOA18Ffh9h8Px3xwOx2Ooxf1Hgfc7HI4T4f8dSuLalsUT9hIzyOrdtc5uZ3eO2lvgM/fhD6zty7R1dyEGo5aGC/0E/MHVn5AElMUXe1EztiI79qlppKaLaxPki/F4g3z/eBt/9I3THL/Uj1Gv4ZHDZfyvzx7izz61j7t35OP2Bukd3Zho4KzHz9/9pAGPL8inP1DLwa15t2xTU5pOaoqe+uvjG+KVdb51hFcv9JGfaeZ3P7QD/RIlDLuq1YjEletjyV4ekzNe/vm5ZjQaic8/sZ0M29Kd1Lurs1GUjRE7T73SxuC4m/vvKuLQtlvPMUCGzURJbirXeqfxeJP7u9UzPMuPX2/HZtbznx/ftmxXbOQ812/ALODvvuJkYsbHY3eXLwjXm5Ekid3V2fgCIVp7hMWFQCVp+bFwYf5nb/rztUX/3lw+D2Hmw0X2Br+ajtHZ7VTYS9ErKSjpo4y5PBRmxZ5nNpp07NhbxMVTPTReHOCuw6VxXfdaWBiHhLIpTV83E+U1WaTajFxrGmb/veUY15j2qm8f57svX2N6zk+W3cSjh8s4tC3vhgtNXVk6bzcO0dI9RVlecq0jQrLMPz5zlXGXl8eOlHFke/6S22k1Eruqs3i7cYiuwRkqC+1JW+PIlIdvvXQNo17L55/YjnmZtH9RtoUsu4nGjomkWh8oisI3X2xlbj7Ax95bs6L9x+7qLH56ooP66+PLvteJ4Mr1Mc40D1Oeb+XX3lO14ra7q7PpHZnjatdE0iwugiGZb77YSkhW+K1H60hfoX6yINNMTloKTV2TST3PF6+Ncsk5RlWRnUcPl6247a7qLF4+30v99XF2VIpUpUCYvq7KvNuPJIF2Xm2J19rsaCQN+doKJF2AxpEV+wxWZOf+YowmHfXn+vBtAkfmhYiGoghri1XQaDRs21NIMCDT2hC7UaesKDz9Zgd//9NG5uYDPH5POX/5mYPcs7Pgljv+2lI1XdS6AR1iL5zuoa1vmrtqsnns7vIVt90TjkhcTmLkSVYUvvXiNXz+EJ94yEFBlmXZbSMRCa8/hLM3eRGJt8JCekdl5i3p3ZvJz7SQm2GmqWtizVH3WPF4A3z3FSc6rcSnH65b1Ydtd7jO6UoSI08vne2hb3SOe3bks6186YhTBElSbwx8/hDXknSe5+YDPPWqE71Ow6c/UItGc2v6dDGVhTY1kty+MZFkweZDiLFV8Lj9pFgMyLOqGNPZ1Tv+GpvaBdk63brmfRuMOnYfLMHvC1K/CXzHIr8JEgoEN0fqdDNTuzMfvUFNNQdjqP2QZYVvvdjKL073kJ1m4k8+sY/HjpQvewdvtxgoyk7F2edi3pe889Ix6OK5U91k2Ix88gNblqzPWUxdWTpGvZbLbeNJsz44cWWAtr5p9tRkc6Bu9SjNnhpVSFxOkpBwzfn40evXSTFq+fj7HAtdfyuxuzoLf0CmJUkprGfe7MI15+fRI+UUriBmIxTnpJJpM9LYMUEwCbWlkzNefnGmB3uqYdWoXYSIYExWqvL5U93MegJ88O5y8jLMq26v1WjYWZnJ9JyfnuHN2SglSC5CjK3CvCeA2WwgOKO6cmjDYqw2qwoloKfP146srP0HadtdhZgtBhov9G+471ikQFxCQRFibFVMKXq27i7AM+fnWmP0JrD/8Vobp5qGwx5d+yheosj3ZnZXZxEMyUnzyQqGZL714jUUReG3Hq6LqvvQoNeyrSKDkUkPgxOJ7XYGcLn9/PREB2ajjo+9tyYqoVNVZCc1Rc+V62NJiUg881Yn874QTx6tXLZO7GYWIk9tiY8wDk24eePKALnpKbw/PK5tNSRJYldVNvO+IM7exI/B+tnJDgJBmSePVmKOshygqsiOxaRLynkemfTw+uV+stNMPBhDh+mu8Hm+nITzLNj8CDG2AgF/kIA/RIpFT8ilijGdTa33KMhMJTSdgx8P3WGD2LWg12vZc7iEYEDmytm17yceKIvSlHJg49OmtwM79xej1Wm4craXUBRRgjcbBjlxZYDinNQlO+qWY0+4HT5ZP9zHL/YzOO7m6O5CtpSmR/28iK3F+ZaRRC1tgWff7sLrD/GhoxWkRTmaSqvRsLs6C9ecn7YEC4n+0TneahyiIMvCvTFYVVQW2ElLNXC5bWxhTFKi+PHr7ciKwq8cq4ppTNTeLep5Ptea2PPcNTTDmeYRSnOtyzYVLIVWo2FPTTbTc36u9yX2PP/0ZAchWeFX7quKqT5tW0UmRr2W860jG26iK9h4hBhbgblZNVKVYjEQdLlAktCmqsX6VrMe3az6A3tltGnZfURD3c4CrHYTVy8P4JraOFfmGwr4RWQsKswWA3U785mb8XG9eeUL0+j0PN97rQ2LScfvfmh71Hf5ACW5qWTaTDR2jCc8NTQ16+PZU12kpuj50L0VMT13d3U2RoOWM83DCb3ADIy7OVk/QH6mmXtjnOV4MJzOPNuS2JFWPznRgaLArx6rRKuJ/qdWo5HYX5uL2xvkamfirELaB1w0dExQU5y2EI2LluriNNKtRi45EysYn3tbHe7+q8cqV02T38yBhfOcOMHYPzbHJecYFQW2mP31jHotu2uyGJv20jm0OcZ0CTYOIcZWwD2neoyZLWqaUpuaiqRTO7UkSSLPUIIS0lI/dnVdFx6tTsPB+yqQQwpn3uiIy9rXwg1iTETGombXgWI0GonLZ1aOjn3v1TYCQZmPPlhDdlpsZrGSJLGnJpt5XyjhY4eeO9WFzx/iw0crYjbNNOq17KnOZtzlpWMgcReYn7/ZiaLAk/dVxjz421GSTlqqgYvXEickuodnaOqcwFGctqzFwUoc3Jp4IfHC6W4AnrinPKoU72I0ksSBulzmfUGaEiQYe0dmaeiYoLrITm3ZrZ5nq7GlJB27xcDFa6MJu4F58UwPAI8cLov5PQQ4EO5GPZeESLJgcyPE2AosGL6aDYRcLrS2G9v18zOshKZzmPRO0Te3PvPWyi3Z5BfZ6WobZ2CDvGcWfMYUERmLhVSbidqd+bim5rnWuHRnpbN3iqbOCWpL06MqNF+K/XWqserpq4mL6IxOz/N24xC5GWbu3rE2a4VD29TXd6Y5MescHHdzqW2M8nwru6pitwXQaFQh4fEFaexIjJB4YZ0X6dJcK7kZZurbxxPStBEROlVFdmqK09a0j4iQSJRgfPGs+h4+fKhsTc+/IcKYgFrLkSkP51pHKMpOZWdl7IIbYGt5Bqkpei60jm6oqbNg4xFibAXmwmLMZNQgz88vdFJGyMswE5pUf5AaRq+u61iSJHHkAbVT6NTx9g35Yi7uphRiLDb2HilFp9dw4e3uJU18fxG+OD9xb8WaLs4AFfk28jPNXG4bx5MgK5TnT3URkhU+eKQsptTaYmpL1YjE+daRhEQkXlp0kV7re3mwTq0/OpsAwTg04eayc4yyPCt1ZdHX2y1GkiQO1eUSCMoJqRN8+Zxan/rIodI1v4cluankZ5ppaB+PuwHs+PQ8F66NUpKTyvaK2KNiERZSlQk4z8cv9qMo8PA63kOdVsNeRzYut18YwL7LEWJsBSLdjUbU/9feJMbyM83Iriw06Lgytj4xBpCdZ2XL9jwmxty01A+ue3+xItKUa8ecamTn/mLm3QEaLtzoyj847qa5a5ItJWlUrcMMVZIkDm/LIxiSOd86ut4l38KEy8uZqyMUZFnWZeap1Wg4uFWNSFxyxldITLi8nG1R17hrHXP9SnJTKcy2UN8+jivOXczHL/ajAB84uPaLNLyTqnyzIb6/BS63nwvXRinIsqwphRoh8nkMBGXOxbn+7mTDIIoCD+4rXtd7WJ5vJS9DvYGJ54gpnz/E6avD2FMN657FenibGoF+qzH5v/mCzYMQYysw71F/pHUBtaheZ7s1MoaswxYqYMQzyrB7/eH6A0fLMRi1nDvZuZAmTRbvdFOCEhRiLFZ27S8mxayn/lzfwmcH4NRVNXV53yqGn9FweFs+kgRvN8VuNLsav7zcj6woPLS/ZFXTytU4ukt9rW/EefbqyYYBQrLC+/YXx1zQvRhJkrhvVyEhWeHtOF4E531BTjcPk241srtmfc7qOelm6srSud7vYmBsbvUnRMnbjYOEZIVjuwvXJXQA7t6ej1YjcaJ+MG4NG8GQzFsNg1hMOvZtyVnXviRJ4uiuAoIhmdNx/M6cax1h3hfk3h23mjTHSmWhjcIsC5ecY8x4NtbeSLBxCDG2At7wXEqdX/VM0tpuHGOSk25WLwgu9c7myjpTlaBGWA4dq8TvC/H28evr3l8siDTl+jAYdew8UEzAH6KnXa1FkhWFs80jpBh1MXesLUW61ciOikw6B2foHIxfgbzPH+LN+kFsZj0H6tZ3AQT1RqW2NJ22vmkGxt1xWGHkIj2E2ahbqFdaD4e25mHQazhZPxg3L6pzLSP4/CGO7ixYc5p3MfeFRe3JOEXKZUXhZP0gBr2GQ0vMGI0Ve6qRXdVZ9I3Oxa0j8HLbGDOeAEe252PQ3zpjNFYOb8tDp42vYDxxZQBJgqMxWJYshyRJ3LurgJCscCoBN1mC2wMhxlYgEt3QelWH5JtrxvQ6DflZZib77WglLQ1j67O4iFC7M5+8IjudznG62pI3cuSd3ymRplwrufmqYI9YlPQMzzI162NPddaSw6vXwgP7VGPJ4xfjN7Xh9NUhPL4g9+0ujNs6I6N/TsQpOlZ/XU0pHt6eF5eLtNmkirpxl5fmOBR4K4rCiSsDaCSJe2K021iOXdVZ2C0GTl8djst4pJauScZdXg7U5i47wzNWIoLk5JX4CMZIWjYekWQAq9nAXkcOw5Me2uLgOdY/Nkf38Cw7K7OiNvJdjUNb89Dr1BsD4Tn27kSIsRWYD9cYaObC7vu2W+t9SnJS8Xk1lKWW0zc3yPj8+ruzJEnivodq0Ggl3nqtLWlzKxcc+BUFJSAiY2vBlqb+OM9MewFoCnfr7VhD199y1JWmU5hl4cK1UabilMo+2TCIViOtOjsxFnZVZ2FPNahCLw4F3ifrVVEXiRbFg8gF/5eX+lfZcnX6RufoHZ1jZ1XmioOsY0Gn1XD3jnw8vmBcuhbPhL3w7tkRH7EIUFeWQXaaifOtI8yuM83mmvPR2jNFZaEtqrFC0RIRjL+8vP4bg4gNRaSmLx6kpujZ68hhdGo+LjcGgtsPIcZWYN4TwGjSEZpR76Z09ltbwItzVBPYXI1qjlkfh0J+gPQsC3sPl+Ke9fPmK23JuVu6wfRVRMbWgjnViEYrMTOtRsaaOifQSBJb19hVtxSSJHH/3iJCssLxS+uPjvWNztE7Msf2ikzsUTrZR4NOq+GBu4qY94UWhNRamZr10dI9RVWhfcVh4LFSnm+jqtBOY8cE/eusy4q40ccj/beYY7sL0WokXj7Xu650qj8Q4vL1MbLsJioLbas/IUo0ksQDdxXjD8q8vk6xc/7aKIpCXNLQi6kpTqMkJ5VLzlFGp9durK0oCudaRjAatOyM4w0WwIP7igB4+fzGTmIRbAxCjK3AvMePKUWvuu9za5oS1K4sAM1sHhIS9XGoG4uw+1AJeYU22lvHVnV3jwcLgk8jiXFIa0SjkbDaTcxMe/H5Q3QNzVKeb43JbT8aDm/Nw55q4PVLA+uORkTqVI5sj6+IAFVImAxaXr3QRyCGYeo3c+HaKAqs2aNtJd5/UJ3JGLF7WAuKonC+ZRSTQcuONXpOLUeGzcTBrbkMT3poWMfg64aOCXz+EAfqctdduH8z9+zMx2LS8ctL/fjWkU493zKCJMG+OIsxSZJ46EAJigKvrkPsdA7OMO7ysqc6C2McUuWLKcuzsaUkjZbuKTE8/F2IEGPLoCjKosiYC7RaNJZb78gjQ55HxkJUpZXTNdPDtM8VlzVoNBruf7QWvUHLW69dX4i2JIqFm26tBsWf3E7OOwlbWgre+QDtfVPIikLlOuwslsOg1/KBg6X4AqF13UmHZJmzLSNYTDp2VMb3Th/AbNJzbHchLrefU+swqz3XMoJGktbdXbcUO6uyKMiycK5lhAmXd0376BiYYWLGy+7q7LjUs93MQ/tVwfjSOgRjJL0W76gTgMmg49ieQubmA2vuWhybnqdjcIa6sE9dvNm7JYdMm5G3G4fWfAMTiX4m4qYA4P0HSwF4RUTH3nUIMbYMwaBMKChjTNETnHahs9uXvJu0mg2kW41qQWf2NgAaxprjtg5bWgp3P1CF3xfitWdbCK4jurAaC5ExvRZ5fm0XpUQSdLmYOX8WT3987RLiTaRu7HqnWvuRCDEGcHRnwUJ0bK0t8W2908y4/eyvzY1pyHEsPLivGJ1Ww4tnetZkAjs65aFraIbasnRsCbhIaySJh/aXEJKVNUfHIhfpeNYRLaYw7PLePuBakzmozx+isWOCgiwLReEbyHhz/13qeX7pXO+azvOVsLltvKNiEXRaDQ/uK8EflHn1QuzpfUVRqL8+TopRR90axjNFw7byDIqyLZxvHWVkypOQYwg2J0KMLYMvXLwfiYwtVS8WobLAxozbT7GxGgmJ88OX47oWx/Y8HNtyGR2a5e3X2uO678Uokd9PnRbZu3EDyyMoioKvr5fxZ5+h53/8OZ3//b8w/C/foOtf/32jl7YiNrs6d7I/bD1RkR+/+pzFGPRaHjlUhi8Q4tm3uta0j8vhbt096zSuXIm0VCP37S5g3OVdk0XDxbBxbCIiOhEObs0lJy2FE/UDjMUYgV58ka4tjV9t4M08eqQcgKdPdsRcQ9rSPUkwJMfFXmU57BbDus5zQ6TZJc5p3sUc3VmA3WLg+MX+mM1+B8bdjLu8bK/IWLe32HJIksSjR8qRFWXN32nB7YkQY8vgC3d/GXWgBIO3uO8vpqpIFWpjowq1mTV0z/QyOBc/R2pJkrj3fTVk5abS2jCUMHd+hUhkTIfs9aLIiRmuuxqB8TEmXnienj/7f+n58z9l8vln8fX1krKlFjQagu7YfasURWF0ep769nFeOd/LD395nX9/oYW//2kjT78Z+8VtJTJz1HT2zIQHm1lPhi1+RfE3c3RXAfmZZk7UD9A3GlsBuqIoXL4+hsWkw7HG+YTR8sjhMkwGLc+d6op51mJ9+ziSxLoc91dDp9Xw+L3lhGSFn7/VGdNzB8bdTMwk9iINUFFg466abDoGZ6hvj612LLJ9vIvOb+aRQ2UY9VqeP92Nzx99FN/jDdLWN015vpW0ODaR3IzRoOWRw+oNzAtnumN6bkPkPUxAOn8xdzmyKclN5VzLSMzfacHtixBjy+ANR8b0qBeOm933FxMZcXN9wMWR/P0AnB46H9f16PRa3vfEVkwpOt569XpChokvjEPS60BRkH3JqxtTZBl3UyP9/+crdH3h/2bimZ8RGB0ldc9d5H/2c1R+9WsU/8EfoTEaozaklRWF1p4p/v2FFv7g66f5wjfO8Pc/beRHr7fz6oU+TjUNU98+zi9O93C9Pz51fgBZ4aYOyR+iMDs17sXSi9FpNXzk/moUBX5wPLau2+6wB9qOyqyEiggAm9nAQ/tLmPUEYqqHmfX46RhwUVVoJzUlvk0QN7O/NpeSnFTONo/QOxJ9AXWyLtIQmW0KT7/ZGfX8WllRaOyYwGrWJyxKG8FmMfDgvmJm3H5ei8EH72rXBCFZSbhYBPUGJstu4sSVAcZd0UdBG9onkCTYnsDIHahp8w/dW4mCGgUVvDsQYmwZIpExvbz0XMrFlOSmYtBpaO+fZltWLVZ9KueHLhOQ4+vVZUtL4b2PbwXg5aevMhHnu6ZImlIKFyDL84lPVSrBIK43T9L9xT9m4O++gudqI6bKKnI/+WkqvvJ3FHzu97Du3Y/WrHoOSTodchQeaE2dE3zpmxf42x9c4VTTMKGQzN4tOXz4aAWfe3wbf/KJvfzPzx7ic4+rdX5n4jhIOMVswGTRYwYKMuNnw7Ac2ysy2VGZybXe6Zhex9VONS2UyIjTYt67vxi7xcBL53qjthe42jmJoiQ+ogPqRfDJ+9SL4PePX49a2CbrIg1QkGXhyPZ8BsbcUY+a6hmexeX2s6Myc91jrqLhof0lpKboeeFMD5Mz0dWe1idR0Oq0Gp64p4JgSOFHr0dX9hG5KahOwk0BwPaKDGqK02jomKCpc/3elYLNjxBjy+ANG63qQuqPyVK2FhF0Wg3l+TYGxtz4/AoH8u/CHfTQMBofR/7FFJamc+zhLfh9IV74SSNzUf7YRcPCaEq96sydyLoxJRTC9dZJur74BUa++y2CkxPYjtxDyZ98iZI//iL2u+9Fa75VyEg63YqRMX8gxLdfusZXf9zAwNgcB+ty+cJH9/DV37ubzz2+jYcPlbF3Sw7l+TZy0lLYU5NNWqqBC62jBILxS8saUo0YkMi2JP6HG+CjD9Zg1Gv5wfHrTM9FF9GMFIInss5pMSaDjo/cX00gKPO9V6OL4jV0RC7SiRc6ANsqMtldnUVb3zSno+j+TGbkLsKHj1aSYtTx9JuduKI418mM3IE62eDJ+yrxBUJRiR1ZUWjumiQt1bBgFZRoDmzNparIziXn2MJNyUq09kyhkBzBDWppykcfrEEjSXzvtbZ12cIIbg+EGFuGVKsJvUGLTVGjTyuJMYDasnQUoLV7iiMF+5GQeL3/7YSYtdZszeXgsQrcs36e+2ED7igvvquxkKY0qh1rodnEeN3MX2+j9398iZHvfIvQ9DRp9z9I+f/8W/I+9ZuYSstWfK6k1SEvI8bmfUG+8uMG3mwYpDgnlT/71D5++7Gt1BSnLZsq1GgkDtbl4fEFuXJ9bJ2v7B1CJjW6qI+hbmY9ZKel8OR9lbi9Qb77snPVz50/EKJ9YIaSnNSkiQiA/bU51JWl09Q5wSXnyu93SJZp6pwky26Kq9Hravz6A9UY9Bp+/EY7c/Mr++1FLtKJLDq/GbvFwIePVjDvC/KjN1YXO9d6ppAkqIuj8fBq3L0jn8pCGxeuja7qKD845mbWE6C2NCOhKf3FaCSJjz1YgyQRldiJ3LgkqotyKYpzUnnPXYWMTs2vy9JEcHsgxNgylFRk8Ed/+X6sAdV9X7tCNyWoqSKAxs4JcszZbM+qo2emjw5Xd0LWt2t/MbsPFuOanOfZ79fjjsdYnMigcIvaDRicjm9dmuzzMfLUt+n7X3+Fr68P25F7KPvrvyXn1z+KLi3KC4VOu2RkLBiS+drPGmnrm2bvlhy++PG7KMm1RrXLe3aqg96PX1z/SJwIkQo0z0TyulKP7SnEUZxGffs4J1Zp8ugYcBEMyWxJUlQsgiRJ/Kf3OtBpNTz1qnPFjrae4TnmfUG2lifvIg2QZU/hsSPlzHoCfP+1thW3vRa+SCf7fbxvVyGleVbONo+seBPh9QXpGJyhNDf+xsMroQmfZ0mCb790bcWmjWRHaCOU5Fq5/64iRqbmeWaVzsXWnilSjNqkRe4iPH53BXaLgV+c7qZfFPPf0QgxtgIajURwenn3/cWU5llJTdFztXMCRVG4v+ReAF7vfTMha5MkiQNHK94RZD+oZ3aNhpURFqIpETE2GT8xFpyepvev/gLXyRMYCoso/uMvkvep30SfHtsPsKRdOk354zfaudY7ze7qLH7nsbqYhl3nZ1rYEfZw6gzbUayXcY8fPzDUO03An5w5nxpJ4jOP1pGaoucHx9voHl7+tXSEX2dNgrsolyI3w8yv3FfJrCfAN19oXTaKd603LHRKknuRBnjf/mIqCmycbRnh7Ap1eK09U5gMWsryohP+8UKjkfjNh2vRaTV8+6Vry4ralu5JQrKSdLEIqth55FAZEzPeFUVt64KgTf5n8cP3VpKTnsIr53px9i79ezc542V0ah5HcTpaTXIvmWaTjk+8fwvBkMK/PN8S11IKweZCiLFVCM1EhoSv3IWkkSS2VWQwPeenb3SOSnsZpdZiGsdbGHaPJmRt7wiyElyT8zzz1GXGR9Z+96QsRMbUYvnAZHwKR2Wfj/6v/A3+gX7sR++j5It/Skpl1Zr2JeluTVM6e6c4frGf/Ewzv/VI3Zp+MB/cVwzAS2d71rSumxmf8eFL0RHwh2hvjV/6czUybCY+82gdoZDC15+5umyarSfcLZhsERHh/r1FbC3PoKlzYtmI5ELUqST5F2mtRsNnHq3DqNfy1KttjC/RcDA542Vkah5HcVrSL9IARdmpPHm0gllPgG+92Lrk3MqmcL1YsqNOER49UkZpnpVTV4e5cO3W30FZVnD2TZOTlkJW2J8vmRgNWj7zSB1I8G+/aMXtvfX70rpB0c8Iu6qyOLqrgP6xOZ55MzbbFcHtgxBjqxB0udCYLWj0qzt/7wp3fF10jiFJEu8tO4aCwgtdryZsfaogK+fweypxz/n5+feu0LdKjcZyLEQo7KkgSfgH4pO2m/j50/gHB7Efu5+cj30iqvdyOW4u4A+GZL7zshMJ+PTDtaQYdWvab11pOhUFNi61jdE1tL7omMcbYN4XxJBtRqORuHymh9AaHMlBPScBfyhqGwNQU+aPHilj3OXl73/WiH+JWYE9w7NYzXrSrYnzdFoJjaRGdqxmvRrVvMmqJRiSaeufJj/THNfh5bGQm27mNx6oZt4X5GtPN93im7UQudugizTAA/uKqStLp7Fjgl+c6r7l8cb2MbQaieqixEyBWA2dVsNvP1qHQa/cbBBiAAAgAElEQVThmy+2MnDTMPaekVnmfcENiYpFqCy08+hhNYL3L8+13PJd28ibggi/9p4qctNTePl8LxeXELWC2x8hxlYh6JpeNUUZYWdlFga9hvMtIyiKws6srZRYi7g82kjfbOJG+EiSxM79xTz4wTpCIZkXftzIlXO9sTcPRCJjRgOGgkK83V3rHhgemptj+uQb6DIyyf7Vj6y79kfSqjVjkdd25uoww5Meju4upLJg7RccSZJ48mglAD95o31djRfj4XRxZpaFrbsLmJn2cuqXK+9TURSmJz10XBvl4qluXnu2hR9/8wL/+uW3+LevvMV3vnaaS6d7CEWZpnjs7nIO1OXS3u/i337RckPUZG4+wLjLS2muNam1WDeTlmrk809sB+DrP796g/N999As/kDya9pu5p6dBRzdVUDf6Bz//uKNKdVrPWo96UakUSNoJInfeWwrmTYTP3+7a2GkEKgNLe1905Tn2zAZ1naTEg/yMy18+gO1+PwhvvZ00w3RJ2ev+h46NvA9BHjsSDnbKtRI7TM3mf629U9jMekSNkYqGkwGHZ//0HaMei3//kIr/WOifuxOY+O+obcBciCA7HajLS6JanujQcvu6mzOtYzQPTxLeb6Nxyof4h/q/41nO17i8zt/M6EXv6raHFKtRl79eTNn3+hkZGCGYx/YgtEU3WlWL9gSEmDZto2pV15m9vw57EfuvmW7jgEXVzsn6RmZZXRqnlmPH19AJhSS0Wo16HUSeq2Gfd4u9vn9pB17Dxr9+guIJV34tYRCBCUNz5/uRqeVePRw2br3vaU0nR2VmTR2THCmeZjD2/LXtJ/IsOksu4n9uwoZ6J2m+fIgM1PzbNmRjznVQMAfYtblZXLczcTIHOOjcwQDNwotnV5DeqYZs8XAyOAM59/soq15hAcfqyVrleYEjSTx6Q/UMj3r46JzjG+92Mqn3l+LRiMxNKFOMCjK3riLS4Sa4jQ++mAN333FyVd/3MAXProHm8VA+4BaHrBREZ3FfPTBGgbG3Vy8NspP00z8yn1qir1j0IXJoKV4Ay/SoM7H/b0Pb+evnrrEvzzfwh98ZBeVhXa6h2eRFajaBO/h/tpcekfmePFsD//wsyZ+/1d3YtBr6RwMn+cEzW+NFo1GFbV/8e2LvHCmh+y0FO7dWcCM28/YtJftFZloNvDGBdTv628+XMvXf36Vv/tJI3/8sT1k2EwbuiZB/BBibAUCURbvL+ZAXS7nWkZ4q3GI8nwbW9KrcaRX0TrZRsN4M7vCw8QTRV6RnSc/tZfXnm2hq22c8eELHHt4C4XRRBgikTGNhP2+9zD9xuuMfu+7hFzT2A7fjdZu52zzCM+f7mZ48p0hthaTDnuqEaNeg06rIRhSCARl/IEQ1g619d60Y1dcXl9EjCnBAPXds4y7vBzbUxi3dNvHHqzhT3rP84Pj19lalrGmFNn4TESMpWAw6nj0Izv55fOt9HVN0dd1a5GwJEF6loWsnFQycyykZ1lIzzRjtZsWxLvfF+TcyS6uXh7g6e9e5u4Hq6ndmb+iuNfrNPzeh7fz5R/Vc6ppGFmG33y4ltEpNQKVk5H8Gp2luG93IWPTavv+l39Uzx/+xu6F5oOKdUQ744VOq+F3n9jOX3/vMi+d7cVs1HFsdyFDEx5qS9OTYqS6GiW5Vn77sa384zNN/J+fNPBHv7FnQehUFiTWdT9aPnRvBSNTHi45x/jGs8187oltdAzOYLMYyLRvvKiwmPT8l1/ZwV//x2W+8/I1zEYd2vC5rSzcHO/h3i05PHFvBc+82cmXf1TPFz66B6t57WUfgs2D9ktf+tJGr2FNeDz+LyX6GNLUOCOvHcdStxXLtu1RPSc7zcTppiE6Bl28Z08hBr2WMlsxpwbP0T7dxZGCA+g0idXAeoOWmq05oEBPxwTOphG8ngAFxWloVxh70+TswjcmUeCwUFpWjrG4GHf9FdyNDUy9+jKdJ8/S3nAd2eNmW00ej97n4GPv28Lj91Twnj1FHN1VyD071LTOsT2FPLC3GP9Lz+CVNbRWHVkYG7Ue3A31+AcGsB97gB+f7mdsep7fergubj9IZpOeFKOOS21jdA/NcnBrbswX20vXxmgfcPG+/cVkhP3qarbmUlqVSXqmmZx8G8XlGWzZkc+eQyUceaCKHXuLqHBkk1dkJy3DjNGkv0FoaXUaSiszyc5Lpadjko5rY7im5ikqS0erW/6c6nVa9m3J5VrvFE2dE/SOzBEMKXQOzvDefSVkp8VHkFksRjye2AYvL6auLJ0ZT4DGjglauqfoGJwhxaDjw0crNjSVGsFo0LK7KotLbaNcbhunZ3iWsWkvB+pyk+o9tRL5mRay7CbOtY5yqW2M/tE5PL4gH7m/es21lPFEkiR2V2fTNTRDU+cELT2TDI57qCtL50Bd3kYvD1CjjFtK0znbMsKF1lG6h2eZmw/w8OEycjbJd6WmyI7XH6KhfYKW7kn21GRjNETfPS5YmvWelyiP8efLPSbE2ArIg72Mv32K1Lv2klJVHdVzNJJEMKTQ1DmJ1WxQnbkNFkJKiKsTrcwHvWzLqk3wytUfvsLSdEoqMxjud9HbOYnz6ggWq5H0LPOSF7jGa534xjUUbLFQlJ2DITcP+9FjBM1WegamyJgZptg7hmOuh6LOSxjPn8B3+QLulqvMX2/D19ONf3iI4NQkwZkZAuNjeE6+zpAll2dcmeyvzcWyTq8jT2srvu4u5B37+OHZYaqK7HzgYOm69nkz5flWBsbcNHVN4vEG2V4Rm8/V6atD9I+5efRw2YK3kyRJWKxG8grtFJamU1CcRmZOKharEU0MnXhpGWaqanMYHnTR1zlFZ9sY+UVpmFOXF6N6nYZ9W3LoGprhatfkgn3Hh+6tiNtFer0/ZJIksb0yk8lZH00dE4RCChUFNu7eURCX9cUDs0nHzsosLreN0Rv2fHroQAn5SRh5FS0luVasZj0Xro3i8QXRaCR+7T1r61xOBFqNxF012Vzvd9Eengd79458qos2rjj+ZtKtRqoL7ZxvHV2wDPnYgzXoV7jpiYV4fFe2lmcw4/bT0DHBlfZxdlVlJtVH7k5EiLE1kgwxFmh3MnXpMra778FYWBT18/IzLfzycj89I7Mc212ojkuyldI43kLzxDWKUvPJs+QkcOXvYLEa2bJTvevs756io3WMgd5pMrItWG5K7TVd68I3rqGw1kJhlrq+IBq+dnaGE0oh1v+/vTuPjqO6Ez3+rapepVZrsWTZ8r5ejME2xsZAwBgIa0hgQiCZSchKBmZgkrzkTXJeJnkzmUnyXiYhk2QyycxLMslkmZCFQDYIEGywwWA7YION8bUl27ItL7LWllqtXqrq/VElqW1Lsmykbhn9PufUqeqq7tKv66qqf3Xr1q1rrmfFW1YTrqvDipWBZZFrbSFz6CDp/ftI6V0kX3mZri2bSWx8lsTG5wAIL1nOc8k4+492cfkFU15X24vehnpSe3bTMGk+r7Q43PqmOcwa5e4ZDMPgwrmTeLm+hZcbWnHcM+saYN3WJo539PKONfPGpMuDcCSAumAKuZxDY30revsRIiVBaqYM3SA/GDBZdX4tnclMf7cWd14zf9RqnUbjQGYYBsvmV9OdyrLvSBdXLasrSj9ow4lFg1ysati44yiu6/KuaxcQDo6vWok5U+NMikfYuqeFqy+ewZK5hXs6wEgELJOViyaz52AHrYk0t181r2h39Q6lujzKotmVPL/jKHOmxrnm4pEf/09ntPaVJfMmkbNdtu1pYfOuZubXlUsbsteh2MmYMRaP6ymE48e7xjzw5GO/pumhh5nxqU8TXbDwjD770DMN/P75Ru5YM4+b/Jqbw91H+ec//StBM8AnV3yEmpLCHiQ721NsXFvP/j1e/2Gz5k1i5ZWzqfGTmR8//Ee6dICVfzaZFep8AB58ag9PbDnI5RdM4UNvWXTKj7frujjd3eQ6OsglOrETCXKJTpxUD67/iJGKa97M9zYcZvNrzVx90TTec/3Cs04C2h5/jJZf/Ixnz7+ZjdkavvaRK8bscT7tXWm+9N8v0dye4roVM3jnNfNx073kWlrIdXb0f9f+cWcCuzvBkdYefjPtar74qVvGJK58jfWtrP39a/SmcsycW8XqGxZSNkz7G9d1efaVI7jA6qWjV+tUU1PG8eOj9/isppYktZVRAsNcVi+mVDpHR3d6XNWKnaylI8WcmVV0JQr3FIgz4TguLYneUbv8Nxa6U1ks0xjVy7yjva88seUgP1u7B9MweNe1C7hm+bRxcWn/XDPa5TLE3xiyYIrfkGAcy7R6SUvgDHuJB+/yxbqXmvjd841cungKlWVh6mJTuHPhbfxk1y/41svf4xMr7iMWLNzBvLwyyk23X0hTYzubN+ynsaGVxoZWZs6t4sIV03H9m/kMw/sB3Hs4wZNbDlJbGeWuG9SgO7hhGFhlZVhlZYSZMeTfft+NcQ639LBuaxOVZWFuOcu7H/vKIt3cjFq6YEwSMae3l/Thw5iHm7ivtJFdDQ2EH+5k5y+ThHKneeyUZVFl28yxz66vtzM1a/4k7vjACtY9qjmwt40Hv7uZS1bP4YLl0wZtH2gYBleOYhI2VqYV8FmUZyMaDoyLdljDqa6IEgkHGNufl7Nnmsa4TsSAgj639Wxdv3IGM2pK+favX+UnT+5m+95W3nfjeeOutlEMb3wfTYos3eInYyN9bmKe0kiQ29fM40ePa370uOZvbr8QwzC4vG4lzT3HefLA03z75e9z39IPURIs7AFp2qxKbptZQVNjO1uebeTA3jYO7G2j79/B8BusP/RMAy7w/pvOe92XYqLhAB+7Ywlf/PGL/Gr9XnozNm+/au4ZX7Lse5D4oq79WGUXYXd3Y8VeX9cCuc4OurZsJlW/h/T+/WRbTuwxfxZgmwHarVK6YzVMnTedujl1BOLlWPE4VrycQHkcKx6n8/mNHP/hDygNF+7SVSwe4ZZ3LmH3jmM891Q9G59q4NWXDrPqqjnMVTVylizEG9yi2VX8wwdW8v1HX+OVhlY+891NvH31XK5aVjdua5fFiSQZG0amrQ0rHh/o2+oMXbWsjs07j7GtvoW1LzVxrd/u4G3zbiSR6WLT0Rf5xtb/4L5ld1MWKmxfRYZhMH12FdNnV3H8aBfb/3QIveMYLi6hsMVrje281tjOBXOrRq1Dxqp4hP/17ov5yoNbefSFRg4d7/Z7YR/5nZDBybW0186l7theeOQ7NDwCZixGqGYyVnk5gfJyAuUVWPE4ZrQEMxrBjESxolHMSBQjHMawLAzLAssk19LCgf/zeZwer6sOq6yMkkXnE6qbRmjaNMJ10wjWTMYoK2Pti0089EwD2aTD/PZy3rF03iltmtJZ7+p5NFTYA6BhGKgLpzBzXhV/enY/O7cd4YlHdlJdG2PZqhnMO6/mjG4UEEKcW6riET7+zmWsf/kwP19Xz0+e3M1TLx7iHWvmsWxBddH7SRPDkzZjQ3Bdl4b77iE4ZSqz/veQbe5Oq7Wzl3/6ry10p3J89I4lXOg3pnVchwf1wzx3eBOV4QruvvA9zI6PrHPZsfLQjkfZuP9F7rvyvTy2NsGL+jifvuviUemSIl+iJ8N3f7uTHfvaiEWDvGPNPK5YMnVEBwvXdfnkNzcwq2M/dy4IkG0+RubYMa82yz71sT8jVXXLW4m/6UqC1cPXJB1r6+EXTzfwkt/T+dy6OG9eMZ2LF9YQDFjse3wt2V/8kH0rbuKGe9951vG8Xp3tPWxev5+GXc24LpTFwyxePo2Fi2tPuXFjNBSivYU4c1Iu408hyiSRzPDrZ/fxzLbDOK5LXXUpN14yk0sX10pN2RCK3WZMkrEh2N3dNHzsfkqXXcS0+z/6uta151AHX/7pNgD++rYLWLbAe4al67o83riW3+19AsswuXH2tVw3a82Y90M2lEfqH+XJA0/zV4vv4ev/dYApk0r4xw9eMiaXuRzX5cktB3lkwz7SWZspVSXcuGoml55fS2iYS6JHWpP83Xc2ccmiydx760AHuq7r4iST5Do7sROd/k0EKZxUL05vypvuTeGk07i2DY6Da9u4tk2wuoba974f4wxqjuoPdfLYpka27WnBBSIhi4sX1jDtqGbuC7/m6OVvYfUH73g9m2hUdLaneGXLQXa9cpRczsEwYPqcKi5aNWNkHQEPwrYderoz9KayZDM2mXSOaDQEJlTXxobty04UliRj408hy+RwS5JHX2hk085j2I5LLBrk0vNruWLJVGae5ikeE40kY2dprJOx9MGDNH7us5RffQ21737v617fq/vb+NdfvkI25/C2K+Zwy+Wz+rs9eK11Nz967Wd0ZrqoLZnMLXOvZ1nNBZhGYX/UflX/O546sJ7V0Tt4/Jku3nXtAq5fOXSj/NHQ3pXm4Q17eX7HUWzH9ZIaVcPK8yajZlSe0pnh09ua+OEfNHfdoLj6omljGttINLf38MzLh9m8s5nWRC+qez9/dnQ9iaveyoq7bi92eP3SvVnqX2tm1/ajNB/uIhiyeN/9lxMcorPIVE+G9pYeOtp66GjtoaM9RbIrTbIrTapn6OeVBkMWdTPKmTa7kumzK6mqLpU2a0Ukydj4U4wyaUv08scXD7Fx+xES/v5bWxll6fxqLlpQzbxp5RO+xkySsbM01slYcvsrNH39q1S//R1U3Tw6XRTsPZzg249spzWRZnpNKXdeM5/Fs70ORVO5FL9p+AMbml7AxWVKaS1X1K1i5ZSLCnbH5a/2/I6nDq4n3nQ1LUcifPX+NxXsURttiV7WbW3i+VeP0pbw7li0TIP508qZP72cWbVlzJpSxsPr9/LCzmN8/u5V1I2jO+5c16WhKUHDUxuYs+GXlN52J9NuubnYYQ3q+XUNbNt0kKtuWsj5S+twHIfW5iTHmhIcPdzJsaYEiY7eUz4XCJqUxsKUxkKUloWJlAQJhQIEwxYVFSUc2N/G4QMddLQOPCorWhpk+qxKps6ooLo2xqSaUgKj2C9XLmvT3ZWmO5Em1ZMhncqRSmVJp7KkUllyGRvbdrBtF1yXyupSaqaUUV0bo7K6hEBgfPURNtokGRt/ilkmOdth+95WNu44yo69baSzXtOOcNBibl2cBdPLWTC9gpm1sQn3mKUJk4wppUzgW8BSIA3crbWuz1v+YeAeIAd8Xmv9u+HWN9bJWOf6Zzj2w+8z5UMfJn7Zm0Ztvd2pLD9fV89zfl9P02tKuXJpHcsX1DCpPMKxnuM8vn8tW45txXEdLMNiQcVcFk9SLKycz9TSWixzbH5AHtrzW9Ye3EDvjsu4eOYC/vq2sX2O5mAc16X+UCcvN7Swc187B451cXJBV8TCPHDf5eOyxqV760sc/rdvUHPnu6i8/sZihzOors5eHvzuZhzbpbK6hM721AkPKQ9HAkyui1M9uZSKqhIqJpX4j2gKDLnN8w9k3YleDjV20LS/nUP72+lJDnSkaBhQMamEqupSyiujlFdGiVdGKa+IUhILDbn+ZFealuZu2lqStDUn6Wjroauzd9hautMxDO+JBpP87xmviBKviFBWEaV0mFhOx7Ydujp7SXT0kurJ0NuTJdWTIdWTpbcnSzZrk834Q940QCgSIF4e6d8m8YrIiLbPUEbyA2PbDr09WTKZnH/ZuS+mHNmsjWO7uK53wuG6LuFIkNKyMLF4mFhZmFARuvhwXRfHcXH8RLt/7DgYhkG0NERwnHXG22e8JMjZnM2uAx1sq29h94EOmlqSJywvj4WYURNjek2MmsooNeURaiqiTCqPvCFr0YqdjBVyL7oNiGitL1NKXQo8ANwKoJSaAnwEWAFEgGeVUk9qrU/TqdPYicybR8XyiyhZtHhU1xuLBvngzYu4dvl0HtvUyJ92Heenf9zDT/+4h6mTSpg7Nc7sqW/ivOmXcdjezc7Edna172FX+x4AgmaQ6bE66mK1VEcmUV0yiUmRSmLBGLFQKWHr7M9m3Ly0Z/WSqa/7u54N0zBYOKPCu0txjZe8Nh7r4sDRLvYf7aKtq5c3XzJrXCZiAPgHqVxXF3ZPEgzTizVvOPk1fSdE3i+eP+kCrvfw9vzluP5o8OVu/rwh3h/B5bq3nsfGp/fR2Z4iXhFlyrQ4tXVxaqeVU1EVfV3bNxaPcN6FUzjvwim4rkt7Sw/HDidobe7m+LFuWpu7aW/pOeVzhgElpSGipSFKY9443Zuj+UiCZNeJPWOblkFZPEJVTSll8QixeJhoaYhINOgPASLRIMGQhWmZWJaJ67i0tSRpPtJF6/FuWpuTtB3vpr311Fgsy+iPJVoa8qeDRKMhgiHLG4IWhmnQ1Zmisy1FR1uKjvYeujp6cZzhzxUtyyAYChAMWZSWhfsTh95UlpZj3TQfOfVHIRAwKSuPeDGVBImWBImUeNN93zUYsgj1xRcKkC7LkehI9dce9l1qzn+dnyyfjVDYIhaPECsLezWm0SDhSCBv8F6HwhZWwCIQMAkETKyAt2+4rksmnaO7K01Pd6Y/xmR3hp5kpr9dYiZjk+0bZ05/s04wZHnbyd9eJaUhoiV+OfrbLRQOEAhaBIMmwZBFIGgN2eYxl7XpSXox9bWZzKRzZNID8WXSXkLr1cY6ODkX23GwcwNJI3gnPCWlIUpiYS+OSICIv53C0YFtFgxaBIImgaC33UbzuBcMWFw4d1L/TWXdqSz1hzppONzJoeZuDh3vZse+NnbsO7HPRAMvUYuXekN5SYiy0hDxkhCxaJBo2CISDhANBYiGLa9PvlCAUHB043+jKWTN2FeBzVrrB/3XTVrraf7024Cbtdb3+q8fBr6otd4y1PoK0QN/ITLlju40W3cfZ+ueFuqbOukd5CATLc0SqmrDKuvADreTCXSCMfjXNwkQIoJFEMsIYhEgYASwjCABI4CJhYEBrgkYGK4BmOAatDqH6DKaiexfw5fff9MZPyC7UMbLmeVgkq/uoOlfvlLsMEYkfuXqERwcR/Y/EImG6O093Y+6ty4Xg5RtkbQDJHMBunNBkrkAvbZFr2PRa5vY7sAPYsSyqQxlqAjnKA/lKA9liQXs/v7wvNX2TZ8U72m+n+tCT86kO2ORzFkksxbJrEky68XRmzNxR7gNAEKmQ1nIJhayiQVtIpZDOOASthx/cAmYLsPtWn0xJbMW3RmT7qxFd3YgtqwzOrUSpuESDTiUBBwiAYeg6cUWMN0Tpk0jr+QMyNgmPVmTnpxJKjcwnTuLuEz/OOa4w29j0xiIqW9sGd5886Sx4xr02gbpnEmvbZLOGWdUhqbhYhkuAZP+01Pb4Yy3+2CxWYbX2W0q623HM+XF5WKZ3nT++g0DTPKmB5mPX5Z9W8M4abcZmO9tMcd1SWdssjmHnO2NM1mbnO2Qs11OzR/cvKlTt7n3v2Rgml53PN75qOHHaQyco/qfNYyB3fvEeUb/d8mPP39/z//rJxwiBtmupmmw5s+vpGxy9SBLR894qRmLA515r22lVEBrnRtkWRcwbH8KlZWFae9RUzO2d5zU1JSxYE41d94AtuNyqLmLfU2dHG3r4UhLkpaOFIlkhq5EGV1HppLJOWDYGOEURqQHI9yDGU5BIIsRyOAEM+QCGQwzDZaNYTqccp1vKP6/ydsuX0RtbXzMvvNoGOtyOVtVq5aRu+E6sh2dgIvr+DVUruNVTjmOdwDzh76DWX9tmfdiIEnqr00D7+gzcHQyTjhSnfyevOWGccL7jz+9HoDEhvWj9r07T/+WU4T84eR7Ol3ANoKkA1Es1yacS55wAM0BHWcZ59nEkjNDZKwoGStC1opgmwFsI4BtBnEMk0ium5JMFyXZBEFn6Mr8rD+MVACo8Id8DgZZK9IfjzcOYxtBbDNAzgxi9w1GgKCdJpxLEskl/XEP4VySoJM+gxRlQATvgH2ynBEkHSgha4XImeG8cZisGfK3l4VtBHBMq38aIGz3EM6lCOeS3tj2YgzZvQTsDCbOIH9xZE4uw0zAH1uR/m3klam3/frmOabln/O6hF2bMjtFKJciZKcI2ymCuV4CToaAk/WHDJaTxXJzmK592m3rYJCxomStyEnbLETWCpMzQ/2xOf0xBrztZgbJmgEcTFxjYCgIyx/OwhmXonvSeAytffA57v6n94z9HxpCIZOxBJD/C2r6idhgy8o4zfG2vf3USwujrRg1MCWWweKZFSyeOfgDkrM5xxtsx2tvYjtksg624/2wO+5A+w7HcbFdl6yTIWNnyTgZHMfBNfzLWYYDfWN/XnVpnCXT5ozbmicY3zVjAOV3vLvYIQyr4i/eT/Z48wgOcCM8AroulVWltLclh3vLievMf+0OMg9OPOse4j34/+8nzxtTY7X+MVhtRWUJHe09Y7NyGKRmZJwby3hHuO6KihI6Okb/98v19wXH9Z776bjgOi6O49VwOc5AmC7esv5mDHnz6R+7p7SSGJg+9buOeMsO8kbXX6fjurhO3mtnIIa+Note7F6zGtdxcfLicU9af/98d5BlDCwzLZM333lbIdqMDbmskMnYc8BbgZ/7bca25y3bDHxBKRUBwsAiYEcBYztnBAMmwcAbr/GkKBzDNAnVThnVdZbUlJEMj98EeaKK15SRHscnLhNReU0ZGSmTcae8MlbUk/xCJmMPA9cppTbiXTz5gFLq40C91vo3SqlvABsAE/g7rfWp99YLIYQQQrzBFCwZ01o7wL0nzd6Vt/w7wHcKFY8QQgghxHgg17uEEEIIIYpIkjEhhBBCiCKSZEwIIYQQoogkGRNCCCGEKCJJxoQQQgghikiSMSGEEEKIIpJkTAghhBCiiCQZE0IIIYQoIuOce7aYEEIIIcQbiNSMCSGEEEIUkSRjQgghhBBFJMmYEEIIIUQRSTImhBBCCFFEkowJIYQQQhSRJGNCCCGEEEUUKHYA45FSygS+BSwF0sDdWuv64kY1cSmltgKd/st9wH8AXwdywBNa688VK7aJRim1CviS1oX/epoAAAYJSURBVHqNUmo+8APABXYA92mtHaXU3wNvwSufj2mtNxct4AnipHJZDvwW2OMv/rbW+mdSLoWhlAoC/wnMBsLA54GdyL5SVEOUyyHGyb4iydjgbgMiWuvLlFKXAg8AtxY5pglJKRUB0FqvyZu3Dbgd2Av8Xim1XGv9UnEinDiUUp8E7gKS/qyvAp/RWj+tlPp34FalVCNwFbAKmAE8BKwsRrwTxSDlshz4qtb6gbz3LEfKpVDeA7Rqre9SSk0CtgLbkH2l2AYrl39knOwrcplycFcAfwDQWr8ArChuOBPaUqBEKfWEUmqtUmo1ENZaN2itXeBx4NrihjhhNABvz3t9MfCMP/0Y8Ga8fecJrbWrtT4ABJRSNYUNc8IZrFzeopRar5T6nlKqDCmXQvoF8Nm81zlkXxkPhiqXcbGvSDI2uDgDl8UAbKWU1CIWRw/wFeAG4F7g+/68Pl1AeRHimnC01g8B2bxZhp8Qw0A5nLzvSPmMsUHKZTPwt1rr1Xi1x3+PlEvBaK27tdZd/g/7L4HPIPtK0Q1RLuNmX5FkbHAJoCzvtam1zhUrmAluN/Bj/yxlN95OUpW3vAzoKEpkwsmb7iuHk/cdKZ/Ce1hr/WLfNHARUi4FpZSaAawDfqS1/m9kXxkXBimXcbOvSDI2uOeAmwH8NmPbixvOhPZBvDZ7KKXqgBIgqZSap5Qy8GrMNhQxvolsq1JqjT99E145PAfcoJQylVIz8U5kWooV4AT1uFLqEn/6WuBFpFwKRilVCzwBfEpr/Z/+bNlXimyIchk3+4pcehvcw8B1SqmNgAF8oMjxTGTfA36glHoW706kD+KdZf4EsPCu7W8qYnwT2SeA7yilQsBrwC+11rZSagPwPN7J3n3FDHCC+ivgm0qpDHAU+EutdULKpWA+DVQCn1VK9bVR+ijwDdlXimqwcvk48LXxsK8Yruue/l1CCCGEEGJMyGVKIYQQQogikmRMCCGEEKKIJBkTQgghhCgiScaEEEIIIYpIkjEhhBBCiCKSri2EEOccpdQPgPcN85ZGf/xdrfXnxziW9wDXa63f6z+p436t9df8ZYvxHhB9mXQcLYQYitSMCSHORR8FpvpDX6eNt+bNW+kP/zKWQfjPrPsSXh9GAO/M/5ta61fxHhL9qbGMQwhxbpOaMSHEOUdr3Yn//DilVMSf3aa1PlrgUD4BPKW1PuS/NgZ5zwPAZqXUN/24hRDiBJKMCSHekJRS+/EvUyql/gG4FFgL/E8gCvwI+L/AvwNrgEPAR7TWf/A/Hwa+CPwFUApsxXuUygv+8ihwD3CX/3qNv06UUi7wAa31D7TWu5RSjcCH8R56L4QQJ5DLlEKIieJqYAlwBfARvMcGbcJ7tNbFgMZr39Xnh8Bq4E5gBV4it04ptdBfvhqIAU/6rzcC9/vTU4Gf5a3r98DbRvXbCCHeMKRmTAgxURjAPVrrJLBbKfXPeM82/QmAUupbwKN+O7ByvCTsAr/dF8DnlFJX4F2avAdYBdRrrdMAWuuMUqrTnz75cukO4H+M7dcTQpyrJBkTQkwUR/xErE8SaMh7nfLHYeAif3qTUip/HWF/AKgFWkb4t48DIaVUpda6/YyiFkK84UkyJoSYKLKDzHOGeG/GH1/GQJLWJ+2PXUbe1KPvWGuP8P1CiAlEkjEhhDhV36XJWq31H/tmKqX+DXgN+CZwBLj+pM+5Q6yvGujRWidGO1AhxLlPGvALIcRJtNb1eA3w/59S6ial1Dyl1BeAe/GSMYDNwDylVDzvo10ASqkVSqlY3vzleDcLCCHEKSQZE0KIwd0NPAp8H68B/k3A27XWT/nLnwG68e6q7LPOn78R+Mu8+WuAX49xvEKIc5ThukPVqgshhBiOf0fmHK31HcO8ZxlegjZLa91RsOCEEOcMqRkTQoiz92XgKqXU7GHe81HgAUnEhBBDkWRMCCHOktb6OPC3wBcGW66UugBYitfTvxBCDEouUwohhBBCFJHUjAkhhBBCFJEkY0IIIYQQRSTJmBBCCCFEEUkyJoQQQghRRJKMCSGEEEIUkSRjQgghhBBF9P8Bfzy95yuX7OgAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plt.figure(figsize = [10,5])\n", "plt.plot(tnow,F1,label = \"F1\")\n", "plt.plot(tnow,F2,label = \"F2\")\n", "plt.plot(tnow,F3,label=\"F3\")\n", "plt.plot(tnow,F4,label=\"F4\")\n", "plt.xlabel(\"Time(t)\",fontsize=15)\n", "plt.ylabel(\"states\",fontsize=15)\n", "#plt.xticks([i*30 for i in range(10)])\n", "plt.legend(fontsize = 8)\n", "plt.savefig(\"plot/states.png\")\n", "plt.savefig(\"plot/states.pdf\")\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.6.5" } }, "nbformat": 4, "nbformat_minor": 2 }