{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# charging of capacitor" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "\\begin{equation}\n", "\\frac{Q}{C}+R\\frac{dQ}{dt}=V\n", "\\end{equation}" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "\\begin{equation}\n", "\\frac{dQ}{dt}=\\frac{V}{R}-\\frac{Q}{CR}\n", "\\end{equation}" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "from scipy.integrate import odeint\n", "import matplotlib.pyplot as plt" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "# function that returns dy/dt\n", "def model(y,t,C,R,V):\n", " dydt = V/R-y/(C*R)\n", " return dydt" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "V=10 #V\n", "C=1e-6 #F\n", "R1=100 #ohm\n", "R2 = 200\n", "R3=400" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "# time points\n", "t = np.arange(0,0.002,0.0001) #" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAaEAAAEgCAYAAAAdRMy/AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAgAElEQVR4nOydd3xVRfbAvycdQmih19BEigQhgNhogsiigCCLFdcCKq67WFbd36qLHRdXXVGRVdcusDZQREQRV1GBICBNpENAagiEkvrm98fcvBTeS30tyfl+Pvfz7p07c+bcm3LezJw5R4wxKIqiKEowCAu2AoqiKEr1RY2QoiiKEjTUCCmKoihBQ42QoiiKEjTUCCmKoihBIyLYClQmGjRoYBISEoKthqIoSqVi5cqVh4wxDT3dUyNUBhISEkhOTg62GoqiKJUKEdnp7Z5OxymKoihBQ42QoiiKEjTUCCmKoihBQ42QoiiKEjTUCCmKoihBw69GSESGisgmEdkiIvd5uB8tIrOd+8tEJKHAvfud8k0icnFJMkWkjSNjsyMzyim/UER+EpEcERlTpP/xTv3NIjLeH+9AURRF8Y7fjJCIhAMvAJcAnYErRaRzkWo3AkeMMe2BZ4CpTtvOwDigCzAUeFFEwkuQORV4xhjTATjiyAbYBVwPvFtEv/rAQ0AfoDfwkIjU883TK4qiKKXBn/uEegNbjDHbAERkFjAC2FCgzgjg7875+8B0ERGnfJYxJhPYLiJbHHl4kikiG4GBwFVOnTccuS8ZY3Y4dV1F9LsYWGSMSXXuL8IavPcq/ORK4MnOhuPH84/09NPPs7JsvaJHTo7n8qL3XK7ChzFlKzOm8AFlv87Dz9fGZcghgkwTRQYxZLiiyCaS1mG7CRNbd1NOO1JddckhgmwiyDaRZBNJjgmnafh++kSuAuCoK45ZGSPJIYJcwsg14eSSd4RxVcxHJITvBmBB5kCWZJ3rvpdfN4z4sCM8VutJt463HpvKcVMTF2G4TBgGseeEcXXMB4yKWQDAt1l9eOLEHbgQjIfj4zrXExd2AoDJ6VNIzk4sdN9FGMYIg6K+5fG4JwBIyW3KyLTXMYBB7Dtz6gPMiPsL50T9BMAzJybw6qmrQArUNfazRfhvfFF/nPuZeh/+jGOuOPJ+GgVl3xP7EhNqvu1+T3889mj+z8upl8dP8UOoE5YOwDVp0/kuqzdFqgAwLPorXqzzVwB25TbngsMfnV4JaBh2mORO18GaNR7vVwR/GqHmwO4C1ynYUYfHOsaYHBE5CsQ75T8WadvcOfckMx5IM8bkeKhfFv1OayMiE4AJAK1atSpBpOIz0tJgy5b8Y9s2W+bNwGRmBlvjgHOKGFKpTzpxpBPHMWq7zyPI4UpmuetOZAZp1CWDGDKJtobFOf8jz3MD/wHgAy7nOt4kgxhchJ/WZyr1qEcaALfxDosZ5FG3K5jDHD4H4BCx3MI/vD5Hr+Nfk4DdBL6ERJ7ido/12rCNx05Mdl/P4jLS8Dx5kZT5HRzdD8A+oljgRU+A7INHsJMnsJZ2fMc5Huu1ytkKJ63MTGJZSaJXmcePZAF5/ceynjM91svIjYB9+9zXm0nw+kxHjoXBMVv3OFlspY3X/l0HDoLzc9pHHXbS0mO9gydj4eRvAByLjGJXTBRkxkFODVuh5kGI/5UTp1ywf7/X/iqCP42QB7tL0Qx63up4K/c0fVhc/eIoVRtjzExgJkBSUpJmAPQlqan5Rmbz5sJG59ChsskKC4O4OKhVyx5Fz2NjIToaIiMLHxERp5d5ux8ebvsRsZ9FD0/leWV5n2DP844i18dPCPO/qcX+wxH2OBTO/sPhHDgcztH0cF558hDn97IGd8rUekydUcfj62jVPIcrf3jKff1h92YcSj3dqADsvedZ+NMjIELYghhO3hQLQHi4ISbaEBMN0dGGqEiD6/NfoL79M+j8f3GcWJ1FZARERBgiIyAy0n4m9boEJtl/bnVThZsfP0lEhJUZHmZfZd55wpXvQttcAC75Pop6K9OL1LN168TVh8t/c+v94kdCTk6ah1dv6NrxLujwZwDO3x/Gp2uPECam0KvPO+L6boBIK/OfGyJIO3rYfS9PtgjE1+sPbW3/LTJhxcZDhX58Bc/bJbwDtex7mnwgjGsOHSxUJ+8zMqIWtN3rfqZlW7NxuQ647+e4ckjNOUSOK5vE5pOhnn2miF+X8PShlxjUcCixEfbntfjgQpYf+YGBDYdQO3EDhMPKgz9z6tuHuL7mbB7qnD9yOm9JT9Jz0nl62BpovAeAB7+4GXY25sXu/2ZY00sB+GDPHO76+Q76NR8Dw1d7/P2pKP40QilQyPy2APZ6qZMiIhFAHSC1hLaeyg8BdUUkwhkNeerLk379i8haUkIbpTykp8P8+bBxY2Gjc+SI9zY1akD79vbo0AHatoUGDU43MnmfMTH5f90hyBdfwLp19stk0aNvX3j/fVsv/TcYN9m7nIMRTSHBnjdoB02a2FdQu3bhzyZNIqBl/p/KCy/ZGcGYGGuLY2Lyj+bNa0Oj2gAMvw6Oj7N1IiKE07+rNXafPf9qcU8cA8QBEN8EZr5dXN1a7rP+l9ujNFx5a+nqNW0Cv/M+aClEtyalqxcNJLU+vdwYQ7Yrm6jwKACyc7P5KX0Rp5qcYnTn0e56zy97np/3/8zd594NDToC8J9V/+GJVU9wfffr+esFdopsw8ENXPhiIh3jO/LL4F/c7e99/yE2p25m0+2baB/fCIAZX6zh1dUvceYZbQlvMRKAE1mb+f7Q94S3CiehT/7PLn3xMdJz0qnTqR7E2C8y8fWb0eBgAxp1qEfrzrbuWVvPpG9qX7q162p/2fyAP43QCqCDiLQB9mAdDa4qUmceMB74ARgDLDbGGBGZB7wrIv8EmgEdgOXYv4jTZDptvnZkzHJkzi1Bv4XA4wWcEYYA91fkgZUCGAM//givvAKzZ8OJE6fXiY3NNzJ5BifvumnTkDYqBTlxwtrVX3+1tnXzZnv+3/9Cs2a2ziuv2GtPpKTknzdoAKNHQ+PGhY9GjaBevUJ2hbvvtkdpGDu2dPXyBn1KPp9v+Zx9x/dxZdcriY6IBmBG8gwWbl3I7b1uZ1BbO9X34cYPGfvfsYzuPJrZY2YDdiTzu3d/R1R4FJmd86eMP/n1ExZtW8TozqPp6Bihk9kn2Zy6mZRj+b8QdaLr0KRWExrUbFBIp2EdhnHgxAFiImLcZb/r8Dua1GpCv9b93GXdGndj8XWLT2u/+Y+biQqPonZ0bXfZvy/7N//m34XqDW43mMHtBpf9pZUBvxkhZ43nduw/+3DgNWPMehF5GEg2xswDXgXechwPUrFGBafeHKwTQw4wyRiTC+BJptPlvcAsEXkUWOXIRkR6AR8B9YBLRWSKMaaLMSZVRB7BGkuAh/OcFJQKcOgQvPWW/a+7oYAPygUX2KOg0WncuNIYmsxM2LrVnnd2/DHXr4eLL4Y9ezy32bQp3wiNGmXtalHjkmdg8oiMzB8VKb4jx5WDIISH2SnJn377if/t/B89mvbgwtYXArB2/1pGzh5Jm7pt+PK6L91tx388ngMnDjCk3RCaxdkf6Op9q/n4l48Z3Haw2wjViKhBrsnlZPZJd9uYiBiGdRhGraha5Lpy3f1P6jWJ0Z1G07lhvsPwuK7jGNR2EI1i838hmtduzm935U9B5vHs0GdPKxvQZgAD2gwoVFY3pu5pZQCNazU+rSxYiCnqJaN4JSkpyWgUbQ+4XPDVV9bwfPyx9UID+9/1+uvhhhugY8egqlhWjh2DpUvhm2/skZxsHeVGjLCPCHDwoH3EqCg7W3jGGda+duhgz3v2tNNjSmDIzMlk9vrZpGWkcUefO9zll7xzCQu3LOS7G77j3JbnAvDU0qe498t7uavvXUwbMg2AzYc3c8b0M2hTtw3b/rTN3X7iJxM5mXOSaYOnuf95J+9NZmfaTno07UGbetZBIMeVg8u43FNxSj4istIYk+TpnqZyUMpPSgq8/jq8+irs2GHLwsJg2DC46SYYPrxSzu3ceSc895y1rXmEhdnBW97IBuzU2bZt0KqVXTxX/MPOtJ2sO7COjg060r5+ewC+3v41kz6bxDktzuG1Ea8BICKM/3g84RLObb1uIyLM/nuLDrdTaEdO5a9B9mnehz/2/iMXtLrAXZZQN4Ffb/+VhrGF0968fOnLp+mU1CyJpGaF/6fm9aeUDX1rStnIzrZOBq+8AgsW5P+nbt0abrzRjnxaenYHDSUOHoT//S9/pPPcc9C/v73XrJk1On36wIUXQr9+cN55p49qRKCNdy9ZpRy8sPwF1uxfw+ODHnevYzy19CleTH6Rfw75J5P7Wq+NiLAINh7aSN2Yuu62UeFRTOw5kdrRtcnKzXIbhTdGvkHNyJpEhud/IeqX0I9+Cf0oSGR4JB3iO/j7EZUiqBFSSsfmzXbE8/rr+fsFIiPtKvrNN8OgQfkuyCGIy2UdA/KMTsHlKoAlS/KN0IQJcOut1m9C8Q2pp1I5mX2SFrVbAHDgxAFGzR5FVm4WK25e4a73+prXSd6bzPXdr3cboV7NezE4dTDNa+dv4+vRtAdrblnjlpfHjOEzTuu7ToxnN3YlNNA1oTJQbdeEZs6EW27J31XfqZOdbrv2WmjoMWNvSGBMYb+H1q1h1y57XqOGdY3OG+n06WPLlIqRmZPJxkMbaRTbyL2I/+7ad7n6w6u5ovMVzLliDgBZuVnUfKwmLuPi1P+dcnudvbnmTY5lHmPUmaMKGR2lcqNrQkr5eeedfAN03XUwcaL97x3CXm27dtnZwnfegWXL7NoNwB13WC+3fv2gVy/rUKCUn11Hd7Fm3xr6J/QnLtruCbpz4Z28mPwi0wZP465z7wLgjPgzqBlZs9CaSVR4FN/+4Vta1G5RaCH/usTrAvsQStBRI6R4Z+5cGD/eGqAnn4R77w22Rl7JzYXPP4cZM+Czz/KXqj74wNpNgLvuCp5+lZkcVw7LUpaxJXUL47vnB5sfM2cMK/auYMn4Je71le5NunNG/BnukQ3YqbP0+9MJk8LTtX1b9g3MAyghjRohxTNffml3OObmwl//GrIGyBh4/HE7Y5g31RYZCVdcYQdw/foV3145nb3pe9l3fB89mvYA7K7//m/0J8eVw+jOo6kVZSMc9Gvdjzoxddx7XwBu7nkzN/e8uZC8osZHUQqiRkg5nR9+sBtisrLg9tvh0UdLbhNAXK7Csb+WLLEGqG1b61Twhz8U3gCqeCfHlUN6Zjr1atjAISv2rKD3K73p3LAz62+z+8BrRNZgXNdxRIZFciLrhNsI/WOI96CkilJa9CuKUpjVq+0+n5Mn7RrQc8+FzPrPgQMwdardDPrDD/nlU6bY2GybN9sBmxog7xR0RPpgwwfUfbIukxfmB6vr1rgbDWs2pGXtluS4ctzlb416i9dGvBZSO+2VqoGOhJR8Nm2CIUNsyoTLL7cu2SHgdr19u50R/OADu00JYNYsONdufnd/Kt6ZkTyD55c/z+29bufWXjbyZ+u6rTmRfYK96fmxfqMjotl/934kRL54KFWf4P+HUUKDnTvhoovsLs4hQ+Ddd20agyDicsHzz0PXrtbo5ObCpZfCp5/CM88EVbWQZtOhTTz9/dNsTd3qLhOEDQc38ENK/hCye5Pu7L97P19c+0Wh9mqAlECiIyEFfvvNbjZNSYHzz4cPP7Sx/IPMI4/A3/9uz8eNg6eeqhTBGALOyeyTRIdHux0EHv/ucd5c8yYGY1MFAGM6jyGxSSJnNznb3S4iLKJQsExFCQY6EqrupKbakc/WrdCjhx1mhEiogNtug7POgo8+gvfeUwPkiYmfTKT+1Pos37PcXTam0xiu7Xat27sNIL5mPOe0OKeQ67SihAJqhKoz6ekwdKjNttapk91oUyd4IU7WrrWh5/KCcDdsaFPajxwZNJVCikVbFzH588mcyj7lLouOiCYrN4tV+1a5yy7teClvjnqTgW0GBkNNRSkTaoSqK6dO2QWWFStsFM5Fi4IWgicrCx5+2KY+eOMNuw6UR3VenigaUuv+r+7n2WXPsmTHEnfZvefdy76793Fbr9sCrJ2i+AZdE6qOZGXZ3ZzffGMzrX35JTQPTpyulSttuqGff7bXt95q46FWZ4wxPPD1A8xZP4flNy93R4qe0HMCO9N2utMZABpfTan0qBGqbuTm2sCj8+dDfLw1QG3bBlyNjAw7+nnqKatS27bWIzwvknV141jmMXeqZRHhx5Qf2Zy6mbm/zHWHypnQc0IwVVQUv6DTcdUJY2wgtTlzIC4OFi7Mz1UdYD74AJ54wrphT55sR0LV0QAZY7jqg6to+I+G7Ejb4S5/eMDDfHP9N1ybeG3wlFOUAKAjoeqCMTZl6Kuv2pwF8+fbRZgAq5C3xnPVVfDttzY+at9qFMfSGMPK31bSs2lPRMR95LhyWLprKQl1EwDcaagVpaqjI6HqwpQp8OyzNrrnhx/CBReU3MaHfP01JCZaT3CwxmjGjOplgAD6vd6PXv/uRfLe/LxUjw54lD137uHqblcHUTNFCQ5qhKoD8+ZZIxQWZjfcDB0a0O5nzoSBA60L9j+qWczLn/f/jMu43Nc9m/akSa0m7Enf4y5rU68NTWo1CYZ6ihJ01AhVB/79b/v5yCM2HXcA+e9/bUoFgAcegH/9K6DdB5Ub5t5A4oxE5v863102ZcAUdk/ezcgzdfOTooCuCVV9UlOtA0JYGNx4Y0C7/uILuPpquxb02GM2CGlVxhiDy7jc4XPOanQWsZGxpBxLcdfJ84BTFMWiI6GqTl7o6UGDoHHgwvD/+COMGmW7njwZ7r8/YF0HhQWbF5A4I5E31rzhLpuYNJGdf97pjlqtKMrpqBGq6rz3nv288sqAdvvTTzYl0fjxMG1a1Y98cPjUYdYeWMvbP7/tLqsZWZP4mvFB1EpRQh+djqvK7Nlj045GRdlhSQC57TabfG7AgJBISeRT0jPTeXnly9SOru3eQDqu6ziMMYztMjbI2ilK5UKNUFVmzhy7IDNsGNSt6/fuDhyAY8egvRNVZvBgv3cZFL7f/T33LLqHxrGNubbbtdSIrEFEWIRuLFWUclDFvqMqhQjgVNzRo9bz+7zzrCt2VeLwycN8sukT9/WQdkO4ucfN/GfEf4iJiAmiZopS+dGRUFVlyxYbIbtWLRg+3K9dnToFl10Gq1bZUVCjKpQnbW/6XhJnJHIq+xQ7/7yT+JrxiAgzL50ZbNUUpUqgRqiqMmuW/Rw5EmrW9Fs3OTk26+n//gfNmtmMEAF0wvM7zeKacV7L8ziWeYyT2SeJRx0NFMWX6HRcVcSYgEzFuVxw0002IEO9enZfUEKC37oLCOmZ6fxl0V8K7e15+/K3+eq6r2hZR1O7Koqv8asREpGhIrJJRLaIyH0e7keLyGzn/jIRSShw736nfJOIXFySTBFp48jY7MiMKq4PEYkUkTdEZK2IbBSRqrOTZe1a2LDBpmrwo3fAX/5ik9DVrAmffQZduvitq4Bx1xd38Y/v/8GdC+90l9WKqoVUdR9zRQkSfjNCIhIOvABcAnQGrhSRonkDbgSOGGPaA88AU522nYFxQBdgKPCiiISXIHMq8IwxpgNwxJHttQ/gCiDaGHMW0BOYWNAIVmrefdd+jhljA5b6iXbtIDoaPvoIzjnHb934nYIZTB+48AEGJAzg3vPuDaJGilJ98OdIqDewxRizzRiTBcwCRhSpMwLI22L+PjBI7FfOEcAsY0ymMWY7sMWR51Gm02agIwNH5sgS+jBArIhEADWALOCY7x4/SBiTvx7kZ6+4W2+1UbGHDPFrN37jVPYpHvz6Qa768Cp3Wcs6LVk8fjE9mwU2zYWiVFf8aYSaA7sLXKc4ZR7rGGNygKNAfDFtvZXHA2mOjKJ9eevjfeAE8BuwC5hmjEkt+hAiMkFEkkUk+eDBg6V99uDxww+wc6dN1+2HdA1z59qZvjyClBXcJ6SeSuWZH59h1rpZrNm3JtjqKEq1xJ/ecZ4m0U0p63gr92Q0i6tfXB+9gVygGVAP+FZEvjTGbCtU0ZiZwEyApKSkovqHHnkOCePG+TxUwaJFcMUVNinrzz9XTgN08MRBGsY2BKB57ebM+N0MWtdtTWKTxCBrpijVE3+OhFKAgu5ELYC93uo402J1gNRi2norPwTUdWQU7ctbH1cBnxtjso0xB4ClQFI5nzU0yMmxURLA51Nxy5blByQdP966Y1c2Xlj+AgnPJTBv0zx32dXdrub8VucHUStFqd740witADo4XmtRWEeDeUXqzAPGO+djgMXGrhLPA8Y5nm1tgA7Acm8ynTZfOzJwZM4toY9dwECxxALnAL/48PkDz9df29g5HTpAjx4+E7tli438c+JE5Q5ImuPK4WT2SRZvXxxsVRRFcfDbdJwxJkdEbgcWAuHAa8aY9SLyMJBsjJkHvAq8JSJbsKOTcU7b9SIyB9gA5ACTjDG5AJ5kOl3eC8wSkUeBVY5svPWB9bL7D7AOO2X3H2PMz356HYGh4N4gH1qJP//ZpiUaPhxeeaVyBSQ9mnGUOjF1AJjUexKJTRLpn9A/uEopiuJGCrqnKsWTlJRkkpOTg62GZzIybKiCY8dg40Y480yfiP3qK7joIhv9Z8uWyhMNISs3i7u/uJtPf/2U5AnJ1K9RP9gqKUq1RURWGmM8LndUou+0SrEsWGAN0Nln+8wAAezaBbGxNildZTFAAC7j4vvd35NyLIUfdv8QbHUURfGCxo6rKvgpTM8f/mCjYwcgE4RPiYmI4YOxH/Db8d84p0Ul3kmrKFUcNUJVgfR0+MRJNfD73/tcfNOmPhfpc4wxPLfsOX5L/42pg21QjNZ1W9O6busga6YoSnGoEaoKzJ1r14TOPx9atfKJyOeesyF5broJIirBb8mmw5u4Z9E95LhyuLrb1XRr3C3YKimKUgoqwb8XpUR8PBW3dy/89a9w8iQkJkLfvj4R61fObHAm0y+ZTnzNeDVAilKJUCNU2Tl0yOZQCA+34Qx8wN/+Zg3QqFGhbYDm/jKX5rWbk9TMOt1MTJoYZI0URSkr6h1X2Xn/fRspYfBgaNiwwuJWr4bXX7dTcFOnllg9aMzbNI+Rs0cyes5o0jLSgq2OoijlREdClR0fTsUZA3ffbT8nTbKBF0KVi9tdzLktz2Vkx5HUia4TbHUURSknulm1DITcZtWUFOuIEB0N+/dD7doVEjd/vo2KULeuTdFQP8T2d647sI4z4s8gKjwKsGF4IsL0e5SihDq6WbWqMnu2Hbb87ncVNkAA//mP/XzwwdAzQLPWzSJpZhL3fHGPu0wNkKJUfvSvuDLjY6+4WbNsUtZx40quG2gS6ibgMi4yczNxGRdhot+fFKUqoNNxZSCkpuN+/RU6drQjoH37oEaNYGvkd349/CtnxJ8RbDUURSkjOh1XFckbBY0aVWED9MEHdm9QqDHt+2ms3LvSfa0GSFGqHjodVxkxxmdTcTt2wNVXW5fsrVtDJ0jpf9f/l3sW3UN8jXi23rHVnY5BUZSqhRqhysjq1bBpk90XNGhQhUT99a+QmQmXXx46BghgxJkjGNN5DCM6jlADpChVGDVClZG8UdAVV1QosNvy5VZUdDQ88YSPdKsgxhhEhKjwKOaMmYNUxhSuiqKUGl0Tqmy4XNaNDSo0FWcM3HmnPZ88GVqHQLDpqd9N5Q9z/0CuKxdADZCiVAPUCFU2vv8edu+2m1TPPbfcYj78EJYutTN699/vQ/3KScqxFB7+38O8ueZNlu5eGmx1FEUJEDodV9l49137OW4chJXvO0Rubr7hmTLFJ/tcK0yL2i347KrP2Hl0Jxe2vjDY6iiKEiDUCFUmsrPhv/+15xWYigsPtzN6L7wAN9/sI93KyYETB2gU2wiAfgn9gquMoigBR6fjKhNffWVTN5x5pk30UwF69IBXXw1uwrqnv3+ajtM7krw3RDYAK4oScNQIVSYK7g0q56L95s0+1KcCuIyL71O+Jy0jjXUH1gVbHUVRgoQaocrCqVPw0Uf2vJxTcZs3Q+fOcOmlNgVRMAmTMGaNnsXCaxZyfffrg6uMoihBQ41QZeGzzyA9HZKSyp3o5957rfFp1Ch403ALtyx0u2BHhkcypN2Q4CiiKEpIoEaoslDBMD3ffGMHUjVrwiOP+FCvMvDC8hcY+s5Qbph3Axo4V1EUUCNUOTh6FD791K4D/f73ZW7ucsFdd9nzv/wFmjXzsX6lpFvjbtSKqsW5Lc7VjaiKogDqol05+PRTG+CtXz9o3rzMzT/6CFautMbn7rv9oF8puaD1BWy9Y6vbJVtRFEVHQpWBtWvt58CB5Wr+5pv28557IDbWRzqVknUH1rHh4Ab3tRogRVEKokaoMrB9u/1s27bMTY2BhATrjFCOmbwKkZaRxqjZo+j97978sPuHwHauKEqlQI1QZWDHDvvZpk2Zm4rAc8/ZpHVNm/pWrZKIDIukT/M+dIjvQPcm3QPbuaIolQK/GiERGSoim0Rki4jc5+F+tIjMdu4vE5GEAvfud8o3icjFJckUkTaOjM2OzKhS9NFNRH4QkfUislZEYvzzJipI3kioHEYoj/BwH+lSBmKjYnlr1FssGb+EGpFVP/24oihlx29GSETCgReAS4DOwJUi0rlItRuBI8aY9sAzwFSnbWdgHNAFGAq8KCLhJcicCjxjjOkAHHFkF9dHBPA2cIsxpgvQH8j26UvwBcePw8GDNulPkyZlapqaateDjh71k25e+PXwr4XSMWhSOkVRvOHPkVBvYIsxZpsxJguYBYwoUmcE8IZz/j4wSKzv7ghgljEm0xizHdjiyPMo02kz0JGBI3NkCX0MAX42xqwBMMYcNsbk+vD5fUPeVFzr1mWOmv3BBzB+vE3fHSh2H93Nea+dxyXvXMKxzGOB61hRlEqJP41Qc2B3gesUp8xjHWNMDnAUiC+mrbfyeCDNkVG0L299nAEYEVkoIj+JyF88PYSITBCRZBFJPnjwYCkf3YdUYCouL/fd6NE+1KcEUo6lIAgiQmxkgF3xFEWpdPhzn5Cn3YhFt8l7q+Ot3JPRLK5+cX1EAOcDvYCTwFcistIY81WhisbMBGYCJCUlBX6bfzmdEiPLDzcAACAASURBVPbvhyVLIDISRo3yuVZe6duyLz9N/IkaETUIDwvCQpSiKJUKf46EUoCWBa5bAHu91XHWaOoAqcW09VZ+CKjryCjaV3F9fGOMOWSMOQl8BvQo57P6j3KOhN5/30ZKGDoU6tb1g15FOJF1wn3eonYL4mvG+79TRVEqPf40QiuADo7XWhTW0WBekTrzgPHO+RhgsbFBxeYB4xzPtjZAB2C5N5lOm68dGTgy55bQx0Kgm4jUdIxTPyB/V2WoUE4jlDcVF4i9QWv2rSHhuQTe+fkd/3emKEqVwm9GyFl/uR37z34jMMcYs15EHhaRy5xqrwLxIrIFuBO4z2m7HpiDNQqfA5OMMbneZDqy7gXudGTFO7KL6+MI8E+sYVsN/GSMme+ft1EB8oxQQkKpm6SkwHffQUwMXHZZyfUrypz1czh08hCLty/2f2eKolQpRKMZl56kpCSTnBzALKDG2Lm0Y8esm3aDBqVqtmoVTJpkw8zlZQP3J8YY5qyfw2UdL9P9QIqinIaz3p7k6Z4GMA1ljhyxBqhWLYgv/RrL2WfD999Dtp93PRljELGecL/vGuCYQIqiVAk0bE8oU9AzrhypDyIjfatOQb7Y+gVD3h7CgRMH/NeJoihVHjVCoUw5nBKWLLHTcf6cZXUZF3cuvJMvt33JG6vfKLmBoiiKF9QIhTLlMEKTJ0OPHrBwoZ90AsIkjC+u/YL7zruPu869y38dKYpS5dE1oVCmjJ5xv/4Kq1dDnTowYID/1AJoFteMJy56wr+dKIpS5dGRUChTxpHQ7Nn2c+RIG+/U17y15i2dflMUxaeUaiQkIknABUAz4BSwDvjSGJPqR92UMhqhvA2q48b5XpVdR3cx4dMJZORk0L5+e85rdZ7vO1EUpdpRrBESkeuBO4DtwEpgExCDjbl2r4isAx4wxuzys57VD2PKFDdu3TrYsMF6cg8a5Ht1WtVpxfRLprN632o1QIqi+IySRkKxwHnGmFOebopId2xIHTVCvmb/fsjIsFYlLq7E6nlTcaNH+881+8YeN5ZcSVEUpQwUa4SMMS+UcH+1b9VR3JRxKi48HOrV832suF8P/0qDmg2oX6O+bwUriqJQgmOCiDwlIrd4KJ8sIlP9p5ZSVs+4v/8d9u2Dfv18p0KOK4ffv/97Ok7vSPLeAIYrUhSl2lCSd9xwnFw6RXgO+J3v1VHclGOPUFSUHRH5iiOnjlA3pi6xkbF0blg0M7uiKErFKWlNyBhjXB4KXU6KbMVflNIIGWNzBw0dWqqlozLRMLYhi69bzN70vdSMrOlb4YqiKJQ8EjopIh2KFjplHp0VFB9RSs+4H3+EsWNtlAR/hOoREZrXLpqVXVEUxTeUZIQeBBaIyPUicpZz/AGY79xT/EUpR0J5XnGXXlquGKceWbJjCTfNu4lDJw/5RqCiKIoXSvKOWyAiI4F7gD86xeuA0caYtf5WrtqSmwu7HK/31q29VnO58vMF+corzhjD5IWTWb1vNe3qteP+C+73jWBFURQPlBgxwRizjvz02EogSEmBnBxo2tSmR/XCd9/B3r3Wga53b990LSK8N/o9nlr6FJP7TvaNUEVRFC+U5KI9U0TO8nIvVkRuEJGr/aNaNaaUU3F5YXp+/3vfTcUBnNngTF4b8RoxEd4NoKIoii8oaST0IvCAiHQD1gIHsWF7OgC1gdeAd/yqYXWkFEYoJ8d6xYHvYsVtOLhBXbEVRQkoxY6EjDGrjTFjgSRgKTY8z0LgJmNMojHmOWNMZgD0rF6UwjNu3z5o2xY6doTExIp3uWDzArq82IXbP7u94sIURVFKSUkBTCOBp4DrgJ2AAI2A54EnReRsY8wqv2tZ3SjFSKhFC+ueffy4b6bith3ZRlR4FAl1EyouTFEUpZSUNB33NFATSDDGpAOISG1gmoi8BAwFSr+lXykdZYiWUKuWb7qc1HsSQ9sPpVWdVr4RqCiKUgpKMkLDgA7G5G+DNMYcE5FbgUPAJf5UrtpSQty4jRshOxvOOsu3Dgnt6rfznTBFUZRSUNJmVVdBA5SHMSYXOGiM+dE/alVjMjOt33V4OLRs6bHKE0/YdaDp0yvWlTGGPy34E9/s+KZighRFUcpJSUZog4hcV7RQRK4BNvpHpWrOzp02/k7LlhBx+kA1IwM+/tieDx1asa7mb57Pv5b/i5GzR5KemV4xYYqiKOWgpOm4ScCHInIDNrOqAXoBNYBRftatelKCZ9yCBZCebmPFdTgtql/ZuKjtRUzpP4WmtZoSF+3j6KeKoiiloKSwPXuAPiIyEOiC9Y5bYIz5KhDKVUtKcErIixXnizA9MRExPNhPQwAqihI8SgzbA2CMWQws9rMuChRrhE6cgE8+sedjx5a/i33H91G/Rn2iwqPKL0RRFMUHlLQmpASaYjzj5s+HkyfhnHNKnXD1NIwxXPnBlSTOSGTdgXXlVlNRFMUXlGokpASQYkZChw5BnToVm4rbf2I/e9P3knoqlWZxzcovSFEUxQeIBw9sxQtJSUkmOTnZv500bGitzZ490Ox0I5GZaePGxcaWv4vMnEx+OfQLiU18EO9HURSlBERkpTEmydM9v07HichQEdkkIltE5D4P96NFZLZzf5mIJBS4d79TvklELi5Jpoi0cWRsdmRGldSHc7+ViBwXkbt9/wbKyPHj1gBFR0OTJh6rREdXzAABREdEqwFSFCUk8JsREpFw4AVsVIXOwJUiUjRE843AEWNMe+AZYKrTtjMwDuuRNxR4UUTCS5A5FXjGGNMBOOLI9tpHAZ4BFvjmqStInnt2QgKEFf7RfPONdUwoL6v3rebR/z1KZo7Gm1UUJXTw50ioN7DFGLPNGJMFzAJGFKkzAnjDOX8fGCQi4pTPMsZkGmO2A1sceR5lOm0GOjJwZI4soQ+crLHbgPU+fO7y42U96MgRGDzY5rg7frzsYo0x3Db/Nh74+gGeWvqUDxRVFEXxDf40Qs2B3QWuU5wyj3WMMTnAUSC+mLbeyuOBNEdG0b489iEiscC9wJTiHkJEJohIsogkHzx4sIRHriBePOPmzbOx4vr0KV/AUhHhsYGPMajNIP50zp8qrqeiKIqP8Kd3nKfQmkW9ILzV8VbuyWgWV7+4PqZgp++OSzFRQI0xM4GZYB0TvFb0BV5GQt9/bz8vqUC42AFtBjCgzYDyC1AURfED/jRCKUDBCJwtgL1e6qSISARQB0gtoa2n8kNAXRGJcEY7Bet766MPMEZEngLqAi4RyTDGVDAsaAXwYoRWrLCfvXqVXaQxhuKMrKIoSjDx53TcCqCD47UWhXU0mFekzjxgvHM+BljsRO2eB4xzPNvaYNOJL/cm02nztSMDR+bc4vowxlxgjEkwxiQAzwKPB9UAgce4cRkZsHat9VM4++yyiTPGMPSdody76F6OZhz1nZ6Koig+wm9GyBmR3I5NB74RmGOMWS8iD4vIZU61V7HrM1uAO4H7nLbrgTnABuBzYJIxJtebTEfWvcCdjqx4R7bXPkIOYzyOhNassfuCOnUq+3rQ6n2r+WLrF7y66lUiwnRfsqIooYduVi0Dft2smpoK8fEQFwdHj7qz1b36Ktx0E1x/PfznP2UXu2LPCnYe3cmYzmNKrqwoiuIHitusql+PQ4WCnnEF1nBuvBGGD7eREspDr+a96NW8HItJiqIoAUCNUKhQTMy4xo3LLu7IqSPUq1GvgkopiqL4F42iHSqUkEeoLGxJ3UKzfzZj4icT0elWRVFCGTVCoYIHz7hvv4UzzoAHy5h3bsmOJWTlZpHjylH3bEVRQhqdjgsVPIyEli+HzZvhwIGyibqpx00MSBhAZHikDxVUFEXxPWqEQgUPRihvk2qSR5+S4mlXv50PlFIURfEvOh0XChhTOIK2Q543eGkjJaRlpPHz/p99qpqiKIo/USMUCuzbZ0Mj5O0Twm4b2roVYmKgS5fSiXlpxUskzkjkwa/LuIikKIoSJNQIhQIepuLyRkFnnw0RpZw0zXZlUyOiBhe2vtDHCiqKovgHXRMKBTx4xpV1Kg7gwX4PMqnXJOrXqO873RRFUfyIGqFQwMNIaMgQm0PowjIOauJrxvtQMUVRFP+iRigU8GCEkpJK7xW3ePtiAAYkDNB9QYqiVCp0TSgUqEC0BGMMf/78zwx6cxAfbvzQx4opiqL4Fx0JhQJF0nr/9BN88w0MGgTduhXfNNuVzdguYwmTMIafMdy/eiqKovgYHQkFm5wc2LXLnrduDcDcuXDnnfD22yU3jwqP4m8X/o1VE1cRHRHtR0UVRVF8jxqhYLNnD+TmQrNmdlMQ+Z5xZYmUoGtBiqJURtQIBZsi60HG5IfrKck9+54v7mH2utnkuHL8qKCiKIr/0DWhYFPECO3aBQcP2uAJBSL4nMba/WuZ9sM0akTUYGCbgTSMbeh/XRVFUXyMGqFgU8QIFQxaWtwMW4f4Drw8/GVST6WqAVIUpdKiRijYFPGMK22khJiIGCb0nOA/vRRFUQKArgkFmyIjoYgIaNSoeKcEzZaqKEpVQY1QsCkSN+7RR21Q7Usv9Vx9//H9dH6xM9OXT1djpChKpUeNUDDJzIS9eyE8HFq0cBeLQJiXn8yba97kl0O/sGjbInXLVhSl0qNrQsFk507rk92qFUREcOgQ1KgBsbHem9x17l20r9+ehLoJAVNTURTFX6gRCiZF1oMeewz+9S94/nm47TbPTcIkjFGdRgVIQUWpumRnZ5OSkkJGRkawVakyxMTE0KJFCyIjI0vdRo1QMCniGbdiBbhc0Lbt6VUzczLJys0iLjoucPopShUmJSWFuLg4EhISdGrbBxhjOHz4MCkpKbQpQzBmXRMKJgVGQjk5NnApePaMe2PNG7R6thUzV84MnH6KUoXJyMggPj5eDZCPEBHi4+PLPLJUIxRMCnjGbdwIp07ZmbkGDU6v+t2u70jLSKN2dO2AqqgoVRk1QL6lPO9Tp+OCSYGRUMFICZ54c9Sb3Jp0K72alyHft6IoSoijI6FgUsAIlSZSQt+WfYkI0+8NiqJUHfxqhERkqIhsEpEtInKfh/vRIjLbub9MRBIK3LvfKd8kIheXJFNE2jgyNjsyo4rrQ0QGi8hKEVnrfA7035vwwPHjcOiQTd/QpInXyNnHMo+Reio1oKopiqIECr8ZIREJB14ALgE6A1eKSOci1W4Ejhhj2gPPAFOdtp2BcUAXYCjwooiElyBzKvCMMaYDcMSR7bUP4BBwqTHmLGA88JYvn79E8kZBrVuDCB9+CB9+ePp03Ks/vUrTp5sy9bupp8tQFKXSEx4eTvfu3enatSuXXnopaWlpANxwww00atSIrl27lllmcW0///xzOnbsSPv27XnyySdLLN+xYwfDhg2jY8eOnHHGGTzxxBPleErv+HMk1BvYYozZZozJAmYBI4rUGQG84Zy/DwwSu7I1AphljMk0xmwHtjjyPMp02gx0ZODIHFlcH8aYVcaYvU75eiBGRAKXmrTIHqGWLWHUKKhVq3C1HWk7yHHl0L5++4CppihK4KhRowarV69m3bp11K9fnxdeeAGA66+/ns8//7xcMr21zc3NZdKkSSxYsIANGzbw3nvvsWHDBq/lLpeL0aNHc8stt7Bp0ybWrl1LcnIyM2f6zkvXn0aoObC7wHWKU+axjjEmBzgKxBfT1lt5PJDmyCjal7c+CjIaWGWMySz6ECIyQUSSRST54MGDJTxyGSgSM84bz13yHLsn72b4GcN917eiKIUR8c9RRvr27cuePXsAuPDCC6lfv365Hsdb2+XLl9O+fXvatm1LVFQU48aNY+7cuV7LFy5cSEJCApdddhkA0dHRTJ8+nWnTppVLL0/4c5Xb00+gaMRNb3W8lXsymsXVL1EPEemCnaIb4qEexpiZwEyApKQk30UMLTASmjLF2qTJk6Fbt9OrNotr5rNuFUUJTXJzc/nqq6+48cYbS6x7wQUXkJ6eflr5tGnTuOiii7y227NnDy1btnRft2jRgmXLlnktj46OJjExsZCMpk2bcuzYMbKysoiKiirNoxWLP41QCtCywHULYK+XOikiEgHUAVJLaOup/BBQV0QinNFOwfre+kBEWgAfAdcZY7aW/1HLQQEj9MHDsHYtTCiQHijXlcvOoztpW89D+ARFUXxLECPSnzp1iu7du7Njxw569uzJ4MGDS2zz7bfflqsvT5H3RcRreXh4OMePHz9NxsmTJ4mI8I358Od03Aqgg+O1FoV1NJhXpM48rFMAwBhgsbFvYx4wzvFsawN0AJZ7k+m0+dqRgSNzbnF9iEhdYD5wvzFmqU+fvDQ4RuhEk3asX28DaXfvnn978fbFtPtXO6758JqAq6YoSuDIWxPauXMnWVlZ7jWh4rjgggvo3r37aceXX35ZbLsWLVqwe3f+ikZKSgrNmjXzWt6/f38+++yzQkZq0aJF9OjRgzBvof7LijHGbwcwDPgV2Ar8n1P2MHCZcx4D/BfreLAcaFug7f857TYBlxQn0ylv68jY4siMLq4P4G/ACWB1gaNRcc/Ts2dP4xNcLmPi4owB8+2naQaMSUwsXGXGihmm5mM1zd+//rtv+lQUpRAbNmwItgrGGGNiY2Pd5z/99JNp2bKlycrKMsYYs337dtOlS5dyyfXUNjs727Rp08Zs27bNZGZmmm7dupl169Z5LTfGmGuuucb87W9/M8YYs2/fPtOlSxezdOlSr/16eq9AsvFmJ7zd0MOPRujQIfvq4+LMM/90GTDmpptOr5aemW6OnDrimz4VRSlEKBohY4wZPny4efPNN824ceNMkyZNTEREhGnevLl55ZVXSi2zuLbz5883HTp0MG3btjWPPvposeVPPPGEOfPMMw1gvvrqKzNq1ChTq1Ytk5iYaLZv3+6x77IaIbH3ldKQlJRkkvNCG1SElSvthqBu3bi66xrefRdmzICJEysuWlGU0rFx40Y6deoUbDWqHJ7eq4isNMZ4DEqmYXuCgYeYcQUjJWw4uMHjQqGiKEpVQ41QMHCMkElowxVXwNChkLexOeVYCl1f7Eq3Gd3IdeUGUUlFURT/o9Ewg4FjhKRtGx67o/CtzYc307hWY85scCbhYeFBUE5RFCVwqBEKBkUyqhZkQJsB7J68m7SMtMDqpCiKEgR0Oi4YOEbo8z1nsXw5ZGUVvh0RFkGDmh4y2ymKolQx1AgFGmNg504AJj7Rmj59YPNme2v9gfVk5WYV01hRFKVqoUYo0OzbBxkZHKh/Jrt2hxEbC2eeCTmuHAa+OZBmTzcj5VhKsLVUFEUJCLomFGicqbjk+IshFXr0sCF7dqSl0Ci2Edm52TSPKxpsXFEUpWqiI6FAk2eEos4F8vcHJdRN4OdbfmbpDUuRcoSAVxSlcuIpqd3u3bsZMGAAnTp1okuXLjz33HOllldS21BLahf0UDiV6fBJ2J5HHjEGzPB26w0Y8+67FRepKErZCcWwPdddd5159NFHzd69e83KlSuNMcYcO3bMdOjQwaxfv75U8oprm5OTY9q2bWu2bt3qjhG3fv16r+W5ubmmR48eZu7cucYYYzIyMszll19uXn75Za/9lzVsj46EAs327RhgxcEEwI6ENh3apC7ZihJkistNVzCR6MyZPs1jV4i8pHZNmzalR48eAMTFxdGpUyd3sruSKK5tKCa1UyMUaHbsIJX6ZLkiqVsX2rWDiZ9OpMm0JizevjjY2imKEiTyktrl/cPPY8eOHaxatYo+ffq4y0qbyqFoW0/J6/bs2eO1fOPGjcUmtfMF6pgQaLZvJ55UDidvZ1/tM8h2ZREZHklEWARJzTzG91MUJQCUNlzjhAmFE1BWlOKS2h0/fpzRo0fz7LPPUrt2bXd5aZLaeWprPDxkVU5qpxQlJwd27QIRJKE1TZtCVHgUi65dRMqdKdSOrl2iCEVRqhbektplZ2czevRorr76ai6//PJCbUoaCXlrW+2S2lW1o8KOCdu3GwMmp1nLislRFKXChKJjQsGkdtdee63505/+VGZ5LpfLa9tQTGqnI6FA4jglND+4ip49YX3KLnak7Qi2VoqihAhnn302iYmJPPnkk7z11lssXrzYPcr57LPPSiVj6dKlXttGREQwffp0Lr74Yjp16sTYsWPp0qWL1/Inn3yS5ORkHn30URYvXsytt97Kzp07ue2229ixY4dPnlmT2pWBCie1e+01tt74GO3ZSuPGMPq123lxxQv8a+i/+GOfP/pOUUVRSkST2vkHTWoXyuzYwQrs7tSkJAhDqBFRgwtaXxBkxRRFUYKDGqFAsn07ydgvA716wfPDnmf/3ftJbJxYQkNFUZSqiRqhQLJ9u3sklBeuJy46TsP0KIpSbVEjFEByt+3kJ3pATBrhrVZ49M1XFEWpTqgRChQZGWz6LY7jxFHvgvcY+kFvbv7k5mBrpSiKElTUCAWKnTtpxl7ea3gHw4bn0ii2Ef0T+gdbK0VRlKCiYXsCxY4d1OUo485az7g7/kV27kQMOh2nKEr1Ro1QoHDyCNGmDQCR4ZFBVEZRFCU00Om4AJG9ZSe38AJ3ZrUiKyc72OooihIieEpql5GRQe/evUlMTKRLly489NBDpZZXUltNaleJj4rEjvvponsMrb41/B3T4+Ue5ZajKIpvCMXYcXlJ7Vwul0lPTzfGGJOVlWV69+5tfvjhh1LJK66tJrWrxqzYXBei0ql9sgOD2w4uuYGiKAFFpggypfCevUvfuxSZInyy6RN32cyVM5EpwoRP8vM57E3fi0wRmj3drEI65CW1ExFq1aoF2IjY2dnZpd5PWFzbUExqp2tCAWLFvpaQeQkP1TyXSf1jgq2OoighRl5SuxtvvNF93bNnT7Zs2cKkSZNOS2qXnp5+moxp06Zx0UUXeW3rKXndsmXLvJZHR0cXm9QuKiqqws+tRigQpKeTnNkVgN4DaxMdoRESFCXUMA+d7q36yZWfnFY2oecEJvQsnNWuWVwzj+1Lg7ekduHh4axevZq0tDRGjRrFunXr6NrV/h8pKamdt7bGwwb5Kp3UTkSGisgmEdkiIvd5uB8tIrOd+8tEJKHAvfud8k0icnFJMkWkjSNjsyMzqrx9+JpTv+zk53b7kOgjnN1DDZCiKPl4S2qXR926denfvz+ff/65u6y06b2Ltq1WSe2AcGAr0BaIAtYAnYvUuQ2Y4ZyPA2Y7552d+tFAG0dOeHEygTnAOOd8BnBrefoo7pnK65gw/8l5hgciTNj9tUx6Znq5ZCiK4ltC0TEhL6ndgQMHzJEjR4wxxpw8edKcf/755pNPPimVvOLahmJSO39Ox/UGthhjtgGIyCxgBLChQJ0RwN+d8/eB6WJX0EYAs4wxmcB2EdniyMOTTBHZCAwErnLqvOHIfakcffzgqxeQx+F9W6iT1pV6UTWpFVXL1+IVRaki5CW1e/HFF/n444/Jzc3F5XIxduxYhg8fXioZv/32G+PHj/fYtmDyutzcXG644Qa6dOkC4LE8L6nd22+/zYABA5g+fbo7qd3HH39MQkJChZ/Zn0aoObC7wHUK0MdbHWNMjogcBeKd8h+LtG3unHuSGQ+kGWNyPNQvTx9uRGQCMAGgVatWxT6wN67t2ZBrZ8SSfe3V5WqvKErVpeiayyef2HWosuwNKki3bt1YtWqV1/vDhg1j2LBhpSq/7777uO++/JWUgQMHlkun4vCnEfK0+FF09ctbHW/lniYhi6tfnj4KFxgzE5gJNrOqhzYlc801cM01aIwERVGUwvjTMSEFaFngugWw11sdEYkA6gCpxbT1Vn4IqOvIKNpXWftQFEVRAoQ/jdAKoIPjtRaFdQqYV6TOPGC8cz4GWOwsYs0DxjmebW2ADsBybzKdNl87MnBkzi1nH4qiVBOMB9dkpfyU5336bTrOWX+5HViI9Wp7zRizXkQexnpKzANeBd5ynAJSsUYFp94crBNDDjDJGJML4Emm0+W9wCwReRRY5cimPH0oilL1iYmJ4fDhw8THx2t2Yx9gjOHw4cPExJRtM77oN4HSk5SUZJKTk4OthqIoPiA7O5uUlBQyMjKCrUqVISYmhhYtWhAZWXgFXERWGmOSPLXRiAmKolRLIiMjaeOkVlGChwYwVRRFUYKGGiFFURQlaKgRUhRFUYKGOiaUARE5COysgIgG2D1NoYrqVzFUv4qh+lWMUNavtTGmoacbaoQCiIgke/MQCQVUv4qh+lUM1a9ihLp+3tDpOEVRFCVoqBFSFEVRgoYaocAyM9gKlIDqVzFUv4qh+lWMUNfPI7ompCiKogQNHQkpiqIoQUONkKIoihI01AiFEKKhfH2GvsuKoe+vYuj7Kz1qhEKLOiLi/pmE2i+yiPQRkYuCrYc3RKSBiNQBMCG42Knvr2KISD8R+X2w9fCGiDQTkZag768sqBEKEUTkMmARMF1EnofQ+kUWkYuBGcCRIuUhYShFZCQ2seG/ReRtEakdbJ0Kou+vYjjv7xlgd7B18YSIXA4sAV4VkfdF5GwRiQ6yWm5C+v0ZY/QI8gGcAawFBgEJwDfAHCDauS9B1m8AcAA427muEex3VkS/lsAPQB8gynl3bwKdgq2bo99AfX8V0q8/cAw4y7muBdQMtl4F9IsHvgB6OtdPYb9wXAJEhYB+If3+dCQUGqQDvwJrjDE7jDH9sLme3obgjohEJBxoiv0GlSMiscDrIvKWiLwrIhFOvWB+o0/H/pGdMsZkGWPGAoeB+0WkZgjo1xx9fxUhFzuCbOyMLmYBb4vIRyLSIYh65ZGJNd4NAYwxfwE2AyOwXzCD/f4MIfz+1AiFBieBo0DPvAJjzOXYX5ppQdPK6pELfAA8CfwL2AYsA6YAtYH/OvUCbihFpLZjJI8DPwI9CqxpTAZigZeCqF9PEWljjHkLeJwQe395GGPSgOVAz1B6fwX0+xa4Bnge2A8sAG4EfgP+GSy9RKSxiNQ0xhzHfmHsISJtAYwxTwOngL8718H4/TtfHkaJrAAADElJREFURLoZY74Brsb+/oXM+3MT7KFYdT2AC4G7gcFAPeBi4Cegd4E6nYF/Bkm/c4E/Ono1dsquA24vUKcmMBeoFQT9hgNv5fUN/M7RZSRQzymLc+rEBEG/IcBe4EHnui5wbQi9v8FYw/gHoC2QCMwHRoXI++sHPIT9Z9nRKRsA/LlIvc+BJkHQ73fA+0AD57oPdgruNqBdgXrzgYZB+vkeAt7AmXpz3umfQuH9FTw0vXcQEJEhwIvYEcalwL3ATcBjwEsiMgVYgTUESSJSwxhzKoD6DceOfL4Gujs6PIn91p5VoOoVWAMaUERkKPAwcKex30IxxswXkRrYf1pNRCQZ6Ig15AH9PXd+vo9jf749RUSMMWki8gF26iaPYL2/fsALwL+xa5D/B1yO1fl+7Ag8mO9vuKPLx0AnoK+I3GuM+VpElhSodxXWuAfsb8Pp9xLs3+rtxphDAMaYZSISif2i1kREVgI1gNYU/psJhH4XY9/fM9j3F4udbfmfc+TVC8r7O41gWsDqegD3YP+BAtQBJgBrsP8QBmOnQD4GkoFuAdatI3ZEluRcXwwsBGoXqCPYP7a1QJcA69fWeVd5768+9tv7hdhv7mcDfwM+w3ordQ+wfn2w6wG9neuFwNNF6gh2VBTw9+f0fxPwVIHrcdg1oG5AV+CBIL6/1sBi4Bznugv2y0+LAnXCnPe3Pgi/f02xjkNPO9f1gOuBsVgHj6bO9QJgHo4zSgD16+b83+jrXL8PvFekTtDen6dDR0LB4QTQA8AYcxSY6axbvghcCSzFfouKMMbsD7Buu7Bzx+sc/RaKyB+x0zXfOnXisf8sxhpjNgZYv71YV/Y4EbkC+AuwCfuHVRP77fRRxwEgyhhzxLsov7AOGGOMWeNcPwWMFZGGxpiDTlkToBUBfn/OiMxg1wX65JUbY2Y5v3+zgUuMMY8E8f39hv0StsrRbb2zqH8+dkEd7N9GDex7DvTv3yHse2osIndh/15/wjpP/BW41hjzuoi8h/37PRFg/X4FrjbGbHKu7wOmiEh3Y8xqp6wOEE1w3t/pBNsKVpcD6yHVyjmvif02/3iB+/HYOeW+QdSvTZGyPBfxucAQ5zwJu6AeHiz9sJ5Ij2BHErc4ZY2xU0yjQuX9OeWNsF8q/likPGiuu9jpmR+BZ4qUPw6MDpZeBfTIC6wc7ny+DFzlnPcnyGsYzu/fLdiprUkFyh8teB0Kh/O3Ohv4vyLlAf37Le5Q77gAICJjgA+B90XkCewoaChwrnONMeYwdu797CDqN1tEHnE23mGMyVu/SAH2ORtqp2IXOnODpN9jwAXAg8BtxpgZjq77gXDsP/2AUkS/h0VkVN49Y8wB7JrLFSLSsUB5wNYJRGSYiLzprFlg7LfzocBwEXmmQNUY7BpCQCmqH6d77e4G9hdYK4oKpn7Oz+5V4G5jzAsFqkZhp4QDSlH9HI9RAIwxx7Dru1eLyLkFygP291siwbaCVf3ADn1/wBqe5ljvmZeBq7Dzx2uwv9DTsNNKHUJAv5eAGwvUeQTrwrscZ8NbkPWbiZ32KFhvLLAaaB8C+hV9f/Wwc/PDg/D71wf7T/x/wCdAZIF7dbFTrK/w/+2dfYwdVRmHn19bAbVf1hpDIyBWGiNVWkywlHSbGOIfaraGxmCVNKZqqAiIxCCIHzVEY1JjQSMhEQFJmxAxabVEEJFiCkmDlkqLRW1NtQoqiWTdNfJV9/WP92z2ara4ez9m5m5/TzLZnbkzZ59u7953zpz3vAduAw5S8QTV4/nRcqdOZpH+gUzWWdoEvwnOuxjYCyxpgh/jvcgZ5eu3gEvH9pu0uSfUe2YCx4ChiHiK7Bo/QA6knwacTw7EHiEfJR1qgN+D5JyHwXLOPOBtwLqIONAAv58CK8qdMWVs6Ivks/DDDfAb+/29DyByXOVeciC4amaRNxGrybGg7S139ENk4smdZCbk+6P6MYIJ/SLi3y09oxfJR9jrIuKJJvjBeI+jZGt+ClgfEb9rgl/5/c2MiNFy3s+AH7fsNwYvalcBkjaRM6eviohnJC0ks1NeGxGfr1WOl/WbFxGblHXETo3xwc6m+M2NiC+XWf0LI+JoE/3qcGpF0pyIGCnJBjcCi8gbnhclLYiIZxvsNwd4iXz/HWmg32vIihMLI+KvDfRbFBFP1+E1WRyEekC5A15J9iC+RKYRryOfF2+OiL9JOgPYRmaoVPrmnaLfByPizw32u7j0QJrqtzYqznBs8ZtLztgfjjIGVT7Ut5D1w35OpuR/LiL+1VC/s8l0/CrH0Kbi91ZybOiFiVur3W8J+f9b71ygl8GP47qMpHeQWW57yDTSm8g36h7yjunbZYB6Rbnk+Yb7/bPhfpWmwLbhV9mH0wR+ryLHAt6j8XI8IxHxMfLD6SvA7RUHoKn6fafiADRVv1srDkBT9bujyQEIHIR6wRLg/oj4UURsAB4mS7icBNwB7Ce7zBuBK8tzefvZr1d+D5ElZgY0Xix1kLyLXhXj85nsZ796qDszYrpt5Iz++4CVLccuIyfazS37s6mhHpf9Tli/T5DLM8wr+6uoOAvTfieGXzubx4S6gKRl5GMXRcTBMpdlGPhhRPymnLMV+FNEXGc/+9Xk98eIuN5+9msSfhzXIcpihjuBTwLfl7SWnPfzJmCNpIFy6qPkmIH97FeXX9UlZOw3zf26Qt1dsX7dyCKUs8lCj4Pl2PnA78mJa2eQmSu7yEc1R6lwoqf97Gc/+/XDVrtAv2/kkgKXMD5T+TxydvdFZf8N5HINp9vPfvaz33Ty68q/sW6Bft/IQenb+e+lDlaR5dQX1+VlP/vZz379sHlMqE2krH0fETeT+fq3SJpXSmbsJlN1aysSaD/72c9+/YCz46ZAmYS4gLwLGY2WSrSS7iJXKNxD1nO6GlgdFVYbsJ/97Ge/fsNBaJIolzf4KvBU2X5JzkYebjlnA1m36RxgU0RUVrDSfvazn/36EQehSaCsmrsV+GZEPFLSJFeQufubI1dHbT3/5Ki2lIf97Gc/+/UlHhOaPHOBs8r324F7yFIt6wAknSfp3PJ6ZbWuWrBfZ9ivM+zXGU336xkOQpMgIl4CvgFcJGlV5JocD5OLqA1IeiVwAfB0Ob/S7qX97Gc/+/Ut7aTUnYgbufTx5eSqngMtxx+iAamS9rOf/ezXj9us9sPXiUVEPC9pGxDAdZLeQj6zfR0VL3cwEfbrDPt1hv06o+l+vcSJCVNE0klk1/hSci2gmyJiX71W49ivM+zXGfbrjKb79QIHoTZRri8f0cA128F+nWK/zrBfZzTdr5s4CBljjKkNZ8cZY4ypDQchY4wxteEgZIwxpjYchIwxxtSGg5AxxpjacBAypiIkzZd0Wfl+kaQfdLHtqySt72J7d0k66/+faUxnOEXbmIqQ9EbgnohY2uV2ZwGPAedGxLEutbkauCQiPt6N9ow5Hu4JGVMdXwMWS/qVpLslPQEg6SOSdkjaKemIpMslXS1pn6Q9khaU8xZLuk/SXkm7S2kXgHcBj40FIElXSjooab9yMTQkvVrSbZJ+UdpdU47PlPR1SQfK+VeUNncDF5YAZ0zP8BvMmOq4FlgaEcvGekUtry0FlpOFLA8Dn42I5ZK2AOuBG8nilhsj4pCkdwI3kwHoAmDv//ycMyPiBUnzy7HrgQcjYkM59qikB0rbZwLLI+LYWMCLiFFJh8kF1FrbNqarOAgZ0wx2RcQIMCLpH8DOcvwA8HZJs4GVwN2Sxq45uXw9FXiypa39wDZJO4Ad5di7gUFJnyn7pwCnAxcCt4z1oiLi2ZZ2niFX8nQQMj3DQciYZtC6UuZoy/4o+Xc6AxiKiGUTXPscGVTGeC8wAAwCX5B0NiBgbUT8tvVCZUQ73sDwKaVtY3qGx4SMqY4RYE47F0bEMHBE0gcgg4ekc8rLTwJvLsdnAKdFxC7gGmA+MBv4CXBFCTpIWl6uvR/YODb2M/Y4rrAE+HU7vsZMFgchYyoiIv4OPFISEja30cSHgY9KepwMDmvK8XvJng/ATGCrpAPAPmBLRAwBNwCvAPaXn39DOf9W4Gg5/jjwIQBJrweei4i/tOFpzKRxirYx0wBJ24FrIuJQl9r7NDAcEd/tRnvGHA/3hIyZHlxLJih0iyHge11sz5gJcU/IGGNMbbgnZIwxpjYchIwxxtSGg5AxxpjacBAyxhhTGw5CxhhjauM/FU7GK2UNdhYAAAAASUVORK5CYII=\n", "text/plain": [ "<Figure size 432x288 with 1 Axes>" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "# Charging\n", "# initial condition\n", "y0 = 0\n", "y1 = odeint(model,y0,t,args=(C,R1,V))\n", "y2 = odeint(model,y0,t,args=(C,R2,V))\n", "y3 = odeint(model,y0,t,args=(C,R3,V,))\n", "\n", "# plot results\n", "plt.plot(t,y1,'r-',linewidth=2,label='R1=100$\\Omega$')\n", "plt.plot(t,y2,'b--',linewidth=2,label='R2=200$\\Omega$')\n", "plt.plot(t,y3,'g:',linewidth=2,label='R3=300$\\Omega$')\n", "plt.xlabel('time(sec)')\n", "plt.xticks(rotation=45)\n", "plt.ylabel('Q(C)')\n", "plt.legend()\n", "plt.show()\n" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAaEAAAEgCAYAAAAdRMy/AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAgAElEQVR4nOzdd3hUVfrA8e+bTg9Eeg3SUUAIICpWEHApKijYXdvatrn7s6xtVdwFVxdxrVhW1wIqugKLgAULVgiI0kRpQui9hUDK+/vj3EwKM6nTQt7P88yTmTPnnvPOhcybe++554iqYowxxkRCTKQDMMYYU31ZEjLGGBMxloSMMcZEjCUhY4wxEWNJyBhjTMTERTqAquS4447TNm3aRDoMY4ypUhYuXLhDVRv6e8+SUDm0adOG9PT0SIdhjDFVioj8Eug9Ox1njDEmYiwJGWOMiRhLQsYYYyLGrgkZY6ql7OxsMjIyyMrKinQox4ykpCRatGhBfHx8mbexJGSMqZYyMjKoU6cObdq0QUQiHU6Vp6rs3LmTjIwMUlNTy7xdSE/HichgEVkpIqtE5E4/7yeKyJve+9+KSJtC793lla8UkUGltSkiqV4bP3ttJnjlp4vIIhHJEZFRxfq/yqv/s4hcFYp9YIyJTllZWaSkpFgCChIRISUlpdxHliFLQiISCzwFDAG6AJeISJdi1a4FdqtqO2ACMN7btgswBugKDAaeFpHYUtocD0xQ1fbAbq9tgPXA1cAbxeJrANwP9AX6APeLSP3gfPoAbMZyY6KKJaDgqsj+DOWRUB9glaquUdUjwBRgRLE6I4BXvOdTgXPEfYoRwBRVPayqa4FVXnt+2/S2OdtrA6/N8wFUdZ2q/gDkFet7EPChqu5S1d3Ah7iEF3xPPEFOo4ase+j+kDRvjDFVVSiTUHNgQ6HXGV6Z3zqqmgPsBVJK2DZQeQqwx2sjUF8ViQ8RuUFE0kUkffv27aU06d+Tq/NIvDyePlv+V6HtjTHmWBXKJOTvuKz4+ahAdYJVXpIybaOqk1Q1TVXTGjb0O+tEqXq060pezT3sSdjJgSMHKtSGMcYci0KZhDKAloVetwA2BaojInFAPWBXCdsGKt8BJHttBOqrIvEFRd8zWxL/7DdkP7GWnMzaoejCGFNFxcbG0qNHD0444QSGDRvGnj17ALjmmmto1KgRJ5xwQrnbLGnb2bNn07FjR9q1a8e4ceNKLV+3bh3nnXceHTt2pEOHDvz973+vwKcMLJRJaAHQ3hu1loAbaDC9WJ3pQP6otFHAXHXrjU8Hxnij51KB9sD8QG1623zitYHX5rRS4psDnCsi9b0BCed6ZUEX3641PXdmgsawcH5uKLowxlRRNWrUYPHixSxdupQGDRrw1FNPAXD11Vcze/bsCrUZaNvc3FxuueUWZs2axfLly5k8eTLLly8PWJ6Xl8fIkSO58cYbWblyJUuWLCE9PZ1JkyZV6jMXFrIk5F2fuRX3xb4CeEtVl4nIgyIy3Kv2IpAiIquA24A7vW2XAW8By4HZwC2qmhuoTa+tO4DbvLZSvLYRkd4ikgFcBDwnIsu8PnYBD+ES2wLgQa8s+GrUoHdNF+aXc7ey/WDFri0ZY45t/fr1Y+PGjQCcfvrpNGjQoELtBNp2/vz5tGvXjrZt25KQkMCYMWOYNm1awPI5c+bQpk0bhg93X9mJiYk8+eSTPProoxX/kMWE9GZVVX0feL9Y2X2FnmfhkoO/bR8GHi5Lm175GtzoueLlC3Cn2vz18RLwUokfIkh6N9sIiW/yUPx1bP30Sp761VPh6NYYUxahGqpdjtsycnNz+fjjj7n22mtLrdu/f3/2799/VPmjjz7KgAEDAm63ceNGWrYsuArRokULvv3224DliYmJdO/evUgbTZs2Zd++fRw5coSEhISyfLQS2YwJYZLW8QCkdyIn7gBr96yNdDjGmChx6NAhevTowbp16+jVqxcDBw4sdZt58+ZVqC/1kxRFJGB5bGwsBw4UHUylqmRmZhIXF5z0YROYhknHbok8tvUV3tgwjpmXHnUgZ4yJJNXQPMog/5rQL7/8wpEjR3zXhErSv39/evTocdTjo48+KnG7Fi1asGFDwZ0pGRkZNGvWLGD5mWeeyfvvv18kSX344Yf07NmTmJggpQ9VtUcZH7169dIKe+EF99/yyisr3oYxJmiWL18e6RBUVbVWrVq+54sWLdKWLVvqkSNHVFV17dq12rVr1wq162/b7OxsTU1N1TVr1ujhw4e1W7duunTp0oDlqqqXX3653nPPPaqqumXLFu3atat++eWXAfv1t1+BdA3wvWpHQuGSvyz4WncqbtvBbew/fPQ5XWNM9XXSSSfRvXt3pkyZwiWXXEK/fv1YuXIlLVq04MUXXyxzO4G2jYuL48knn2TQoEF07tyZiy++mK5duwYsHzduHOnp6YwdO5a5c+dy00038csvv3DzzTezbt26oHxmUZvPrMzS0tK0wst7r1nD5uNP5fm6f+aT2/YzT8byryH/4qbeNwU3SGNMmaxYsYLOnTtHOoxjjr/9KiILVTXNX307EgqXli05JLW4f9+fWPRRO0SE9XvXRzoqY4yJKBsdFy7x8aS2zKHB+p3s+mYk3704iB4dKjYNkDHGHCvsSCiMpG0qaaRDTg3WLLUEZIwxloTCKdVLQsCCBa5o5Y6V5GnxVSaMMaZ6sCQUTm3a0BuXfdLT4YI3L6DTU534cv2XEQ7MGGMiw5JQOBU6EkpPh44pnaibWNdmUDDGVFs2MCGcUlNpzkb61f6BVoO6cXO3O7j/jPuoEV8j0pEZY0xEWBIKp9RUBPiq5kCYshVIjnRExhgTUXY6LpyaNoXERNi2DQ4e9BXnaR5Lti6JYGDGmEjxt6jdhg0bOOuss+jcuTNdu3Zl4sSJZW6vtG2jbVG7iM/HVpUelZo7Ll+HDqqg2z5foV99pXo457C2e6Kdxj4Qq1sPbK18+8aYMonGueOuvPJKHTt2rG7atEkXLlyoqqr79u3T9u3b67Jly8rUXknb5uTkaNu2bXX16tW+OeKWLVsWsDw3N1d79uyp06ZNU1XVrKwsvfDCC/W5554L2L/NHRft2rRhJw1odHonBgyAWBLodFwnmtVpxqpdqyIdnTEmgvIXtWvatCk9e/YEoE6dOnTu3Nm32F1pSto2Ghe1syQUbqmppLCL1g32kZkJK1bAyyNeZt0f1nFKy1MiHZ0x1ZZI4Efh1awnTSq5bkXlL2qX/4Wfb926dXz33Xf07dvXV1bWpRyKb+tv8bqNGzcGLF+xYkWJi9oFgw1MCLfUVAB6p6zll13dSU+Hq09IiXBQxphIKWlRuwMHDjBy5Egef/xx6tat6ysvy6J2/rZVPxNW26J21Y2XhNISfgAKZk4AyMrJYum2pZGIyphqT0tYm+6GGwrq3XBDpdexKyLQonbZ2dmMHDmSyy67jAsvvLDINqUdCQXa1ha1q+KPoAxMmD9fFfTj469XUO3TxxWv2rlK6/29nrb4ZwvNzcutfD/GmBJF48CEwovaXXHFFfr73/++3O3l5eUF3NYWtTO+I6Ge2+cAsHgxHDkCbeu3pUGNBjSp3YQtB7ZEMkJjTITkL2o3btw4Xn31VebOnes7ynn//ffL1MaXX34ZcFtb1K6Kq9SidvlUoU4dOHiQju1y+GlVLN9/D926we5Du6lfo35wgjXGlMgWtQsNW9Qu2on4joamPvwTO3e6BARYAjLGVDuWhCLBS0Inxq2gQYOj396RuYOVO1aGOShjjAk/S0KR4CUh1h49e/b7P79P08ea8vvZvw9zUMYYE36WhCLBS0K6Zi1jxsDxx0Nmpnvr5BYnkxCbQEJsArl5uREM0hhjQs9uVo0ELwnJurUs3wBr1rhRcqecAg1qNGDLn7ZQJ7FOhIM0xpjQsyOhSGjTxv1cu5bevd3TwoPuLAEZY6oLS0KRkH9NaN060nq5IfKFZ07I99POn9iwd8PRbxhjzDHCklAk1K0LDRrAoUP0Pn4XUPRICGD8F+Pp+GRHHv/m8QgEaIwx4WFJKFLyh2nXWEVCAqxcCfv2Fbx9RpszqJ1Qm9iY2AgFaIwJB3+L2mVlZdGnTx+6d+9O165duf/++8vcXmnbVqtF7YDBwEpgFXCnn/cTgTe9978F2hR67y6vfCUwqLQ2gVSvjZ+9NhNK6gOIB14BlgArgLtK+zxBmTsu36hRqqD6xhualuaefvJJwdt5eXl68MjB4PVnjCkiGueOy1/ULi8vT/fv36+qqkeOHNE+ffro119/Xab2Stq2Wi1qJyKxwFPAEKALcImIdClW7Vpgt6q2AyYA471tuwBjgK64pPO0iMSW0uZ4YIKqtgd2e20H7AO4CEhU1ROBXsBvRKRN8PZAKQrdK3T99fDww9C6dcHbIkLN+JphC8cYE3n5i9qJCLVr1wbcjNjZ2dlIGRcrKmnb6raoXR9glaquUdUjwBRgRLE6I3BHIwBTgXPE7a0RwBRVPayqa3FHMX0Cteltc7bXBl6b55fShwK1RCQOqAEcAQqdEAuxQiPkbrgB/vKXgrxUmKryxfov2H94f9hCM6Y6kgcEeaDoF/2wycOQB4QZK2f4yiYtnIQ8INwwo2B9h037NyEPCM0ea1bh/osvapebm0uPHj1o1KgRAwcOLNeidoG2jcZF7UKZhJoDhYd2ZXhlfuuoag6wF0gpYdtA5SnAHq+N4n0F6mMqcBDYDKwHHlXVXcU/hIjcICLpIpK+ffv2sn720hUaIVeSS9+9lP7/7s/U5VNLrGeMqZryF7VLSUlh165dvkXtYmNjWbx4MRkZGcyfP5+lSwvWGps3bx6LFy8+6jFgwIASt1U/E1Yfy4va+Tt2LP5JA9UJVnlJffQBcoFmuOtJfxKRtkdVVJ2kqmmqmtawYUM/TVVQsal7Pv8cxo+HXcXS4IDUATSv05w8zQte38aYo+j9it5f9CtqxiUz0PuVYR2H+cpu6HUDer8yaVjBmt/N6jRD71c2/WlTufsNtKhdvuTkZM4880xmz57tKyvr8t7Ft61Wi9oB/YA5hV7fRbGL/8AcoJ/3PA7YgUsaRerm1wvUprfNDiCueN8l9PEUcEWhtl4CLi7pMwV1YMKhQ240Qny8ak6OnnaaezlrVtFqh3MOa05uTvD6NcaoanQOTMhf1G7btm26e/duVVXNzMzU0047TWfMmFGm9kratrotarcAaC8iqSKSgBtoML1YnenAVd7zUcBcL+DpwBgRSRSRVKA9MD9Qm942n3ht4LU5rZQ+1gNni1MLOBn4MYifv2RJSdC0KWRnw8aNpHkrbRS/XyghNsGGaRtTTeQvavf0009z1lln0a1bN3r37s3AgQMZOnRomdrYvHlzwG2jcVG7UA/RPg/4CVgN3O2VPQgM954nAW/jBh7MB9oW2vZub7uVwJCS2vTK23ptrPLaTCypD6C2V74MWA78X2mfJ6hHQqqqp5ziDn8+/VRff909HT7cf9XDOYf1lcWvaOaRzODGYEw1FS1HQsea8h4JhXQCU1V9H3i/WNl9hZ5n4YZK+9v2YeDhsrTpla/BXecpXu63D1U9EKjvsGnTBr76CtauJe2UMwD/0/cAjHprFDN+msHBIwe5qfdN4YvRGGNCyGZMiKRCI+TatYN69WDzZtjk59rmFd2uoEvDLjSrU/EhoMYYE20sCUVSoRFyMTHQq5d76e9oaGSXkSy5aQkjOhW/1coYU1HqZ2iyqbiK7E9LQpFUbJh2797QogUUG5YPQIzEECP2z2VMsCQlJbFz505LREGiquzcuZOkpKRybSf2D1B2aWlpml58+FplrF0Lbdu6zLNhAzk5UNr9X/sO7+PZ9Gfp1bQX57Q9J3ixGFPNZGdnk5GRQVZWVqRDOWYkJSXRokUL4uPji5SLyEJVTfO3ja2sGkktWkBMDGzcCIcPE5eYWOomLyx6gTs+uoP+rfpbEjKmEuLj40n1N1eWCSs7vxNJ8fHQsiWowvr1vuI9eyAz0/8m1/e8nqEdhnJ3/7vtNIIxpsqzJBRpxeaQu/56qF8f/vc//9XrJNZhxiUzGNRuUJln1TXGmGhlSSjSig1OaOaNwA50v5AxxhxLLAlFmp8RcnD09D3FbT+4nf/74P+46X9246oxpuqyJBRpxZJQ/hxyCxdCXgkTZx/MPsiEbybw/KLn2XJgS4iDNMaY0LAkFGmFFrcDaNLEDZrbvx9Wrixhs+Q2TBw8kW+v+5YmtZuEPk5jjAkBS0KRVuxICMp+Su6WPrfQq1mvEAVmjDGhZ0ko0po2hcRE2L4dDh4ECpJQeQYnHDxyMATBGWNMaFkSirSYGGjd2j33hmlfdBHMmAH33FP65odzDnPVe1eROjGVfYf3hS5OY4wJAUtC0aDYKbl27WDoUGjUqPRNE+MSWbt7LTsP7eSzdZ+FMEhjjAk+m7YnGvi5LlQeT573JHUS6pBa36YgMcZULXYkFA2KjZADmDULLrsM3nmn9M27Ne5mCcgYUyVZEooGfo6Eli+HN96ADz4oX1M/bP2B3LzcIAZnjDGhY0koGhSbPw7KPky7sJv+dxPdn+3OuyveDV5sxhgTQpaEooGfI6GTTgIR+OEHKOtyJ92bdKdmfE02H9gcgiCNMSb4LAlFg5QUqF0b9u6F3bsBqFMHOneGnByXiMri6h5Xs/4P6/ld39+FMFhjjAkeS0LRQKTEmRPKetNqUlwSKTVTghycMcaEjiWhaOFnhFxFZk4At9b7R2s+YsnWJcGJzRhjQsSSULTwcyR0yilw3nlw8snla2ritxMZ+OpA7v/0/iAGaIwxwWdJKFoEGJwwcybceGP5mhpzwhha1WtFn+Z9bAlwY0xUsxkTooWfYdoV1aR2E9b8bg2xMbGVbssYY0LJjoSiRYCpezIz4YsvYNGi8jVnCcgYUxVYEooW+QMT1q2DQqfQ3ngD+veHf/yj/E3m5OXwxpI3mPD1hKCEaIwxwWZJKFrUrQsNGsChQ7B1q684f7nv8syckG/59uVc9u5l3PPJPezM3BmkQI0xJnjsmlA0SU2FXbvcKbkmbsnurl0hKQlWrXL3sdavX/bmujXuxs1pN3NS05OonVA7REEbY0zFhfRISEQGi8hKEVklInf6eT9RRN703v9WRNoUeu8ur3yliAwqrU0RSfXa+NlrM6EMfXQTka9FZJmILBGRpNDsiTLyc10oPh569HDPFy4sf5NP/eoprut5HYlxiUEI0BhjgitkSUhEYoGngCFAF+ASEelSrNq1wG5VbQdMAMZ723YBxgBdgcHA0yISW0qb44EJqtoe2O21XVIfccBrwI2q2hU4E8gO6k4orwAj5PJvWp0/v3LN2+zaxphoE8ojoT7AKlVdo6pHgCnAiGJ1RgCveM+nAueIiHjlU1T1sKquBVZ57flt09vmbK8NvDbPL6WPc4EfVPV7AFXdqaqR/ZYOMELu9NPdzxkzKtasqvL8wudp/6/2bN5vk5saY6JHKJNQc2BDodcZXpnfOqqaA+wFUkrYNlB5CrDHa6N4X4H66ACoiMwRkUUicru/DyEiN4hIuoikb9++vYwfvYICJKEhQ6BmTdixAw4cKH+zIsKMn2awds9aXl78cuXjNMaYIAnlwATxU1b89v1AdQKV+0uaJdUvqY844DSgN5AJfCwiC1X14yIVVScBkwDS0tJCO/2An/njAGrVgmXLoHVrN9dpRUwcPJHLTryMi7teXLkYjTEmiEKZhDKAloVetwA2BaiT4V2jqQfsKmVbf+U7gGQRifOOdgrXL6mPz1R1B4CIvA/0BIokobDKT0Lr10NuLsTGHvVWRaXWT7UlwI0xUSeUp+MWAO29UWsJuIEG04vVmQ5c5T0fBcxVN9nZdGCMN7ItFWgPzA/UprfNJ14beG1OK6WPOUA3EanpJaczgOVB/Pzll5QETZu6RYQyMvxW2bABKntWcNvBbcz8aWblGjHGmCAIWRLyjkhuxX3ZrwDeUtVlIvKgiAz3qr0IpIjIKuA24E5v22XAW7ikMBu4RVVzA7XptXUHcJvXVorXdkl97Ab+iUtsi4FFqhr5b+YS5pAbOxZatYLnnqt481sObKHjkx25eOrFrN+7vuINGWNMEIjNslx2aWlpml6RqQvK4/LL4fXX4d//hquvLvLW++/Dr37lbmBdurTiXYyeOpq9WXt5buhztE5uXbl4jTGmFN719jR/79mMCdEmwAg5gAED3Mw+y5a5R9euFevi5REvkxSXhFR0lIMxxgSJzR0XbQKMkANISIALL3TP33yz4l3UiK9RJAHZ0bAxJlIsCUWbEo6EAEaPdj+nTCky2XaFbDmwhSv+ewUPfPZA5RoyxpgKsiQUbUpJQmeeCQ0bws8/w+LFletq1a5VvPbDazz+zePsP7y/co0ZY0wFWBKKNi1buvuDNm2Cw4ePejsuDi66yE1s+v33levqtFanMXHwRBbesJA6iXUq15gxxlSAjY4rh7CMjgN3NLRuHfz0E7Rvf9TbmzdDjRqQnBz6UIwxprIqPTpORNKA/kAz4BCwFPhIVXcFLUpTID8JrV3rNwk1bRqabj9d9ylpzdJs7SFjTNiUeDpORK4WkUXAXUANYCWwDTfn2oci8oqItAp9mNVMCSPkCsvJcdeGguGBTx/grFfO4sHPHgxOg8YYUwalHQnVAk5V1UP+3hSRHrgpdezW+2AqZXACuOnl0tLcsO316yGmklf3zmt/Hv/46h8cV/O4yjVkjDHlUGISUtWnSnm/kuOzjF9lSEItW7rlHX75Bb78Evr3r1yXvZv3JuO2DJKT7EKTMSZ8Sjsd94iI3Oin/I8iMj50YVVzJcwfl0+k4J6hyty4WpglIGNMuJV2Emco3lo6xUwEfhX8cAxQpiMhKEhCb7/trg8FywerP+D0f5/Onqw9wWvUGGP8KC0Jqarm+SnMw/9icSYYmjSBxES3ZkMJS6medBK0awfbtsFnnwWna1Xloc8fYt76eTzx7RPBadQYYwIoLQllishRY4S9Mr+DFUwQxMS4ZVSh1FNyY8a458E6JSciPHXeU4w7Zxx3nnZncBo1xpgASktC9wGzvKHaJ3qPXwMzvfdMqJTzlNx331V+Lrl83Rp3447T7iAhNiE4DRpjTAAlJiFVnQWcD5wFvOw9zgRGqur7IY6teitjEjrhBJeA5s93R0bBdij7EAs3LQx+w8YYQxlmTFDVpRQsj23CpQwj5PL16BGaELYc2MKpL53Knqw9rLx1pd1DZIwJutKGaE8SkRMDvFdLRK4RkctCE1o1V8YjocICzHlaYY1rNaZt/bY0rd2UrQe2Bq9hY4zxlHYk9DRwr4h0A5YA24Ek3CwJdYGXgNdDGmF1Vc4kdOONMGkSvPceDB8enBBEhDcufIPkpGTiY+OD06gxxhRS2owJi4GLRaQ2cB0uAa0FxqvqyjDEV30Vnj9OtdQLPqmprtqUKcFLQgANazUs8jo3L5fYmNjgdWCMqdZKOx0XLyITgF+AK4HRwD+BC7z3Twp5hNVVSgrUrg379sHu3aVWv/hi93P6dMjMDH44OXk53PHhHQyfMpy8o28dM8aYCiltiPZjQB2gjar2VNWTgM5AWxF5Bng31AFWWyLlOiWXmgp9+8LBgzBzZvDD2XpgKy989wJzVs3h24xvg9+BMaZaKi0JnQdcr6q+tZ9VdR9wEzAGuCSEsZlyjJCD4M8lV1jzus15a9RbzL1qLv1a9gt+B8aYaqm0JJSnfpZeVdVcYLuqfhOasAxQ7sEJF13kfs6cCfv3l1y3Is5pew6ntz49+A0bY6qt0pLQchG5snihiFwOrAhNSMannEmoRQu3pENuLixYEMK4gMVbFjPqrVFk5WSFtiNjzDGttCHatwDvisg1wEJAgd64VVYvCHFspowrrBb23HPQuDE0aBCakADyNI/L372cZduX0f3L7tx7xr2h68wYc0wrbYj2RqCviJwNdMXNnD1LVT8OR3DVXgVuWO3cOUSxFBIjMbx24Ws8l/4c/3fq/4W+Q2PMMUv8XPIxAaSlpWl6enr4Oty/H+rWhaQkN+66HJPD5ebCjh3uqMgYYyJJRBaqapq/90q7JmQiqU4dd79QVhZsLfu0OV99Bc2bw5VHXc0LjZy8HMZ+PpadmTvD06Ex5phhSSjaVeCUXMeOsHMnfPyxWxcv1G6bcxv3fnIvv57269B3Zow5plgSinYVSEIpKTBggDsl924Ybif+8yl/dmsQnXpH6DszxhxTQpqERGSwiKwUkVUictQynSKSKCJveu9/KyJtCr13l1e+UkQGldamiKR6bfzstZlQWh/e+61E5ICI/Dn4eyAIKjBCDoK/4mpJWtVrxXe/+Y5TW50a+s6MMceUkCUhEYkFngKGAF2AS0SkS7Fq1wK7VbUdMAEY723bBTcjQ1dgMPC0iMSW0uZ4YIKqtgd2e20H7KOQCcCs4HzqEKjAkRDA+edDQgJ8+ils3hz8sIqLkYL/Sgs3LWTN7jWh79QYU+WF8kioD7BKVdeo6hFgCjCiWJ0RwCve86nAOSIiXvkUVT2sqmuBVV57ftv0tjnbawOvzfNL6QMROR9YAywL4ucOrgomoXr1YMgQN7P21Kml1w+WD1Z/wKkvncrIt0ZyKPtQ+Do2xlRJoUxCzYENhV5neGV+66hqDrAXSClh20DlKcAer43iffntQ0RqAXcAD5T0IUTkBhFJF5H07eG4yl9cOeePKyx/LrlPPw1aNKXq07wPLeq2oG/zvkWOjowxxp9Sl/euBH83tRS/KSlQnUDl/r7VSqpfUh8P4E7fHZAS7r9R1UnAJHD3CQWsGCqtW7uf69e7kQaxZV/LZ/hwN1y7b98QxeZHclIy6Tekk5yUHL5OjTFVVij/VM0AWhZ63QLYFKiOiMQB9YBdJWwbqHwHkOy1UbyvQH30BR4RkXXAH4C/iMitFfuoIZSUBM2aQU4OZGSUa9NataBfP4gJ8wFJ4QSUnZtt14eMMQGF8utpAdDeG7WWgBtoML1YnenAVd7zUcBcb9bu6cAYb2RbKm458fmB2vS2+cRrA6/NaSX1oar9VbWNqrYBHgf+pqpPBnMHBE0FR8gVtjMC95HuyNzBWa+cxVmvnMWOzB3hD8AYE/VCloS866dJmIwAACAASURBVC+3AnNwM26/parLRORBEclfgPpF3PWZVcBtwJ3etsuAt4DlwGzgFlXNDdSm19YdwG1eWyle2wH7qFIqODgB3MCEESOgUSNYvTrIcZWiTkIdcvJyyM3LZdP+4gfBxhgT2mtCqOr7wPvFyu4r9DwLuCjAtg8DD5elTa98DW70XPHygH0UqvPXkt6PuPwkVIEsIuJm/8nLg9deg/vvD3JsJUiMS+Sdi98hPjaeRrUaha9jY0yVYcOXqoJevdzP2bMrtPmvvdl0/vnP8EzjU1jzus2LJKDvNn+HTZprjMlnSagqGDzYzaa9cCH89FO5Nz/nHDj3XNi3Dx4ocUB6aE1eMpm059O4/cPbLREZYwBLQlVDUhJc4K0hOHlyhZp49FE3Su7ZZ+HHH4MYWznUiK9BjMRQK6EWJQ2LN8ZUH5aEqopLLnE/J092ow3K6cQT4dpr3a1Gt98e5NjK6PxO57PkpiX89cy/RiYAY0zUsSRUVZxzDjRsCCtXwuLFFWriwQehfn1o29bddhQJnY7r5Hu+N2svzyx4xk7NGVONWRKqKuLi4OKL3fMKnpJr0gR++QUef9w1F0l5msfQyUO5+f2befSrRyMbjDEmYiwJVSX5p+SmTHFjriugTp0gxlMJMRLDH/r+gXYN2nFR1xJH0BtjjmGWhKqSfv2gVSvYsAG+/LLCzajC22/D6afDwYNBjK+cRnYZybKbl9EmuU3kgjDGRJQloaokJqZgtboKnpIDl4TGj4d58+Cxx4IUWwUlxCb4nr/+w+vc/8n9do3ImGrEklBVk39K7u23ITu7Qk3ExBQkn0ceCc+id6XJ2JfBNdOv4cHPH+SzXz6LdDjGmDCxJFTVdO8OnTrBjh3w8ccVbuaMM9yccgcPwr33BjG+CmpRtwX/Of8/PHTWQ5zZ5sxIh2OMCRNLQlWNSNF7hirhkUfcKLmXXoIffghCbJU0+oTR3HP6Pb7XWTlZEYzGGBMOloSqovwk9N//wqGKL6HdoQPcdJO7RvTnP1foHtiQ2ZO1hzNePoMHP3sw0qEYY0LIklBV1L49pKXB/v0wc2almrrvPqhXz03ls21bkOILgi/Xf0n6pnReXvwye7P2RjocY0yIWBKqqoJ0Su6442DOHDcRQ+PGQYgrSH7V4Ve8NeotPrnqE+ol1Yt0OMaYELEkVFWNHu2uD82cCXsrd6TQty/UqBGkuIJoZJeRtE5u7Xv94eoPyc3LjWBExphgsyRUVTVv7u42PXwY3nsvKE3u2eNGyu3bF5Tmgmraj9M497VzGTZ5GHlasdkijDHRx5JQVRakU3L5rrgCxo6FceOC0lxQ1UuqR8OaDTmrzVnEiP23NeZYIXZ3etmlpaVpenp6pMMosHOnm5VUFTZtgkaVW0L7m2/czEBJSe4aUatWQYozSLYd3EbDmg19axHtPrSb+jXqRzgqY0xpRGShqqb5e8/+pKzKUlJg0CC3SNDbb1e6uZNPdpeasrLgL38JQnxB1qhWI18C2pO1h97P9+aGGTfY/UTGVGGWhKq6IJ+SGzcOEhPh9ddhwYKgNBkS6ZvSydiXwYJNC2yuOWOqMEtCVd2IEW5o25dfwvr1lW6uTRv4/e/d89tui64bWAsb0HYAX137Fe9e/C414qNwaJ8xpkwsCVV1tWvDsGHu+ZQpQWnyL39x9w998UWlVowIuZ5Ne5JaP9X3+i8f/4VxX4yzIyNjqhBLQseCIJ+Sq1cPnnkGPvoITjstKE2G3IrtKxj3xTjumXsPS7ctjXQ4xpgyivAizyYohgxxmWPxYlixAjp3rnSTo0YFIa4w6tywM9MvmU7GvgxObHxipMMxxpSRHQkdCxITYeRI9zxIR0OFpafDrl1BbzbohnYYyo1pN/pef7f5O95a9lYEIzLGlMaS0LGi8Cm5IF4TefRR6N3b3cRalew/vJ+Rb41k9NTRvL2s8sPXjTGhYUnoWHHWWW4G0lWrYOHCoDU7YICbou7JJ13TVUXthNr88eQ/ckrLUxjWcVikwzHGBGBJ6FgRGwsXX+yeB/GUXI8ecNVVbiXxSy5xq0dUBSLCb/v+ls+v/pykuCQAMrMzeWHRCzYJqjFRxJLQsST/lNyUKW4WhSAZPx7atnXXhs4/382ZWlXExsT6nj/8+cNcP+N6rpl+TQQjMsYUFtIkJCKDRWSliKwSkTv9vJ8oIm96738rIm0KvXeXV75SRAaV1qaIpHpt/Oy1mVBSHyIyUEQWisgS7+fZodsTYXLyye5u002bYN68oDXbqBF8+KGbpm7uXLj00qDmuLDp07wPreu15qa0myIdijHGE7IkJCKxwFPAEKALcImIdClW7Vpgt6q2AyYA471tuwBjgK7AYOBpEYktpc3xwARVbQ/s9toO2AewAximqicCVwGvBvPzR4RI0O8Zyte2rVv8rl49+OEH2L49qM2HxYhOI/j5tz9zcouTfWXjvhjHy4tfthtcjYmQUB4J9QFWqeoaVT0CTAFGFKszAnjFez4VOEfcDJUjgCmqelhV1wKrvPb8tultc7bXBl6b55fUh6p+p6qbvPJlQJKIJAbt00dKfhKaOhWOHAlq0926uSOiL75wR0VVUXxsvO/56l2rufeTe7lm2jUs2bYkglEZU32FMgk1BzYUep3hlfmto6o5wF4gpYRtA5WnAHu8Nor3FaiPwkYC36nqUVc7ROQGEUkXkfTtVeHP/xNPhK5d3Y09H34Y9OZ79y5YBlwVvv026F2ETdv6bfn3iH9z12l30a1xN195Tl5OCVsZY4IplElI/JQVP+cRqE6wykuNQ0S64k7R/cZPPVR1kqqmqWpaw4YN/VWJPiE6JVfcHXe4y1AvvRTSbkJGRLi82+U8fM7DvrKl25bS/l/teWf5OxGMzJjqI5RJKANoWeh1C2BToDoiEgfUA3aVsG2g8h1AstdG8b4C9YGItAD+C1ypqqsr+Dmjz5gx7ud770FmZsi6yV/07vrr4b//DVk3YfXMgmdYt2cdn6z7JNKhGFMthDIJLQDae6PWEnADDaYXqzMdNygAYBQwV90V4unAGG9kWyrQHpgfqE1vm0+8NvDanFZSHyKSDMwE7lLVKJ4rugKOPx769IGDB2HGjJB1c+ut8Ne/Ql6ey3ufHAPf208MeYJJQycx9uyCKSIy9mXYwnnGhEjIkpB3/eVWYA6wAnhLVZeJyIMiMtyr9iKQIiKrgNuAO71tlwFvAcuB2cAtqpobqE2vrTuA27y2Ury2A/bhtdMOuFdEFnuPyq2PHU3CdEruvvtcMjpyBIYPd/cSVWWxMbFc3+t6kpOSAcjTPEZPHc0JT5/Azzt/jnB0xhx7xIamll1aWpqmV5Vv2c2boXlziI+HrVshOTlkXeXlwRVXwBtvuHWI5s2DTp1C1l1Ybd6/mYGvDmTrwa388odfqBlfM9IhGVPliMhCVU3z957NmHCsatrUzSd35Ai8+25Iu4qJgZdfditKNGzo1tk7VjSt05TvfvMdcy6f40tAOXk5nP3K2Uz8ZiJHcoM7DN6Y6saS0LEsTKfkwB1wTZ3qjoJatAh5d2EVHxtPz6Y9fa9n/jSTT9Z9wr/m/4sYsV8hYyrDfoOOZSNHuuwwdy5s2RLy7mrWhBTvDixVtzprVZnwtDyGdRzGe6Pf45+D/klcjBuQmZmdyb1z72Xz/s0Rjs6YqsWS0LGsfn0YPNhdtHkrvIu7PfQQ3Hxz1ZvwtCxiJIYRnUYwvONwX9mLi15k7LyxjJ46OoKRGVP1WBI61oXxlFxhl11WdMLTnGN8EoJTW53KBZ0u4LZ+t/nKdmbuZMlWmw7ImJJYEjrWDR/uzpN98w2sXRu2bo8/3k14mpzsxkXceGNQF3yNOj2b9uTd0e9yfqfzfWUTv51It2e78fDnD5ewpTHVmyWhY12tWjDCmzd24sSwdt2tG/zvf1CjBrz4ojs62rkzrCFEVJ7mkRSXxJltzvSVbdi7gR2ZOyIXlDFRxpJQdXDzzW4c9cSJ8M9/hrXrU0+Fd95xiWjyZLj77rB2H1Fjzx7Lxts2cmqrU31lD3z2AE0fa8p/vv9PBCMzJnpYEqoOTjutYJbRP/0Jnn8+rN0PGeLWILrwQjdgId+xfHouX4MaDYq8Pph9EFWlT/M+vrKFmxby9YavbU0jUy1ZEqourroKnnjCPf/Nb9wS4GHUrp07IsqfiPzIERgwAF5/vXoko3yTR05m65+30um4giklHvjsAU556RReWPRCBCMzJjIsCVUnv/0tjB3rvvWvuMJdsImQKVPcyLnLL3djJzZujFgoYZdSs2A5K1XlhEYn0LJuS0Z0KljzcdqP03hx0YvszdobiRCNCRubO64cqtTccYGowp13wiOPQGIizJrlpveJQBgvvwx//CPs3Qt168Jjj8G117pVyqsbVUUKffC+L/Rl/sb5TB45mTEnjIlgZMZUns0dZwqIwLhxbsz04cMwbFhElkcVgV//GpYvdyHs2+fWJRo4MKwjyaNG4QSkqtzY60bOa38eIzoWHB098OkDXPT2RSzYuCASIRoTEpaEqiMReOopN2b64MGCkQMR0KwZTJvmZuBOSYGPP4aPPopIKFFDRPj1Sb9m5qUzqRFfw1f+5rI3mbp8KpnZBQsVLt22lK83fE1uXm4kQjWm0iwJVVcxMfDvf7t7iHbvhnPPhZ8js16OiJvYYflyd8nquusK3jt4MCIhRaWZl87k6fOe5pSWp/jKJnw9gVNeOoV/fh3eoffGBIsloeosPt6NEDjnHLfm0IABsH59xMJp1MjdR5R/Zmr1amjdGv7xj2N/2p+ySK2fyk29byI+Nt5X1jq5NW3rt+Xc48/1lb2y+BX6vdiPN5e+GYkwjSkXS0LVXVISvPce9OvnEtCAAS4hRYH33nMzLNx+O5xyCnz2WfUazl0W951xH6t+u4pujbv5ymb+PJNvMr5h16FdvrINezfwn+//w7o96yIQpTGBWRIybhW699+H7t3dKblzz3Wn6CLsT39yYbVoAQsWwJlnQpcu8PjjsGtXqZtXGyJSZGDDSyNeYtqYaVzQ+QJf2dvL3+aq965i7OdjfWVZOVl8uf5LsnKywhqvMYVZEjJOcjJ88AF07OgGKQwZEhWLAQ0ZAsuWwb33usVif/zRDev+618jHVn0qp1Qm+Edh9OkdhNfWet6rRnRcQQD2w70laVvSue0f5/GqS+dWmR7m9vOhJMlIVOgUSP48EN3Iebbb91iQFmR/yu5bl148EH45Rc368K558INNxS8P2MGPP20u9/I+Deyy0jeG/Meo08oWO/o4JGDnNjoRE5ufrKv7MCRAzR5tAkd/tWBnLyCC3F2P6EJFbtZtRyOiZtVy2LVKujf363GOmyY++aPjy99uwjp2xfmz3crVlx6qZuVKM3vbXHGnzzN8y1TvnjLYvr/uz8dUjqw8IaFvjqnvnQqMRLDlJFTaF63OQC5ebnExsRGJGZTtdjNqqZ82rVzR0QNGrjDjKuugtzovA9FFW67zU36kJkJL7wAvXtDr15untYDByIdYfTLT0AAPZr0YM8de/jfJQVTOu0/vJ9vMr7h6w1fF5ly6NfTfk3rx1sze9VsX1lmdibZudnhCdwcEywJGf9OOAFmz3aDFiZPdstBROFRswiMHu3mofvxR5eQGjSARYvcKbswz9N6TIiNiaVpnaa+13US67D5T5v59OpPSYpL8pX/uONH1u9dT3JSsq/s+YXPU/vvtXng0wd8ZYdzDrNx30Y7pWf8siRkAuvd2x0JJSXBpElucaCXX47aO0g7dnTzz23cCK++CoMGwZhC06798Y9uqqCXX3ZTA9l3Ytk1qtWI01qdVqTsq2u/YsUtK+jRpIevbNP+TRzJPVJkCYuFmxfSYkILznql6ByFr/3wGnPXzi1y7clUP3ZNqByqzTWh4mbNct/m+/a513Xruosv110HPXtWiRlHVaFJE9i2raCsZUs4/XQ44ww32KF168jFdyzZf9iNqqyTWAdwM4JfM/0aBrcbzOsXvg7AoexD1PxbTeJi4jh09yHiYuIA+Pu8v7PlwBZu6XMLHVI6ROYDmKCza0KmcoYMcYcXL77obmrdtw+efdZd/e/Z081DFwX3FZVm1iy3sOyIEe6U3YYNbj2j4qftNm50w8Lt77OKqZNYx5eAAEZ0GsHO23fy0vCXfGWZ2Zlc1f0qLux8oS8BgZsf74n5TxRZwuKxrx4jdWIqz6Y/6yvbdWgX7yx/h8VbFof405iQU1V7lPHRq1cvNaq6ZInqH/6g2qCBqvuuVk1KUr38ctVPP1XNy4t0hKXKzVX94QfVf/1LddQo1fT0gvcefNB9pOOOU73gAtXHH1f97jvVw4cjF291MWPlDH3ki0d0b9ZeX9n1069X/oo+8c0TvrK5a+Yqf0X7v9S/yPYD/jNAz59yvh44fMBX9uP2H/XH7T/qoexDof8Axi8gXQN8r0b8i70qPSwJFZOVpTpliuqAAQXJCFTbt1cdP151y5ZIR1ghf/+7arNmRT8SqMbGqg4bVlAvL0/1gw9U16xRzcmJXLzHusM5h3XljpW67cA2X9m3Gd/q8MnD9e6P7/aVZR7JVP6Kxj8Yr7l5ub7yX73+K+Wv6PQfp/vK3lvxng56dZBOSp9UZPv/rfyfLtq0KMSfqPopKQnFlXakZExAiYluaNro0bBmjZuV+6WX3NQ/d9zhZiMdNsxdOxo0CGKrxj0ld97pwl+9Gj7/3M1Z98UXbjBDYmJBvc2b3bUkgIQEOP54aN/ePTp0gKFD3VIVpnISYhOOuj7Up3kfpo2ZVqQsPjaer675ip2HdhYZdt6kdhPa1m9bZMTf8u3LmbN6Dt0bd/eV/bL3F4ZOHkq7Bu34+bcFM8oPfm0wew/v5e2L3qZF3RYAzPp5Ft9v/Z5Bxw/ipKYnAbDv8D4y9mXQsGZDGtZqGLwdcIyzJGSCo21beOghuP9+N7T7hRfc8uH//a97NG/uZiFt1859S7dr5x5NmkTlwAaRghCvucaVZWUVve/owAE3n91PP8GmTbBihXvkmzevIAlNnAhz5kDjxkc/WrZ0SctUTlxMHP1a9juq/IXhLxxVdumJl9KtcTdaJxeMRhGEQccPolmdon85pG9KZ+ehnUWuXf33x//y/KLnSU5K9iWheb/MY+jkoQxuN5hZl80C3JmmTk91onZCbeZfN993c+8zC55h1a5VXNfzOjo37AzAmt1rWLZtGe1T2tPpuE6Au5E4MzuTmvE1iyTWY4klIRNccXHuEGDoUHeo8MorLiGtXg1vv310/Vq1Cr7t8x/5SappU7fuUZRISnKPfB06wCefuOcHD7qJJn76yR0I/vwzdOpUUPfTT93ACH9OPdUdaYFb7LZPHzeDUvFklZzsZodo4k0Jl5kJeXluF0ZhHo9qrZNbF0lAAB2P68jsy2cfVfeb675h64GtHFfzOF/ZkHZDSE5KplfTXr6yuJg4Oh3XidTkVF9ZZnYmP+38iaS4pCKzS7yz4h0+Xvsxg9oN8iWhmT/N5Hezf8ctvW/hyfOeBCBjXwatH29Ny7otWf/HgmVWLnjzAjbs3cCUUVNo16AdAFOXT+XD1R8ysstI39IeWw5sYfrK6bSs25Ih7Yf4tv96w9fExsTSs2lPX3Ldd3gfeZpH7YTaRRJuqIW0JxEZDEwEYoEXVHVcsfcTgf8AvYCdwGhVXee9dxdwLZAL/E5V55TUpoikAlOABsAi4ApVPVKRPkyQNG3qzm3dfjt89507TFi1quhj5074/nv3KK5GjYJzXMcf75ZerV3bPerUKfqz8POkpLB/K9eq5SYh797d//tjx7ojqq1bj34U3mb79pIXuX3vPTe6D9w9Uffd5z5qnTruUbeu+9mihZttKd8jj7gEl5TkTinmJ9SkJDjxROjsvgfZswfWrStaJ/95QoI7o1rdEl67Bu18X/T5Luh8QZFZygEGtRvEinYripQlxSXx060/ceBI0ak7bul9C4OOH+Q74gFoWa8l57U/r8iyHIeyD1Ezvia1E2oX2X7J1iWs3r3aXdj3fLXhKyYtmkSHlA6+JPTjjh/5zf9+w+mtTy+ShM597VwOHDnA3jv3UjexLgA3zbyJN5a8wWsXvMZl3S4D4J3l73DjzBsZ1XkUzwx9pmw7rJxCloREJBZ4ChgIZAALRGS6qi4vVO1aYLeqthORMcB4YLSIdAHGAF2BZsBHIpJ/wiJQm+OBCao6RUSe9dp+prx9qGp0zk9TlcXEuHl0evU6+r3du49OTD//7H5u3w5Ll7pHecTGFiSmwskpIcHNgRcf747Y8p+XVFa4PDbWfZb8h0jR18Ufhd7vGhNDVxFIAY4TOMH7Jhdxj/+5141yYlj0RF227k5g694ktu5JZOueBLbtSWRvZhwtM1bDB/tBBF3VlpqJbcg8HMu+fW7k/MaNrtlt6w/Bx1/5Msb4sf3Ztd///H9jr13D3ZevBxHmfn4cI+/rGnDXbp/+Ncclu5tLL7y7M599X4/4OCUuVomPVeLj8oiPU87tvYcJv10DwJad8Vz6YCdiY5XYGCU2Bvcz1j2/58r19GjvboB+59MU3v+mge+92BglxtumYXI2f7ligy+Wh//TkqwjMcTEKDFCkZ9n9thLvxPc/Uo/b0hi5tcNfLs6RtT3XAR+PWQLSYnuy3zWN/XZtCOh0PuK4P4ZWzfO4vQe7l65A5kxzPgqpeCfEC3yz3l69700bpBNLHB4TU3WbkhiHcsK1e9AB+nAym/20Kr3OgDOpwFJG8eR/Ysw44OlvmT/Vuy35O3PY/Xk+Rzf3E0o/HKjv7EwM4/lL+xllSwBoGVeT34bfy815rUmO2se8XFKo4PruaDmhaRsbs0Hjy3x7bu2uZ3Jisnkp9eWk9b+iPu12bKHmtTm51nb+Ogj95fQN9nL2ZG1g7VLMiDla3eLRrAFGrFQ2QfQD5hT6PVdwF3F6swB+nnP44AdgBSvm18vUJveNjuAuOJ9l7ePkj6TjY4Lsz17VBcuVH3zTdW//U319ttVb75Z9Yor3NjpgQNVTz5Z9YQTVNu0UU1JUU1MPHpY2zH+yCZWd1NP19NCl9JFv6avfkm/InXGcbvezUP6J/6ht/KEXsckvZz/6EW8qe9yvq/ebM7V7nynHVmhrVmrTdikyezSJDI1lmzdTT1f3bP5KGBYF/Gm78Uq2pb4ET6gYHTl7YwLWK8tq4oUJLMrYN2/c4fvxVuMKrH/XST7XpzDhwHrjeKtMn+mDzmnTJ8pldVl/kzjuD2kn2k1qUfXic1Sam3VhjWWqzZuXOFfZSI0Oq45sKHQ6wygb6A6qpojIntxfys2B74ptm1z77m/NlOAPaqa46d+RfrwEZEbgBsAWrVqVeIHNkFWr567GbZnz/Jtl53tRg3kP/bvdz+PHHHv5T9ycoq+Lq0sL+/oh6r/8uLv508Am//7Xfh58dclPc9/XUicKslAcvH3tGCanDt0ftF9VKwN9HQABpHFIH7v5/381yf6iqbljicrbwI5Gku2xpGtsWTnxZGjsdSOPQRJbpqfprmJfLTv9+QSS67GFPzUGHKJoVvdGpDg6o7av4oOB8eRSwy5WrheLHVjD0KTgqmD7s54g0N5ieQh5GlMoZ8xnJx8EJJd3XYH6vC7rW+jiPeAPI3xPU9M7QWxhwEYnLGS1oeyiryvKuQRQ986G6GZa7P2kfqMWfORrw7gnqt3RNuqBdRydbtuzeb8nZ+DuDp4bQI0TtgF7fv7PtOgFYs4mJd0VD1FaNs4ERq6uk32NmbwhoKvr/z6+eI69YG4QwD0WruNmAML8Kdbvb3QyrVZ40gK56z0PyNM/bj9cPLJft+rrFAmIX9njrWMdQKV+7tKXVL9ivRRtEB1EjAJ3LQ9frYx0SY+HurXdw8TMrW9R2lqAueUsc3e3iOwO3zP/lxivct9z07yHoGdX8Y2zwX+AEBjYHKJdQv+ALjSewQ23Pes5Pl2C5JAfyDAOBdPweXt8SXW6w1cB0BT4KMS655V4rsVFcqhRxlAy0KvWwCbAtURkTigHrCrhG0Dle8Akr02ivdV3j6MMcaESSiT0AKgvYikikgCbhDA9GJ1pgNXec9HAXO984fTgTEikuiNemsPzA/UprfNJ14beG1Oq2AfxhhjwiRkp+O86y+34o4LY4GXVHWZiDyIu0g1HXgReFVEVuGOTsZ42y4TkbeA5UAOcIt6o9b8tel1eQcwRUTGAt95bVORPowxxoSHLeVQDtV2KQdjjKkEW8rBGGNMVLIkZIwxJmIsCRljjIkYS0LGGGMixgYmlIOIbAd+qUQTx+HuaYpWFl/lWHyVY/FVTjTH11pV/S6yZEkojEQkPdAIkWhg8VWOxVc5Fl/lRHt8gdjpOGOMMRFjScgYY0zEWBIKr0mRDqAUFl/lWHyVY/FVTrTH55ddEzLGGBMxdiRkjDEmYiwJGWOMiRhLQlFERPwttGcqwPZl5dj+qxzbf2VnSSi61BMR379JtP1HFpG+IjIg0nEEIiLHiUg9AI3Ci522/ypHRM4QkdGRjiMQEWkmIi3B9l95WBKKEiIyHPgQeFJE/gXR9R9ZRAYBzwK7i5VHRaIUkfNxCxs+LyKviUjdSMdUmO2/yvH23wRgQ6Rj8UdELgQ+BV4UkakicpKIJEY4LJ+o3n+qao8IP4AOwBLgHKAN8BnwFpDovS8Rju8sYBtwkve6RqT3WbH4WgJfA32BBG/f/QfoHOnYvPjOtv1XqfjOBPYBJ3qvawM1Ix1XofhSgA+AXt7rR3B/cAwBEqIgvqjef3YkFB32Az8B36vqOlU9A7fq7WsQ2SMiEYkFmuL+gsoRkVrAyyLyqoi8ISJxXr1I/kW/H/dLdkhVj6jqxcBO4C4RqRkF8TXH9l9l5OKOIBt7RxdTgNdE5L8i0j6CceU7jEveDQFU9XbgZ2AE7g/MSO8/JYr3nyWh6JAJ7AV65Reo6oW4/zSPRiwq3cDYBgAADMFJREFUF0cu8A4wDngCWAN8CzwA1AXe9uqFPVGKSF0vSR4AvgF6Frqm8UegFvBMBOPrJSKpqvoq8DeibP/lU9U9wHygVzTtv0LxzQMuB/4FbAVmAdcCm4F/RiouEWksIjVV9QDuD8aeItIWQFUfAw4Bf/VeR+L/32ki0k1VPwMuw/3/i5r95xPpQ7Hq+gBOB/4MDATqA4OARUCfQnW6AP+MUHynAL/14mrslV0J3FqoTk1gGlA7AvENBV7N7xv4lRfL+UB9r6yOVycpAvGdC2wC7vNeJwNXRNH+G4hLjL8G2gLdgZnABVGy/84A7sd9WXb0ys4C/lCs3mygSQTi+xUwFTjOe90XdwruZuD4QvVmAg0j9O+7A3gF79Sbt09/Hw37r/Ajrgx5ygSZiJwLPI07whgG3AFcBzwMPCMiDwALcIkgTURqqOqhMMY3FHfk8wnQw4thHO6v9iOFql6ES6BhJSKDgQeB29T9FYqqzhSRGrgvrSYikg50xCXysP4/9/59/4b79+0lIqKqe0TkHdypm3yR2n9nAE8Bz+OuQd4NXIiL+S7cEXgk999QL5b3gM5APxG5Q1U/EZFPC9W7FJfcw/a74fU7BPe7equq7gBQ1W9FJB73h1oTEVkI1ABaU/R3JhzxDcLtvwm4/VcLd7blc++RXy8i++8okcyA1fUB/B/uCxSgHnAD8D3uC2Eg7hTIe0A60C3MsXXEHZGlea8HAXOAuoXqCO6XbQnQNczxtfX2Vf7+a4D76/103F/uJwH3AO/jRiv1CHN8fXHXA/p4r+cAjxWrI7ijorDvP6//64BHCr0eg7sG1A04Abg3gvuvNTAXONl73RX3x0+LQnVivP23LAL//5riBg495r2uD1wNXIwb4NHUez0LmI43GCWM8XXzvjf6ea+nApOL1YnY/vP3sCOhyDgI9ARQ1b3AJO+65dPAJcCXuL+i4lR1a5hjW487d7zUi2+OiPwWd7pmnlcnBfdlcbGqrghzfJtwQ9nriMhFwO3AStwvVk3cX6djvQEACaq6O3BTIbEUGKWq33uvHwEuFpGGqrrdK2sCtCLM+887IlPcdYG++eWqOsX7//cmMERVH4rg/tuM+yPsOy+2Zd5F/dNwF9TB/W7UwO3ncP//24HbT41F5E+439dFuMETfwGuUNWXRWQy7vf3YJjj+wm4TFVXeq/vBB4QkR6qutgrqwckEpn9d7RIZ8Hq8sCNkGrlPa+J+2v+b4XeT8GdU+4XwfhSi5XlDxGfBpzrPU/DXVCPjVR8uJFID+GOJG70yhrjTjFdEC37zytvhPuj4rfFyiM2dBd3euYbYEKx8r8BIyMVV6E48idWjvV+Pgdc6j0/kwhfw/D+/92IO7V1S6HysYVfR8PD+119E7i7WHlYf39LetjouDAQkVHAu8BUEfk77ihoMHCK9xpV3Yk7935SBON7U0Qe8m68Q1Xzr19kAFu8G2rH4y505kYovoeB/sB9wM2q+qwX61YgFvelH1bF4ntQRC7If09Vt+GuuVz0/+2dfYxdRRmHnx8tBbW0pdYYqoBYqEaqtJjwUdI2McQYbFpDo1ggjQENiIBICIL4UdNoTEosaEDU8iFpEyImRUsEESlaiE21lLYIamuqRVBIJHVr5Kv29Y93NnusW9jde/ecs9vfk5zsnnPnzj579+59z8y8MyPpXZXrtY0TSDpL0p1lzILIu/MPAfMlragUPZwcQ6iV/f34/6zdp4HnKmNF45r0K3+7W4GrIuKmStFxZJdwrezvVzJGAYiIHnJ89zxJsyvXa/v/fV2ajoKj/SCbvr8mA8/byOyZ7wLnkv3HW8g39PVkt9IJLfD7DnBhpcwyMoV3I2XCW8N+3yO7ParlPgY8DhzfAr/9X78jyb75+Q28/04lP8R/BawFDq08NonsYl0J3AY8Sc0TVA/kR+VOncwi/TOZrDOjDX79lDsH2ARMb4Mffa3IQ8rXbwMX9Z636XBLaPgZA+wFdkfEM2TT+EFyIP1o4HRyIHYn2ZW0vQV+D5FzHhaUMhOB9wKLI2JbC/x+DpxW7owpY0NfJvvCd7TAr/f1mw8QOa5yHzkQXDdjyZuIeeRY0JrKHf1uMvHkTjIT8iNR/xhBv34R8Z9Ky+gVsgt7cUQ80QY/6GtxlGzNzwJLIuKPbfArr9+YiNhXyv0C+GnlvDV4U7sakLSUnDl9RUQ8L2kKmZ3y5oj4YqNyvKbfxIhYqlxH7KjoG+xsi9+EiPhqmdU/JSJ2tdGvCacqko6IiD0l2eAGYCp5w/OKpMkR8UKL/Y4AXiXffztb6HckueLElIj4ewv9pkbEs014DRQHoWGg3AHPJlsQXyHTiBeT/cXLI+I5SccCq8kMlVrfvIP0+3hE/LXFfueUFkhb/RZFzRmOFb8J5Iz9nihjUOVDfQW5ftgvyZT8L0TEv1vqdyKZjl/nGNpg/N5Djg293H9tjftNJ/++zc4Feg3cHddlJL2fzHLbQKaR3ki+UTeQd0w3lQHq08pTXmq5379a7ldrCuwQ/Gr7cOrH743kWMBZ6luOZ09EfJL8cPoacHvNAWiwft+vOQAN1m9lzQFosH53tDkAgYPQcDAdeCAifhIRFwCPkEu4jAPuALaSTeaLgctLv7z97Ddcfg+TS8zMVd9iqQvIu+g50TefyX72a4amMyNG20HO6L8fmF25dgk50W5COR9PA+tx2e+g9fs0uT3DxHI+h5qzMO13cPgN5fCYUBeQNJPsdlFEPFnmsvQAP46I35cyq4CnI+Ja+9mvIb+/RMR19rNfm3B3XIcoFzNcC3wG+KGkReS8n3cCCyXNLUU3kmMG9rNfU351LyFjv1Hu1xWaboqN1INchHI8udDjgnLtdOBP5MS1Y8nMlXVkV80uapzoaT/72c9+I+FoXGCkH+SWAufTN1P5FHJ299nl/O3kdg3H2M9+9rPfaPLryu/YtMBIP8hB6dv5360O5pDLqU9ryst+9rOf/UbC4TGhISLl2vcRcTOZr3+LpIllyYz1ZKpuY4sE2s9+9rPfSMDZcYOgTEKcTN6F7IvKSrSS7iJ3KNxArud0JTAvalxtwH72s5/9RhoOQgNEub3B14FnyvFbcjZyT6XMBeS6TScBSyOitgUr7Wc/+9lvJOIgNACUq+auAr4VEY+WNMnTyNz95ZG7o1bLHxb1LuVhP/vZz34jEo8JDZwJwAnl+zXAveRSLYsBJJ0i6eTyeG1rXVWwX2fYrzPs1xlt9xs2HIQGQES8CnwTOFvSnMg9OR4hN1GbK+kNwBnAs6V8rc1L+9nPfvYbsQwlpe5gPMitjy8ld/WcW7n+MC1IlbSf/exnv5F4jB16+Dq4iIiXJK0GArhW0rvJPtu3UPN2B/1hv86wX2fYrzPa7jecODFhkEgaRzaNLyL3AroxIjY3a9WH/TrDfp1hv85ou99w4CA0RJT7y0e0cM92sF+n2K8z7NcZbffrJg5CxhhjGsPZccYYYxrDQcgYY0xjOAgZY4xpDAchY4wxjeEgZIwxpjEchIypCUmTJF1Svp8q6UddrPsKSUu6WN9dkk54/ZLGdIZTtI2pCUnvAO6NiBldrncs8BhwckTs7VKd84DzI+JT3ajPmAPhlpAx9fENYJqkxyXdLekJAEmfkHSPpLWSdkq6VNKVkjZL2iBpcik3TdL9kjZJWl+WdgH4APBYbwCSdLmkJyVtVW6GhqQ3SbpN0m9KvQvL9TGSrpe0rZS/rNS5HjizBDhjhg2/wYypj2uAGRExs7dVVHlsBjCLXMhyB/D5iJglaQWwBLiBXNzy4ojYLulU4GYyAJ0BbNrv5xwXES9LmlSuXQc8FBEXlGsbJT1Y6j4OmBURe3sDXkTsk7SD3ECtWrcxXcVByJh2sC4i9gB7JP0TWFuubwPeJ2k8MBu4W1Lvcw4rX48CnqrUtRVYLeke4J5y7YPAAklXlfPDgWOAM4FbeltREfFCpZ7nyZ08HYTMsOEgZEw7qO6Uua9yvo/8Pz0E2B0RM/t57otkUOnlw8BcYAHwJUknAgIWRcQfqk9URrQDDQwfXuo2ZtjwmJAx9bEHOGIoT4yIHmCnpI9CBg9JJ5WHnwKOL9cPAY6OiHXA1cAkYDzwM+CyEnSQNKs89wHg4t6xn97uuMJ04HdD8TVmoDgIGVMTEfEP4NGSkLB8CFWcB1woaQsZHBaW6/eRLR+AMcAqSduAzcCKiNgNLAMOBbaWn7+slF8J7CrXtwDnAkh6K/BiRPxtCJ7GDBinaBszCpC0Brg6IrZ3qb7PAT0RcWs36jPmQLglZMzo4BoyQaFb7AZ+0MX6jOkXt4SMMcY0hltCxhhjGsNByBhjTGM4CBljjGkMByFjjDGN4SBkjDGmMf4LtI7zGNjy7/YAAAAASUVORK5CYII=\n", "text/plain": [ "<Figure size 432x288 with 1 Axes>" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "#discharging\n", "\n", "y0 = V*C\n", "y11 = odeint(model,y0,t,args=(C,R1,0))\n", "y12 = odeint(model,y0,t,args=(C,R2,0))\n", "y13 = odeint(model,y0,t,args=(C,R3,0,))\n", "\n", "# plot results\n", "plt.plot(t,y11,'r-',linewidth=2,label='R1=100$\\Omega$')\n", "plt.plot(t,y12,'b--',linewidth=2,label='R2=200$\\Omega$')\n", "plt.plot(t,y13,'g:',linewidth=2,label='R3=300$\\Omega$')\n", "plt.xlabel('time(sec)')\n", "plt.xticks(rotation=45)\n", "plt.ylabel('Q(C)')\n", "plt.legend()\n", "plt.show()" ] }, { "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": 4 }