{"nbformat":4,"nbformat_minor":0,"metadata":{"colab":{"name":"fiche python Résolution d’un système d’équations différentielles avec odeint.ipynb","provenance":[],"collapsed_sections":[],"authorship_tag":"ABX9TyMnpMy3A+MUu8t7W1XMAE/w"},"kernelspec":{"name":"python3","display_name":"Python 3"},"language_info":{"name":"python"}},"cells":[{"cell_type":"code","source":["# mécanisme réactionnel en 2 actes élémentaires :\n","# A => B (k1)\n","# B => C (k2)\n","\n","#système d'équations différentielles\n","#da/dt=-k1*a\n","#db/dt=k1*a - k2*b\n","#dc/dt=k2*b\n","\n","import numpy as np #Pour faire divers calculs\n","import matplotlib.pyplot as plt #Pour tracer des graphiques\n","import scipy.integrate as integr #Importe la bibliothèque pour faire la résolution d'équation\n","\n","## Définition des données pour la réaction étudiée et de l’état initial\n","k1=10\t\t#valeur de k1 en s-1\n","k2=1\t\t#valeur de k2 en s-1\n","ABC_0=np.array([1,0,0]) \t#valeur des concentrations initiales en mol/L\n","\n","## Définition de la fonction F image du système d’équation différentielle\n","def F(Y,t):\t\t#renvoie la valeur de la fonction F pour les 3 concentrations \n","\treturn ([-k1*Y[0],-k2*Y[1]+k1*Y[0],k2*Y[1]])\n","#Attention localement : [A]=Y[0], [B]=Y[1], [C]=Y[2]\n","\n","## Définition des np.array du temps\n","tmax = 2\t#tmax en s\n","nbpoints = 150\n","T=np.linspace(0,tmax,nbpoints)\t#t est pris entre 0 et tmax (min), par pas de 2s\n","\n","ABC=integr.odeint(F,ABC_0,T)\n","#résolution du système [A]=ABC[:,0], [B]=ABC[:,1], [C]=ABC[:,2]\n","\n","A=ABC[:,0]\t#extraction de la première colonne de abc pour avoir [A]\n","B=ABC[:,1]\t#extraction de la première colonne de abc pour avoir [B]\n","C=ABC[:,2]\t#extraction de la première colonne de abc pour avoir [C]\n","\n","plt.figure(1)\n","plt.xlabel('t en s') \t#légende de l’axe des abscisses\n","plt.ylabel('concentrations en mol/L') \t#légende de l’axe des ordonnées\n","plt.plot(T,A,'k-',label='[A](t)') \n","plt.plot(T,B,'r-',label='[B](t)') \n","plt.plot(T,C,'b-',label='[C](t)') \n","plt.title('évolution des concentrations en fonction du temps') #affiche l’étiquette de la courbe\n","plt.legend() #affiche le titre de la courbe\n","plt.grid() #affiche le quadrillage\n","plt.show() #affiche le graphique"],"metadata":{"colab":{"base_uri":"https://localhost:8080/","height":295},"id":"ZdW7dDf_k86m","executionInfo":{"status":"ok","timestamp":1641476337275,"user_tz":-60,"elapsed":470,"user":{"displayName":"Lionnel Malara","photoUrl":"https://lh3.googleusercontent.com/a-/AOh14Gg9uzWMkTid1qlQikWp7IMi29R7n5oSh3ECmY5O=s64","userId":"16342666118884215780"}},"outputId":"d3770b6e-58d1-424d-98aa-8f71694f9cf9"},"execution_count":null,"outputs":[{"output_type":"display_data","data":{"image/png":"iVBORw0KGgoAAAANSUhEUgAAAYIAAAEWCAYAAABrDZDcAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nOydd3gU1frHPy8phBJ6DS2RpnQJAiogWAC5YkFURLCgooA/Ua8N9KqIil4bKnotiCIqoGIBBbEROtKkCQYpoReRXgIkOb8/3llcQsom2dndZM/neebZnbIz352dnXfOOW8RYwwWi8ViCV+KBVuAxWKxWIKLNQQWi8US5lhDYLFYLGGONQQWi8US5lhDYLFYLGGONQQWi8US5oSlIRCREiIyV0T+5Yd9xYuIEZHIfH7+JhH5oaA6fDyWEZF6gThWuCAitUXksIhEBFtLfhCRqiIyS0QOicjLATyuve5DiLA0BMA7wCvGmO8CedCsjIYx5hNjTOdA6ghHCmqwvfaTIiKXeuaNMZuNMaWNMekFVxkU+gN7gDLGmH+7cYDCet3765opDBT5L5gVxpibg63BEnqISKQxJi3YOgJMHWC1sZGl4Y0xJmwmIA6YBPwFbATu9Vp+DKjgte256JNSFNpyehzYBOwGPgLKOtvFAwaIdOZTgEu99vMU8LHzfrOz7WFnOh+4FZjjtf0FwCLggPN6gde6JGA4MBc4BPwAVMrh+z4E7AC2A/2cY9dz1hUHXnI07QLeBko46yoB3wL7gb3AbKBYNsdoDPzobLcLGOq1/5HOsbc774s76zoCW4F/O+dzB3Cb1z5LAC875/sAMMdLW1tgnqNtOdDRl/OTw7mfC7wK/A08A9QFfnHm9wCfAOWcfYwDMtBr5TDwcBa/fxww2Tkf64A7M10Ln6HXzyHgd6CV1/pHgG3OumTgkmzOeU6/XY7nNtN+PgROAiec73OpG79bDuc+GNd9EnCH17an6ci0nzN0O8v7AWuAfcB0oI7XZwwwEPjT0TocvabmAQed3z860/kcil5rKcBNXvvqBqx29rMNeNC1e2Ogb8bBmtCb+RLgCSAaOAvYAHRx1v/C6X/aF4G3vX74dc5nSgNfAuOcdfH4bghO2zbzhQhUcC6uvmhr7UZnvqLXRbweaID+wZKA57P5vl3Rm0QToBTwaaY/xKvoDasCEAtMAUY460agN5coZ2oPSBbHiEX/cP8GYpz5Ns66p4EFQBWgsvNHGO71B0hztolyLvijQHln/ZvOd6sBRKA3ieLO/N/O9sWAy5z5yrmdnxzOfRrwf875LgHUc/Zb3NE9Cxjp9ZnMv+9p+3W2f8s5Hy3Qh46Lva6FVEd/hHOeFzjrGgJbgDiv/dbN5rfN6bfL8dxmsa8PgWe85t343bI798G47pPw3RBkpfsq9F5wjqP1cWCe13oDfAOUQR+SjgM/o/eOsuiN/ZZM5/MV5zxdBBwBGjrrdwDtnfflgZau3R8DeTMO5gS0ATZnWjYE+MB5fwfwi/Ne0D9lB2f+Z2Cg1+caok9SkZkvFgpmCPoCCzNpnA/c6nURP+61biDwfTbfd4z3n8X5Exn0RifOBVfXa/35wEbn/dPOxVwvl3N6I/BbNuvWA9285rsAKeafP8CxTOdhN/q0X8xZ1zyLfT6CY4C9lk33+mNle35yOPebc/mOV3t/xyx+31P7BWoB6UCs1/oRwIde18JPXusaAcec9/Wcc3ApEJWDntx+u2zPbTb7+5DTDYEbv1t25z6g173XvgpiCKYBt3vNF0ONYR1n3gAXeq1fAjziNf8yzoMF/xiCUl7rPwP+47zfDNyFjt+4en8Mp8HiOkCciOz3TGiTrKqzfhJwvohUBzqgXQCznXVxaHPXwyb0j18V/5L5OJ5j1fCa3+n1/ijaQsluX1sy7cdDZaAksMTrXHzvLAdtDa0DfhCRDSLyaDbHqIXeOHz5LpucZR7+Nqf3x3u+SyX0aTqr/dYBrsv0G7YDqntt4+v58eB9jjxeNBNEZJuIHAQ+djT5Qhyw1xhzyGtZbr9fjDM2sQ64DzUWux0N3ufLQ26/HWR/bn39Dv7+3fJ6TM9x/X3d+4M6wGte534vapy9te7yen8si3lv7fuMMUcy6fWc72vRVtcmEZkpIuf76TucQTgZgi3oU1M5rynWGNMNwBizD+17vAHoDUwwjllG+xrreO2rNmrJvX9gD0fQP6qHal7vDTmT+TieY23L5XNZsQO9UXvvx8Me9IJs7HUuyhpjSgMYYw4ZY/5tjDkLuBJ4QEQuyeIYW9Amry/fpbazLDf2oN0ndbM53rhMv2EpY8zzPuw3u3OfeflzzrKmxpgyQB/0j57bfkC/XwURifVa5vPvZ4z51BjTDj1vBnghi81y/O38gBu/W6hc95Dz/zMzWeneAtyV6RosYYyZlw+tAOVFpFQmvdsBjDGLjDFXod10X6OtBVcIJ0OwEDgkIo84cQQRItJERM7z2uZT4Gagp/Pew3jgfhFJEJHS6M1iosnaw2QZ0EtEokSklbMvD3+hLY3sbp5TgQYi0ltEIkXkBrT74Nt8fN/PgFtFpJGIlASe9KwwxmQA7wGvikgVABGpISJdnPdXiEg9ERF08C7d0Z2Zb4HqInKfiBQXkVgRaeOsGw88LiKVRaQSOjbzcW6iHW1jgFdEJM75nc4XkeLO57uLSBdneYyIdBSRmj6cj9zOvYdYdGDwgIjUQAcevdmV3T6MMVvQPvURjrZmwO348L1FpKGIXOx8z1T0Zn/GOc/tt/MDbvxuIXHdOywDeohISSe24PYc9pWV7reBISLSGEBEyorIdfnQ6c0wEYkWkfbAFcDnzvxNIlLWGHMSHWjO6j/oF8LGEBj1874CHcDbiD7BjEYHcDxMBuoDO40xy72Wj0E9RmY5n01FBxiz4j/oU9E+YBheBsUYcxR4FpjrNC3bZtL4t6Px3+gg6MPAFcaYPfn4vtNQj49f0G6eXzJt8oizfIHTBfITOvYBeg5+Qm+I84G3jDEzsjjGIXRgtTvadP8T6OSsfgZYDKwAVgJLnWW+8KDzmUVo0/sF1GtpCzpYNxT9k25Bb9S5Xse5nXsvhgEtUQP4HeoY4M0I9Ea5X0QezOLzN6J9y9uBr4AnjTE/5aYPHSx8Hr0ud6JPgUOy2Tan366guPG7hdJ1/yrqJbULGIt6hWW3rzN0G2O+cr7XBOfcrwIuz6tOL3ai94rtjpa7jTF/OOv6AinOce4GbirAcXJE/un9sFgsFkugEJGOqCOJLy1aVwmbFoHFYrFYssYaAovFYglzbNeQxWKxhDm2RWCxWCxhTqFLOlepUiUTHx+fr88eOXKEUqVK5b5hgLG68obVlXdCVZvVlTcKomvJkiV7jDGVs1zpduiyv6fExESTX2bMmJHvz7qJ1ZU3rK68E6rarK68URBdwGJjU0xYLBaLJSusIbBYLJYwxxoCi8ViCXOsIbBYLJYwxxoCi8ViCXNcMwQiMkZEdovIqmzWi4i8LiLrRGSFiLR0S4vFYrFYssfNFsGHaNm47LgczXJZH+gP/M9FLRaLxWLJBtcMgTFmFpqKNjuuAj5yXFwXAOWc6mCuMGfOHN577z1POTiLxWKxOLiaa0hE4oFvjTFNslj3LVpbdI4z/zNa23NxFtv2R1sNVK1aNXHChAl51vL555/z1ltvMXnyZGJjY3P/QAA5fPgwpUv7q8CU/7C68kao6oLQ1WZ15Y2C6OrUqdMSY0yrLFdmF2nmjwkt0LEqm3XfAu285n8GWuW2z/xGFo8bN84AJjk5OV+fd5OiGMXoJlZX3glVbVZX3iiKkcXbOL22aE3yV6PUJypX1hQbf/31l1uHsFgslkJJMA3BZOBmx3uoLXDAGLPDrYN5DMHu3bvdOoTFYrEUSlzLPioi44GOQCUR2YoWkY4CMMa8jRas7obWFT0K3OaWFrAtAovFYskO1wyBMebGXNYbYJBbx8+MNQQWi8WSNWETWRwTE0OJEiWsIbBYLJZMhI0hAChXrpw1BBaLxZKJsDIEZcuWtYbAYrFYMhFWhsC2CCwWi+VMrCGwWCyWMCcsDYGx+YYsFovlFGFnCE6cOMGhQ4eCLcVisVhChrAyBGXLlgVsLIHFYrF4E1aGoFy5coA1BBaLxeKNNQQWi8US5lhDYLFYLGFOWBkCO0ZgsVgsZxJWhqBEiRI235DFYrFkIqwMAWgWUmsILBaL5R+sIbBYLJYwxxoCi8ViCXOsIbBYLJYwxxoCi8ViCXPCzhBUqVKFo0ePcvTo0WBLsVgslpAgLA0BwK5du4KsxGKxWEKDsDME1apVA2Dnzp1BVmKxWCyhgTUEFovFEuZYQ2CxWCxhTmSwBQSaypUrIyLWEFgslpAmPR3WrIHFi2HJEli6FLp1q0jHjv4/VtgZgsjISKpUqWINgcViCRmMgfXrYdEinRYv1hv/kSO6vlQpOPdcEHHn+GFnCEC7h6whsFgswcAY2Lbtnxu+53XfPl0fEwMtWkC/fnDeedCqFTRoABERkJT0tyuarCGwWCwWFzl2TG/08+fr9OuvsGOHrouIgKZNoWdPvemfdx40bgxRUYHVGLaGYM2aNcGWYbFYihjGwJYtMG/ePzf+336DtDRdX68eXHIJtG6tN/3mzaFEieBqhjA2BDt37sQYg7jV6WaxWIo8x4/rQK7npj9/PmzfrutKltSb/YMPwgUXQNu2ULlycPVmR9gaghMnTrB//37Kly8fbDkWi6WQcOSI3uxnzdJpwQI1BgAJCdCxI5x/vt74mzYNfBdPfglbQwCwY8cOawgsFku27NsHc+bA7Nl641+yRLt5ihWDli1h0CBo105v/s5tpVDiqiEQka7Aa0AEMNoY83ym9bWBsUA5Z5tHjTFT3dQEpweVNWrUyO3DWSyWQsKuXXrDHz++HoMHw8qV2u8fHa39+g8/DB066I2/TJlgq/UfrhkCEYkA3gQuA7YCi0RksjFmtddmjwOfGWP+JyKNgKlAvFuaPNjoYovFAnDwIMycCT//rNOqVbo8JqY67drBsGF642/dOjQGdd3CzRZBa2CdMWYDgIhMAK4CvA2BATx2tSyw3UU9p7CGwGIJT44f1z5+z41/4UKN4I2J0S6em27Sfv7Dh+dw6aUXBVtuwBBjTN4/JHKfMWZkLtv0BLoaY+5w5vsCbYwx93htUx34ASgPlAIuNcYsyWJf/YH+AFWrVk2cMGFCnjUDHD58mNKlS2OMoUuXLvTo0YO77747X/vyJx5doYbVlTdCVReErja3daWnw7p1pVm6tDxLl5Zn5cqyHD8eQbFihoYND9Gy5T5attxHkyYHiY7OCJiu/FIQXZ06dVpijGmV5UpjTJ4nYLMP2/RExwU8832BUZm2eQD4t/P+fLS1UCyn/SYmJpr8MmPGjFPv69SpY/r27ZvvffkTb12hhNWVN0JVlzGhq80NXTt2GPPhh8b06mVMhQrGaC+/MY0bG3PvvcZ8840x+/cHXpc/KIguYLHJ5r6a364hX5zvtwG1vOZrOsu8uR3oCmCMmS8iMUAlYHc+dfmMjS62WIoGJ09qd8/33+v022+6vFo16N4dLrsMLr4YqlcPrs5QJr+GwJf+pEVAfRFJQA1AL6B3pm02A5cAH4rIOUAMEJCCwtWrV2f9+vWBOJTFYvEzW7bA9OkwbRr89JMO+kZGqv/+iBHQtSs0a6ZunpbcydYQiMgh9IYvnH7jFyDX8XNjTJqI3ANMR11DxxhjfheRp9EmymTg38B7InK/c4xbnSaM61SrVo25c+cG4lBnsmOHXsVz5kByMq03bdLIkzJloGpVaNJEnZQvuUTnLZYwJyND8/VMmQKTJ8OKFbq8Vi244Qa4/HJ96i9bNrg6CyvZGgJjTGxBd240JmBqpmVPeL1fDVxY0OPkh2rVqrFnzx5OnjxJVKDC/xYtgpdegi+/1KiU8uWhWTMOn3UWJWvW1MearVth1Kh/whUTE9WV4cYbC3fEisWSR44e1af9KVPg229h5059wm/XDl58Ebp1g3POcS81cziRU4tgCTAHmAYkGWNSA6YqAFSrVg1jDLt376ZGjRruHuzAARgyBN5+G8qVg8GD4eab9cm/WDFWJyVRxbvaRHq6dnT+8IMajQce0EiW66+H++7TBCYWSxFk+3a96U+ZokYgNVUbyl27wpVX6pN/hQrBVln0yGmMoA3QDh3MHSYif6PdPNOMMWsDIc5N4uLiANi2bZu7huCPP/TRZdMmNQBPPw2xuTS2IiI0CXmrVjB0qJYpeu89GD0aPv1UR8CefVaTmVgshZz162HSJJ0WLtRlCQnQv79e6h06aGSvxT1y6hpKA5KcCRGJQ43CMyJSD1hgjBkYAI2u4Ln5b9uW2ZHJj8yeDVddpf3/c+ZoXHp+OOcceOUVeOop7Tb67381f23v3mpYzjrLr7ItFjcxBlav1hv/2LGt2LBBl7dqpc83V16pOfltl0/g8NlryBizHRgDjBGRYqjff6GlZs2aAGzdutWdAyxZou3ZWrXUtSEhoeD7LFNGWwh3363G4PXX4bPP4KGH4LHHNO+txRKCGKOlFydN0t7O5GS90TdpksYrr0CPHlCnTrBVhi85jRFMIQc3UWPMla4oChCVKlUiKirKnRbB5s1wxRVQqRLMmOF/B+YKFeD557Wr6dFH4bnnYPx4eOstNT4WSwhgjD4PTZigBiAlRXs9O3bUS/fqqyE5eRkd3ajGbskTObUIXgqYiiBQrFgxatSo4f8WwfHj2rb1uDy4GcVSvTqMHQu33QYDBuhI2nXXwRtvWLdTS9D4/Xe9+U+YAOvWac9o587wxBP616hY8Z9tk5ODp9PyDzmNEcz0vBeRaKCBM5tsjDnptrBAULNmTf+3CJ58EpYvV7eHxo39u+/s6NhRj/niizB8uLZC3n4brr02MMe3hD0bNvxz81+5Ut08L75YG6w9eqintCV0yTXuTkQ6An+iKaXfAtaKSAeXdQUEv7cIFizQm/Htt2vXUCCJjtZxgqVLtbO1Z0+NP9i3L7A6LGHD9u0wcqSWYKxbVy+/2FhtkG7bBj/+qH8FawRCH18CsF8GOhtjLjLGdAC6AK+6KysweFoEfglmPnkSbr0VatZUD59g0aiRJl4ZNkwHkps00Xy7FosfOHoUPvkEunRRP4j774cTJ+CFF3QMYO5cuOceG/tY2PDFEEQZY0715DkxBIWkEmfO1KhRg2PHjrHPH0/N772nHZ6jRgW/dFFUlHbILligWi67TB/X0tKCq8tSKMnI0OItt9+uN/g+ffRSf+wxDZNZulTjHa3XT+HFF/fRxSIyGvjYmb8JWOyepMDh7UJaoSDhiocO/VPKKNBdQjmRmKgJWu69Vz2LkpI0IM3+Yy0+sG4djBsHH32kT/ulS6svwi23QPv2NqFbUcKXn3IAWifgXmda7Swr9PgtqOyVV2D3bvXtD7UomFKl4P331QCsXAktWqgjt8WSBQcPauO2XTuoX199D+rXh48/1nq+Y8bARRdZI1DUyLVFYIw5DrziTEUKvwSV7d+vieSuvRbatPGTMhe48UYtvNqrl2odOBBefllr9FnCGmM0tcO778LEiXDkiAazP/+8dgO5nYrLEnxyNQQicgUwHKjjbC+AMcYEuSO84FSvXh0RKViLYPRoOHwYHn/cf8Lcom5dHc0bOlSNwPz58MUXNkVFmLJvnz7pv/eeNhZLldLnhTvv1LyGoda4tbiHLw28kcAtQEVjTBljTGxRMAIAUVFRVK1aNf8tgrQ09ZXr2FG7XAoD0dHagpkyBTZu1HGEKVOCrcoSIIyBWbOgb1+Ii9Pho+LFtTWwY4cahdatrREIN3wxBFuAVYEqGBNoatSokf8WwVdfaTqJ++7zr6hAcMUVGv+fkKDhnkOHWq+iIsyBA5G8/LJ2+Vx0kRZ36ddPs50vWqStgNyS4lqKLr54DT0MTBWRmcBxz0JjTJEYM6hZs2b+S1a++qp2t4SSp1BeOOssmDcP/u//tL7fggWas8impygyLF2qHs2ffHI+J05oKccPPlDvn1Klgq3OEir40iJ4FjiK1hOO9ZqKBPmOLl61SvvYBw3STFqFlZgY7Q8YM0a/T8uWOo5gKbScOKH2/MILtefvs8+ga9edrFypP+2tt1ojYDkdX1oEccaYJq4rCRI1a9Zk//79HDlyhFJ5+XeMG6cG4Kab3BMXSG67Dc49V1NTdOyorrCFZdzDAmjKh3fe0WnXLnX7HDlS/f6XLfuTJk2s+48la3xpEUwVkc6uKwkSHhfSPI0TZGRonH3XrlClikvKgkCLFhqA9q9/wQMP0GjYMHUst4QsxmjNoxtu0DjB4cO1wMu0aRr1O3iwVke1WHLC14Cy70XkmIgcFJFDIlJk7g61a9cGYPPmzb5/KClJs2r16eOOqGBSrpwOgr/wApVnz1Y/wlWrgq3KkomTJ/VZJDFRo3x/+EE9gP78U2v+du1qg74svpPrpeK4ixYzxpQoau6jAHWcdAspKSm+f+jjj9XF4qqr3BEVbETg4YdZ9vLLcOCABsp9/HHun7O4zv792muXkKDPIceOaVfQ1q0aGlK3brAVWgojYf/MULNmTSIiIti0aZNvHzh2TIOwevaEEiXcFRdkDrRoof6FiYnqeD5woBbesQSclBT1Uq5VCx55BBo2hO++0yIw/fvbwV9LwQh7QxAZGUmNGjV8bxH88osmmbv+eld1hQzVq2sa6wcfhP/9T/shfDWalgKzaJH2/9etC2++qeUdly7Vn6RbN9v9Y/EP9jJCu4d8bhFMmaKPX+FUZzUqSgvueKqOt2ypo5EWVzBGq5xecolG+U6frnZ440Z1Vjv33GArtBQ1fDIEIhIhInEiUtszuS0skMTHx/vWIjBGR+I6dw7PZG3XXKNeRTVrqmfRE09AenqwVRUZMjJ0nL5NGy0hsWaNZgPZskULvzgObhaL3/GlVOX/AbuAH4HvnOlbl3UFlDp16rBt2zZOnsylFPNvv6m3UPfugREWitSvr4FnN9+svoqXXw5//RVsVYWakydh7FgtJtejB+zdqwPAGzfCv/9tUz9Y3MeXFsFgoKExprExpqkzNXNbWCCJj48nIyMj9wjjKVPUo+Zf/wqMsFClZEnNU/Duu5rBrGVLTU9hyROpqdrvX6+eRvtGRWlE8B9/6ABw8eLBVmgJF3xNOnfAbSHBxONCmus4wZQp2m4vSkFk+UVEM5XNmweRkVqd7Y03tPvMkiOpqZr/p25dre9bs6b2OC5bpuUiIn2J97dY/Igvl9wGIElEvqMIJp0DbRFALrEEO3dqts5nnw2IpkJDy5bqxnLzzRrRNHeu1mgoXTrYykKO1FQ9NSNGaDqI9u01PKNjR5v22S1OnjzJ1q1bSU1NzdPnypYty5o1a1xSlX980RUTE0PNmjWJivK9tLwvhmCzM0U7k8+ISFfgNSACGG2MeT6Lba4HngIMsNwY0zsvx/AHtWrVAnJpEcyYoa+di2y2jfxTvjx8842OaD7+OCxfDpMmQaNGwVYWElgDEDy2bt1KbGws8fHxSB5O9qFDh4gNwcGZ3HQZY/j777/ZunUrCQkJPu/Xl1KVwwBEpKQx5qivOxaRCOBN4DJgK7BIRCYbY1Z7bVMfGAJcaIzZJyJB6XMpXrw4cXFxObcIkpKgTBmbiC07ihWDIUO066xXL/V7fO89LXkVphw/rqfAGoDgkZqammcjUJgRESpWrMhfeXTg8MVr6HwRWQ384cw3F5G3fNh3a2CdMWaDMeYEMAHInJPhTuBNY8w+AGPM7jyp9yO5xhIkJek/2Xbg5szFF6t3VYsW0Lu3jnoe9fn5oUiQng4ffggNGmiph7p1NQ5x5kzo1MkagUATLkbAQ36+r+RWeExEfgV6ApONMec6y1bllppaRHoCXY0xdzjzfYE2xph7vLb5GlgLXIh2Hz1ljPk+i331B/oDVK1aNXHChAm+f0MvDh8+TOls+q6HDx/OmjVr+PTTT89YF71nDxdcdx3rBgxgqwsRxTnpCiYF0SVpaSSMGUPt8eM5Eh/P7088wdE8NFXd0uUmhw4dZunSeD74IIFNm0rRsOFB7rhjI4mJ+4J+8w/Vc+a2rrJly1KvXr08fy49PZ2IEKwz4quudevWceDA6T4+nTp1WmKMaZXlB4wxOU7Ar87rb17LlvvwuZ7ouIBnvi8wKtM23wJfAVFAAuqhVC6n/SYmJpr8MmPGjGzXPfrooyYyMtKkpaWdufKTT4wBYxYvzvex86srmPhF1/ffG1O5sjElShjz3nvGZGSEhi4/8+OPxjRseMCAMeecY8ykSX75qn4jFM+ZMe7rWr16db4+d/DgQT8r8Q++6srqewOLTTb3VZ/cR0XkAsCISJSIPAj4Mpy+DajlNV/TWebNVrSlcdIYsxFtHdT3Yd9+Jz4+nrS0NLZv337myqQkKFvWjg/khy5ddPD4ggvU3bR37yJV4+DXXzUVxGWXwf790XzwAaxcqYFhwW4FWIJPSkoKJUqUoIXXvePrr79GRPjjjz9OLVu/fj0tWrQ4rXX022+/cfvttwOQlJTEvHnzTq0bNWoUY8aM8ZtOXwzB3cAgoAZ6I2/hzOfGIqC+iCSISDTQC5icaZuvgY4AIlIJaIC6qwYcjwvpxo0bz1w5Y4b6yYdgU7FQUL26Jsx55hmtm9iypbriFmLWr9e6v23b6o3/tdfgo49+5dZb7WViOZ26deuybNmyU/Pjx4+nXbt2jB8/PtttAJ577jnuvfde4ExD0K9fP9544w2/afTFa2gPkOd6jMaYNBG5B5iO9v+PMcb8LiJPo02Uyc66zs5gdDrwkDHm77weyx94+hHXr19Phw4d/lmxbRusWwcDBgRDVtEhIgIeewwuukg9ic4/XxPrDx5cqB6d9+5VezZqFERHw1NPaRqI0qUhKckG04Uy99133xk32+zwtS++RYsWjBw50mcNhw8fZs6cOcyYMYPu3bszbNiwLLc7dOgQK1asoHnz5qSkpPD2228TERHBRx99xJtvvkn79u2Jj49n4cKFtG7d2ufjZ4erLjDGmKnA1EzLnvB6b4AHnCmo1KlTh8jISNatW3f6Ck/qhHbtAi+qKNKunYbQ3nYb3H+/utN88E6cWt0AACAASURBVAFUrBhsZTly/Di89ZamVzpwAPr1g6ef1saOxeIr33zzDV27dqVBgwZUrFiRJUuWkJiYeMZ2ixcvpkkT9ceJj4/n7rvvpnTp0tx1112n4ghatWrF7NmzQ98QFCYiIyOJj48/0xAsXKiPfs2bB0dYUaRiRQ1Ae/11eOghHXsZNy4kU3sbo7FxjzwCGzbokMeLL0LTpsFWZskreXlydyugbPz48QwePBiAXr16MX78+CwNwY4dO6hcuXKO+6pSpcpp4wwFwRoCL+rVq5e1IWjRwmYA8zci2i3Urp12FV18sSbdHz48ZM71woXaaJk3T2/833+vhsBiyQ979+7ll19+YeXKlYgI6enpiAgvvvjiGb7/JUqUyDUtRmpqKiX8VCXRl4CywSJSRpT3RWSpiBTJPAseQ2A8sRXp6Zp/3w9NL0s2JCZqANqdd+qjdtu2sHp17p9zkV27tOeqTRttBYwerRKtEbAUhC+++IK+ffuyadMmUlJS2LJlCwkJCcyePfuMbc8555zTHkpjY2M5dOjQadusXbv2VPdRQfHFa6ifMeYg0Bkoj8YDnJEzqChQr149Dh48yJ49e3TBmjVw+LA1BG5TqpQm4P/mG63Cnpioo7EBzmR68iS88opGBH/yiXYHrV0Lt99uPYEsBWf8+PFcc801py279tprT/Me8nD22Wdz4MCBUzf/7t2789VXX3HhhReeMhxz587lsssu84s2X7qGPG2WbsA4x/On8Lh55AGP59C6deu0f27hQl3Rpk0QVYURV16pvpi33665Gb77DsaMCciI7I8/avLUP/7QWjsjR6pBsFj8xQxP4kovPO6hWdGvXz8mTpzIHXfcQYMGDVixYsWpsYvffvuNxo0bU9FPTha+tAiWiMgPqCGYLiKxQIZfjh5ieBsCQKOFypXTyiGWwFCtmibnf/NNDeRr1kxbCi6xYYNW4OzcGdLS9NBTp1ojYPEPERERHDhw4LSAsqzwBJRVrVr11LIBAwZQPJvxsj179jB8+HC/6fSlRXA7GkS2wRhzVEQqArf5TUEIER8fT7Fixf4xBAsXwnnnaWZNS+AQgYEDNUNbnz5w9dVawmvkSI3w9gOpqfD88zpFRmqG0PvvD5lxaksRoVatWmzZsiXX7bIKKIuJiaFv375Zbu+vLiEPud7hjDEZaM3iRiLSAWgMlPOrihChePHi1K5dmz///FMzZq5cabuFgsk552h95MceU/fSJk00QrmA/PSTegENG6atgeRkePRRawQs4UuuLQIReQG4AfBE/4IWkZnloq6gUb9+fW0RLFumXkPnnRdsSeFNdLSG8l51FdxyC3TtSoN//UvTVJQpk6dd7dypUcCffgr16+u4wKWXuqTbYilE+NLncTVavL6bMaa7M13ptrBgcSqWYMUKXWATzYUG552nJTEffpjq06bpI/3PP/v00YwMePttOPts+OILePJJ/XmtEbBYFF8MwQY0TXRYUK9ePfbt20fqwoXaH12rVu4fsgSGmBh44QV+e/11fX/ppTBokLr4ZoMn8emAAeqVumKF5geKiQmcbIsl1PHFEBwFlonIOyLyumdyW1iw8HgOnVyyRD1WiqanbKHmYOPG2nV3//3wv//p75TJNe/oUQ1UTkxUz6Bx43RsoGHDIIm2hCWZ01BHRETQokULmjdvTsuWLU9lFC0MaagnA8OBecASr6lI0sDxGyy+dq3eYCyhSYkSGv01c6Z6dV18sUYn79vHzJn60738siaHS05W5yNr0y3BwNsjqESJEixbtozly5czYsQIhgwZcsY2HkItDfVYp56Ax7M62Rhz0m8KQoy6detSNyKC6NRUawgKA+3ba3/PsGEcfOldHvmkHW8fu4W6dQ0zZkgo5rGzBIv77tOWpA+USE/3LZy8RQt1a84HBw8epHz58lmuC7k01CLSERgLpKBRxrVE5BZjTJH0GoqKiqJL9eqa6sAagsJByZJM6/gCd40dzrZdETzAyww/+1dK1n8VradksYQGx44do0WLFqSmprJjxw5++eWXLLcLxTTULwOdjTHJACLSABgPnJk7tYhwocct0U8JnSzusXevPuiNGweNGkUzb1IabeZlwBNToNF0eOEF6N/fBgVa8vTkfsylNNSeriGA+fPnc/PNN7Nq1aozso8GOg21L/+OKI8RADDGrKWIexE1M4Z1wPGoIv01Cz1ffgmNGsH48fD44+pd2ubCSK1xsHIltGql7kIXXaRJhCyWEOL8889nz549/PXXX2esC7k01MBiERktIh2d6T1gsV+OHqLU2r+fFaARxpaQ49ChSPr0gWuvhRo1NFP4GWUM6tVTN6ExY+D337Wb77HH1J3IYgkB/vjjD9LT07NMHBeKaagHoFHF9zrTamdZ0eToUcrs2sUKYM2aNcFWY8nEjz9Cv37nMWGCpohYsCCH4nEiWlhgzRotfvPcc9qEmDw5oJotFg+eMYIWLVpwww03MHbs2CxrI4dcGmpjzHHgFWcq+qxejWRksAooZg1ByHD0qNYHGDUKatdOY9q04rRq5eOHq1aFsWM1vfXAgZquont3LZUZH++mbIvlNNLT03PfyCEk0lCLyGfO60oRWZF58svRQxHn5r8vLs62CEKEhQvh3HPVCNx3H7z77hLfjYA3HTpoqbEXX4RfftHWwXPPaWV6i8UFikIa6sHO6xV+O1phIDkZIiIo1bQpq4NcMjHcOXlS+/6few7i4vTe3akTJCUVoBxGVJSGHN9wg0YmP/aYthbefNMmH7L4nUKfhtoYs8N5O9AYs8l7Agb6VUUokZwMZ51FgyZNSE5OzlNTzuI/kpPh/PPVEPTpo05AnTr58QC1amkGumnTNMvsZZdpTur16/14EIulcODLYHFWpudyfwsJGf74Axo25JxzzuH48eOkpKQEW1FYYQx88IFmmU5JURfRDz/0Wz2aM+naFVat0mbHjz9qd9GQIZDJQ8NiKcrkNEYwQERWAg0zjQ9sBIrmGEF6Ovz55ylDANjuoQBy4AD07q35gdq00cyhmWp9u0NMjN78166FXr20bFmDBtpllFEkq7JaLKeRU4vgU6A7mnSuu9eUaIzpEwBtgWfzZh04bNiQpk2bArB8+fIgiwoPFizQtC2ff/7Pw3mNQGeHiIvTm/+CBVCnjpbHbNtWq6RZLEWYnMYIDhhjUowxNzrjAsfQymSlRaR2wBQGkmQngPrss4mNjeWss86yhsBl0tO1XnC7djo/Z44+nPuS78s12rSBefPgo48059QFF8BNN8GmTUEUZSmMZE5DvXPnTnr16kXdunVJTEykW7durF27NvTTUItIdxH5E9gIzESTz03zm4JQwmMInKT1zZs3t4bARbZvh86dYehQ6NlTE0O2bRtsVQ7FikHfvtpdNHQoTJqk18XDD8O+fcFWZylEeDyCjDFcc801dOzYkfXr17NkyRJGjBjBrl27Qj8NNfAM0Bb4yRhzroh0Aopm11ByMpQrB06yp+bNm/P1119z5MgRSpUqFWRxRYtp0/Q+e+wYvP++BgCHZL2A0qXh2WfhrrvgiSfgpZdU8OOPa3CarXhfaMhDFmrS00v4PQv1jBkziIqK4u677z61rHk2YfGBTkPti9fQSWPM30AxESlmjJkB5CecJ/RJTtanPueO1Lx5c4wxrFy5MsjCig7p6fCf/0C3bjoGsGSJDg6HpBHwpnZtdV9aulST2T3wAJxzDkyYYAeULT6xatUqEhN9S9qcVRrq+++/n7lz59K+fXvgnzTU/sCXFsF+ESkNzAI+EZHdwBG/HD3USE6GSy45Nevp11u+fDltQ6bPovCye7em/PnlF735jxqlhcYKFS1awPTp8MMP2k10441aCu3FF4OtzJILeakfc+jQMVfSUPtKKKahvgqtW3w/8D2wHvUeyhUR6SoiySKyTkQezWG7a0XEiEjwWhqHDsG2bacVta1Tpw5ly5a14wR+YM4cTRMxb572rLz/fiE0At507qzNmbFjYdcu6NSJZg89BIsWBVuZJURp3LgxS5b4VuU3pNJQi0gE8K0xJsMYk2aMGWuMed3pKsoR57NvosFnjYAbRaRRFtvFouksfs3XN/AXa9fq69lnn1okIjRr1swaggJgjHard+wIJUuqZ2a/fsFW5SciIuDmm7Ul+eKLxK5dC61ba/DDqlXBVmcJMS6++GKOHz/Ou+++e2rZihUrsuzeCak01MaYdCBDRPIT19kaWGeM2WCMOQFMQFsXmRkOvADkbP7cxnPS69c/bbHHcyjD9gPnmf37oUcPrRNz1VVaNyDblNGFmRIl4MEHWfDpp5ob++eftf7BTTf9c11Zwh4R4auvvuKnn36ibt26NG7cmCFDhlCtWrUztg25NNTAYWCliPyI19iAMebeXD5XA/DOtrQVaOO9gYi0BGoZY74TkYey25GI9Af6A1StWpWkpCQfZJ/J4cOHs/1s7Z9+4ixg9tatpP/9T4MnJiaGI0eOMH78eGq4FOGUk65gUhBd69aV5sknG7NrV3EGDtxAz55b+e234Otyk8PGkNShA5HNm1N7wgRqTJpEsQkT2HH55Wy6+WaOV6kSPG2hes5c1lW2bNkznqR9IT09PV+fy8zhw4fJyMg4ta/Y2Fjef//9M7bzPpbnfe/evRk7diy33HIL1atXZ+7cuaSnpxMREcGcOXNo0KAB0dHRWepMTU3N23k1xuQ4AbdkMd3sw+d6AqO95vsCo7zmiwFJQLwznwS0ym2/iYmJJr/MmDEj+5V33mlM5cpnLF60aJEBzOeff57v4xZIVxDJr67x440pUcKYuDhj5szxryZjCtH52r7dmHvuMSYqypjoaGMGDTJm8+bQ0BYiuK1r9erV+frcwYMH/XL8zZs3m5o1a5rmzZvnuN26detM8+bNzVlnnXVq2bFjx8xHH32Upa4ffvjBbNy4Mdv9ZfW9gcUmm/uqL4PF5YyODZyagPI+fG4bUMtrvqazzEMs0ARIEpEUNFZhctAGjDdsgLPOOmNx06ZNiYqKYvHiIl2d0y+kp2vxmBtv1KRxS5fChRcGW1UQqV4d3nhD81fdfDO8846W0BwwwEYpBxC9BwYHTxrqzMFimfEElK33yn6bWxrq+GyKKuXn+/piCG7JYtmtPnxuEVBfRBJEJBroheYtAk6lsKhkjIk3xsQDC4ArjTHBueNu3AgJCWcsLl68OM2bN2fhwoVBEFV42LcP/vUv+O9/4e671UXUq8ZGeFOnDrz3no4X9OunLlP16sGdd+p1Z3GNmJgY/v7776Aag0BijOHvv/8mJiYmT5/LdoxARG4EegMJIuJd5DUW2OuDoDQRuQeYDkQAY4wxv4vI02gTJXQKx6alacK5Xr2yXN26dWvGjRt3qn/Ocjq//66DwZs360Nv//7BVhSi1KkD//ufpqx44QU1Dh98oK2FoUPVOFj8Ss2aNdm6dSt//fVXnj6Xmpqa55tpIPBFV0xMDDVr1szTfnMaLJ4H7AAqAS97LT+Ej2mojTFTgamZlj2RzbYdfdmnK2zdqsYgixYBqCF46623SE5OplGjMzxgw5qvv9ZUEaVKwYwZYd4V5Cu1amk03ZAhGoj2zjsaj3DjjRqk1qxZsBUWGaKiokjI5n+dE0lJSZx77rkuKCoYbunKKfvoJmNMkjHmfGPMTK9pqTEmze9KgomneZ7FGAFwKpeH7R76h4wMePJJdZlv1Ehjq6wRyCM1ami468aNWjbzm2/Uv/Zf/4LZszUIw2IJAL5kH+0hIn+KyAEROSgih0TkYCDEBYwNG/Q1myeHhg0bEhsbaw2Bw5EjcO218PTTmrJ/5swg1A4oSlSrplF3mzZpbc5Fi6BDB7Ws33xjcxlZXMeXweL/ooO4ZY0xZYwxscaYMm4LCygbN2qUaK1aWa4uVqwY5513njUEaBaO9u1h8mR9mB0zRgt8WfxAhQqa1TQlRbuOduyAq6+GJk004d2JE8FWaCmi+GIIdhlj1riuJJhs2KDZJSOzHzJp3bo1y5cvzzX/R1Fm6VLNoPDnnzBlCgweXAiyhhZGSpaEQYP0RH/yCURFaZ7uhASt4vN3rhleLJY84YshWCwiE0XkRqebqIeI9HBdWSDZuDHb8QEPrVu3Ji0tLWzzDn39tbYEIiM1cVy3bsFWFAZERmoR52XLYOpUHYwZOlRbrnffDX7KPGmx+GIIyqDZRzvzT93iK9wUFXA2bMh2fMCDZ8B4fpjVrzVGHVt69ICmTeHXX/XVEkBE4PLLtZDzihVqHD78UOshdOumy+3AsqUA5GoIjDG3ZTEVlfyROvK5e3euLYIaNWqQkJDgt0IQhYGTJ4U77lCPxuuuU/fQLPJjWQJJ06YwerQGbQwbpv11nTv/s/zYsWArtBRCfPEaaiAiP4vIKme+mYg87r60AOFxHfXB17hDhw7Mnj07LKIU9+6FRx5pxpgxOn45fnwhrx9Q1KhSRUtnbtqkrYPISI1Url1bf7AtW3LdhcXiwZeuofeAIcBJAGPMCjRdRNEgD4agffv2/PXXXyR7itwXUTZuhAsugFWryjJunHo0FvPlSrEEnuLF4ZZb4LfftMl2wQXw3HMQH68eR9OnW/dTS6748vcuaYzJ7DdZdALKNm/W1zp1ct20Q4cOAMyaNctNRUFl6VI4/3ztLXvppeX06RNsRRafENHqP998o2Nejzyio/pdu0KDBtSaONF6G1myxRdDsEdE6gIGQER6oqknigZbtqh7ng+54uvVq0e1atWKrCH44Qe46CJ9yJw7F5o1OxBsSZb8EB+vrYItW+DTTyEujrpvv61RfzffrGXiwqB70+I7vhiCQcA7wNkisg24DxjgqqpAsnWr/kF86PsQEdq3b18kB4w/+kgzG9StC/Pnq0OKpZBTvLjmL5o1i0VjxsDtt6sf8Pnna57wd9+Fg0UrSYAlf/jiNbTBGHMpUBk42xjTzhiT4rqyQLFlS7YRxVnRoUMHNm/ezKYikk/eGH14vOUWbQ3MmgVxccFWZfE3RxIS4M03NTT8f//TcYO77tKaCbfdZnMbhTm+eA09JyLljDFHjDGHRKS8iDwTCHEBIR+GAIrGOEF6ugawPvaYltedOhXKFK3kIZbMxMZqMNqyZdr0u+kmmDRJcxs1bAjPPw/btwdbpSXA+NI1dLkxZr9nxhizDygacaUZGfqElIfc3U2aNKFSpUr8+OOPLgpzn2PHoGdPfTh85BHtGoqODrYqS8AQgbZttXtoxw51Qa1eXVNj164N3bvDV1/ByZPBVmoJAL4YgggRKe6ZEZESQPEcti88/PWXJvLKQ4ugWLFiXHbZZfzwww+FNp5g71649FJ1MHn9dX0ItO6hYUypUto3OHMmrF0LDz2kecV79NCHpAcfhJUrg63S4iK+/P0/AX4WkdtF5HbgR2Csu7IChCfoJg+GAKBz587s2rWLFSt8qs8TUuzYoWMBixfD55/D//1fsBVZQor69TWx3ebN8O23mgr7tde0WE7z5ppvZOvWYKu0+BlfBotfAJ4FznGm4caY/7otLCB4Lug8lnXr3LkzAD/88IO/FbnKhg3Qrp0GjE2dqjUFLJYsiYxUN7Ivv9Qxg1GjNCvqww9r19Ell2h3kvU6KhL41CFgjJlmjHnQmaa7LSpg5LNFEBcXR5MmTZg+vfCcipUr1Qjs36+F5S+5JNiKLIWGypXVq2D+fO068qS2uO02qFpVa31/+60dTyjEhHeFsi1bdIS0cuU8f7Rz587Mnj2bo0ePuiDMvyxYoN1BIuol6CRStVjyTv368NRTWith/nyNTfjpJx1cjouDgQMhKUld0iyFhvCuULZ1q3YL5aO6SpcuXThx4gQzZ850QZj/+PFHffqvWFGjhRs1CrYiS5HA43XkqaQ2ZYpeaGPHQqdOGqR5zz0amGKNQsgT3hXK8hhD4E379u0pUaIE3333nZ9F+Y8vvtBu3nr1tCUQHx9sRZYiSVQUXHEFTJigSao++0yrGI0Zo03RWrXg3nv1ScQmwAtJwrtCmadFkA9KlChBly5d+Prrr0PSjXT0aLjhBu0GmjnT1hGwBIhSpbR4xeefq1GYMEFTWrz3ng5S1a4N999vjUKIEb4VyjzBZPlsEQBcffXVbNu2jcWLF/tRWMF5/XVNTd+liyaSK1cu2IosYUnp0vo0MmmSGoVPPoFWreCtt9QoxMXR4KWXYNo0OH482GrDmuyrtTsYY24LhJCAs2uXejkUwBB0796diIgIvv76a8477zw/iss/r74KDzwA11yjD2M2WtgSEsTGaonN3r3V5XTqVPjqK6pMmQLffafru3XTC/fyy22ukwDji9dQTRH5SkR2O9MkEclff0ookc8YAm8qVKjARRddxFdffeUnUQXjpZfUCPTsCRMnWiNgCVHKlFGX04kTmfv112oIevXSwjq9eqkXX7du2p20c2ew1YYFvnQNfQBMBuKcaYqzrHCTzxiCzFx99dWsWbMm6FXLXnhBMwNcf72moI+KCqoci8UnTHS03vTffVcD12bP1nD35GTo31/zH7VuDU8/rWkv7LiCK/hiCCobYz4wxqQ504doSurCjSfDYo0aBdrN1VdfDRDUVsFzz8Gjj2rq+U8+sUbAUkiJiNCxg5degnXrYMUKrZMaEaGxC61a6f/1jjs0Id6hQ8FWXGTwxRD8LSJ9RCTCmfoAhb/m3Y4deoFVqlSg3dSqVYsLLriATz75JCjeQ888808a6Y8+0swAFkuhRwSaNoXHH9fAtV279ALv0EH9onv00OCYyy7TXEjr1gVbcaHGF0PQD7ge2ImWqOwJFP4B5B07NDzeD2k3+/Tpw6pVqwKehG7YMPjPf6BvX43jsUbAUmSpXFkv9IkTNWtwUhLcd596/t13n0Y8160LAwZoa2H//lx3afkHX5LObTLGXGmMqWyMqWKMudoYs9mXnYtIVxFJFpF1IvJoFusfEJHVIrJCRH4WkdwryPuLnTu1/9EPXH/99URGRvLxxx/7ZX+5YQw8+aS2lm+9FT74QBs3FktYEBWlgWr//S+sXq3ZFEeNgiZN4OOPtbVQqZJmTh02TFsUaWnBVh3S+OI1NFZEynnNlxeRMT58LgJ4E7gcaATcKCKZExz8BrQyxjQDvkDTWQSGHTv8ZggqVqxIt27d+PTTT0l3OZzeGG0FPP009OsH779vjYAlzElI0KR433yjxTZmzdICO2lpagguuEANQ48e8Pbbajgsp+FLv0izLCqUnevD51oD65yaxyeACcBV3hsYY2YYYzxZ2xYAgXNL3bHDr+G2ffr0Yfv27SQlJfltn5kxBoYOhWef1YCx996zBWUsltOIitL0FsOHw6+/wp49mvLiuuvU62jAAO1CqldPazaPH6/3gjBHchvgFJHlQEfHACAiFYCZxpimuXyuJ9DVGHOHM98XaGOMuSeb7UcBO40xZ9RDFpH+QH+AqlWrJk6YMCHXL5YVhw8fpnTp0pCezkWdO7PppptI6dcvX/vKzIkTJ+jRowdt27bl8ccfz5+uHDAG3nnnLCZOrM2VV25j8OA/XTcCvugKBlZX3glVbQHVZQwltmyhwuLFlF+8mHIrVhB55AgAR2vVYt+557K/eXP2t2jBvujoIne+OnXqtMQY0yrLlcaYHCfgZuAPYLgz/QH09eFzPYHRXvN9gVHZbNsHbREUz22/iYmJJr/MmDFD32zfbgwY89Zb+d5XVvzf//2fiY6ONrt3786frmzIyDDmgQdU8sCBOh8IctMVLKyuvBOq2oKq6+RJYxYtMua//zWmWzdjYmP1TwbmcJ06+mf77DNj8vh/dpOCnC9gscnmvurLYPFHQA9glzP1MMaM88EAbQO8o7VqOstOQ0QuBR5DU10HJuGIJ1rRz5nYBgwYwIkTJxgzJtchFJ8xRnN0vfKKxtmMGpWvrNkWiyUzkZEam/DQQxrdvHevFu94/nmOV6mirnjXXw9Vqqgr66BB2pXkCUYtQvjkcGiMWQ2szuO+FwH1RSQBNQC9gN7eG4jIucA7aBfS7jzuP/94+gT9NFjs4ZxzzqFjx468/fbbPPjgg0QUcBTXGBg8GN54Q19ffdUaAYvFNSIjoU0baNOGFW3a0PHCC7W4d1KSpr8YO1YT5oFmJGjXTj2T2rVTj6VC7LXhWi+zMSYNuAeYDqwBPjPG/C4iT4vIlc5mLwKlgc9FZJmITHZLz2l4DIELuZkHDBhASkpKgctYGqMtgDfe0PxB1ghYLAEmKkpTaA8Zoml89+/XAefXXtPlM2dq8Z0WLaBCBejaVSM8Z8yAQlC50BtXQ5CMMVOBqZmWPeH1/lI3j58tLnUNgaacqFatGq+++irdunXL1z4yMrQV+vbb2mp94QVrBCyWoBMZCS1b6nTvvfq0tmkTzJmj9RXmzNF6zsbotueee6qFQZs26qkUon/k8IxF3bEDypeHmBi/7zo6Opr77ruPRx99lCVLlpCYmJinz2dkwN13q2voo49qHqEQvXYslvBGRMv+xcdDnz66bN8+DWCbOxfmzdNoz1GjdF358ppAz2MYWrcucIobfxGehmDnTldLdg0YMIARI0YwYsQIvvjiC58/l5GhCRfff1/zBw0fbo2AxVKoKF9es6l6egPS0zX6+ddf/5meeeafLKp1657eamje3JUH1NwIT0Pgx6jirChTpgyDBg1ixIgRJCcn07Bhw1w/k56uSRU//FBbl089ZY2AxVLoiYhQj6OmTfUPDnD4sA5CewxDUpLmjgftUmrcWL2ZEhP1tWlT141DeMalumwIAAYPHkzx4sUZMWJErtump8Ntt6kReOopjYq3RsBiKaKULg0dO8Ijj8CXX2rivK1btaTnww9rb8XXX8PAgdp9FBur4w133kmZVatckRR+LQJjXO8aAqhSpQqDBg3i1Vdf5cEHH6RJkyZZbpeeDs8/fw4//aT5g/7zH1dlWSyWUKRGDc2F1KOHzhsDmzerl9Lixfr6lkMCbgAAEBFJREFU5ZeUqFLFlcOHX4vg4EE4dsz1FgHAkCFDiI2NZejQoVmuT0vTzLo//VSVZ5+1RsBisTiIQJ06ahieew6mT4c9e9h98cWuHC78DIFLwWRZUbFiRR599FGmTJnC7NmzT1uXlqaOBuPHQ//+68nGVlgsFosignEpaC38DIGLMQRZce+99xIXF8cDDzxwKkX1yZNaVnLiRHjxRbjxxqIXsm6xWAoP4WcIAtgiAChZsiQvv/wyixcv5p133uHECejVS6vtvfIKPPhgQGRYLBZLtoSfIdi1S1+rVg3YIW+44QYuvfRSHn30Ca68MpUvv4SRIzWZnMVisQSb8DMEu3err2758gE7pIgwcuT/OHLkQ6ZPj+H11w2DBwfs8BaLxZIj4WcI/vpLw7oD6Kh//Dg88kg9MjKuAAZSvvwnATu2xWKx5EZ4GoLKlQN2uGPH4JprNN35W29l0K7dSgYNGkRKSkrANFgsFktOhKchcCkoIzNHjsAVV8D332sSuQEDijFunNb06d27NydOnAiIDovFYsmJ8DQEAWgRHDyo6cmTkrSehSfNSHx8PKNHj2b+/Pnce++9ruuwWCyW3Ai/FBO7d7tuCPbtUyOwdKkGjF1//enrr7vuOoYMGcKIESNo0aIFZ599tqt6LBaLJSfCqkUgJ0/CgQOuGoI9e+CSS2DZMo0VyGwEPAwfPpxu3bpxzz33MGfOHNf0WCwWS26ElSGIOnBA37hkCHbtgk6dYM0a+OYbuOqq7LeNiIhg4sSJJCYm8vTTT5OUlOSKJovFYsmN8DIE+/frGxcMwaZN0KEDbNigHkJdu+b+mdKlSzN16lTi4uK44oor+OWXX/yuy2KxWHIjrAxBtMcQ+Nlr6Pff4cILtUUwfTrkJUFgxYoVefnll0lISKBbt25MmTLFr9osFoslN8LKELjRIpg/H9q318pzs2ZBu3Z530fFihVJSkqiadOmXH311bz++usYY/ym0WKxWHLCGoIC8P33cOmlULGi1qpu1iz/+6pYsSIzZsyge/fuDB48mP79+3Ps2DG/6LRYLJacCCtDEH3ggNYQ9UOeoXff1WCxs89WI5CQUHB9pUuX5ssvv2To0KGMHj2a1q1b8/vvvxd8xxaLxZIDYWUIovbv18f3Yvn/2unp8O9/w113QefOMGOGf4ccihUrxrPPPsv333/P7t27SUxM5JlnnrFRyBaLxTXCzxAUoFvo8GHNG/TKK3DvvTB5MpQp40eBXnTp0oUVK1Zw1VVX8Z///IcWLVowbdo0O3ZgsVj8TlgZguj9+/P9+L5+vQ4ET50Kb74Jr72m2azdpGrVqkycOJEpU6Zw4sQJunXrxmWXXWYD0CwWi18JK0MQlc+o4s8/h5YtYfNmjREYONAFcTlwxRVXsHr1akaOHMmKFSto3749F110EZ999pntMrJYLAUmvAxBHruGUlNh0CBNE9GoEfz2G3Tp4qLAHIiOjmbw4MGkpKQwcuRINm3axA033EDNmjV5+OGHWbt2bXCEWSyWQk/4GIKTJ4k6dMhnQ7BqFVxwAbz1ltYVnjUL6tRxWaMPlCxZksGDB7N+/XqmTZtG+/btefXVV2nYsCHNmjXj8ccf59dffyUjIyPYUi0WSyEhfAzBnj36moshSE2FYcO0K2jLFh0QfvFFiIoKgMY8EBERQdeuXZk0aRJbtmzhlVdeoWLFijz//PO0bduWuLg4rr/+el577TWWLFlCWlpasCVbLJYQJXzSUP/1l75mM1hsDEyaBA8/DBs3Qu/eWmA+gMXM8k21atW4//77uf/++9m7dy/Tpk1j6tSpzJkzh88//xyAUqVK0axZM5o2bXrqtX79+lSrVg0JYNlOi8USerhqCESkK/AaEAGMNsY8n2l9ceAjIBH4G7jBGJPiihiPIch0Zz95Ug3A88/D8uXQpAn8+KNGDBdGKlSowE033cRNN90EwJYtW5g7dy7z5s1j+fLlfP7557z77runto+JiaFOnTokJCQQHx9PQkICcXFxVKlShc2bN9OgQQMqVapEdHR0sL6SxWJxGdcMgYhEAG8ClwFbgUUiMtkYs9prs9uBfcaYeiLSC3gBuMEVQV6GYO9emDdPE8RNmKC9RmefDR98AH36uO8WGkhq1apFr1696NWrFwDGGLZv387KlStZv349KSkpbNy4kZSUFBYuXMjevXuz3E+5cuWoUKECZcqUoUyZMsTGxp5675kvXbo0xYsXJyYm5tSr93vv14iIiFNTZGRkjvPFihWzrRaLxUXcvOW1BtYZYzYAiMgE4CrA2xBcBTzlvP8CGCUiYlyImnr/m0o8z1p2t63LwUO6rHhxrRlw881w+eUFCjguNIgINWrUoEaNGlmuP3jwIDt37mT37t388ssvVKlShd27d7N792727dvHoUOHOHjwINu3byc5OZmDBw9y8OBBUlNTXdXtbRgAoqKiThkHETk1ZZ7P7za5fS4rjh07RokSJVw6A1njq4E8evQoJUuWDAkt3ripqyAPD/7W5a8HmWuvvZaOHTv6ZV/euGkIagBbvOa3Am2y28YYkyYiB4CKwB7vjUSkP9AfNMgqP0Vcjkb8TZOKGUj7GCpXOUGjRgc5++xDFC+u3jWzZuV5l37j8OHDIVmYpmXLlpQuXdqnUpppaWmkpqZy4sQJTpw4wcmTJ7N89UwZGRmnpvT0dNLT089YltVrRkYGqampRDrNNmPMadHWnnl/L/O8z+0cRAawOZmX5yW3teX32c0tXQV9lgz0b+krkZGR7twrvC9+f05AT3RcwDPfFxiVaZtVQE2v+fVApZz2m5iYaPLLjBkz8v1ZN7G68obVlXdCVZvVlTcKogtYbLK5r7rZGbINqOU1X9NZluU2IhIJlEUHjS0Wi8USINw0BIuA+iKSICLRQC9gcqZtJgO3OO97Ar84lstisVgsAcK1TjCjff73ANNR99ExxpjfReRptIkyGXgfGCci64C9qLGwWCwWSwBxdTTEGDMVmJpp2RNe71OB69zUYLFYLJacCQOHSYvFYrHkhDUEFovFEuZYQ2CxWCxhjjUEFovFEuZIYfPWFJG/gE35/HglMkUthwhWV96wuvJOqGqzuvJGQXTVMcZkmU+50BmCgiAii40xrYKtIzNWV96wuvJOqGqzuvKGW7ps15DFYrGEOdYQWCwWS5gTbobg3dw3CQpWV96wuvJOqGqzuvKGK7rCaozAYrFYLGcSbi0Ci8VisWTCGgKLxWIJc4qMIRCRriKSLCLrROTRLNYXF5GJzvpfRSTea90QZ3myiHQJsK4HRGS1iKwQkZ9FpM7/t3f+oXYUVxz/fFGrJg0mqdQGa1sjQmpATZSi0fqjtqgRG6WFpigYjdY0UZo/FCwpIoViRGhLodCCFBQkalOV1jbQaFKEhGeMNsmLv/LjJdiE0JSoMaGSaDj9Y86tk9u97+3TO3vTd88Hljc7Z2b3e8+evefuzmMmsx2RtMG39im8S+uaJ+lf2flvz2y3SNrq2y3tfQvr+kWmaYuk9zJbSX/9TtJeSZs72CXpV657k6SZma2Iv2pousm1DEpaK+m8zLbT6zdIWt8tTaPQdoWk/dn1uj+zDRsDhXXdm2na7DE12W1FfCbpDEmr/XvgNUk/qmhTNr46rVjz/7SRprneDkwFPgNsBM5pa7MQ+I2X5wJPevkcb38icKYf57gGdV0JjPPyD1u6fP9gD/01j7YV5bx+MjDkfyd5eVJTutra302a3ryov/zYlwEzgc0d7LOBFYCAi4CXGvDXSJpmtc4FXNvS5Ps7GWE1wMLargCe+7Qx0G1dbW2vJ62RUtRnwBRgppcnAFsq7sei8TVWngi+BmwzsyEzOww8AcxpazMHeNTLy4GrJMnrnzCzQ2a2A9jmx2tEl5mtNrN/++4AaSW30tTxVyeuBlaa2Ttm9i6wErimR7q+Dyzr0rmHxcxeJK2Z0Yk5wGOWGAAmSppCQX+NpMnM1vo5obnYap17JH914tPEZrd1NRJfZrbHzF718gHgDdJ67jlF42usJILTgX9k+7v4X0f+t42ZfQTsBz5Xs29JXTnzSVm/xUmS1ksakHRDlzSNRtd3/DF0uaTWsqPHhL/8FdqZwKqsupS/6tBJe0l/jYb22DLgr5JekfSDHugBuFjSRkkrJE33umPCX5LGkb5Q/5BVF/eZ0ivrGcBLbaai8VV0YZqgPpJuBi4ELs+qv2xmuyVNBVZJGjSz7Q1J+hOwzMwOSbqT9DT1jYbOXYe5wHIzO5LV9dJfxyySriQlgkuz6kvdV58HVkp6038tN8WrpOt1UNJs4Fng7AbPPxLXA2vMLH96KOozSZ8lJZ7FZvZ+t45bh7HyRLAbOCPb/6LXVbaRdDxwCrCvZt+SupD0TWAJ8G0zO9SqN7Pd/ncI+Bvpl0IjusxsX6blEeCCun1L6sqYS9tje0F/1aGT9pL+GhFJ55Ku3xwz29eqz3y1F3iG7r0OrYWZvW9mB738F+AESafSY39lDBdfXfeZpBNISeBxM3u6oknZ+Or2wEcvNtKTzRDpVUFrgGl6W5tFHD1Y/JSXp3P0YPEQ3RssrqNrBmlw7Oy2+knAiV4+FdhKlwbNauqakpVvBAbs48GpHa5vkpcnN6XL200jDdypCX9l5/gKnQc/r+Powbx1pf1VQ9OXSGNes9rqxwMTsvJa4Jpu+qqGti+0rh/pC/Vt912tGCily+2nkMYRxjfhM//cjwG/HKZN0fjq6oXv5UYaVd9C+lJd4nU/Jf3KBjgJ+L3fGOuAqVnfJd7vLeDahnU9D/wT2ODbH71+FjDoN8IgML9hXQ8Cr/n5VwPTsr63uR+3Abc2qcv3HwCWtvUr7a9lwB7gQ9J72PnAAmCB2wX82nUPAheW9lcNTY8A72axtd7rp7qfNvo1XtJNX9XUdlcWXwNkyaoqBprS5W3mkf6BJO9XzGekV3YGbMqu1ewm4yummAiCIOhzxsoYQRAEQfAJiUQQBEHQ50QiCIIg6HMiEQRBEPQ5kQiCIAj6nEgEQVCBpImSFvZaRxA0QSSCIKhmImnG2iAY80QiCIJqlgJn+dzzD7cbJd0saZ3bfyvpOK8/KOlnPpnagKTTKvpens15/3dJExr4PEHQkUgEQVDNfcB2MzvfzO7NDZK+CnwPuMTMzgeOADe5eTxpOo7zgBeBOyqOfQ+wyPt+Hfig0GcIglrE7KNBMHquIk3C93Ja0oKTgb1uOww85+VXgG9V9F8D/FzS48DTZrarrNwgGJ5IBEEwegQ8amY/rrB9aB/P23KEinvMzJZK+jNpPpk1kq42szfLyQ2C4YlXQ0FQzQHSsoFVvAB81+elR9JkZWtNj4Sks8xs0MweAl4mzaYaBD0jEkEQVGBp7v41voD5w22214GfkFar2kRaHnDKKA6/2I+7iTQL5oqROgRBSWL20SAIgj4nngiCIAj6nEgEQRAEfU4kgiAIgj4nEkEQBEGfE4kgCIKgz4lEEARB0OdEIgiCIOhz/gMVm7OK6oXeigAAAABJRU5ErkJggg==\n","text/plain":["<Figure size 432x288 with 1 Axes>"]},"metadata":{"needs_background":"light"}}]}]}