{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "## Solving and plotting 2D dynamical systems" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Example: $\\dot{x} = y$, $\\dot{y} = -x$" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "# importing the needed python libraries\n", "import matplotlib.pyplot as plt\n", "import numpy as np\n", "from typing import Tuple\n", "from scipy.integrate import solve_ivp # used in numerical integration" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "def simulate(n_timesteps: int, state0: list=[1.0, 0.0], h: float=0.02) -> Tuple[np.ndarray, np.ndarray]:\n", "\n", " def diff(t, state):\n", " x, y = state \n", " return y, -x # dx/dt = y, dy/dt = -x\n", "\n", " t = np.arange(0, n_timesteps * h, h)\n", " \n", " sol = solve_ivp(diff, \n", " y0=state0, \n", " t_span=(0.0, n_timesteps*h),\n", " dense_output=True)\n", " \n", " return sol.sol(t).T, t" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Trajectories starting from to different initial states" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "[]" ] }, "execution_count": 3, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAD4CAYAAADhNOGaAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8vihELAAAACXBIWXMAAAsTAAALEwEAmpwYAABBHUlEQVR4nO3dd3hUVfrA8e9JL4T0EAiEhN5r6IKgohQVRVFZxS66yq6uusr+XPuqWHZdO6KLYsUOqBRFadKL9F4SSCEhBNL7nN8fd4IhJKTNzM3MvJ/nmWdunftmZjLvveeeorTWCCGEcF8eZgcghBDCXJIIhBDCzUkiEEIINyeJQAgh3JwkAiGEcHNeZgfQEBERETouLs7sMIQQwqls3rw5U2sdWXW5UyaCuLg4Nm3aZHYYQgjhVJRSSdUtl6IhIYRwc5IIhBDCzUkiEEIINyeJQAgh3JwkAiGEcHM2SQRKqdlKqQyl1M4a1iul1OtKqYNKqe1KqX6V1o1RSu2zrptui3iEEELUna2uCD4Expxn/Vigo/UxFXgHQCnlCbxlXd8NmKyU6majmIQQQtSBTdoRaK1XKqXizrPJBOAjbfR5vU4pFaKUagnEAQe11ocBlFJzrdvutkVcQjREUWk5OYWlZBeWkl9STmFJOUWl5RSWllNQYjyXlFmwWDTlWlNuMR5bjp5i+b4T+Hh64OPlgbenwsfLmPb39iQkwIewAB9CA30IC/QmNMCHiGa+tA71JzYsgMggX5RSZv/5wg05qkFZDHCs0nyydVl1ywdV9wJKqakYVxPExsbaJ0rhkrTWZOWXkHq6iPScIjLzijmRW2w8W6dPFRg//DmFpRSXWRp1vJJyCyXljXuNquIjAukSHUT3Vs3p1zaUvm1C8ffxtOkxhPtyVCKo7jRHn2f5uQu1ngXMAkhISJDRdMRZ8ovLSDyZz5HMfI6cyCf5VCGp2YWknC4k9XQhRaXn/jAH+3sT0cw4K+8Y1Yxgf2+C/b1pXuk5yNcLP29P/H088fc2Hn4+Hvh4euDpofDy8MDDA+NZgdZGIigtt1Barikps1BSZqGgtIzTBaWcyi8hq6DEeM4v5UReMSmnCjh2qpATucU1/n1HMo2/bdHO4zVu07t1MCM7RzGqSxS9YoLx8JCrC1E3jkoEyUCbSvOtgVTAp4blQlTrVH4J+9Jz2Xc8l33puRzKyCPxZD7pOWf/iEY08yUm1J8u0UFc3CWKViH+tArxJ7q5H5FBvoQ388HXy/Zn1EqBn4cnft62e+2ycguJJ/PZnZbLrpRsNiedYlPSqXO225aczbbkbF775cA5667o3YorerXkws6Rdvm7hXNzVCJYAEyz3gMYBGRrrdOUUieAjkqpeCAFuAH4k4NiEk2Y1pqjWQVsT85mR0o2e9Jy2Hs896yz5mB/bzpENWN4x0jiIwKJjwgkLjyQuIgAAnycshutanl5etAhKogOUUFc2btVtduUllvYnHSKZfsyWL73BPvSc89a//22VL7fdvY5VutQfyYPjOW6hDZEBvnaLX7R9ClbjFmslPocGAlEAOnAk4A3gNZ6pjLugL2JUbOoALhNa73Juu844L+AJzBba/1cbcdLSEjQ0umcazmZV8ympFP8fvQ0O1JOsyM5m5yiMgB8PD3oFN2Mzi2a0zm6GZ2jm9MlOogoublaJ6fyS1iy6zjfb09l9cGT5902LjyAqSPaM7FfjE2vakTToJTarLVOOGe5Mw5eL4nAuWmtOZyZz6bELDYlnmJz0ikOZ+YD4OWh6NIyiJ4xIfSMCaZX62A6tQjCx0vaPtrazpRsPt9wlM83HMVynp+BW4fGMXVEO1qF+DsuOGEXkgiEqVJOF7L6YCZrD51k9cFMMqxFPKEB3vRvG0r/tmEkxIXSMyZYzkRNtDkpi3dXHOan3ek1bvPnke25e0Q7QgJ8HBiZsAVJBMKhCkvKWXMok2X7MvjtQCaJJwsAiGjmw5D2EQxtH86AuDDaRwZK8U4TlpFbxJw1iby17FC164P9vXlmQneu7N1KPkcnIIlA2N2xrAJ+3ZvBsn0ZrD10kuIyCwE+ngxpF87QDhEM6xBO5xZB8oPhxDJyinhr2UHmrK12fBMm9W/NI2O6yM3nJkoSgbA5rTUHMvJYtOM4i3amsfe4UVMlLjyAUV2iuKhLFAPjw6S6ogvbduw0T3+/iy1HT5+zbki7cGZc05O24YGOD0xUSxKBsAmtNbvTcli4I41FO49z+EQ+SkFC21Au6x7NxV1bEB8h//juqKCkjJnLD/H6rwfPWdclOog3JvelY4sgEyITFSQRiEY5llXA/K0pzNuaysGMPDw9FIPbhTGmR0su69aCqOZ+ZocomhCLRfPphqM8Pu/cDolHdo7k5Wt7S/GRCSQRiHrLLijl++2pzPs95UxL1oFxYUzo24qxPVoSFii1RkTttNbM35rKA19sPWfd7cPi+ce4Lnh7SvVgR5BEIOpEa82GI1nM3XiMhTvSKC6z0DGqGVf1jWFCn1a0Dg0wO0ThxMotmpkrDvHykn3nrPv8rsEMaR9uQlTuQxKBOK+TecV8vTmZLzYe43BmPkG+XlzZpxU3DIilR0xzqekjbC6/uIxHv9nOD9vTzlo+vmdL/n1db2lPYgeSCES1dqZk8+GaRBZsS6WkzEJC21BuGBjLuJ7RLtVfj2jaNiZmcf27a89q4ezj6cGCvwyjS3Rz8wJzMZIIxBml5RaW7DrOnDWJbEw8RYCPJxP7xXDzkDg6Sa0OYaLScgtPLtjFZ+uPnrX8mQnduXlInDlBuRBJBIKCkjK+3HiM91YdIeV0IbFhAdw8pC2TEtoQ7O9tdnhCnOXn3enc9dHZ/+e3DGnLk1d0l7EWGkgSgRs7lV/CnLWJzFmTyKmCUgbEhXL3iPaM6hKFp/xDiSbueHYRV7+9mrTsojPLhneMYNaUBBmlrZ4kEbihzLxi3l1xiE/WHaWwtJxLukZxz4XtSYgLMzs0IeqtoKSMuz7adFZX2r1bB/P51MFyP6uOJBG4kaz8EmatPMycNYkUl5UzoU8Mfx7ZXsr/hUsot2ieXLCTT9b9cR+hb2wIn905WK4QaiGJwA1kF5Yya+UhPlydSEFpOVf2bsVfL+5I+8hmZocmhM1prXlqwa6zOsAb0i6cj+4YKA3UamDXRKCUGgO8hjHK2Pta6xlV1v8duNE66wV0BSK11llKqUQgFygHyqoLsipJBGcrLivnk3VHeePXA2QXljK+Z0vuv7ij9Osi3ILFonmiyhXCbcPieOLybtL+pQq7JQKllCewHxiNMUj9RmCy1np3DdtfAfxNa32RdT4RSNBaZ9b1mJIIDBaL5ocdaby8ZC/HsgoZ3jGC6WO70L1VsNmhCeFw5RbNfZ9uYfGu42eWzZjYkxsGxpoYVdNSUyKwxfXTQOCg1vqw1roEmAtMOM/2k4HPbXBct7b12GmufmcNf/38dwJ9vPjo9oF8fMcgSQLCbXl6KGZO6c+upy+jbbjRFcr0b3cQN/1H9h7PMTm6ps0WVwTXAmO01nda56cAg7TW06rZNgDjqqGD1jrLuuwIcArQwLta61k1HGcqMBUgNja2f1JS9QNjuLqTecW8vGQfX2w6RkQzXx4d04Wr+8ZINVAhqkg5XciwGb+emW8T5s/Pf7vQrbuusOcVQXW/QDVllyuA1RVJwGqY1rofMBa4Tyk1orodtdaztNYJWuuEyMjIxkXshMotmo/WJjLqleV8vTmZOy+I59eHLuTa/q0lCQhRjZgQfxJnjGfWlP4AHMsqpMvji3lt6QGTI2t6bJEIkoE2leZbA6k1bHsDVYqFtNap1ucM4DuMoiZRyZ60HCa+vZon5u+iR0wwi+4fzmPjuxHkJ62BhajNpd2jOfLCOCb2iwHg1aX7iZv+I8eyCkyOrOmwRdGQF8bN4ouBFIybxX/SWu+qsl0wcARoo7XOty4LBDy01rnW6Z+BZ7TWi893THe5WVxUWs6bvx5k5opDBPt78+SV3bmiV0upCSFEA53KL6Hvsz+fmb+ydyteu6GP2/xP1VQ01OjmeFrrMqXUNGAJRvXR2VrrXUqpe6zrZ1o3vRr4qSIJWLUAvrN+CF7AZ7UlAXexKTGLR77ZzuET+UzsF8Pj47sRKgPBCNEooYE+JM4Yz9wNR5n+7Q4WbEtlwbZUFt0/nK4t3beXU2lQ1sSUlFl4del+Zq44RKtgf56f2JMLO7nfPREh7K2otJzBL/zC6YJSAK7t35pXJvU2OSr7kpbFTmB/ei4PzN3K7rQcbhjQhn9e3o1mvtKHihD29OP2NO77bMuZ+XX/uJjoYNccg9uetYZEI2mt+WD1ES5/4zfSc4qYNaU/M67pJUlACAcY36slu5+57Mz84Bd+4X+/HTExIseTKwKTZReW8sjX21iyK52Lu0Qx45peRAb5mh2WEG7p3RWHeGHRXgBaBvux+tGLXGrsA7kiaIK2J5/m8jdW8cueDP45vivv35IgSUAIE919YXtWPTIKgLTsItr930KOnnT9aqaSCEygtebjtYlc+85ayss1X94zhDuHt3ObKmxCNGVtwgI4/Pw4YkL8ARjx8jK+2Hi0lr2cmyQCBysuK2f6Nzt4fP4uhnUI58e/DqdfbKjZYQkhKvHwUKyefhFPXtENgEe/2cE9H282OSr7kUTgQJl5xdz43nq+2HSMaaM68L9bBkjbACGasNuGxbPo/uEALN51nLjpP1JWbjE5KtuTROAgu1KzufKN39iZms0bk/vy8GWdXeomlBCuqmvL5ux8+o9aRR0eW0RmXrGJEdmeJAIHWLn/BNfNXIsGvr5nKFf0bmV2SEKIemjm68WRF8bRLjIQgIR/LWXf8VyTo7IdSQR29s3mZG7/cCOx4YHMu28YPWJkvAAhnJFSil8fGslk60A3l/13Jcv2ZZgclW1IIrATrTVvLz/IQ19tY1C7ML64ezAtmrtma0Uh3MkLE3vyxOXGTeTbPtjInDWJ5gZkA5II7EBrzQuL9vLS4n1M6NOKD24dSHPpMloIl3H7BfF8cOsAAJ5csIs3f3XuMQ4kEdiYxaJ5asEuZq08zM1D2vLqdX3w8ZK3WQhXM6pLFAumDQPglZ/28+LivSZH1HDyC2VDFovmsXk7mLM2ibuGx/P0ld2lZpAQLqxX6xCWPGAMqvjO8kM8MX+nyRE1jCQCG7FYNNO/3c7nG4w2Av83rqu0FBbCDXSODmLZwyMB+GhtEs98v9vcgBpAEoENaK155ofdfLkpmb9e1IGHL+ssSUAINxIfEchvjxp9FM1efcTp7hnYJBEopcYopfYppQ4qpaZXs36kUipbKbXV+niirvs6g1d/3s+HaxK544J4/ja6k9nhCCFM0Do0gJ//ZhQTvfLTfj5Zl2RyRHXX6A7vlVKewFvAaIyB7DcqpRZorateH63SWl/ewH2brPdWHub1Xw8aA8mMl+Igt6I1lBUZD62hONd4eHiBtz/4BoFPM/CSbkTcRccWQcy7bxhXvbWaf87bSUQzH8b0aGl2WLWyxcgnA4GDWuvDAEqpucAEoC4/5o3Z13Tfb0vluYV7GN+rJc9d3VOSgCsoK4GcZDiVBKePQm4a5GVA/gnIz4SCTOPHviQfSvJAN7LfmaCWENMf2gyC+BEQ3Qs8pMTWmfVpE8Indwzipv+t555PtjjFeMi2SAQxwLFK88nAoGq2G6KU2gakAg9rrXfVY1+UUlOBqQCxsbE2CLtxNidl8dBX2xgQF8p/ruuNp9QOci4lBZCxB07sMZ4z9kDmfshOBqoM1uQfBoGRxiOq2x9n+j6BxsPLDzw8jdc4th4sZVBaaE0YecZ8TXLTYO8PxqM6PSdBnz9B/IXGMYRTuKBjBM9e1YPH5+1k7Gur2PL4aMKacAeTtkgE1f0CVh32bAvQVmudp5QaB8wDOtZxX2Oh1rOAWWCMUNbgaG0g6WQ+d320mZgQf2ZNScDXS/5BmzSLBU7sheSNkLIZUrdA+m7Q5cZ6Lz+I6ASxQyCsHYTEGo/QtsYZu6edGgNqDaeTjJiOrocjK4w4K9vxlfGorM9NMPQvENXFPnEJm5gyuC3bjp3m683J9Hv2Zw48NxZvz6Z5tWeLRJAMtKk03xrjrP8MrXVOpemFSqm3lVIRddm3qckrLuP2Dzdi0ZrZt0o30k2S1pB5wPhhTVwFib9BwUljnV+wURQz/EFo2ds4ww+NM+dsWynj2KFx0OOac9cXZMGub2Hr55BSaWjWrZ8Yjwo9J8FFjxuJSzQpr0zqzbrDJ0k+VUjHxxaROGO82SFVq9FjFiulvID9wMVACrAR+JO16Kdim2ggXWutlVIDga+BtoBnbftWx6wxi7XW3PvpFn7anc7HdwxkaPsIh8cgalBWAkmrYf9i2LfIONMGaN4a4odD3HCIHWyc8TvzvZzC07D5Q1j9GhRmnbveOxCuehu6TXDuv9OFaK2J/8dCACYPjOWFiT1Ni6WmMYttMni9tbjnvxg/7LO11s8ppe4B0FrPVEpNA/4MlAGFwINa6zU17Vvb8cxKBLNWHuL5hXv5v3FdmDqivcOPL6ooL4VDy2DnN7D3RyjJNYp54i+ETpdB+1EQGu/aP4ilhbDuHfjl6erXXzUTet/g2u+BE8guKKX3Mz8B8OFtAxjZOcqUOOyaCBzNjESw9tBJbnx/HWN6RPPWn/pJDSGzaA0pW4yikV3zjLNiv2DoegV0Hg/tRoJPgNlRmif1d5h3H2RUuagOjIKbvjaKw4Qp1h46yeT31gGw7YlLCQ5wfEeUkgga4XRBCWP+u4oAX08WTLuAZr62uLUi6qUo27hpuvlDOL4DvAOg81jocS10uBi8fM2OsOk5sR++uvXcpDD4Prj0WamFZIJnvt/N7NVHAEy5XyCJoIG01kz77HeW7DouA8uYIesIrHsbfv8ESgsguif0v824QerXtOtmNyl7F8LcyWcvC+8Aty2CZuYUU7iruOk/AnDr0DieurK7Q48tiaCBvt6czMNfbeORMZ25d2QHhxxTYBT/rH4N9iwA5Qk9r4WBd0GrflLe3RjFufD17XDgp7OX/2ULhMt9L0fILiyl99PG+//T30bQqUWQw44tiaABUk4XctmrK+nWqjmf3zVYGo05Qtp2WPY87F8Evs0h4TYYdA80l3GebUpr4ybzkn+cvXzaZoiQEx57m781hfvnbgXgyAvjHHbPsaZE0DRbNzQBWmuemLeTcovm35Ok5bDdZR6AL2+Gd4fD0TVw0T/hb7tg9DOSBOxBKRhyLzyVDdf874/lb/aHF2KNNgzCbib0iSHQx7hH8+g3202ORhJBjRbtPM4vezN4cHQn2oS5cS0UeyvKhiWPwduD4eAvMOIRuH87jPi73ANwlJ7XGgnh2g+M+eJseCkePrnGaJUt7GLjPy8B4MtNyaRlF5oaiySCamQXlvLUgl10b9Wc24bFmR2Oa9LauAH8Rn9Y+xb0ngx//R0uegz8Q8yOzj31mGgkhFGPGfMHl8IzobB7gblxuagAHy+evKIbAENe+NXUWCQRVOO1pQfIzCtmxsReeDXRvkGc2qkk+GgCzL/PaPA1dRlMeFNqrzQVFz4Cj580bswDfDkFngqGopzz7yfq7bZh8Wem529NMS0O+ZWrIjEzn4/XJXJdQht6tpaqojZlscCG9+DtIUZHa+P/A7cvgVZ9zY5MVOXpZSToe9f/sWxGG9j2hXkxuajl1mEu75+7FbMq70giqGLGor14e3rw4KUy0phN5Z+Ez6+HhQ8bff7cuw4G3CF97zd1UV2M4qKEO4z576bCK53AUm5uXC4kLiKQNmH+ADw+f6cpMch/YSUbjmSxeNdx/nxhe6KC/MwOx3UkroaZF8Dh5TDuFbjpGwhpU+tuogm5/D9w30ZjOi8dngmDU4mmhuRKljxgDHH5ybqjFJc5PslKIqjk3z/tIyrIlzuHtzM7FNegNax/F+ZcbgzdeOdSo1GYNAhzTpGd4MnTxtgNAK/1lqIiGwnw8WJCH6Oa9N0fb3b48SURWK0/fJL1R7K458L2+PtIHyyNVl4KPz4Iix6BTmPh7hXS4ZkrUAqmbYRLrZ0EfzcVFvzF3JhcxL8nGf8fy/edoKTMsdV2JRFYvfHrQSKa+TB5oPnDYDq94jz4dBJsmg3DHoDrPzGGdxSuY+g0mLrCmN7ykVH0JxrFy9ODy3sZA93f+6ljrwokEQC/Hz3FbwczmTqinVwNNFbhKfj4KjiyEia8BaOflhvCrqpVH3j4oDF9fAc8HWoUB4oG++/1fQBYuicDi8Vx76VN/kOVUmOUUvuUUgeVUtOrWX+jUmq79bFGKdW70rpEpdQOpdRWpZTjR5sBZq9OJMjPixsHyVB/jZKfCXOugLRtcN0c6HuT2REJe2sWCY+lG9PaAk+HSDJoBC9PD4a0Cwfg5Z/2Oey4jU4ESilP4C1gLNANmKyU6lZlsyPAhVrrXsCzWAehr2SU1rpPdZ0h2VtGThGLdqQxqX8bAmWcgYYryoaProLMgzD5c2OgGOEevP3giUp9Ez0dYlooruDdm/sD8M7yQw47pi2uCAYCB7XWh7XWJcBcYELlDbTWa7TWp6yz6zAGqW8SPttwlDKL5uYhcjXQYKWF8PlkOLEXbvgEOlxidkTC0Tw8jdbIFf5r3ri8zq653x8jl/12INMhx7RFIogBjlWaT7Yuq8kdwKJK8xr4SSm1WSk1taadlFJTlVKblFKbTpw40aiAK5SVW/hs/VFGdo4kLiLQJq/pdiwW+OZOSFoDV8+UJODOPL3gn9b/zdNH4du7zY3Hic27bxgAN/1vfS1b2oYtEkF1lcKrLSRUSo3CSASPVlo8TGvdD6No6T6l1Ijq9tVaz9JaJ2itEyIjIxsbMwCrD50kI7eYGwZI46YGWzED9v4AY14werEU7s3LBx4xhmJk+1zYPMfceJxUnzYhZ6ZLy+1fldQWiSAZqPxL2hpIrbqRUqoX8D4wQWt95hpSa51qfc4AvsMoanKI+b+nEOTnxcjO0tlZg+z9EVa8CH1uNAaPEQIgIAz+vNaY/v6vxnCjot4u7dYCgNd/OWD3Y9kiEWwEOiql4pVSPsANwFn91iqlYoFvgSla6/2VlgcqpYIqpoFLAYd0tlFQUsbiXccZ37Mlft5SZbTeTh+F7+4xOowb/x9pLSzO1qKb0Z0IwOt9ZFyDBnjxml6A0cbJ3hqdCLTWZcA0YAmwB/hSa71LKXWPUqriNPEJIBx4u0o10RbAb0qpbcAG4Eet9eLGxlQXy/aeoKCknCv7yOhX9WaxwLx7jeqCkz40ao0IUdXAuyDY2kBTGpzVW2igz5npwhL79j9kk3YEWuuFWutOWuv2WuvnrMtmaq1nWqfv1FqHWquInqkmaq1p1Nv66F6xryP8sjedkABvBsWHO+qQrmPDLEhcBZc9D6FxZkcjmrL7txnPGbsgaa25sTiha/oZFSzfWW7fqwK3bPJpsWhW7DvBhZ0iZSzi+spJg1+egY6XQr+bzY5GNHUeHnDLD8b0B2OksVk9TR/bBYDX7Vw85JaJYHtKNifzSxglN4nr75dnwFIKY1+S+wKibuKHQ5DRhw4/PmhuLE4mMsj3zLQ9B61xy0Swcv8JlIIRnWxTDdVtpGyGbZ/B4HshLL727YWo8JctxvOm2UbPtKLO4sIDAGO8FHtxy0SwKekUnaKCCKt0M0bUwYqXwD8Mhj9kdiTC2fgEQDdrhwOfTjI3Fidz36gOALyzwn5dTrhdIii3aH5POkX/uFCzQ3Eu6btg/2IY/Gfwa252NMIZXfuB8Xx4mQx1WQ8T+hgdNSzfZ5seFarjdongQEYuucVlJLSVRFAvq18Dn2Yw4E6zIxHOysMTOl5mTC95zNxYnIiPl/1/pt0uEWw7dhqAvrGSCOqs8BTsmgd9/mS0GhWioa6dbTyvf8fcOJzUqfwSu7yu2yWC/el5+Hl70DYswOxQnMeu76C82OhKQojG8G32x3Tq7+bF4WSu7G00fF2087hdXt/tEsGBjDw6RDXDQ9oP1N22uRDVTcYcFrZx3UfG87c1djYsqqhIBN9vO6cbN5twv0SQnkvHKBk/t84KsuDYBqPGh7QbELbQ9UrjOXP/+bcTZwzvFAHA2sMna9myYdwqERSVlpOWXUQ7GXug7g79Cmhof7HZkQhXUfmEIi/DvDiciK+XfTvGdKtEkJ5TBEB0sHSSVmeHl4FfCMT0MzsS4UoG32s8r5Obxk2BWyWC49mSCOotdRvE9Deq/glhKxWJYPVr5sYhADdLBOm5xQBEN5dEUCdlJcY4xNE9zI5EuJoQ61hWWhqW1VdBSZnNX9OtEsHJPCMRhDfzrWVLAUDWIaODuRaSCIQwW6cWRtXbo1kFNn9tt0oEBdbBHQJ8pJijTnKsVdVCYs2NQ7g2Gb2sTqKCjJKMjJxim7+2TRKBUmqMUmqfUuqgUmp6NeuVUup16/rtSql+dd3XlgpLyvFQ4OuAJtsuoaJGR6D00irsoOJ7lZdubhxOIsraJXVGbhNMBEopT+AtYCzQDZislOpWZbOxQEfrYyrwTj32tZnC0nICfLxQUh++bgoyjefACHPjEK6puXWY2Fz7NJJyNZFnEkGRzV/bFqfGA4GD1mEnS4C5wIQq20wAPtKGdUCIUqplHfe1mdJyC16ekgTqzGK9KeUp3XULOwi0DgyVn2luHE6ioiSjpMz2RWm2SAQxwLFK88nWZXXZpi77AqCUmqqU2qSU2nTiRMO6Y/Xy8KC8XIbKE6JJqLgyl+Er66TE+ttlj95IbfGK1Z1iV/1ka9qmLvsaC7WepbVO0FonREY2rMza21NRKjem6sH68Ujf8cIeCqzdJQSEmxuHk6i4EvDxtH0i8LLBayQDbSrNtwaqFvrVtI1PHfa1GS9PRZlcEdRdRZfThVln9xophC3kWnvSbCZjh9dFSblxQtZUrwg2Ah2VUvFKKR/gBmBBlW0WADdbaw8NBrK11ml13Ndm/L09KbNoisvkDLdOmrUwnvPsNzKScGM5KcZzcJvzbycASLdWG40Ksn07qEZfEWity5RS04AlgCcwW2u9Syl1j3X9TGAhMA44CBQAt51v38bGVJNQ6xjFpwtKadFc2hLUqiIR5KQA/U0NRbgwD6nOXRdHMvMBiI+w/dW5LYqG0FovxPixr7xsZqVpDdxX133tJdyaCE7mldBCupmoXURHQEHGbuh2pdnRCOHWDp/IA6BNmL/NX9utUnFogJEIThXYZ7g3l+MTCOHt4fgOsyMRrqYgy+wInI7FenszwMcm5+9ncatEUHEVkHq60ORInEh0L2NIQaniJ2xp4/vGc8Lt5sYhADdLBDGh/nh6KJJO2r7TJpcVP8K4RyCjSQlbWvGi8TxkmrlxCMDNEoG3pwetQ/1JPJlvdijOo4N1ZLKDS82NQ7iWilbr4e3NjcNJlJXbt/2TWyUCgLbhgZII6iMkFiK7wJ4fzI5EuIqktcazX7C5cTiR3Wk5AHSMsk97HrdLBB2jmnEgPc/uGdal9JwER9dA1mGzIxGu4IsbjedrPzA3DifyzeZkAK7s3cour+92iaBX62CKyywcyMgzOxTn0XsyKA/Y+rnZkQhnV172R9cSFcWOolafbTgKwKQE+zS+c7tE0DPGuBzdkZxtciROJDgG2l8Mmz+EUtt3gSvcyGLrkCNdLjc3DidTau0ax17jrbtdIogLDyTI14ttyafNDsW5DLsf8jPg94/NjkQ4K61h43vG9MRZ5sbiRLQDqm67XSLw8FD0bRvKhiPSoKVe4i6A1gNh9evGoPZC1NevzxrPIW2NxoqiTn7da4wU2K1lc7sdw+0SAcAFHcI5kJHH8Wwp5qgzpeDCRyH7KGyQszlRTxYLrPq3MX3PKnNjcTIvL9kHwN8v62y3Y7hlIhjWwRh6cfVBGRmpXjpcDB0uMRoDVYxnLERdfDLReG47TKqN1tPe47kAjOxsv7HD3TIRdI1uTnigD6sOSPfK9aIUXPYClBbA0qfMjkY4i8yDcHiZMX2z3XqZd0lFpX90mW/PsdbdMhF4eChGdo7i170Zdhn/06VFdoKhf4Wtn8K+RWZHI5zBm9YuzCe+D5627zDNlf136QEArunX2q7HcctEADCuZzQ5RWWsPiTFQ/U2cjq06AEL/iIDj4vz+2LKH9O9JpkXh5OaueIQAP8c39Wux3HbRHBBxwia+XqxaEea2aE4Hy9fuPpdKMqG7+6RMY1F9Q7+AnusRUH/J/9n9VW5tKJiUC17aVQiUEqFKaV+VkodsD6HVrNNG6XUMqXUHqXULqXU/ZXWPaWUSlFKbbU+xjUmnvrw9fLkkq5RLNmVLkNXNkR0DxjzAhz8GZY9Z3Y0oqnJz/zjBvFti8AnwNx4nNCrS40ef8f1jLb7sRp7RTAd+EVr3RH4xTpfVRnwkNa6KzAYuE8p1a3S+le11n2sD4eMVFbh6n6tyS4s5add6Y48rOtIuAP63WxUC9z5rdnRiKaivBRetvYqOugeaDvU3Hic1DvLjWKhFyb2svuxGpsIJgBzrNNzgKuqbqC1TtNab7FO5wJ7gJhGHtcmhneIICbEn7kbj5odinNSCsa9Am0GwXd3w5GVZkckzKY1PGtUzya8A4x90dx4nFTFsJQAwf7edj9eYxNBC611Ghg/+EDU+TZWSsUBfYH1lRZPU0ptV0rNrq5oqdK+U5VSm5RSm06csE21Tw8PxfUD2rD64EmSpGvqhvHyhclzIaw9fD4ZUraYHZEw07OV6rr/ZbN5cTi5Wz7YAMCr1/d2yPFqTQRKqaVKqZ3VPCbU50BKqWbAN8ADWusc6+J3gPZAHyAN+HdN+2utZ2mtE7TWCZGRtmtYcV1CGzwUfLZergoaLCAMpnxnPH9yDaRuNTsiYYbnW4Ol1Jh+/KS5sTix0nILx7KM4XSv7mvfaqMVak0EWutLtNY9qnnMB9KVUi0BrM/VNjdVSnljJIFPtdbfVnrtdK11udbaArwHDLTFH1Uf0cF+jOvZks/WHyWnqNTRh3cdzVvCzfONPmTmXAFH15kdkXAUreHpMCgxWsDyzxPSXqARHv1mO+CYm8QVGls0tAC4xTp9CzC/6gbKaA73P2CP1vo/Vda1rDR7NbCzkfE0yN0j2pNbXMbnclXQOGHt4PbF0CwKPr5ahrd0B2Ul8HQIaGvNu39mgJd9qzq6Mq01325JAeDV6/s47LiNTQQzgNFKqQPAaOs8SqlWSqmKGkDDgCnARdVUE31JKbVDKbUdGAX8rZHxNEjP1sEM6xDO7NVHpKVxYwW3NqoLhrWHT6+DDe+ZHZGwl7wT8K9KxbRPZBn3jESDvbvSGAWwfWQgvl6eDjuuckRf17aWkJCgN23aZNPXXLn/BDfP3sCzV/VgyuC2Nn1tt1ScC9/cBfsXGdVMx74Invav/SAc5NAy+PgqYzq8g9wYtgGtNfH/MM6fN//zEsKb2T6pKqU2a60Tqi5325bFVQ3vGMHAuDBe/+UAhSXSwKzRfIPghk+NAW02/Q8+HA+nj5kdlbCFhX//IwkMvk+SgI1UdDcdGxZglyRwPpIIrJRSPDKmMydyi/lgzRGzw3ENHp4w+hm45n+QvhtmDoM9P5gdlWioohx4KviP8Shung9jnjc3JhdhsWjetjYg+/GvFzj8+JIIKkmIC+OiLlHMXH6I0wUyCpfN9LwW7llp3Ez+4kajs7oiGTPaqWz7AmZUGjj90URoN9KsaFzO377cCkBC21CC/BxfhCqJoIpHxnQmr7iMf/+03+xQXEtYO7j9Jxj2APz+Cbw1GPYtNjsqUZuibOMq4LupxvyAu+CpbPCvse2nqKdT+SXM35oKwKd3DTIlBkkEVXSJbs6UwW35dH0SO1PkrNWmvHxg9NNw51LwD4HPr4evboPsZLMjE9X55VmYEfvH/H0bYfwr5sXjovo++zMAD47u5NCaQpVJIqjGg5d2JjTAhyfm78Ricb5aVU1eTH+YugJG/h/sWwhvJMDyF6GkwOzIBBg1gp4KhlXWH/0h04yrgMhO5sblgiqPkvjXizuaFockgmoE+3vz6NgubDl6mi83SU0Xu/DygZGPwrSN0HkMLH8e3hwAWz6G8jKzo3NPJ/YbCaCiRpCHt3Ev4DLpZtweLBbNlP8ZfQrNv2+YqbFIIqjBtf1aMyg+jOd+3EPq6UKzw3FdIbEw6UO49UejRfKCafDWANj+pQx44ygVCeCtAX8sm7ocnsiUewF2dM3MNQB0iGpG7zYhpsYiiaAGHh6Kl6/tTbnWTP92B87Y8M6pxF0Ad/1q9GTqHQjf3mVcIWyaDaWSiO0iefO5CeC6j4xioFZ9zYvLDWxKzOL3o6cBWPLACHODQRLBecWGB/CPsV1Yuf8EczdKEZHdKQWdx8LdK2HSHKNR2g9/g1d7GPcQ8qrt01DUh9aw9m0jAbx/0R/LKxJAt3p1KiwaoNyiuXbmWgA+vXMQnh7K5IhAugisxY2D2rJo53Ge/WE3A+PDaB/ZzOyQXJ+HB3S/yvhRSloNa94w7iEstzZeGvYAXPS49HBZH9kp8MVNkFplvIhbvod4889I3cmwGb8C0CU6iGEdIkyOxiB9DdXB8ewixr2+iqggX+bdNww/b3OqeLm1zAPw5jldpBj3F7pf7fBwnEJZMfzyDKx98+zlkV2NVsFBLcyJy419uPoIT32/G4DDz4/Dw8FXAzX1NSSJoI6W7cvgtg828qdBsTx/dU+HHltUUpQN8+6FvVW6qvDyh6veNpKCMv9S2zQlBbDsuXN//AHGvgQDp7r3+2OiY1kFDH9pGQCrHhlFm7AAh8dQUyKQa+s6GtU5irsvbMe7Kw4zKD6MCX2axLDL7scv2OjMDiB9F3wxBbIOQVkhfH2b8QAY+hcY/pB71HpJ2QI/PQ5Jv527rtf1MO5l430Tpikrt5xJAk9c3s2UJHA+ckVQD6XlFm58bz3bkk/z1T1D6NU6xOExiBqkbIbvH4Dj289d5+kLI/4OA+90jcSQuhVWvwa7vq1+fb+bYfSzRutt0STETf8RMHoWXfnIKNPisEvRkFIqDPgCiAMSgeu01qeq2S4RyAXKgbKKQOq6f1VmJQKAzLxiJry5mnKLZsG0YUQ19zMlDnEe+SdhxYw/esmsTtth0Pcm6DK+aZ8tnzwEWz+FrZ9Bblr12wSEw6XPQa/rjB5fRZNy76ebWbjjOABHXhiHMrFozl6J4CUgS2s9Qyk1HQjVWj9azXaJQILWOrMh+1dlZiIA2J2aw7Uz19CpRRBzpw6Wm8dN3clDsPq/sOWj2rdte4FRi6bdhdCqn2OGXczPhGPr4fAKOPATnKqlG/TmMcY4D31vMsaIFk3W3A1Hmf7tDgB2Pn0ZzXzNLY23VyLYB4zUWqdZxx9errXuXM12iVSfCOq0f1VmJwKAxTuPc88nmxnXM5o3JvdrEnWBRR2Vl8KBn40z7ao3nevCwwuCWhqPwAijvYOnt1EE5elttIguLYTiHGOktsIsyE2HvOP1P1ZgJPS5EXpPhsjOcqPXiaw+mMmN768HYNH9w+nasrnJEdkvEZzWWodUmj+ltT6nEFYpdQQ4BWjgXa31rPrsb103FZgKEBsb2z8pKanBcdvK+6sO868f9zBlcFuemdDd1Es+YQMl+XB0HRxZYZydp211zHF9g6F1AsQPh46XQlQ3+cF3cgfScxn96koAZt7UjzE9WpockaHBtYaUUkuB6GpWPVaP4w/TWqcqpaKAn5VSe7XWK+uxP9bkMQuMK4L67Gsvdw5vx4ncYt5deZioIF/+YmLvgcIGfAKhw8XGoyZaQ1kR5KVDThoUnDTO+stLjIelDJQHePuDb3PwbQZ+oRAUbTykDN/lZeQWnUkC/zeuS5NJAudTayLQWl9S0zqlVLpSqmWlop1q+wDQWqdanzOUUt8BA4GVQJ32b8oeHdOFE7nF/Pvn/YQEeDNlSJzZIQl7Usr4kQ+NMx5CVHK6oISBz/0CwA0D2jB1RHuTI6qbxvY1tAC4xTp9CzC/6gZKqUClVFDFNHApsLOu+zd1Hh6KF6/txSVdo3h8/i4+W3/U7JCEECbILiylzzPGIDPDOoQz45peJkdUd41NBDOA0UqpA8Bo6zxKqVZKqYXWbVoAvymltgEbgB+11ovPt7+z8fb04K0b+3FRlyj+77sdzN0gyUAId5JbVErvp38CYEBcKJ/eOdjkiOpHGpTZUHFZOXd/vJnl+04wY2JPbhgYW/tOQginll1QSu9njCTQu3Uw86ddYHJENavpZrF0Q21Dvl6ezLypPxd2imT6tzt4b+Vhs0MSQthRWnbhmSTQrWXzJp0EzkcSgY35eXsy6+b+jO/ZkucW7uHlJXtlUBshXNDBjFyGvGB0KX1RlygW3j/c5IgaTjqdswNfL09en9yX5v7evLXsEKcKSnl2Qg9pdCaEi9iYmMUk6+AyNw2O5V9XOXePxJII7MTTQ/H81T0IDfDm7eWHSM8u4vXJfQk0uYm5EKJxPl6XxOPzjIqPD43u5BLth6RoyI6UUjwypgvPTujOsn0ZTJq5lrRsGX9XCGf14BdbzySBt/7UzyWSAEgicIgpQ+KYfesAjmYVMOHN1exIzjY7JCFEPWitSfjXUr79PQWAxQ8MZ3yvpt9iuK4kETjIyM5RfPPnoXh7enDtzDV8uyXZ7JCEEHVwMq+Y+H8sJDOvGICtT4ymS7T5HcjZkiQCB+ocHcS8+4bRp00ID365jSfm76SkzGJ2WEKIGqzcf4L+/1oKgL+3J4eeH0dIgAO6JncwSQQOFhnky6d3DuLOC+L5aG0Sk99bR3pOkdlhCSGqeHzeTm6evQGA24fFs+fZMS5b809aFpvoh+2pPPL1dvy9PXllUm9GdYkyOyQh3F5+cRndn1xyZv7TOwcxrEOEiRHZjrQsboIu79WKBdOGERnky20fbuTp73dRXFZudlhCuK1l+zLOSgJbnxjtMkngfCQRmKxDlHHf4NahcXywOpGr3lrDwYxcs8MSwq1orbnx/XXc9sFGAK7p15rEGeNd8n5AdaRoqAn5dW86D3+1nfziMh6+tDO3XxDvsmWSQjQVu1NzGPf6qjPz3/x5KP3bVjtQotOToiEncFGXFix+YDgjOkXy3MI9XPfuWg6dyDM7LCFcktaaO+dsOpMEooJ82f+vsS6bBM5HrgiaIK0187em8uSCXRSVlvPQpZ24fVg8Xp6St4WwhZ0p2Vz+xm9n5t+d0p/Lulc3Iq9rafCYxcLxlFJc1TeGoR3Ceey7nTy/cC/zfk/luat70DfW/c5WhLCV0nILY19bxcEM40o72N+bDY9djK+Xe48l3ahTTKVUmFLqZ6XUAevzOb9SSqnOSqmtlR45SqkHrOueUkqlVFo3rjHxuJqoID9mTenPOzf242R+MRPfWcNj3+0gu6DU7NCEcDpfbjxGx8cWnUkCH9w6gG1PXur2SQAaWTSklHoJyNJaz1BKTQdCtdaPnmd7TyAFGKS1TlJKPQXkaa1fqc9xXb1oqDp5xWX856f9fLjmCGGBPjw6pgvX9GuNh9xMFuK8jmTmM+qV5WfmL+4Sxfu3JKCU+/3v2KtoaAIw0jo9B1gO1JgIgIuBQ1rrpEYe1+008/XiiSu6MbFfDI/P38nfv97OnLWJPD6+G4PahZsdnhBNTmFJOWNeW0nSyYIzy1Y9Moo2YQEmRtU0NfaK4LTWOqTS/CmtdY2F2Eqp2cAWrfWb1vmngFuBHGAT8JDW+lQN+04FpgLExsb2T0py31xisWi+357Ki4v2kppdxNge0fxjbFdiw+ULLoTWmqe/382HaxLPLHt9cl+u7N3KvKCaiJquCGpNBEqppUB1t9MfA+bUNREopXyAVKC71jrduqwFkAlo4Fmgpdb69tr+GHcsGqpOYUk57606zDvLD1FmsXDjoLbcO6o9UUF+ZocmhCneX3WYf/2458z85IFteP7qnm5ZDFSdBhcNaa0vOc+LpiulWmqt05RSLYGM87zUWIyrgfRKr31mWin1HvBDbfGIP/j7ePLXizty/YA2/Hfpfj5el8QXG49x27A47h7RnuAAb7NDFMIh5m9N4f65W8/Mx0cE8sNfLpARAeuose/SAuAWYIb1ef55tp0MfF55QUUSsc5eDexsZDxuqUVzP16Y2IupI9rz6s/7eXv5IT5el8TU4e24eWgcwf6SEIRrWrwzjXs+2XJmPsDHkxV/H0VkkK+JUTmfxt4jCAe+BGKBo8AkrXWWUqoV8L7Wepx1uwDgGNBOa51daf+PgT4YRUOJwN2VEkONpGjo/Han5vCfn/exdE8GQb5e3DI0jtsviCcs0D36TRGu7+vNyTz81bazlq34+0jahgeaFJFzaPA9gqZIEkHd7EzJ5q1lB1m86zh+Xp78aVAsdw1vR3Sw3EMQzqnqPQCApQ9eSIeoZiZF5FwkEbixA+m5vL38EAu2peKh4Iperbj9gnh6xASbHZoQtSouK+fhr7bz/bbUM8t8vDz45cELpSpoPUkiEBw9WcDs1Uf4atMx8kvKGdwujDsuaMfFXaKkYZpoctJzirjmnTUknyo8s6xdRCBz7x4sNeMaSBKBOCO7sJQvNh7lw9WJpGYXERsWwJ8GxTKxX4z8gwnTffd7Mn/74uzy/4l9Y3jx2l54S8eLjSKJQJyjtNzC4p3H+XhdEhuOZJ1Z/tDoTtw3qoNcJQiHyS4sZdpnW1h1IPOs5c9M6M7NQ+LMCcoFSSIQ57Xt2GkmvLX6rGVBfl58dPtA6fFU2IXWmk/WJfH4/F1nLW8V7Mdndw0mLkJqANmaJAJRJ1n5Jfz5k82sr3SFADCycyT/ua6PVEEVjbY9+TS3zN7AqSq96N47sj0PX9pZrkTtSBKBqLeV+09w8+wN5yyfPLAN/zeuK0F+0lBN1M3BjFwe+GIrO1Nyzlreu3Uwb9/Un5gQf5Micy+SCESDlVs0by87yL9/3n/OuimD2zJ9bBdpyi/OkZiZz8NfbWNT0rn9SH5w6wBGdYkyISr3JolA2ERhSTkvLt57Vs+OFUZ1juTpK3tIL6hubOX+E0z/Zjup2UXnrPv3pN5M7BcjHcCZSBKBsLm84jJmLNrDJ+uOnrMuJsSfxy/vxmXdW8g/vgsrKCnj/VVH+E81V4sAz13dgz8NjJXvQBMhiUDYVVm5hY/WJvHMD7urXT+2RzR/v6wz7SKlKwBnprVm5YFMnvtxN/vT885Z3zrUnxkTe3FBxwgTohO1kUQgHGrtoZM8v3APO1Kyq10/qX9r7h3VgXipItikaa1ZffAkr/9ygA2JWdVuM7FfDA9d2llu+DoBSQTCNMVl5Xyy7ijPL9xDuaX679slXaO4dWg8wzqESzGCiQpLyvnu9xQ+XHOk2jN+gPaRgfxzfDdGdo6Uz8rJSCIQTUZuUSkfrU3ijV8PUFRqqXYbLw/FrUPjuH5AGzq2CHJwhO6h3KJZeeAEX248xqKdx2vcrnOLIP56cUfG9IjGU+r4OzVJBKLJKiot54ftaXy45sg59cyruqx7Cyb2a82FnSLx8/Z0UITOLyOniCW7jvPNlhS2Hjt93m3H9ojm1qFxDIwPkzN+FyOJQDiVtOxCPl9/lK82J5NWTVXEqkZ0imR0txZc0CGCuPAAt/wBKy23sO3YaVYdyOTn3ensTjt/UgXo3zaUa/u3ZkKfVgT4SFsQV2eXRKCUmgQ8BXQFBmqtq/11VkqNAV4DPDFGLpthXR4GfAHEYYxQdp3W+tzWJ1VIInBPuUWlLNp5nHm/p7Dm0Mk679c+MpDebULo3TqE3m1C6BId5JRXE6fyS9iWfJptx7Ktz6c5mV9Sp30DfTy5rHs0V/WNYViHCCnicVP2SgRdAQvwLvBwdYlAKeUJ7AdGA8nARmCy1nq3UuolIEtrPUMpNR0I1Vo/WttxJRGIyvKLy1h14ARL92Sw5mBmtY2Z6sLH04N2kYHERwTSorkf0cF+RAX50qK5Hy2a+xLZzA8/Hw98PD0adMVRbtGUlFnIKy7jRG4x6TlFHM8pIi27iPRsY/pIZj5HswoaFD9Ar9bBDG0fwSVdo+gbGyo/+OIsNSWCRl0Laq33WF/8fJsNBA5qrQ9bt50LTAB2W59HWrebAywHak0EQlQW6OvFmB4tGdOj5Tnryi2a/em5bE8+zbbkbLYnn67xPkRJuYW9x3PZezzX3iE3SICPJ33ahJy5uunTJoQWzX3dshhM2JYjCgVjMAaur5AMDLJOt6gYrF5rnaaUqrHzEaXUVGAqQGxsrJ1CFa7G00PRtWVzurZszvUDat5Oa01ecRlJJws4kplPek4RGdaz9vScIjJyisnILSavuKxR8fh4eRDk60Wk9UojurkfLYKN5+hgX+LCA2kTFiADsAiHqjURKKWWAtHVrHpMaz2/Dseo7nSl3uVRWutZwCwwiobqu78Q56OUIsjPmx4xwTKWs3A7tSYCrfUljTxGMtCm0nxroGIU6nSlVEvr1UBLIKORxxJCCFFPjrj+3Ah0VErFK6V8gBuABdZ1C4BbrNO3AHW5whBCCGFDjUoESqmrlVLJwBDgR6XUEuvyVkqphQBa6zJgGrAE2AN8qbWuGJtuBjBaKXUAo1bRjMbEI4QQov6kQZkQQriJmqqPStUEIYRwc5IIhBDCzUkiEEIINyeJQAgh3JxT3ixWSp0Akhq4ewSQacNwbEXiqh+Jq34krvpx1bjaaq0jqy50ykTQGEqpTdXdNTebxFU/Elf9SFz1425xSdGQEEK4OUkEQgjh5twxEcwyO4AaSFz1I3HVj8RVP24Vl9vdIxBCCHE2d7wiEEIIUYkkAiGEcHMumQiUUpOUUruUUhalVI1VrZRSY5RS+5RSB61jJlcsD1NK/ayUOmB9DrVRXLW+rlKqs1Jqa6VHjlLqAeu6p5RSKZXWjXNUXNbtEpVSO6zH3lTf/e0Rl1KqjVJqmVJqj/Uzv7/SOpu+XzV9XyqtV0qp163rtyul+tV1XzvHdaM1nu1KqTVKqd6V1lX7mToorpFKqexKn88Tdd3XznH9vVJMO5VS5UqpMOs6u7xfSqnZSqkMpdTOGtbb97ultXa5B9AV6IwxBnJCDdt4AoeAdoAPsA3oZl33EjDdOj0deNFGcdXrda0xHsdoBALwFPCwHd6vOsUFJAIRjf27bBkX0BLoZ50OAvZX+hxt9n6d7/tSaZtxwCKMUfkGA+vruq+d4xoKhFqnx1bEdb7P1EFxjQR+aMi+9oyryvZXAL864P0aAfQDdtaw3q7fLZe8ItBa79Fa76tls4HAQa31Ya11CTAXmGBdNwGYY52eA1xlo9Dq+7oXA4e01g1tRV1Xjf17TXu/tNZpWust1ulcjDEvYmx0/MrO932pHO9H2rAOCFHGyHt12dducWmt12itT1ln12GMEmhvjfmbTX2/qpgMfG6jY9dIa70SyDrPJnb9brlkIqijGOBYpflk/vgBaaG1TgPjhwaIstEx6/u6N3Dul3Ca9dJwtq2KYOoRlwZ+UkptVkpNbcD+9ooLAKVUHNAXWF9psa3er/N9X2rbpi772jOuyu7AOLOsUNNn6qi4hiiltimlFimlutdzX3vGhVIqABgDfFNpsb3er9rY9btV65jFTZVSaikQXc2qx7TWdRnyUlWzrNF1ac8XVz1fxwe4EvhHpcXvAM9ixPks8G/gdgfGNUxrnaqUigJ+VkrttZ7JNJgN369mGP+wD2itc6yLG/x+VXeIapZV/b7UtI1dvmu1HPPcDZUahZEILqi02OafaT3i2oJR7JlnvX8zD+hYx33tGVeFK4DVWuvKZ+r2er9qY9fvltMmAq31JY18iWSgTaX51kCqdTpdKdVSa51mvfzKsEVcSqn6vO5YYIvWOr3Sa5+ZVkq9B/zgyLi01qnW5wyl1HcYl6UrMfn9Ukp5YySBT7XW31Z67Qa/X9U43/eltm186rCvPeNCKdULeB8Yq7U+WbH8PJ+p3eOqlLDRWi9USr2tlIqoy772jKuSc67I7fh+1cau3y13LhraCHRUSsVbz75vABZY1y0AbrFO3wLU5QqjLurzuueUTVp/DCtcDVRbw8AecSmlApVSQRXTwKWVjm/a+6WUUsD/gD1a6/9UWWfL9+t835fK8d5sreExGMi2FmnVZV+7xaWUigW+BaZorfdXWn6+z9QRcUVbPz+UUgMxfo9O1mVfe8ZljScYuJBK3zk7v1+1se93y9Z3v5vCA+OfPhkoBtKBJdblrYCFlbYbh1HL5BBGkVLF8nDgF+CA9TnMRnFV+7rVxBWA8Q8RXGX/j4EdwHbrh93SUXFh1ErYZn3sairvF0Yxh7a+J1utj3H2eL+q+74A9wD3WKcV8JZ1/Q4q1Vir6btmo/eptrjeB05Ven821faZOiiuadbjbsO4iT20Kbxf1vlbgblV9rPb+4Vx0pcGlGL8dt3hyO+WdDEhhBBuzp2LhoQQQiCJQAgh3J4kAiGEcHOSCIQQws1JIhBCCDcniUAIIdycJAIhhHBz/w+/R8vEPKmIaQAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "data1, data_t = simulate(500, state0=[0,1], h= 0.02)\n", "data2, data_t = simulate(500, state0=[0,0.5], h= 0.02)\n", "fig = plt.figure()\n", "ax = fig.add_subplot(1, 1, 1)\n", "ax.plot(data1[:, 0], data1[:, 1])\n", "ax.plot(data2[:, 0], data2[:, 1])" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Add the vector field" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYoAAAEGCAYAAAB7DNKzAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8vihELAAAACXBIWXMAAAsTAAALEwEAmpwYAACUhklEQVR4nOyddVhUWxeHf4duQVBBRQSxEBtbMa7dddVrd93PuBbY3d117cZWDEzExEAwQJRUujum1vfHyFyGZjhnMM77POcZ5sRee4Yze5299gqGiMDDw8Pzq8EwTCYAjWy7CIAYgAYVY+BjGGYhgJEAvgGIBuAJYDARNfp+vDqAs0TUmGEY5vvxBgDKAjgDoByAlwC6AmgMwBLAIQBNAah+PzaYiD4o+FE5R620O8DDw/PrwjDMj/QkykA65kmk4zkkAHyJyCbfCximMYAhABp+v9YDwBsAiQzDNCAiTwBjABz9fklDAF5ERAzDLAXwhIhWMAzTA8BEACCiVwzDXAOwCoA2gJM/spIAeEXBw8Pze6NbyPE2AC4TURoAfB/gAeBfAGMYhpkFYDCkswNAOmu49f1vewD9AYCIbjAME5+t3RUAXgHIADC9pB+Ca3hFwcPDwxlExJSWbIZhMgBoZu8OFDA9fb8uJxcBLAXwAMAbIor9vr8zgAGFXAtIzVJ6ANQBaAFILUZ/lI5KaXeAh4eHhyPUIR2o04mIISIVIlIvppJwA9CPYRhthmH0AfQCACLKAOACYC+AIwDAMEwZAGrZlIYbgGHfj3UDYJSt3QMAFgM4BWC9oh9QWfAzCh4enl8SIlJloQ0PhmHOQbpAHQzgcbbDpyA1Ld35/r4TgHvZji8HcIZhGA8AjwB8BQCGYUYCEBHRaYZhVAE8YximAxE9KGl/uYL5Fb2eTExMqGrVqqXdDR4enl+YyMhIiMViVKxYEQAQHBwMExMT6OoWtuzxY/LmzZsYIiqX17FfckZRtWpVvH79urS7wcPD84vSr18/CAQCPHjwACYmJqXdHVZgGCY4v2O/pKLg4eHh4ZLLly+XdheUCr+YzcPDw8NTILyi4JiAgABcvnwZv+JaEA8Pz+8BryhYhojg7e2NlStXomHDhqhXrx6sra3xPRKUh4eH56eDX6NgASLCmzdvcOnSJVy6dAm+vr6yY05OTqhbt24p9o6Hh4enZJSqomAY5jCAngCiiMg2j+MMgO0AugNIAzCaiDyU28uCISKsXbsWixcvhkQikTs2f/58/Pnnn6XUMx4eHh52KG3T01FIc6PkRzcA1b9vEyGNgvyhYBgGs2bNQt++feX2d+vWDStXrlRaP4gI/v7+OH36NGbMmIGxY8ciISFBafJ5eHh+XUp1RkFEbgzDVC3glD4Ajn8PuX/BMIwhwzBmRBSunB4Wzps3bzBy5Eh4e3tDS0sLGRkZsLa2xunTp6GqWuLA0HyJi4vDy5cv4e7uDnd3d7x8+RKxsdLMATVq1ICrqysMDQ05k/+7IhRLkJwhQlK6EMkZIqQLxdJNIEamSPqaIRRDJCGIJQQxEcRi6esp96+IScmEhqqKdFNTgbqqCtTVGOhqqMFIRwNGuuow1NFAWR0NGOqoo7yBFiobaaOyoTZM9DShosKvdfEonx99jaISpDngswj5vi+XomAYZiK+p/GtUqUK5x0TCoVYu3YtVq5cCZFIhD///BPz5s1D+/btceXKFc4H6QcPHmDSpEmIi4uT21+tWjU8ePAAZmZmnMr/VSAixKUKEJ6YgYjEDMSkZH7fBIhJyURsigBxqQIkpguRlCFEmkBcYpmZIgkyRRIgk4UP8B1TAy3UMNWHbUUDNKpihEYWRiirq1H4hTw8ReBHVxR5PT7l6WdKRAcgTbQFOzs7Tn1RfX19MWLECLx69QqGhobYs2cPhgwZArFYjNOnT6NOnTpcigcA1KtXDw0bNsT9+/dl+ywtLfHw4UNUqlSJc/kSiQQxMTEIDw9HWFiY7LVXr16oX78+5/KLQ0KaAEGxaQiOTUVwbBqCY9MQlpCO8MR0hCdmSAftHOhrqsFEXxPGuhqwMNaBoY46DLTUYaCtDgMtNRhoq0NfSx06GqrQUleFlroKtNWz/laFmioDVYaBqsr3jWHAMIBYQhBJCAKxBEKRBEIxQSCSIDlTiMQ0IeLThIhLEyAhVYD4NCEikzMQGp+OkPh0xKTkr1kikjIQkZQBt8/R+Z5TsYwW2lQvB/sa5dDK2hiGOrwi4SkaP7qiCAFgnu19ZQBhpdQXSCQS7Nq1Cw4ODsjIyEDnzp1x+PBh2cCspqaGXr16cdqHqKgorFixAvv374dIJIK2tjbS09NhYWGBhw8fwtzcvPBGFOTGjRtYsWIFwsLCEBERAZFIJDumrq6OvXv3lpqSkEgIIfHp+BKVjM+RKfgSlQz/qBQExaYhMV0od66pgRYqGmqhTqUy6GRTAWZltFHRUAumZbRR7rty0FLnxmyopspATRWstk9ECE1Ih29EMt6HJsLjawI8guORkimSOy8sMQPnXn/DudffcrVRvbweutc1Q/e6ZqhRQY935+aRo9STAn5fo3DOx+upB4D/Qer11AzADiJqmvO8nNjZ2RHbuZ6+fv2KMWPG4MGDB9DR0cGmTZswefJkpf2g0tLSsHXrVqxfvx7JyckoX748li1bBj09PSxYsACPHj2ClZUVp30gIowZMwbHjh2T21+hQgVcunQJLVu25FR+FmkCEXzCk+EdloiPYUn4GJaEL1HJyBD+NzMwNdCCdXk9WJrowsJYBxbGuqhqrAPzsjqcKYEfGSJCUGwaHn+JhtvnGLh9iYYgj5lUTlpbm+BPu8roUsf0t/zeficYhnlDRHZ5HitNRcEwzBkA7QCYAIiEtBCIOgAQ0b7v7rG7IPWMSgMwhogK1QBsKgoiwokTJzBt2jQkJSWhefPmOH78OKpXr85K+4UhFotx/PhxLF68GKGhodDR0cHs2bMxd+5c6Ovrw93dHcbGxrC2tuasD2lpaTh79iz27NmDN2/eyB2zs7PD5cuXUblyZU5kh4VHYOjf85CsYQyBfiUIdCsgU9MQYKQOe+qSTFTSJXRoVAs1KuijegV9WJfXQxltdU768ytCRHgXkoibH8Jx630EvsalFXh+K2tjjGlpiQ61yvOL678QP6yi4Aq2FEV0dDQmTZqEy5cvQ11dHcuXL8fcuXOhpqYci52Liwvmzp2L9+/fQ0VFBWPGjMGKFStkaY25xtfXF/v27cPRo0dlrrZZs4Znz55h+PDhOHDgALS1tVmTmZgmxMugOLwOjsPb4AS8C02QzRREKXEQhH+BINIfgkh/qCaFYcmc6Zg5cwbU1dlVDEFBQcjIyEC1atVYb/tnIVMkxgOfKJx/E4IHn6LyPY9hgP+1t8a41pb8usdPDK8oFODq1auYOHEioqKiYGtrixMnTqBBgwbsdLAQPD09MW/ePNy9excA0L17d6xfvx62trmsc6wjEolw7do17NmzR7ZQrquri+HDh2PKlCmoX78+xo8fDxsbG/zzzz8lNr0lpAnwMjAOLwLi8CIgFj4RSSACVBmCnjABiQGeiPZ5icywTxAn/bdQ261bN+zZswdc1R2JiIiAra0tEhMTUa1aNdSqVUtuq1mzJoyMjApvqJhIJBIwDPPDrhEIxRK4fIzA4SeB8PiakO95k+ytMMHeCiZ6mvmew/NjwSuKYpCUlISZM2fiyJEjYBgGc+fOxYoVK6Cpyf0N//XrVyxevBgnTpwAEaFhw4bYtGkTOnTowLnssLAwHDx4EAcPHkRoaCgAwMbGBlOmTMGIESNQpkwZuXMVndWIxBK8/ZYAt8/RcPscjXehiSAC1FUAYyQhNdAT/s9vIT3kEyAWQk1NDa1atUJsbCw+fPgAU1NT7NixAwMHDmR1MCUixMTEwN/fX7ZdvHgR7969y/P81q1bY+PGjWjevDlrfQCkpsaePXvC0tISvXr1Qvv27aGlpcWqDLYJiU/Dvkf+OPnia57HDbTUsKinDQY0qgxV3lT1w8IriiLi5uaGkSNHIjg4GJaWljh+/Dhat27NQQ/lSUtLw4oVK7Bt2zZkZmaiSpUqWL16NYYOHQoVFW6D51+8eIFNmzbhypUrEIvFUFNTw4ABAzBlyhTY29uzMhhHJ2fivk8kXH2j8dQ/BskZIqgwgJlGJiTh3vjy2BmxX94AYqmXTo0aNdC5c2d07twZ7dq1g76+Pv744w/UrFkTa9asKXGMysePH/Hs2TM5peDv74+kpKQCr2MYBv3798ecOXMUUhBisRhPnjwp9Lxr165hy5YtAAAdHR106tQJvXr1Qo8ePWBqalpsuQBkyt/MzKzAe2rbtm3o27dviWZqUUkZ2O8WgENPAvM83rF2BSztZQPzsjoKyygKKSkp0NXVLZXZmUAggIaG8s1wYrFY4UDfghQFiOiX2xo3bkyKcOHCBQJAEydOpKSkJIXaUASBQEDVq1enMmXK0IYNGyg9PV1psnfv3k0AyNzcnFatWkXh4eGstPslMpn2PPSjfrufUFVHZ7JwcKbma+7RvPNedONdGEUnppKBgQEBoDJlylD//v1p//79FBAQkGd779+/Z6VfRERLliwhSONxZJupqSm1atWKRo4cScuXL6dTp07RqlWrCABpaWnRlClT6MuXLyWSm5KSkktucbemTZvS1atXSSKRFEv2pEmTZJ/FxsaGevXqRf/88w/t2rWLbt26RZ8/fyaBQEDz5s0jVVVVGjJkCL1586ZEn5eIyM/Pj9TLW5HV2K1k4eCca7Pf8ICe+cUQkfR3UNzPVRD9+vWjBg0a0IEDByglJYW1dgvjyZMnVKlSJTp8+DCJRCKlyRWJRNS+fXtavnw5CQSCYl8P4DXlM6byM4oceHp6Km0tIqdcc3NzGBsbK1VuUlISHj58iB49epRokZ6I4BuZDGevcNx8H46AmFQAQN1KZdCxdgV0sqmA2mb6ck93Z86cQdWqVdGkSROlOQgA0lnUkydPUK1aNVSrVg1WVlbQ09PLdd6cOXOgr6+PqVOnoly5PEsJF4vMzEzMmjWr0PMCAgJw+/Zt2XsNDQ20b98evXr1Qs+ePWFhYVFs2du3b8e1a9fg5+eHb9++5VkfRUVFBWXKlEF8fLxsX4cOHTB37lx06dJFoSdzd3d3TJ06FR4e33N5qqiibLN+0LcfnetcbXUVmPhewfH1jrC0tCy2rOyIRCJ06dIFDx48AAAYGBhg1KhRmDJlCmrXrl2itgtj06ZNcHBwgEQiQd26dbFhwwaFv7/i4O3tjVatWiEhIQFNmjTB8ePHUatWrSJfz5ueeDjDLyoFzu/C4PwuHH5RKVBhgBbVjNG1jik6fg9m+1kpyTS+JAwcOBCPHz9Gjx490KtXL3Ts2BH6+vqstZ+ZmYmgoCD4+fnJmd78/PwQEBAAoVCY65q6detizpw5GDJkiEImlejoaNy7dw8uLi64c+cOwsOlWXi0LBuhfI9/wOjKOwZoMmIcGNUUbWspZmrL4sOHD9i7dy9OnDiB5ORkAEC7du0wdepU9O3bV86jLWssZGNA//DhA+bNm4dbt24BADp27IiNGzdy/hAaGhqKsWPH4s6dO9DS0sL69evxv//9r0gmbF5RKJHU1FRcu3YNgwYNKpVBRhnEpQpw1TMUFz1C8CE0CQwDNK1aFj3rV0Q3W1Pe06UEiEQieHh4wM7OjvP1qZwkJSWhfv36CAoKAiDNNGBtbS3n7dWmTZsSB3YSET58+IA7d+7gzp07cHNzQ0ZGBjQr1YJJbweoGcjP3upVLoM9wxqhspHiaxrJyck4deoU9uzZg/fv3wMATE1NMWHCBEycOBGVK1dGSkoKZs+ejZ07d7K2vnD//n3MnTsXb9++BcMwGDFiBFauXMlpPjoiwr59+zBnzhykpaWhQ4cOOHLkSKEy+TUKJRAVFUVLliyhsmXL0uHDh5Uun2sEIjG5fAinCcdekfWCG2Th4Ew9drjRoccBFJGovDUVHu64desWrVu3jq5cuUKfPn1SyM6tCGlpaeTi4kKdOnWSrcVomttS5elncq1pLLr8njKEitv9JRIJPXnyhIYOHUoaGhoEgFRUVKhv377k4uJCampq1KpVK4qIiGDt84nFYjpx4gRVqVJF+tk0NcnBwYESEhJYk5EXnz9/pubNmxMAMjAwoGPHjhW4BoQC1ihKfVDnYlOmovDz86MpU6aQlpYWASB7e3tWF+RKm6+xqbT+lg81XnmHLBycqfHKu7TK+SP5hCeWdtd4fiHev39PhoaGuRbv1dTUqGyT3rkURlVHZ/IIjiuRzMjISFq7di1VrVo1l9zKlSvT69evWfp0UtLT02nDhg1UpkwZAkDGxsa0fft2yszMlJ2TnzOHogiFQlqzZg2pq6sTAOrXrx9FRUXleS6vKDjg1atXNGjQIFJRUZHdXOrq6uTt7c25bK4RiSV052MEjTrsTlUdncnS0ZnGH3tF97wjSCgSl3b3eH5BYmJiKCwsjGJiYigpKYkyMjJILJa/11IyhPTP2be5lMbSqx9IUIL7UiQSkbOzM5mYmMgpCy0tLTp9+nRJP1ouYmJiaObMmbLBu1q1auTk5EQSiYTq1atHL1++ZF3m27dvydbWlgBQ+fLl6cqVK7nO4RUFywQHB1Pfvn1JVVVV7sZauHAhp3Lz49u3b7R161by8PAoUTuJ6QI68MifWq69TxYOztRk1V3afMeXQuPTWOopDw87PPkSnUthtN/4kCKTim8GFQqFNH78+Hxdkh0cHDhxc/X396fBgwfL5DRr1ox0dXXJ0NCQFdfknGRkZNC8efOIYRgCQGPGjKHExP8sAwUpCn4xW0GePXuGXr16yQoHWVlZ4cOHD6zmPSqI0NBQXLx4EU5OTnj69CkmTJiAAwcOKNTWt7g0HH4aCKdX35AqEKOZZVmMblkVHW0qQF21tKvl8vDkT2qmCP+c88Qd70i5/ecnt0CTqmWL1AYRITk5Wa6uSva/w8PD0aBBA6xbt46T3/fLly8xZ84cPH78WLbPyMgIDx484MRL6smTJxg5ciQCAwNhYWGBo0ePol27dvxiNttcvHhRtibRpk0bAkC3b9/mVCYRUVhYGO3cuZPatGkjeyoAQHZ2dgoF6b0PSaApJ1+TpaMzVZt/g2aefUvvQxI46DkPD7dIJBL693FArlnGVc/Q0u5aoUgkEjpy5IjMFJW1GRsb07t37ziRmZSURBMnTpTJ+ueff/gZBZvs3LkTM2bMgKqqKg4dOgR7e3s4Ojri7NmznMjLIjY2FmPHjsW1a9fk9hsbG+PNmzfFCsJ6ExyHXQ/88NA3GvqaahjW3AKjW1aFaZni5xRKTEyUywPFw1PavA6Kw8B9z+X2rehTByNbVC2dDhVCWFgY9u7dC3d3d7x8+RKJiYmyY+XKlYOrqytsbGw4kX3z5k2MGzcOERERAMDHUZQUiUQCR0dHbNy4Efr6+rh48SI6deoEAIiPj+ckk2hOXr9+jZ49eyIyUjrNVlFRgYuLCzp27Fik65/7x2LH/S94HhALIx11jGttiREtqha5dkN0dDTevHkDDw8PvHnzBu/evcOaNWvw559/KvyZeHi4IjAmFe03ucrtW9C9FibaVyudDhUBiUSCz58/w93dXbZFRUXh7t27xYqyLg6xsbGYMmUKzp8/z5ueSkJGRgYNGTKEAJCZmRm9ffuW1fYLQyQS0apVq0hNTY0AyLwz1qxZU6TrPb/G0/B/X5CFgzPZrbpLB938KSVDWOh1CQkJtHLlSurbty+Zm5vnclu8cOFCST8aDw/nRCSmU8MVd+RMUudefS3tbhWZtLQ08vPz46z9wMBAGjp0KG96KgkJCQno27cvHj16hNq1a+P27ducRlXmJCgoCCNGjMCTJ09QpkwZ7N27FzExMbh//z4uXbpUYPTu58hkbL7jC5ePkTDSUcff7a0xvLlFsUpaXr9+HcOGDZOlPwCkEbtOTk7o169fiT4bD48yiUnJRPtNrkjO+K+W+P4RjdGlTsnShPzshISE4NKlS5gxYwZvelKEb9++oVu3bvj48SPs7e1x5coVpZiYAOlM79SpU/j777+RlJQEe3t7HD9+HBYWFvj48SMqV66c79pAVHIGNrt8htObb9DTUMP4NlYY27oq9LWKV6nt3bt3WLhwIZydnWX71NXVceHCBfTu3btEny8/hEIhdu3ahejoaKSkpCAlJQXJyclyr23btsXatWuho8NtmmqeX5NvcWlos+Gh3L57s+xhXZ69fFo/I7zXkwJ4eXlRxYoVCQANGjRIqam/4+PjZaYuNTU1Wrt2bZH8uNMFItr14AvZLL5F1gtu0MrrHykuJbPQ63Li7+9Pw4YNk3lW2djYkLGxMWloaJCzs7MiH6lY3Lt3j3R0dHL5s6urq9OWLVtyBWKxQXx8PKWmprLeLs+Pi29Ekpw5ym7VXUpKV07akh8R8AF3xePevXukr69PAGjWrFmcDEz54erqKlsPqFmzZpHSCEgkErrmGSoLlJt4/BUFRhc//35YWBhNnTpVthZiYWFBx44dI5FIRE2bNqWbN28q8pGKTFxcHB08eJDatWsn5/4LgKytrVlPqZCdL1++kIGBAZmYmFCjRo2ob9++NH36dNq0aRM5OTnRixcvWKvVkR2BQEAxMTGst8tTdG69D5NTGPPOe/1SaXiKCq8oisGJEydITU2NGIahrVu3KtxOccnMzCQHBwfZADl58uQiPeEGRKfQsIPShepu29xkRWCKQ3x8PM2fP1/2FF+uXDnasWMHZWRkyM75/PlzsdstCmlpaeTk5ER9+/aVJWnD90JKWX8PHz6cs0JSYrGYvnz5QhcvXqQOHTrkGZnLMAwNGTKEPn78yLp8iURC9vb21KZNG9qwYQN5e3srdZDKHpn7OyORSGj1DW85hfHc//dS4LyiKCLbtm2TZXc8f/68Qm0oQkREBDVq1Eg2SF+7dq3Qa9IFItp8x5eqL7hJtktu07FngSQSF3+AOXz4MBkZGREA0tfXpxUrViilul9cXByNGjVKNnMDQBUrVqTZs2eTh4cHubu7k66uLh07dox12bdu3aJx48ZRkyZN8jRxZd8GDRpEHz58YEWuQCCgnTt35tqGDx8uJ7NatWo0c+ZMun//PmsZXD09PenRo0cUFyefSG/RokU0cuTIElfvyw+JRELu7u5KnZVnERoaKpdwryikZYqo8cq7MmVRZ8ntInkI5iQ2NrbY17BBScynvKIoIo8fP6aKFSvS48ePFbpeUUQiEbVu3Zq6d+9epPTGz/1jqO2GB2Th4EzTTntQZAnSfJ89e5Y0NTVp1qxZFB0drXA7xUUkEpGZmRkZGhrS+PHj6cGDB3LrMJ6enuTr68uJ7NWrV8spp65du9LcuXNp3bp1sv0DBw5ktfwqkWKlUC0tLYv04FAYY8aMkcuM2q1bN5o3b15WRC6pqqrS6NGjWXfDfPv2LQGgKlWqkIODQ65IY2dnZ07SbUskEmrTpg1ZWVnRqVOniq2oXgfFys0ujj4NLPK1Fy9epDJlynCSULAg4uLiqE6dOrRkyRKFZqW8oigGiixas/G0lJCQUOg/NyVDSEuuvCcLB2dqs/4BPf5c8oFdLBZTSEhIidtRhPfv38uZt5RFQEAAubq65nrqO3jwIPXv35+8vLw4kSsQCGjv3r25tpEjR8oph1q1atHcuXPJzc2NhMLiP83mxdmzZ2ns2LFkZ2dH2tra+SomVVVVGjt2LGvprj08PKhnz56ydS8AZGtrS2vXrqWgoCBavXo1WVpakru7OyvyskhNTZVLuFevXj1ydnYu1gAqkUho4eV3cgojIa3wGd727dtln3fYsGGc153I4t27d1SuXDkCQKNHjy72bJRXFBzy5s0bWrduHedynvpFU+v196mqozMtu/aBUjPZGUB4/kMZJrecZK1RdOjQgbZu3cqZCSg7IpGIfH196cKFC1S/fv08FYaamhqNHz+eAgMDWZEZHR1Ne/fupdatW8vJMTU1lclbv3496yaqN2/eUJcuXWTyWrduLWcxCA8Pp/j4+ALbCI1Pk1MWB938C5X78uVLql69uswpxM3NraQfpUj4+fmRtbU1AaAuXboU657mFQUHCAQCWr58OampqdHTp085k5MhFNHK6x/JwsGZ2m54QO4BpWP75OEGoVBY6EDFFa6urqShoUHm5ubUqlUr+uuvv8jBwYF2795N169fJy8vL06ehoOCgmjt2rVUt27dXAqqc+fOrFaXy+Lhw4eyam8AqEePHuTl5UXPnj2jxo0bF2lNIedid7qgYJf1lJQUWeI9FRUVWrBggVKqBkZFRck+a8OGDSksLKxI1/GKgmW8vb3Jzs5ONqXlykvlS2QyddvmJisBmZbJfk58nt+X+Ph41kxbxUUikdCSJUvynM1UqFCBXFxcOJF55coVqlOnjsybLfuAWhQ35fCEdDllURQvwytXrpCxsTHhe6ZnrtbespOamkp9+vSRzWiKUlCtIEXBR2YXA4lEgu3bt2P+/PnIzMwEAOzZswdTpkxhVQ4R4eyrb1h+/SO01VWxcWB9dLSpwKoMHp7SJDQ0FCdOnMg3+j41NRVTp07F2LFjwTAMq7LFYjFOnTqFpUuXIigoSLa/fv36uHfvHkxMTAq8nogw4fhr3POJAgB0sqmAgyPzDmjOIjw8HGPGjIGLiwt0dHSwdetWTJgwAQzD4MmTJ2jdunWJP1dOxGIxpk2bhr1798LIyAjXrl0rUA4fmc0CAQEB1LZtW7knHz09Pdbt2qmZQppxxoMsHJxp2MEXFJGHR5NIJKKjR4/SrVu3WJXNw/O7kJycTKNGjco1m6lbt26+NaVz8tRPvspeYW60YrGYtm/fTpqamgSA+vTpQ5GRkWRkZMRJRTsi6Swqy5uvMLd/8KanknP58mU5F0N8D4pjk4DoFOq85RFVdXSm7fc+kzhHXETW1NnGxoZq1qxZamYDHp6fnWfPntHs2bOpXbt2ZGBgIPe7rlOnTpHXSdIFIjll4REcV+g179+/p3r16hEAKlu2rMw8xGWE/smTJ0ldXb3AQGJeUbCARCKhYcOGyRamAJCnpydr7d/5GEG2S25T/eUu5Oqb+4nm0aNH1KJFC9nNfPnyZdZk8/D8zmRF5589e5bmzp1LHTp0oJYtWxYrZcvUU29kymLzncLXIDIyMmjcuHFyCqpLly6c1ObO4v79+zKl+M8//8g8zLIW8nlFwQKbNm0iQJp/6fz589SiRQtW2pVIJLTrwReycHCmXjsf07c4+cjKt2/fUrdu3eRuqFatWpVqLprfMQ8Oz++FRCKh5OTkYl1z6324TFn02fWkwHMPHz6cayYDgBYvXlySbheKl5cXVapUSZZ1ID09nTp37kwxMTG8oigpLi4upKKiQgYGBvTp0yeSSCT04sWLErebKRTTbCdPWYR1Tne7R48eyf6p2Tcu3XHzIy4ujg4cOEDTpk0rlSA5Hp6fgZwxFwJR3nEhIpGIHj58SFOmTKHy5cvL/b7ZiMQviG/fvpGtrS0BkHl9jR07llcUJcHPz4+MjIyIYRhWU2zHpWTSoH3PyMLBmbbe9c33Kd3NzU2u6Hr//v1Z60NhZGZm0pUrV2jAgAGkoaFBurq69OnTJ6XJ5+H5Gcm5bhGdXPCDlVAopPv379OkSZPIxMSEypQpw3ngZXx8fC7nHF5RKEhSUpLM53r16tWstEkkfepov+khVV9wk668zT99hpubG+nq6hIAaty4MamqqnI+UEskEnr69ClNmTJFttCWtZ08eZJT2Tw8vxLtNz2UKYuipv0XCoV09+5d2rVrF2eJFAUCATk6OlKZMmV4RVFSxGIx9evXjwDQn3/+yZpd/ktkMrVYc49sl9ymFwWkMXZ1dZUpib1795KXlxdNmTKFlT4UxJMnT6h27dq5zF3jxo3jXHZepKWl0Zs3bygoKKhU5PPwlIR/zr2VKQvPr/Gl3R0ZMTExtGLFCjIxMeEVRUlYvny5zK86JaX4RYDywutbPDVccYcar7xD70MS8j3v4cOHsvTX+/fvJyLpk74yUj1IJBLavn27nJKwtbXlvPpbRkYGeXl50enTp2nhwoXUp08fsra2JoZhaNiwYUpJfcDDwwWbXT7JlMUDn8jS7o4cqamptHv3brKysuIVRXG5cuWKzMeZrSyar4PiqM6S29Ry7X0KKGAaeu/ePdLW1iaGYejff/9lRXZRiY+Pp/79+xMA0tHRISMjI9LR0SlS+H9JuH79er4ZTf/+++9SqWXAw8Mmp14Ey5TF3Y/s57IqKUKhsEBFoZZvPLcSYBimK4DtAFQB/EtE63IcbwfgKoDA77suEdEKLvvk7e2N4cOHQ1VVFU5OTrC0tCxxmx5f4zHq8EuY6GngzMTmMCujned5d+/eRe/evZGZmYlDhw5hzJgxJZZdVNzd3TFkyBAEBQXB1tYW586dw7Jly9CjRw/Url2bU9mdO3fGqFGjsG/fPrn9CxcuxMqVK1lP4QAAvr6+cHBwgI6ODvT19aGnp5fr1dzcHG3atGFVvlAoxKtXr9CsWTOoqqqy1i7Pj83QZlVgqKOOqac8MP74axwaZYc/av84aXnU1ApRBflpEK43SJWDPwArABoAvADY5DinHQDn4rat6IwiLi5OlqKXrTKob7/Gk+2S22S/4QGFJaTle97t27dJU1OTGIaho0ePsiK7KIjFYtq4caMsf/6ECRNkZiau6jJkkZmZSfv27SMLC4tcM4lNmzZxKpuIaM6cOXnOYvB9TSYykhszQZ8+fcjY2JhGjBhBTk5OSi1H6unpSRs3bix25TcednD2CvthzVD4EU1PAFoAcMn2fj6A+TnOUaqiyCp0MmrUKFYWrz+EJpDt0tvUZv0DCo3PX0m8f/+eNDU1SUVFhU6cOFFiucUhKw2ynp6eUityBQUFyepiMwxDf/31Fy1YsIBUVFTo0KFDnMretWsXVa1aNU8F0aRJk+IX0ZFIiNLiiKK/EH19SeT3gMjHmcjLiejVYSL3A0Qv9pN4f3uKWmhGwXPL06sJuvRwlA5dHaJNR/rq0ImxtejuhpEU+uQ0UXwwEcvmNldXVzp16hR9+fKF1NTUqFatWnTnzh1WZeSFRCKha9euyaWb8fcvvJ4DG/j5+ZVajZGvX7/me/yaZ6hMWbwKZLdsQGxsrMJjV0GKotSyxzIMMxBAVyIa//39CADNiOh/2c5pB+AigBAAYQDmENHHfNqbCGAiAFSpUqVxcHBwsfv0+fNnrFy5EgcPHoSWllae53h6eqJBgwaFtvU1Ng0D9j2DmgqDC1NaopJh3uYmQJqVdvr06WjZsiWGDh1a7H6XhGfPnuGff/7ByZMnUb16daXJlUgksLOzQ506dbBw4ULUqlUL586dg6qqKgYOHMip7O3bt2Pu3Llo2rQpfHx8EBcXh3LlymHt2rUYM2YMVFRU5C8QpgOxfkCsP5D4DUj4BiR8BRJDgNQoIC0WkIg47TMAoIw5UK0DYP0HYNUO0CpT5EsHDx4MJyenXPv79++PzZs3o2rVquz1Mxt3795F586dUa1aNSxcuBDDhw/HuHHj0LhxY8yYMYMTmQCQkpICOzs7EBHOnTtXpN8sWyxbtgxbtmzB7du30bJlyzzPueQRgllOXgCAB7PbwqqcXonlvn79Gr1798bcuXPxzz//FPv6grLHlqai+BNAlxyKoikRTct2jgEACRGlMAzTHcB2Iip0NOMqzfjz58+xbt06XL16tcDzYlIyMXDvM8SnCXFxSgtYl9dnvS9sQkScrAMURmZmJjQ1NWXvxWKxUuz2ycnJUFVVhaamJkxMTDBy5EgsX74chmXKAPGBQJgnEO4JRH0CYnyB+GBIJxzf0TSQDtqG5oBeeUDHBNA1kb7qlAU0dAF1nf9eVTUARgXijGR4erzCnVs3cOncCehpMDAz0kXbJnXQwsYCNUy1oZUeBcQFAEmhxftQtgOBRiOAqvZATkUHwNXVFbdv38bVq1fx6dMnuWNaWlpwdHTEvHnzoK2d/wONIvj4+GDx4sW4ePEiAMDS0hLa2trw9vbG0qVLsXTpUk7uvZSUFEyZMgUnT56EpqYmtm7dismTJyvlPt+1axemTZsGPT093Lp1K9/U3jvuf8GWu58BAG8WdYSxnmae5xUVHx8fNGvWDKmpqXBxcUHHjh2Ldf0PmWYcRTA95XFNEACTwtrmItfThw8fyMjIiP73v/8VeF5qppB67XxMNRfdpNdBhWeS5Ck94iJDyP/eYaIHq4mO9iJaa0601EC6rTAh2tOSyGk00cO1RO8vEoV5EaXFl0imUCikrl270syZM+n+/ftFXyuQSIiifIme7yU6+ed//cxvW1GO6P4qotT/TBtxcXFUq1atfNdlatSoQY8ePSrR58uPd+/e0Z9//kkMw8jJnD59OmdebRKJhA4fPizzqBs4cCDFx8eTWCymzZs3cyIzi927dxMA0tXVLfA7/efsf3EWmcKSfw/Xrl0jhmGobNmyxTbx4Qc1PakB+AzgDwChAF4BGErZTEsMw5gCiCQiYhimKYALACyokE6zPaMICgpCq1atEBYWht27d2Pq1Kl5nieREP4+7QGXjxE4MMKOLzb0oyEWAaGvgS93gcBHQNhbqcmIUQEq2AKVGgMVGwBmDYDyNoCaButdyLp1WX2yTYsDPl4G3hwFIt7lfY6qJsLqTII71YWuviH09fVzeXlln91xQVhYGLp37w4vLy+5/aNGjcK///5buOeNgnh7e2PQoEH4+PEjLC0tsXbtWvz111+4f/8+2rdvz4lMANi/fz8mT54MHR0d3LhxA+3atcvzvK7b3PApIhkAELSuR4nlrlq1CosXL4atrS2eP38OPb2imbV+yBnF9x9Md0iVhT+Ahd/3TQYw+fvf/wPwEVKPqBcAWhalXTZnFJGRkbIi6QDo4cOH+Z679a5vkYuvF5ewsDA+6EwR0hOIvM4ROY0iWltF+rS9zIjo305Ed5cRfb5LlK48ryOlERdI5Dwr/xnHhfFEScrz55dIJLRjxw7q27cvVa9ePdfMom/fvpSenrtIF1ukpqbS+PHj5WSamppyUp87OwcOHCAApK2tTffv38/zHIlEIptV/HXgeYllSiQSGjBgAAGgAQMGFHlxGz+i1xOXG1uKIjExkRo1aiR3c+V3Y918J3V7m+3kyXoabl9fX+rSpQurbf7SpCcQvT1NdGqw1IS01IBoY3Wiy1OJPlwqsfnopyT6M9Hpv/JWGk6jpF5bSiQtLY08PDzo+PHj5ODgQD179qTx48ezlgUhLzIyMqh79+5yv+dOnTpxHtD577//EsMwpK2tTffu3SMiymWOyp5IkI0HzeTkZKpbty4BoFWrVhXpGl5RKEB6ejq1b99e7qYyMjLKUwn4RSVT7cW3qN/uJ5QhZLfwyMuXL8nExEQpeZ5+asRiIv+H0iflleWlA+BmG6Jb84m+urPubvpTI0gnerAmb6Xhtkm6HlJKcFXrRCwWk6Ojoyw1TvatqANpSThy5AgxDENaWlrk4uJCtWrVymWdCIhOkSmLL5HFq4WRF/7+/lS2bFliGIauX79e6Pm8olCAoKAgunLlCtnb28tuqFatWuU6L10goi5bH1HDFXcoPIHdqfOtW7dkN3ZBtW5/a1KiiR5tINpiKx3o1pgTXf+HVw5FRZhJdHdpboWxrZ7UfPWLkZycTMePH6fOnTvLKlWqqKhwtoifnWPHjhHDMLKyAU2aNMmlGM++/C/Vh0hccqV59+5dUlFRIX19ffLx8SnwXF5RKEhYWBhpamqSsbExTZ8+ncaPH5/rnIWX33ESZXn8+HFZtDQAio6OZrX9opKSkvJj1uYOf090ZarUu2epgdRr6d15IkH+gY08hZAcRXSoa26l8eVeafeME0JDQ2nz5s3UoEEDqlixIkVF5S5BzCZBQUHUsWNHudnMuXPncp3X4Xt6crtVd1mRu2XLFplXW1Zi0bxmbryiUJDZs2fLpqZCoTBXZbms0oerb7CbNG/jxo1yN1ODBg1Ybb+oPHnyhIYPH14qsvPlqzvRiQHSAWyVKdH1mUSRBT8p8SjAm+O5Fcab46XdK854//69bP2AC9LS0mjGjBlyRcgAkJWVVS4X6eyL2ydfBJVYtkQioREjRhAA6t69O4lEIlqxYkUuZcErCgWIjo4mHR0dKlOmDCUkJOQ6HpuSSY1W3KGeOx6z4v+cxbp162RT4qxt1qxZrLVfFNLT02nevHlKzztVIEHPiI71lg5Y6y2JHm2UixHg4YhQj9wK4x1vBlWUwMBAGjlypJzX144dO3Kd9zU2VaYskjNKPqNPS0sjOzs7AkCjR48mFRUVcnV1lTuHVxQKsHDhQgJAixYtyvP4/057kPWCG/QpnP1cMk+ePJFTFGyWYC0MDw8PWVU/IyMjSksrZVNOdm+dDdZET3cQZXLnGcOTD/HBRMvLyiuMwMel3auflnfv3lHv3r0JAJmYmOT5MLro8nuZsmCDb9++yRUqGjx4sNxxXlEUk/j4eDIwMCBdXd081wayTE477n0ukZy8EAqF1KBBA9mTRtmyZZWS2EwoFNLKlSvl1kX++ecfzuXmR1JkMJHzbGnMw+pKUm+cTG6LJ/EUgRi/3DOMlPwrNfIUzNOnT6lNmza0YMGCPI+zZYJ68+YNVatWTe4BVF1dXc7dn1cUxWTFihUEgObOnZvrWHKGkJquvkvdtrmRQMS+V82OHTsIAHXt2pUkEgk9efKEdRk58fHxoSZNmuRyG/T19eVcdnYkEgndu3uXto5pRhkrKkqVhPMs6SIrz4/FV3d5ZXF5aqm61f7MSCQScnV1JZEot2t9dhNUSccbDw8PsrW1lfuNr1mzRnacVxTFICkpicqWLUtaWloUHh6e6/jamz5k4eBMHsHsByiFh4eTgYEBaWpq0pcvX1hvPz9CQkLo1KlTclXm/vjjD6XJT0xMpJ07d1KnxtZ0d4SOdOA5+IfUs4nnx+bRBnmFwTsWsM6ow+5k4eBMfXaV/KExIyNDtv4IgCwsLGQKqiBFkTvN5G/Ovn37EBcXhwkTJsDU1FTuWGBMKg49CcDAxpXRsIoR67Lnzp2LpKQkODo6wtramvX286NSpUp49uwZ0tPTYWxsDACYMmUK53K9vb3x999/o1KlSvA4NAuXOkeiaSU1JLddAYy9A5jact4HnhJiPxdYEA6ofc86u6cZcGli6fbpF+PwqCYAAM9vCfgam1aitjQ1NbF+/Xo8fvwY1apVQ3BwMFxcXAq9rtSSAnKJokkB09PTUbVqVcTHx8Pf3x/m5uZyx8cdfQX3wDg8mNMW5fXzrlehKK6urmjfvj2srKzw4cMH1tM9F8ShQ4cwfvx4mJub4/Hjx+jduzdev34NdXV1zmTu378ff//9N8poSHCgpxYG2KjjQaAIGV23oftfEziTm0VKSgrCw8MRFhaGsLAwhIeHQ0NDA5MnT+YsOd0vz4eLwIWx/713CAK02X+g+h25+CYEs89LkymykTgQkP4G5s2bh5CQEFy7du3HTQrI1aao6SlrfWDChAm5jr0KjCULB2fa89BPobYLQiAQkI2NDQGgGzdusN5+QTx//pw0NDRIS0uLXr9+TURS7whlsHfpVAqcoUeZi/RpdgsNGjJ4EOcyk5OTc6VmAUD29vacB1xJJBLauXMnvXr1ivP8Qnmxb98+ThPvEZHUIy27Kerzf0FjMTH8ondJyFqreMlyVbx79+5RSkoKv0ZRFCQSCdWpU4dUVVXJ39+fXr16RbGx//1Dhux/To1X3qW0THZzOREROTk5yTJoKhOJREINGzYkAHTq1CmlyvY+t5xSF+jTt3/0qYWFJpUvX15p0efLli2TUxJTpkxRSg3pzMxMatGiBQEgMzMzGj9+PF29elVWo5xLXrx4QY0bN6YmTZpQSEgI5/LIebZMWXzaNZiEQiE1b96ce0WVDX9/f1kksjKRSCQUGBjIertZD6sFuctGRUUp/BDCK4oiEh8fT1euXCEior59+8qesJ9+iSYLB2c6/CRAoXYLQyKR0JUrVygoKIiT9gvi27dvtHv3buUJlEikBXWWGpD/4pr09skdatu2LV28eFFpXbhx4wbVrl2b1NTUaN++fUqTm5aWRhYWFjIFZW5uTlOnTs0V+MQF3bp1k0uv/ezZM85lUsCj/5TFrAoEgBwcHLiXS0Senp5kZGRE3bp1y9ObiCtEIhF17NiRypcvz4mSylIU931yZ7E+cOAAaWtry8aw4sIrimLi5eVFAMjJyYmIiP7c94yar7lH6QLl3XC/JCIh0dVp0sHjylQikbS+xosXL5TelYULF5Kbm5tSZUokEpo4cSKtWrWKvLy8OMuUmpO0tDTq16+f3CxKQ0ODDh8+zLnswHfP5ExRKioqSvl/Z49Enj9/vmy/MmaOM2bMIAD0999/s962d1hivrMKNzc3AkB2dnYK3VsFKQp+MTsPhgwZgnPnzmH9+vXoOnQieu96ikU9amN8GysWe/mbIRIAF8YAn5yBNnOADouAUqjTLeuOSPTbLFonJCTg3bt3ssX77K8DBw7krJa0QCDA2rVrcXDvLoRMFsj225wzg8dbT2hpsesQkpNv377Bzs4OUVFROHfuHCpVqgRvb29MmMCts0RiYiJq1aqFqKgovHr1Co0aNWK1/aqONwAAztNaw7ZSGbljbdu2hZubG27fvo0uXboUq11+MbsYfPr0SeZjPGnSJJp+xoPqLLlNSenKqS4nkUiUkvJYqYiEROdGSJ8qnynRzMVT6nz58oVmzZpFRoZl5GYWDvNyB7NywePHj0ldXZ10dHSoadOm1LBhQ6XM5E6dOkUAqGnTpqw7Ljz8FJnvrOLOnTuykgjF/Zzg4yiKzpo1a6Q2OQC+36Jw4104Bjcxh74Wd66i2bl69SquXLmiFFlKQSIBrk4FvK8CnVcDLfKuN87za2JtbY3NmzfjW0goDlbcKNu/Tns/Xrx4wbn8Vq1aYe7cuUhLS8PLly/x9u1bvHz5knO5f/31F9q3b4+XL1/i33//ZbXtdjXLy/6OSxXIHevYsSOaNm2Kp0+fws3NjTWZvKLIRkBAAE6dOiV7H6xiBjERRresqhT5IpEICxYsQEREhFLkKYX7y4B356Smppb/K+3e8JQSurq6mDBxImhJvGyf6fkeSE9P50ymSCTC5MmTsWbNGrn9+/bt40xmFgzDYPfu3VBTU8P8+fMRExMDiUSCmJgYVtqf2bG69PWcZy65ixYtAgCsWrWKFVkAryjkWL9+PcRi8fd3DASVGqGFlTHMy+ooRf7x48fh4+Pz6yiKt6eAp9sBu3HSdQme3x5GRQVYFA0AqKovguDKdM5kqampYd++fTh+/DjKlPnPln/27FnExcVxJjeL2rVrY/bs2YiLi4OjoyP27duHs2fPstL2tA5SReH2OTrXsZ49e6J+/fq4d+8ea7M2XlF8RyAQoH79+ti7dy8AoGmPv6BmaIr2VZUTIZ2eno6lS5cCACIjI5UiMyfu7u7sNfbtJXB9BmDVDui2vlQXrnl+MNQ0gLkBAIAyvk6A333ORDEMgxEjRuD9+/fo1KkTACAjIwPHjh3jTGZ2Fi9eDHNzcxw6dAhz587FgwcPWGlXVYWBhpp0+H7mJz9LYRgGCxcuBACsXr2aFXm8oviOhoYGpk6dKpsKl2/WG7oaKuhYy0Qp8nfv3o2QkBAAKJUZRVxcHGbMmMFOYxmJwIVxgEFF4M+jgKpy1nd4fiJ0jYFR16V/n+wvvWc4xNzcHC4uLti9ezd0dHSwb98+2VokVwQGBmL06NGy92lpaXB1dc1mtSgZZyY0AwCMOpJ7zWXAgAGoXbs2nJ2d8fbt2xLL4hVFDry9vQFGBYECPXSyMYWleSXOZSYkJMjZUePi4iAQCAq4gn0cHR3x9evXkjdEBDjPApJCgQGH+Fw/PPljaQ/Yfc8Nta4K5+IYhsHUqVPh5eUFExMTPHz4kFN5lpaWGDx4sJyFID4+Hl5eXqy039iiLABAKM6t8FRUVLBgwQIA/80qoqKiFJbFK4oc+Pj4QLNSbaQIgU42poVfwAIbNmxAfHy83L6S/FOLy9OnT3Hw4EFkZmaWvDGf68CHC0C7+YB5k5K3x/Nr03Prf38/26UUkdbW1nBzc5NlSuaSgQMH4vbt2zAwMJDtu3+fPVObiZ4mAMA9IDbXsSFDhsDKygoXL17Ew4cP0bt3b4Xl8IoiG0QEb29vlK3bFhqqKmhbs5xSZDZq1AhXr14FADRu3BgVK1ZUmvlJKBRi8uTJAFByRSFIBW7PByrYAq3/KfJlkZGRpbYuw/MDMMtH+npnISDkzgsqO6qqqqhfv75SZLVv3x5ubm6ysgVsrVMAwKY/6wEA5l54J7dfIpFAJBJh7ty5AIDOnTvD3d0dSUlJCsnhFUU2oqKiEB8fDx3rpmhmVRZ6mtxH7jIMg4EDByI8PBwAMGjQIHh4eKBs2bKcywaArVu34sOHDwBYUBRum4CkEKDHZkC1aN+dUCjEoEGDoKqqWjLZPD8vBhWBBsOlfx9oV6pd4Yr69evj2bNnqF69Oh4/fsyaaTkrpuJrnHydColEgr/++ktWV0YkEgEAPn36pJAcXlFkw8fHByo6hhDpmKCVtXIWsbN48+YNAOmMokKFCrCy4j5dSFBQEJYtWyZ7LxKJFF9oS44AXuwB6g0BqjQv8mWzZs3Cs2fPlKYYeX5Q+nw3O0V/AhJDS7cvHGFpaYmnT5/CxsaGk6C/DOF/v101NTWcOXMGbdu2lTvHx8dHobZ5RZENb29vaFW2AQA0qarcgSsrNxXbeWHyg4jwv//9DxKJRG6/wrOKpzsAsRBo51DkS44cOYJdu3ahQoUKUFHhb8XfGoYBemyR/r0r73RDvwLlypXDgwcPWL3fe9QzAwCcdpd3RtHS0sLVq1fRoEED2T5vb2+FZPC/zmz4+PhA07wO1FWAujmSbXFJRkYG3r9/D2traxgZKcdLSCAQYPPmzTh48CAAwN7eHlpaWoopitRY4PVhoN5goGzRZkIvX76UrY3kLDmrLBISEpTqNMBTCE3GSV+FaZy7y5Ymenp6aNmyJWvtTf8efLf9/pdcx8qUKYPbt2/LSivzMwoW8Pb2hqZZTdiY6sqCWZTB+/fvIRKJ0LhxY6XJ1NTURM2aNWUeGIsXL8aDBw8UyyLqdQYQpQOtihZlGxERgf79+8vstBUqVCi+zBJy9+5ddO7cGYaGhkqXzVMAWU4QZ4aWbj9+Imqa6gMAEtOFeR6vUKEC7ty5AzMzM15RKEr2oBsfn09QL2eBBhbKXZ/IMjvZ2Sl3yk1EuHv3LrS0tNC6dWu0aNGi+AMnEeBxDDBvBpSvXejpAoEAQ4cORUJCgmyfMmcUqamp+Pvvv9G5c2d07twZGhoaSpPNUwQ6LJG+Bj8p3X78YlhaWuLOnTuIj49HRkZGsa//7RXFyZMnIRKJkJqaCpGmAVQ0tFHT1KDwC1kkPDwcqqqqSlcU4eHhEAgEaNOmjeK1AUI9gJjPQMMRRTpdTU0N9+/fly2i16hRQ+mK4uPHj1BVVcXEiROVJhcAxGIxNm3aJPMyUyZhYWFwcnJSulyRSIRPnz7h2rVr8PX1LfwCFRVAXVf6d9DTEskODQ1FSkpKidpQBIFAIPNiVBYq3w0B3v5fc607ZmFra4tr167JMkAUi/zyj//MW3HqUQwfPpxOnjxJRER3P0aQhYMzvQ6KK/L1bJGamkoCgXJqXmRHLBaXrFb1/ZVEywyJ0or/nX39+pXCw8Pp2rVristXgPj4eHJ2zr/uMFeIRCI6dOgQffv2TalyfX19qUqVKqSnp0dJSUlKkZmamkrp6enUrl07MjIyIm1tbXrw4EHRLva7L61bsammwvLXrFlDAOjMmTMKt6EI3759Iw0NDerUqZNS5c4650kWDs6k36gnBQQoVrIZfD2K/AkICMDq1ashkUgQ/N0X2cpEV+n90NHRgbq68nMiqaiowMSkBKa2L3ekZicFUnWYm5vD1NQUvXr1Uly+AhgaGqJHjx5KlQlIg7zGjh2LypUrK1WumZkZOnfujJSUFLk0+lySkJAAKysrPHr0CPHx8UhPTy96EGm1DtLXZMWfyrM8fe7evatwG4pQqVIl6Ojo4M2bN5znkspOu+/BwVpV6im8DlEQvKIICICPjw8uX76MiMR0aKmrwFBHOQO2Mm8kTkiPB8K9AOs/SrsnPAWgr6+PgwcPwtnZGVevXlXKfVexYkVcuHBBrtysQtkGBGmFn5MH9vb20NDQwN27dyESiXDp0iWF2ikuDMOgcePGiIuLQ1BQkFJkAkBjC+mDmmbl2gq7wBbEb60o0tLSZDfvqlWrEJ6YAbMy2pzUD85JRkYGLl++zLkcTgn/njagkvK8tXgUp0ePHjh58iSnxYKy07JlS+za9V/+pmIpitrfZ5nvzhVbrp+fHyIiItCiRQt8+/YN9vb2OHeu+O0oStZaY5aTijIwKyNdY1TVNeJnFGyTXeN7enrCOygMFQw0lSI7MjISCxYskIXW/5SEf8+CaaqcnDk8JcfY2Bg6OsopxAUAEydOxKRJkwAUs85KVkqPDxeLLVNTUxP169fHo0ePAADPnz9X3FlDAbLc3LOyLSiD7A+3/IyCZQICAuTefw2PhqG2csxOERER8PX1xdGjR5UijxNivwC65aS1BXh48mHHjh1o2bJl8WYUlvbS16DHxZZnbm6OFStWyO3T1FTOAyBQOjOK7Pj4+LBuXixVRcEwTFeGYXwZhvFjGMYxj+MMwzA7vh9/xzAMq/ktsisKNTU1ZAgJibG5SwtyQdaPZtmyZUozBbBOShSgXzpR1Tw/DxoaGrhw4ULxEj9qlGzWM336dLnssMpUFFWrVkXZsmXx5s0bSCQSpS+oJyYmsp59utQUBcMwqgB2A+gGwAbAXwzD2OQ4rRuA6t+3iQD2stmH4OBgLFmyBLq6uqhQoQLKmVWCnrZyArCypuGhoaFydlxlcOXKlXx9rYtFcgSgxysKnsIxMzPDtm3blCZPTU0N+/fvl5lklKUoXrx4gRs3bqB27dpISEhAhw4dcOjQIaXIzg7b5qfSnFE0BeBHRAFEJABwFkCfHOf0AXD8u5vvCwCGDMOYsdWB+fPnY/ny5ahVqxZCQ0ORKSZUqcRa8wWSXeOvXbtWLlKZa1xdXXHgwIGSNyRMAzSU70rM83NSvXp1pcpr1qyZLJ+YstYoqlWrhqFDh+LpU2mw4KNHj0ollxnbC9qlqSgqAfiW7X3I933FPQcAwDDMRIZhXjMM8zo6umjmo6z4gdq1paknSCKBREkuqznLI65fv14pcgFphlhHR0d2igUpwUOMh0dR1qxZgwoVKihtRlGuXDnMmTNHbl9p5DL7lWYUeY0wOUfpopwj3Ul0gIjsiMiuXLniVabLUhQgMUR51J/lguwzCj09PezYsQNhYWFKkZ2ZmYnExETMmjWr5I397LEgPL80hoaG2Lp1q1LXKGbNmoXy5cvL3pfWjEIgELDmeVWaiiIEgHm295UB5Bwpi3JOibGxkS6NkEgIERu2+yKQmJiIefPmAQA6duyIwMBAxYsGFZOsVOKnT5/GvXv3FG9I06BY6aCJCHfu3FFcHg+PAgwZMgRdu3ZVmjw9PT0sXrxY9l7ZMwptbW14eXmhVatWePy4+F5jeVGaiuIVgOoMw1gyDKMBYAiAaznOuQZg5Hfvp+YAEomI9WxbWTMKsTAT6QLlDNZHjhzBqlWroKWlhdevX6N8+fIwNzcv/EIWyJ49curUqQplkwQg9XhKKbr5imEYzJ49u2jJ4Xh4WIJhGNSrV0+pMidOnCirUqnsGUV6ejri4+Px+vVr1iplFqooGIb5H8MwrFfTISIRgP8BcAHgA8CJiD4yDDOZYZjJ30+7CSAAgB+AgwCmst0PQLoApa6uDmFyHOLS8s7pzjbm5uZQV1dHgwYNEBISws56QRHJXpzoy5cvWLdunWIN6ZsCSaHFMj9pa2ujb9++Chd55/lNSI2RvqopL1COTTQ0NLBy5UoApVeYC4DyFAUAUwCvGIZx+h73wNrqJRHdJKIaRFSNiFZ/37ePiPZ9/5uI6O/vx+sSEScRLGpqaqhRowbSE6IRnaTcmIas4BxlRnFmVxQVK1bE4cOHFXvKL1dLanpKKro1sEKFCvj06RNGjBjBjosuz6/J59vS1+qdSrcfJWDIkCFo2LAhirtmyiZVq1ZlpZ1CFQURLYI0juEQgNEAvjAMs4ZhmGqs9OAHoXbt2hCnJiAqUbEkZIqSFe6vzChOkUgkW8iuV68egoODYWFhUfyGzL4HNGWl8igCWU9X165dkz1xKYPo6Gi8e/dOafJ4SsjrI9LXBsNKtx8lQEVFBfv27VNKVujskdhNmzYFAJQvXx56enqstF+kNYrvucojvm8iAEYALjAMs4GVXvwA1K5dG6KUWMSniyEQKe9JtzTC/Xfv3o1NmzbB1NQUjx49gkAgUMzPvIItwKgCIS+LfEn2afiyZctw7VrOZSluMDExwYQJE7Bu3bqfO7/W70Lo999D9S6l248SkjVoc01QrPQB19RAS7aQzpbZCSjaGsV0hmHeANgA4CmAukQ0BUBjAANY60kpY2NjA1FcKCRg8C1eebOKWrVqyfLXKwsbGxswDIOOHTsiPT0dz549U6whDR2gSgvgS9FTFGRXFGZmZtixY4dS1mcYhsHYsWMxf/58tGnTBv7+/pzL/FEojXT2JfLgy26SVPmt09EVGVffKABA+1rl0KNHD9SvX1+5igKACYD+RNSFiM4TkRAAiEgCoCdrPeEYobDgReo6deqgkoF0ihgYnaqMLgGQro+0b98etWvXRlqacs1eHTt2RO3atZGaWoLPW6MLEPkBSPhW+LmQrlHUq1cPlpaWGDt2LO7du6c098Fhw4ahQoUKqFy5Mmu226JARLhy5Yosi6oyEYvFcHBwwMePH5Uqd+PGjfD09ISTkxMEAkHxLn51UPpatY3C8h88eIAvX74oXUlGRETA1dVV6c4aD32lQcbV9aSz5UWLFsHS0pI9AfmVvvuZt7xKoT5//pycnJxy7c9OQqqALBycaf8jvwLP+1WQSCQlbyT6i7Rs5eOtRTr906dPFBISQkKhsOSyFeDdu3ckFouVLlcgENCrV6+ULvf69esEgBo2bEiRkZFKkztixAjS1tYmAOTg4EBPnz4lkUhUtIuXGki3+GCFZEskEqpcuTJpaGhQamqqQm0oysGDBwkArVmzRqlyLRycycLBmYwqVCKJREJisZiePn1arDbAl0IFatSogeHDhxcYYFZGRx1mZbTwPvT3cN1kxYHNxBowbw54HC+Sm2zNmjVRqVIlucpnyqRu3bpQKQVzhrq6umw9Spns27cPFStWxJo1a+SihbkmIiJClhV5/fr1uHv3btGyx2Ym//e3YRWFZPv6+iIkJAStW7dWau0N4L+1xiwnFWVjY10VDMNARUUFLVu2ZK3d0vm1lgJly5aFkZER+vbti4cPH6JJkyZ5ntfA3BCe3+KV3LufnMajgSuTgUA3wKptafeG5ztBQUEoW7YsPnz4ACMj1kOhCiR7ipo+ffrIRSoXyNmh0tcW/yu2zE+fPsHDwwOxsbEAgE6dlO9am7XWqExFIZb894AmS0fEMr/NjAKAzB7frVs3fPr0Kc9zGpgb4ltcOmJTMvM8rgj0q+dDqtNXWsDoyZbS7glPNsqVK4fjx48rXUkA/ymKWrVq4fjx40WbxUnE0ocNAOhUfNdphmEwbNgwrF27FoD04fDYsWNKi9fJzMyEl5cXLC0tYWysvGJez/2lijEzwk+WjohtfjtFAQCxsbHo1KkTvn3LvQDbsIr0R/UmmL1ZxdmzZxESEsJaez8c6tpAqxlAgCsQ/Ly0e8PzHV3d0kkBLxKJEBMTAwMDA1y5cgUGBgZFu/DSBOmrZVuFvJ2yEv+Fh0uz/EyaNAmRkZFKMzV++PABQqFQ6WanfY+kHnzJb67xMwo2yK5tQ0JC0LlzZ8TExMidU9+8DLTVVfH4S0zOyxWmfPnyaNmyJSe1bH8Y7MZJZxX3V/AZZX9zoqOjQUQ4efIkatasWbSLBKn/1ccedl4huTkzxNra2uKff/5RqC1FyFqfUPZa1BM/6ViV+tGVn1GwQU5tq62tjeXLl8tNTTXVVNGymjHcvrBXErVly5aIiopC69atFY9Z+NHR0AE6LAK+PgPenSvt3vCUIhEREVixYgV69epV9It2fn8KbzIBUFMsJXhORbF//36lREW7ubnh7du3paIoROL/xi5dHW3OEov+loqicuXKAIDx48dj586duaam9jXKITg2DUEx7MRTaGtro2XLloiPj0fHjh1x/fp1Vtr94Wg4EqjcBHBZCKTzDgG/KxYWFli4cGHRL/B/CCR/TwrdfaPCcrNnF5gwYQKrXj8FkZGRgZYtW+LKlSsApCYoZRUiu+QRCgBID/JErVq12PFkzIPfSlGYmZlh+PDhuH//PlRUVLB+/fo8g4Ha1ZQm8brrzV7EcIcOHQBIUwD369cPhw8fZq3tgrh8+TLi45U0aKuoAD22SJWE8yzeBPWbUrZs2aKvC4hFwIm+0r/H3S1RxcSsGUW5cuUUz4isAKampsjIyJCZsWfOnMmZCSgn8y5K85fFP/iXU5m/laJgGAaHDx9GjRo1MGTIEHz9+hUnT57MdZ6FsS7qViqD6+/Yq5GUpSgAabTsuHHjsHr1as49ogwNDWFra4vbt29zKkeGWT2g/QLg4yXA85RyZPL8vGz8nmbCvBlgXrK8SKqqqlBVVcWWLVtQtmxZFjpXNHKmEW/VqhV69uQ+aUX2sUMYHcTZQjbwmykKADKb5YIFCwAAa9euzTNJXK/6ZngXkohAlsxPTZo0kcvk2KZNGxgYGKCo9b0VpV27djAwMEC3bt0wefJkpKSkcCoPAND6H2n6hZtzgSh2i7zz/EI8WPVfhcSxLqw02b17dwwbptyMs8bGxnLBhOvXr+fMBJSdTZelHoZlSRogzM8oOKBOnTro378//Pz84OTklOt4z3oVAQDXPNmZVairq8Pe3h61atWCqqoqPn36hOHDh3MeLcswDCZPltaB2r9/P+rXr89aecR8UVEF+h+Ulko9PQhI4VYZ8vyEBD4G3L6vR8z2LZHJKTt79uxRyiCdHVVVVVnNiV69eqFVq1ZKkbv7pdSkHHVN+j26u7tj+vTpnFgpfltFAUC24LZ69epcQTkVDbXRspoxnF5/k4t8LAldunTB2bNnMX36dERHRxdvwa8EjBw5Etra2gCAgIAAtG3bFnPnzlW8BGpRMDDDx3qLIE6KAJ39CxByKIvn5yLWHzj23TQz5Iy0UiJLZDmqKBtTU1MwDIPVq1crRV5i+n9JTkM/vQUgtY7Ur1+fE0X5WyuKRo0aoUePHvD29sbly5dzHR/R3AKhCel4+CmKFXlTpkxB/fr1sWzZMpiZmWHfvn1KqUNhZGSEIUOGyN4TEZycnLBp0yZO10hqdvgLi96YgAl5hYD1rREfo7xyrzw/KCnRwM5G0r/bLwJqdS/d/rCEqakpRowYgbp16ypF3v9OewAABJ7/eVCam5tjxIgRnMj7rRUF8N+sYtWqVbkGzY42FVDBQBMnXgSzIitrfcTAwABbtmwBEWHq1Kkly91fRKZMmSL3ftSoUVi0aBGn03Q1NTWMWHsR/7sthJXoC9ymW2PKxPHw8ip6RTyeX4iUaGCTtfTveoOBtnNLtz8sUqVKFaxYsUJp8rICgivHe8r2OTg4QENDgxN5v72iaNGiBf744w94enri5s2bcsfUVVXwV9MqePQ5Gn5R7C4CDx48GH/88QdevXqFf//9l9W286JJkyZo3LgxJk2aBBMTE6xcuRKXLl3iXK6NjQ3M+y3F9FsZ6FNDBR0TTqNZ4wawt7fHuXPnOFOSYrEYMTExv36erZ+FpLD/lIR1J6D/gdLtD8s4OjoqVk5YAY4+DQQAVDLURjUrac0JU1NTjB07ljOZv72iAKRFPgBg5cqVuQaW4c0toKWugr2u7FZEYxgGu3fvhrq6OubPn8+59xMAzJkzB2vWrMH58+ehqqqKkSNH4sOHD5zLnT17NtzRANNvZWCAjTpchuvAQINgYWFRtNTTCqCqqopLly7BwMAAFhYWaNGiBfr374///e9/WLt2LRISEjiRWxDJycmFn8QByqggWHAHPgJbvrtu2vQFhl8o1e5wAatFggph2XVpKqDLU1vKqtjNmTNHtg7JCfkVqviZt7wKFxWERCKhVq1aEQC6d+9eruPLr30kq/k36Gss+0VQ5s+fTwBozJgxrLddEDt37iQAZGVlRbGxsZzL8/b2Jk1NTRpiq0aZi/Tp4/8M6IXLec7lPn78mMqVK0cAZNvChQs5l0skva98fHxo48aNZG9vTxMnTlSK3CxEIhE5OjrSsmXL2ClSVURSUlLo9OnTREQkfnf+v0JE1/9RWh9+VS55fCMLB2eymn+DiKSFkoyNjSk5ObnEbaOAwkWlPqhzsRVXURAR3bp1i1q3bk0vXrzIdSw8IZ2qL7hJ8y+9K3a7hZGamkq2tra0adMmpf6YJRIJjR07lv744w+lKAoiog0bNpCKigr9u2A4JTjqk2CVOZG/K+dyg4ODqVGjRjJFMW3aNE7lSSQSOnfuHNWoUUNOQenp6XEqNzsJCQnUvXt3AkAGBgZ09uxZpcmeNm0adenShV4va/Ofknh1SGny09LSqFevXnn+lrnm2rVrtHTpUkpPT+ek/axKdsEx0ofW+/fv0+rVq0ksFtP//ve/ElVR5BVFEShskF54+R1Vm3+DAqJTit12YRS5RCTLZGRkKLUkqUgkoilTphARUcSHx0Q7mxAtMyR6sp2IYyWZmppKf/31Fw0ePJhiYmI4lZVFeHg4/fvvv9SnTx/S0dGhYcOGKUXup0+fqGbNmnJK6siRI0qR7erqSpqq+E9BLDWgoCcFlyBmE4lEQsOHDycA1K9fP6XJJZLOpKpUqUIAOFFSJ18EyRRFFhEREZSQkEDnz58nANS9e3eF2//tFEX9+vVp69atrA7AkUnpVHvxLZp84jVrbf6OyP1PMpKIzg6XDign+hMlhXMqWyKR0LNnzziVkR/p6elKkX3jxg0yMDAgAKSqqkqVKlUiOzs7mjJlCmVmZnIqOzk5mQY1qyynJHZtWqWUByE3NzciItq8eTMBoBo1alBCQgLncrPj6OhIAGjs2LGsty2RSGRKIiwhLdex+vXrE4AS3WO/naJo3Lgx9e7dm9q0aUP+/v4Kf3E52X7vM1k4ONOrQOWYan4LJBIi9wNEK8sTratK5H2ttHv005KWlkY3b94kDw8PCg8PV+5MVSwmvwXVZArC+S9tAkCWlpZ08eJFTkUHBweTiYkJXb58mVRUVMjAwIB8fHw4lZkTHx8fUldXJyMjI4qKimK9/fmX3pGFgzN13Oya69j169cJAP3xxx8lkvFbKopbt24RANLV1aX9+/ezYv9PzRRS09V3qffOxyQSK2894bcgypdo33eb9rkRRIlhpd0jnqIS4CY3i+hgqUpdu3al69evK0VZLVy4kAAQwzDEMAxdv36dc5nZkUgk1KFDBwJA+/fvZ739pHSBbDaRlin/fUokEmrWrBkBoIcPH5ZITkGK4pd1j+3cuTMsLS2RmpqKSZMmoUePHggLK1neJh0NNSzoXhteIYk45c5OEB7Pd8rVAMbdAzosBj67ALuaAO77pWmoeX5MMhKB1WaydBy+ccCsxNHY5+KDW7duoWfPnpy5P2chEAhkcUhEhDJlyuD+/fvw9PTkVG52zp07hwcPHqBp06YYP3486+3XW34HADC+tSW0NeS/z/v378Pd3R2tWrVC27ZtWZedxS+rKFRUVDBp0iTZ+1u3bsHW1hZnzpyRTqUUpHf9imhT3QQbbvsiIjF3/qL09HTs2LFDOVlafzXUNAD7OcDU54B5E+DWPGBfa+DzHb62xY8EEXB5MrCuCiBMAwBcNZ6KymvDsGXrdlSvXl1pXbl69apcnEhKSgoMDQ1Rp04dTuVKJBLs378fSUlJmDVrFhiGwZ49e1ivz33XO1J26y/qmTs77KpVq6THOM6yUOpmIi62LK+nqKgo0tDQkHl+mJiY0JgxY+jz588lmqIFxaRQjYU3adLxvBe2V61aReXLl6edO3dyvoD4yyKREH28QrS9gdSkcbQXUahHaffq90YiIXqwRs7MRHcWl2qX2rdvL/t929jY0OvXynE28fT0JE1NTRo4cCABoKlTp7IuQyT+bwH7Y2hiruNubm4EgOzs7FgxreN3XKPIYujQobIbSVdXl969YycWYvfDL2Th4ExX3obkOpaamkoVK1aULeadPHmSxGIxK3J/O4SZRM/3Eq2zkA5MJ/8k+qa4rziPgrhtklcQ++yJBNzEChQVHx8f2drE7NmzOYtdyIss7yoApKamRtevX6eAgABWZbRad58sHJyp964neR7v3LkzAaArV66wIu+3VhRubm5UqVIl2rJlC6uRyEKRmPrtfkK2S29TSHxaruMHDx6U82OvV68eOTs7KzWo7pciPYHo0Yb/FMbxvkR+9zmPv/itEQmJrkyVVxCbahGlxZV2z4iIaMaMGWRpaUmPHj1SuuwePXrI/b4NDAzowYMHrLV/zTNUNpvIa8xwd3cnAFS3bl3WHkJ/a0UhkUhkA/TYsWMJAHXq1ImVQLOgmBSyWXyLBu17lssLSigUko2NjdzNBIAmTJhAAoGgxLLz48iRI3To0CGl+5ArjYwkosdbiTZ8d8Xc2YTo5b9EmewHQv62xAURbasvryB2NCZKUU6gYlFIS0sjBwcHSkpKUrpsgUBAenp6st+0qakpvX37lrX2kzOEMiXhEy5vcjp48CClpqZS7969CQCrEfe/taLITkZGBjVv3pwA0OzZsxX5LnNx7tVXsnBwpp33c697ZPk3Z23Dhw/nfEaRnJxMtWrVIi0tLRo0aBBdu3aNU8VUaggziN6e/s+ldo35f6ki+FlG8ZFIpBHy2ZXDUgOi03+VuokpL0pzZv7s2TPZb7p69eqsm5yylITjxdxm8g4dOlDHjh1lQYVsuh/ziiIboaGhZGZmRgDo5MmTxf0ucyGRSGjaaQ+q6uhMj3yjch2zt7cnbW1tWWK6RYsWcX6Tv3v3jrS0tOQW8f/3v//Rixcvfj3Tl0RCFPyc6N/O8gPcGnMi/4el3bsfn0+3ciuHpQZEb46Xds9+WFatWkUAqEmTJqwH14058jJXmo7sZI1dAKhFixY0a9YsCg0NZUU2ryhy8OLFC9LQ0CAtLS1WvCRSM4XUecsjarDchb7FyWeYdXd3pz///JM+fPhA5cuXJwC0YMECzgfsnGskAEhFRYXWrFnDmeyEhAQKDQ0tPWUUF0S0t1XuQW9dVaJPN/mZBpH0O/A8Q7S8bO7v6VBXztOo/Ap06NCBunTpwkrG1uxceRsiUxLpgtwzhfj4+Fy/6Y0bN7Imn1cUeXD48GECQJUrV6aIiIiifpf5EhCdQrZLblPPHY9z/ZO/fftGREQfP36kChUqEABycHDgdECVSCRyHl8AaNiwYZyaoUQiEY0fP56MjIyodevWNGnSJNq5cyc9ePCAk7QGBRLg9t86Rs7tpgNRSrRy+1OaxH/NvSidtW22IQrh85cVlbS0NE7WGYNiUmRK4n1IQp7nZDd5AaB58+ax2gdeUeTDtGnTCAC1bt2alXiHOx8jqKqjM00+8ZrE+aT48PHxIVNTUwJAc+fO5VRZJCUlyVJdZy2+tWzZkoKDgzmTKZFIaO7cubmefABQuXLlaMuWLZzJzpegZ0S7muY9UC41ILq9gCjhm/L7xRVRvkTXZuT/eQ/+QRTOfsr834GMjAzWXd2zL17veeiX73mHDh2S/ZbGjh3L+thRkKJgpMeVC8MwZQGcA1AVQBCAQUQUn8d5QQCSAYgBiIjIrijt29nZ0evXrws9TygUonPnznB1dcXkyZOxd+/eIn+G/DjoFoDVN30wyd4K87vXzvMcX19ftG/fHuHh4Zg1axY2bdrEWVSll5cX2rRpg/fv32PMmDF4+PAhjIyMcPToUfTu3ZsTmQCwbt06zJ8/X25fhw4dcPv2bVntcLbx9fVFbGwsUlJSkJycnOfrpJGDYBN7C3i2M/+GyloBDYYBdQcCRlVZ6VtQUBAsLCzY/z9H+QDvzgFe54DkAlLUtF8INJ8CaOqzK5+nREgkBKsF0hLMHWqVx+HRTfI9d86cOdi8eTP69u2L8+fPQ01NjdW+MAzzJr8xtrQUxQYAcUS0jmEYRwBGROSQx3lBAOyIKKY47RdVUQBAdHQ0mjRpguDgYOzfvx8TJ04sjqhcEBGWXP2IEy+CsbqfLYY1s8jzvM+fP6N9+/YICwvDzJkzsWXLFs6UxYMHD9ChQweIxWKsXr0ay5cvh0QiwcyZM7F+/XrOCrLv3bsXf//9N7LfY9WrV8eiRYswdOhQ1m/0z58/Y+bMmbh161aex8eMGYNDhw7Jf8/JkcCL3cDz3YCkkLxSZaoAlvZA1VaARSvAsAqQz/9MLBbj+fPnuH79OpydnVG3bl2cPXu2+B9KIgHiAoCvz4FAN8D/PpAWW/A1OsYQNZmI/W+EmDDdgbP/L0/JqbbgJsQSgroqgy+ruxd4bo8ePZCWloZbt25BS0uL9b78iIrCF0A7IgpnGMYMgCsR1czjvCBwrCgAwNPTEy1btoRIJMLjx4/RrFmz4ojLhUgswcQTb+DqG4VdQxuhe12zPM/z8/ND+/btERISAkdHR6xdu7ZEcouKq6srhg4divDwcNjZ2cHV1RW6urqcyDpz5gxGjhyJ7t27Q0VFBVeuXAEAWFtb4+TJkyX+rnNCRHB2dsbMmTMREBAgd0xDQwPNmjXD2LFjMXr06PwaAMI8AM8zwMdLhQ/KeSAhICIVCE0UIz6DkCogVKpSFXpljGBjYytVSIJUIDMZyEySykgOB0hS/A9sUBmoN0i6lZfOYF1cXDB9+nRoaWnhyZMnePr0KerXrw8zs7zvQ7aQSCRITEyEkZERACAjIwPBwcGoWTPXT5tzIiMjUaFCBaXLFYlEEIlERRrIO255BL8oaU64gDXdoaJS8IPiwIEDcfjwYRgYGOR5XCgUlmi2/iMqigQiMsz2Pp6IjPI4LxBAPKR2uf1EdKCANicCmAgAVapUaRwcXLzsrufOncOFCxdw5MgR6OnpFevavEgTiDDy0Et4hSTgwAg7tK9VPs/z/P390adPHxw8eBAtWrQosdyiEhUVhVGjRqFixYo4dOgQp7Ju3LgBd3d3rFixAl5eXlixYgVu3boFf39/zgavjIwMbNq0CWvWrEF6ejp0dHRgbGyMb9++YfHixVixYkXxGiQCYr4AQY+BoCdA8DMgJYKTvufCsApg3gyo2hqo1kH6Pg8CAwPxzz//4OrVqwAAbW1tCAQCiMVi7Nu3Ty5JJtsIhUKMHTsWw4cPh729PQ4ePIj169fDwMAAHz584DyLbBZEhH///RczZszAuXPn0KtXL6XIBaTfwfDhw5GUlITLly8XqCyG/fsCT/2kDyGfVnaFlnrB3w8RITY2FiYmJnke9/HxQffu3XHw4EF07NhRof4XpCg4W1AGcA/Ahzy2PgAScpwbn08bFb+/lgfgBcC+KLIVKYVKxH4QT2K6gHrscKMaC2/SM7/8o1pLqxSqWCymjIwMpchKS5NPc8KGp1lRCA4OpoEDB5KOjg4JBAIKCAhgze88J2KxmNzd3WnZIkfq2tyGmlZSpc7VVKlvLTU6NrcPPds3g+jtKSKvc0Q+zkQBj4hC3kjdeksQ1JaamkpLliwhTU3NXA4E1tbWNGHCBHr8+DF7HzQH6enp1KtXLwJAQ4cOlTlrqKqq0ujRoykxMXdCOy5ITEykIUOGEADS0tKiY8eOKUUukTRaOytBYJ06dSg6On+vuhlnPGSL14npJfeeio+Pp+rVqxMA2rp1q8Lt4EfzegLgC8Ds+99mAHyLcM0yAHOK0r6iioILYlMyqeNmV7JZfIte+P84KRB+N+7du8eZgsiP4OBg2r17N3Xr1o369u3LevsSiYQuXrwoq9Oc1zZt2jTW5WYnKSmJ2rVrJydTTU2Nxo0bx2p1yfxkZxUpev36NVWrVo0AUK1atVhL/lkUMjMzqX///rLcSwW5gk/PpiSikkr+kCYSiahr164EgEaNGlWih92CFEVpmZ42Aoil/xazyxLRvBzn6AJQIaLk73/fBbCCiG4X1n5x1yi4JiopA38dfIHQhHQcHGmHNtXLlXaXeJRMamoqdHR0WHVYICIkJibm6+GV9Tp06FBUrFiRNblZxMTEoFu3bsj5W1u+fDmWLFnCurzsEBGGDx8OTU1NNGzYEHPmzIFAIMDo0aOxa9cuztbcciIQCDBo0CBcvXoV9erVw/379/M1D407+gr3P0UBAJ46dkAlQ+0Sy3d0dMT69evRpEkTuLm5lWiRu1RMTwVtAIwB3Afw5ftr2e/7KwK4+f1vK0jNTV4APgJYWNT269Wrp7BW5Yro5AzqsvURVV94k+55K8fswsPDFSEhIXkmvTQxMaF27dpxauoi+i/zgKqqKuF7CQFlmpqIpDEVWSa3Bg0aUExM/haDP/c+k80kIpPYyZ115swZAkAVKlSQBfWWBPxoMwquUVFRoTFjxmDOnDmoXTvvWIbSICFNgJGHX8I7LAmbB9VHnwaVSrtLPDzFxs/PD4MGDYKOjg5sbW1Rp04d1KlTB7a2tihfPm+nDTZ5//49mjZtiowMaYVJNTU13L9/H/b29pzL/vjxI+rUqYPMzEwMGDAAN27cQMOGDXHv3j2ULVs21/lEhBZrHyAiSdpXj8WdUFa35O7Kb9++RatWrSASieDq6oqWLVuWuM0fbkbB9YZsTzi9e/emJ0/yLvzBBsW1CSamC2jQPunTxV5Xv18vSR/PL09ycnKp3bdZ2ZGRYyZjbm5OXl5enMoWi8XUsGFDev36NXXr1o0AUOPGjfOtbyMQiWWzCLYWromklTuz1qUOHDjASptEBc8o2I14+gFJSkpCQkICJBIJ6/VsAWlg1dy5c1GuXDn8+eefhdYLNtBSx/FxTTHn/Dusu/UJYQnpWNqrDlQL8aHm4flRYMN9XFH+/vtvfPr0CYA0LqZ3794YPnw4unXrxnlg4ZkzZ/D27Vt07NgRCQkJaNKkCVxcXGRxI9lJTBei/vI7svdfVneDumrJxx+hUIhBgwbh69evmDJlCiZMmFDiNotEfhrkZ96Qwz3Qzy///ClsEBERISt92rBhQ1q7dm2hMsViCa254U0WDs407ugrSs4oeSElHp5fmSNHjhAAsre3pwMHDlB8fLzSZKenp5OFhYVsXNHT0yNPT888z/UOS5TNIpqsusvq7CsrP12bNm1YyU+XHfxo7rFcb+XLl6dy5crR6NGjZcno3N3dWfo68+bZs2ekrq4uNx1u1KgRrVu3rkA3wWPPAslq/g3qtMWVAqP5Km08PHmRmJhI27Zto8DAwFKRn71GdtZmZWWVq7Ld2ZfBMiUx+QQ7WXkPHTpEYrFYlvHa3NycIiMjWWk7O7+domjYsCHNmzePJBIJLVmyhACQjo6OzOeaK3bt2pXrZurcuTN5eHgUeN3TL9FUf7kL1V16O1fxo8L4/Pkzffz4sdSC9nh4fnXi4uLIyMhI9puuUKEC7dq1K9cT/YRjr2RK4uxLdjI0p6WlkZGRES1ZsoTVGjp58dspisaNG8sNnAcPHiRVVVVSUVGhffv2lfT7zBeJREIjRoyQUxSdOnUqUqDX19hU6rL1EVk6SlMN55emPCepqanUoUMH0tXVpVatWtG0adPo6NGj9P79e1bqgvPw/O5kpc03MDCgVatW5SpYlJoplFu0zlnnuiQcPXpUbjw5deoUa23n5LdUFDm5ceMG6ejoEABauHAhZ14bqampVK9ePVJXV6cmTZoQADI2NqZLly4Vem1KhpCmnnpDFg7ONOqwO8UkFy1yMzU1lTp16pRrNqOtrU0tWrSgly9flvRj8fD8lgQHB1OZMmVo9uzZecZJuAfEyimJFJbXGps3by73m27QoAErJZzzglcU33n16pWsHOnIkSNZXwzKws/Pj1q0aEECgYAWLVpEKioqBIDGjRtXaPlEiURCx58HUfWFN6np6rvkHpC3611O0tLSqEuXLrmUxc6dO9n4SHkiEokoPDycUlJSeDdfnl+Su3fv0tevX/M85njRS6Ygpp56w7rst2/fyv2WNTQ0aOPGjZyZmXlFkQ1/f39ZAq1OnTpxlrAsJCRE9vfjx4+patWqBICqVatGL168KPT6D6EJ1G7jQ7J0dKbNd3xJICq8qlZ6ejp179491801Y8YMTha/iIhOnDhBenp6xDAM6evrk5mZGVWvXp0aNWpE9vb2NHjw4Hz9zHl4fkZiUzLlZhGuxVxXLCoTJ06Uc4z58OEDJ3Ky4BVFDqKjo6lFixYEgOrXr6+UZHEJCQk0cuRIWdqBFStWFLqGkJwhpH/OviULB2fqueMxfYlMKlRO9rQCM2fOJGNjY5k735IlSzhRjJ8/f6bGjRvnmZTu5s2brMsrjIyMDPLz86MHDx6Qt7e30uXz/LocfRrISRBdThITE0lXV5dUVVVp6dKlnNa6z4JXFHmQmppKffv2JQBUpUoV+vjxY1G/zxJx7tw5MjQ0JEBavzogIKDQa26+C6MGy12oxsKbdOhxQKEL3ZmZmdSnTx8KDAykxMREWrJkCenq6srWS7Zs2ULp6ezkm8kuc9asWbkUhYaGBvXt25fOnz/PusxXr17Rtm3baNasWTRw4EBq2rSpLMU1AKpcuXKB+Xe4IiUlhW7fvs2b434h0jJFcgpik8snTuXt3r2bateuTa9eveJUTnZ4RZEPIpGI/v77bwJAhoaG9OjRoyJdV1K+fv0qS82sr69Px44dK3RQiUxKp7FHXpKFgzMN3Pu00NlFZmamXK2JyMhImj59uizWw9zcnA4dOsS6Z9SNGzeoXLlyBID++OMPMjMzkw3cBgYGNHr0aLp79y4rdtaMjAzas2cPVa5cOc/ZTK1atWjWrFl05MgRev36da6aGGzy9etXWUpxTU1NmjlzJmey8iMpKYl8fX2VLvdX5/jzIDkl8S0ulXOZR44c4fR+zQteURSARCKhDRs2yJ5+lWWqEIvFtGHDBtnAvXLlykKvkUgk5PTqK9Vb5kLVF9ykrXd9KUNYvAE3MDCQRo0aRQzDyGyfbC+OhYWF0R9//EG7d+8mkUhE9+/fp7Fjx1KZMmVkg7ipqSlrT0sZGRm0a9cuqlSpUp4KI2tTUVGhmjVr0pYtW0osUywW04sXL2jRokVUv379XLI6depEf/31l9w2evRoFj5tbpKTk2ndunVkbGwsy9oqFArJx8eHnJycaPHixZznQcqJWCwmNzc3unDhglLlEkkfAJ2cnEp8X0cnZ8gpiNlOnoU+0IWFhdGnT9zONvKjpP9jXlEUgdOnT9OMGTOUbi7w8PCg5s2bU1BQUJGviU7OoGmnpQVQ/tjsWmTPqOy8f/+e+vTpQw4ODsW+tiiIRKJcP5j09HS6ePEi9e/fn4yMjCgpqfA1l+KQnp5OO3fulKVTOXfuHF29epVWrVpFgwYNotq1a5OqqiqtWrWqRHIkEgmdPn2a2rdvT2pqagUqp+ybjo4OS59USkpKCm3YsIFMTExkMoYPH06NGjUiLS0tOdnbt29nVTaR9HvYu3ev3ADl5eVFDg4OsqR1FStWVFowqEQioStXrsjSnx85ckThdhZfeS+nJELjC3+6v3z5MhkbG5OtrS3rZtaCkEgk5OjoSAzDkJOTk8Lt8IriB0dR5fTAJ5Jarr1PFg7ONP2MB4UnFP/mFIsL96biAi5LsKanp9P27dtp7969eR5jU0HFx8fT2bNnadiwYXLRuxs2bCAvLy+5ja2qa6mpqbRp0yaZq3deW5UqVahHjx7k6OhIp0+fztfFU1Hi4uJkVd3c3d1pzZo1ZGtrK5OvpqZGPXv2pNOnTytlIfbhw4dyMQc9e/ak9+/fF7+dT5FyCmL3wy+FXpOSkkITJkyQzVoXLlyolM9MJDUxDx8+nLIixt+8UdxNl1cUvzCpmULa5PKJqi+8SbUX36JdD74U2xz1q6Ls2aFQKCQ3NzeaO3cujRgxghP5V69eJXNz83wVxNixYzlPlvf48eN8+9C6dWvau3dvgTWjS4JEIqGrV6/K3r9580Yufqh169YKlRWISpI3M9mtukupmYWv3718+VLmbm9hYUFubm7Flq0oCQkJ9McffxAAqlmzZonzYPGK4jcgOCZVlmvGfsMDcvYK471uShmuvn+xWExhYWH04sULcnJyok2bNtH06dOpb9++1KJFC/r8+TMnckUiEa1cuVIWQJq1GRsb09q1a4tlPlUEsVhMU6ZMoY4dO5Kvry8NGjRI1od69erRjRs3iv2dZwrFNGDPUzkl8T4kodDrRCIRrVq1SmZ6HD58OCUkFH4dW4SEhFC9evUIALVq1YqVWCVeUfxGuH2Oos5bHpGFgzP13vmYnvkp3z2U59cjJCRE5qmXc2MYhg4ePMipfLFYLAtA09PTk5VAtbKyotOnTxfbhCqRSGjBpXdyCmL/o6KVIwgMDKTWrVsTACpTpgydOXNGkY+kMO/fv5d5+g0YMIA176iCFMUvWQrVzs6OchZ8/50QSwiXPEKw5e5nhCdmoF3NcpjbpSbqVCxT2l3j+QlJSEjAunXrkJGRAX19fejp6UFPT0/2t76+PvT19dG0aVNOioNJJBJMmjQJ//77r2xf2bJlsWrVKowbN67YBYv+fRyAVTd8ZO971a+I7YMbQCWf4mEvXryAhoYGGjZsiFOnTuHvv/9GUlIS2rZti+PHj6NKlSqKfTAFcHV1Rd++fZGYmIgZM2Zg8+bNUFVVZaXt364UatmyZRVayPrVSBeIaP8jP6q3zIUsHJxp4vFX9CFUedNjHp6SIhKJaMyYMblmMWpqarR69epieVQdehwgN4PotMWV0jILvj4yMpIqVapE27ZtoyFDhhAAUldXp3Xr1ik9tf/p06dJQ0ODANDmzZtZbx+/24zC0tKS/P39OXm6AYA7d+5AW1sbjRs3ho6ODicy2CQxTYjDTwNx+GkgkjNE6GxTATM6VudnGDw/NGKxGGPHjsXx48cBAKqqqrCxsUHjxo3ltsJmFEefBmLZdW/Z+7K6GnCZaY9y+poFXicSidC5c2c8fPhQtq9WrVo4deoUGjVqVIJPVjyICJs2bcK8efOgoaGB48ePY/DgwazLKWhG8UsqCh0dHXrx4gXq1avHSftRUVFo3749fH19UbduXTRt2hTNmjVDs2bNULt2bc4UVElJTBfiyNNAHHoiVRgdapXHJHsrNLUsC4bha3bz/DgQEbZs2QJvb2+ZQqhXrx60tbWLfP2amz44+DhQts9QRx13/2lbqILIwsHBARs2bJC9t7a2xsuXL/Oskc02Hz58gK2tLcRiMf755x/s3LkThoaGuHr1Kuzt7TmR+dspCoZhSF1dHStXrsScOXNYs+FlJzIyEu3atZMVes9CX18fHTt2xKFDh5RyQylCYroQR58G4djzIMSlCtDA3BCT7K3QuY4pVPOx0/Lw/AxkCMX4+5QH7n+Kku0z0dPAzRltUF5fq8jtXLx4EQMHDsy1v0+fPjh58iT09PRY6W9epKenw9bWFnfu3MHcuXNx+fJlVKlSBbdu3YKNjQ1ncn87RWFpaUmxsbFITk5Gy5Ytcfz4cVSrVo11OeHh4Wjbti2+fPki26evr4/nz5+jTp06rMsDgJSUFKSkpKBcuXIlVoDpAjEueITgoFsAvsaloaqxDka3rIoBjStDX0udpR7z8HBPUEwqBu57hpgUgWxfU8uyODy6CfQ01YrV1qdPn9CkSROkpKQAAOrVq4dBgwbhzz//RI0aNVjtd14sXLgQa9asQbly5RAdHY369evj5s2bqFixIqdyfztFYWdnRxcvXsTo0aPh6uoKXV1dbN68GRMnTmTdxBIaGoq2bdvC399fts/U1BTLly/H2LFjoaZWvJu0MEQiERwcHLB9+3aYmprCzMwMFStWlL1WrFgR5ubm6Ny5c5E/q1hCcPkYgQNuAfD8lgBdDVX0b1QZI1tYoHoFfVb7z8PDJhffhGD2eS+5fUOamGN1v7oKzY6Tk5PRrFkzqKqqypRDrVq12OpuoXh7e6NBgwYQCoUAACsrK7x+/Vop1onfzuspK45CLBbT1q1bSVNTkwBQ165dOak98fXrV7K0tKRq1arRlClTZD7etWvXpmvXrnESeHXs2DHZ58q+6ejo0PXr1xVu1/NrPM0650nVF94kCwdn+uvAc7ruFcpHe/P8MMSnZtKIQ+5yHkwWDs502SOk8IsL4fXr10orOZATsVhMbdq0yfWbbtGiBUVFcVMcKTv43byecsZReHt7Y+TIkXjz5g2MjIywd+9e1r0GgoODMX/+fJw+fRqfPn2Co6Mjrl69CgBo27YtNm3aBDu7vJW1orx8+RL9+vVDWFiYbJ+JiQn+/vtvTJgwAZUqVVK47diUTJx7/Q2nXnxFaEK6bP/CDhXRzFQVAoEg361Ro0ZKmaLz/D4QEc6/DsG8i+/k9pfX18SFyS1RxfjH9z4sjMOHD2PcuHGy9zY2Npg3bx7++uuvYseKKMJvO6PIjkAgoKVLl8qe9ocMGcJ6ic6cEZJubm7UtGlT2ZPBkCFDilSoqDiEhYXJkqGpq6vL0parqqpS//796e7duyWa0YjEEnL1jcr19FZx/F5SN6mSZ64fZefR5/l1+RCaQK3W3c91/22540uiQgp4/UxERUVR2bJlCQDZ29uTs7Oz0hN24nefUWTn1atXGDlyJD59+gQzMzMcPnwYXbt25awvRITz589j/vz5CAgIgIaGBv73v/9h4cKFKFu2LCsyMjIyMGXKFNy+fRuenp44dOgQ9u/fj69fvwIAatSogSlTpmDUqFElsnX6RaVgzNGX+BaXLrc/PeANYm5uhSQ1AQ0aNEDfvn3RuXNnNGnShPU1miyICBEREfj06ZPc5ujoiPbt23Mik0c5RCRmYJaTJ575x8rtty6vhyOjm8C87M8/e8jJ2LFjkZiYiLlz56J58+al0gd+RpGDtLQ0mjlzpuwpeNKkSZScnFxUxasQGRkZtG3bNtlTg5GREW3atIm1dNtZtQGyEIlEdO3aNeratavsc2pra9PYsWNZKRh05W1Irqc8CwdnMunjSCq60lKvhoaGNGDAADpw4AArCePevHlDo0aNoqZNm5KBgUGu2Uzz5s0pJiZGqckQJRIJeXh40J49e/gkjCUgLCGNxh19mec9det9eGl3j1PS09N/iMqE4GcUefPgwQOMGTMGX79+hZWVFY4fP45WrVpx2reEhASsXbsW27dvR2ZmJqpWrYrVq1djyJAhnAXq+fv7Y//+/Th06BDi4uIAAE2aNMGUKVMwePDgEkWXiyWEOf/exuUASa5jwpAPiLqxDaKECABAzZo10blzZ3Tp0gVt27ZVyBf9w4cP2Lx5M06dOiXzDMmJoaEhqlWrlmuztrZGxYoVS/w9p6en48GDB7h+/TqcnZ0RGhqK9evXY968eSVqt7gQEUJCQmBubq5UuWwREJ0Cx4vv8TIoLtexRT1qY2wry3zzL/Gwz2/pHlvUpICJiYmYOXMmjh49CnV1dQQGBpZoEbioBAcHY/HixThx4gQAYNOmTZg9ezanMjMyMnD+/Hns2bMHL168ACD1Eff09Cyx23BaWhpWrFqDcvZDsdM1MNdxFUEqom7vQqrPYwCAuro63N3d0bBhQ4XkhYaGYseOHdi3bx+SkpIAAG3atIFIJIK/vz+ioqLyvE5TUxNz587FypUriyUvPDwcN27cwPXr13Hv3j2kpaXJHa9evTq0tAoO6NLW1oa7u3ux5OaFWCzG5cuXsXHjRvz555+YM2eO7FhGRgYCAwPh5+cHf39/2TZnzhx06NChxLKLgkgkwuvXr+Hi4gIAWLp0KQBAIiFc9QrF4isfkZIpynXd/G61MK61JdRUS6bI3717hwMHDmD16tUoU0Z5aWqICM7OzmAYBj179lSaXEAaX3X06FH8/fffCv+WedNTEbh8+TJt3Lix2NeVFA8PD+rXrx/rC+tFkTthwgTOPrNYLKHjzwLzNCVYODhTtSGLKTAyocRyEhMTafPmzVS5cmVydHSU7U9KSiJPT0+6cOECrV+/niZOnEh//PEHVa1alTZs2FAsGRKJhO7du0fTpk2jqlWrFqnsaV5bSUuhpqam0q5du8jKykrW5tSpU2ns2LHUtm1bqly5sqwWes6NrSRyaWlpdPr06Vz7g4KC6MCBAzRgwAAyNDT8r1aFRQ2aedYj3/vg1ItgErOwKJ2RkUGnTp2iVq1ayWTv2rWrxO0WlZcvX1Lbtm1l1QWVVeGOiOjJkyeye0LR8q9EvOmJ5wfA42s8ll/7CK+QxDyP/92+GibaV0MZbcUiwoVCIZ4/f85ZHpwsiAgfP36Es7Mzrl+/jufPn4OIcOTIEYwaNarQ6xV52ouOjsauXbuwe/duxMbG5nuejo5Ovia3KlWqlNix4N69e5g8eTLat2+PrVu3wtXVFXfu3IGLiws+f/4MAFDR0kO5NkOh1ah3nm3YWRhhWe86sK3EzpN+UFCQzKwaHR0NAGjYsCGmTp2Kv/76C7q6uqzIyY/AwEAsWLAAZ8+eBQA0a9YMGzduRJs2bTiVCwCZmZlYunQpNmzYACLC6NGjsW3bNoVnUbzpieeHIiVThINuAdh+/0u+54xsYYEJbax+eA+X6Oho3Lx5E58+fcKaNWtYj/y/d+8eHB0d8e7duzzXZHr16oV58+ahWrVqMDU15SS5Y1RUFGbPno2TJ08CkGYeiI2NhVAohJqhGco0HwC9+vl7Ds7rWhOjW1aFjoZiikooFOLly5ey9UOxWAwXFxfs2bMHN2/eBBFBU1MTgwcPxtSpU9G0aVPOk1zGxcVh9erV2LVrFwQCAapVq4a1a9di4MCBSkmw6eXlhREjRuD9+/coV64cDhw4gL59+5aoTV5R8PzQRCZlYNcDP5x4EZzvOY2qGGJ4cwt0r2sGLXX2kzz+6IhEIgQGBsLX11fOHfjbt2+4f/8+rK2tWZeZNVOaO3cu4uLiwKhpQqe2PfQbdoemWfV8rxvVwgJ/d7AuVhK+/IiOjsaff/4Je3t7TJs2DYcPH8a+ffsQFBQEAKhWrRomT56MMWPGwNjYuMTyCiMjIwO7du3C6tWrkZCQgLJly2LJkiWYMmWKUoLixGIxNm7ciCVLlkAoFKJPnz44cOAAypcvX+K2eUXB81MRm5KJY8+Dse+RPwSi3N5UWbSwMka/hpXQta4pDH7jJIZCoRDq6ux+/peeHzF5zQGEaVSCVuX8E1zqaKhiSttqGNHCAoY67A6UXl5e6NOnD4KDg2FtbY2vX79CIBBARUUFvXr1wpQpU9CpUyfOvAWJSDY7kEgkOHv2LBYsWIDg4GBoampi5syZcHR0hKGhISfyc+Ln54dRo0bh2bNnMDAwwI4dOzBy5EjWZjC8ouD5qRFLCI8+R+HUi69y6aPzQktdBT3qVkTH2uXRpka5YmcO/d1ITBfC1TcK93yi4PIhAgJx/ooZAKqqJaA6E4EKlAAHh3nQ1+cmaeT58+cxevRoOe+yChUqYMKECZgwYQLn5UePHTuGmjVronnz5nj48CHmzp2LN2/egGEYjBgxAitXrlRaCVQiwv79+zF79mykpaWhffv2OHLkCCwsLFiV88MpCoZh/gSwDEBtAE2JKM9RnWGYrgC2A1AF8C8RrStK+7yi+PURiCR48CkK17xCcfN9RJGusS6vh2aWZdHcyhjNrMqyYhr5GYhIzMCLgFg894/Fi8BYBMemFX4RgL4NKqJPw0pobW0C9RK6rBYViUSCpUuXYtWqVbmOjRw5Ev/++y/rs6ecnDp1CiNGjMC2bdtw584d3LhxAwDQsWNHbNiwQWGXbkUIDQ3FuHHj4OLiAi0tLaxbtw7Tpk3jZBb1IyqK2gAkAPYDmJOXomAYRhXAZwCdAIQAeAXgLyLyznluTthWFCEhIRCJRKhcuTJnKSl42CEyKQN3vCNx3ycSrr7RxbpWQ00FdSuVgW1FA9hWKoNapgawMNGBvqbaD1cBkIiQkCZEYGwq/CJT8D40UbaJJcX7TXesXQGdbMqjQ60KRa7+xgVJSUkYPnw4rl+/DgAwMjLKlUa/f//+aNq0KWd9OHfuHIYOHQqJ5L+ZVd26dbFhwwZ06dKF0/sgu6kLAM6cOYOpU6ciISEBdnZ2OH78OGrXrs2Z/B9OUciEM4wr8lcULQAsI6Iu39/PBwAiWltYu2wriuTkZAwePBh3795FlSpVYGVlJbdZWlqiWrVqP2xFOx6p+epzZLL0qTpAuiVl5A76UgRdDVWUN9CCkY46dDXVoKOhCl0NNenfmqrQzPY0nv3XJhQTUjNFSM0UIeX7lpopQlKGCFFJGUgViFnpX1ldDTS3ks6kWlgZw7q83g+n+ABpoFxycjLMzMxgZmZW5LKnbHHhwgUMGTIEYvF/3/vUqVOxY8cOTqpk5mTBggVYs2YNYmNjMXXqVDg5OUFVVRWLFy/GggULOJ9J/ayKYiCArkQ0/vv7EQCaEdH/Cmu3Tp069PHjR1b7KhKJMG3aNOzbty/XsXLlyuHixYtK8Z3m4Y6kDCF8wpJkT+ZfIlMQGJOKdCE7Azbb6GuqwaqcLqqV04NtpTKoV7kMapsZQJdflyk2V65cwZ9//gmRSP7hQVNTE6dOncKAAQM4lX/gwAFMmjQJJ0+exJw5cxAREYFatWrhxIkTrJcnyI+CFAVndxTDMPcAmOZxaCERXS1KE3nsy1erMQwzEcDErPcTJkzAypUrYWqaVxeKj5qaGvbs2QMrK6tcOX2aNm3KWibYnEgkEkRGRkJbWxtaWlrQ1NT8IZ8GfwUMtNTRzMoYzayK7mYplhBSMkSITslAfJoQaQIx0jJFSBWIkSYQITVTjEyRGEy22znr36eqwkBPUzrz0NNU/f6qBn0tdZTT14S+phqf60gJXL9+HYMGDYJIJIKFhQWaNWsm2xo2bFiiXGhF4dmzZ/jf/6TPv8OHDwcAzJw5E2vWrFH6rCpf8gvZVsYGwBWAXT7HWgBwyfZ+PoD5RWk3qyaDrq4urVy5klJTU4sYxF40nJyc8qwu17Vr1xLXf8iL06dPU5kyZQgAMQxD2traZGRkRBUrVqRq1apRhw4dKCYmhlWZPDy/A+Hh4bRixQq6evUqRUREKF1+SEgImZqayo0jU6dOVXotCqKCU3j8yIpCDUAAAEsAGgC8ANQpSrsNGzakFStWkK6uLgGgypUr04kTJ1j98p8+fUrGxsY0f/58OnDgANWqVUv2j65fvz4dO3aMMjMzWZMXFBSUZ5lEAHT+/HnW5PDw8CiH9PR0ucJmWZuJiQkdOHBA6f354RQFgH6QejJlAojMmjkAqAjgZrbzukPq+eQPqcmqSO1nJQUMCwujcePGyRKl2dnZkZubG2tf7JcvX+jQoUNEJK136+zsTO3bt5f9wytWrEjr1q2juLg4VuSJRCJavXo1qamp5bq5bG1tacWKFfTp0ydWZGXx4cMHOn/+PO3fv5/Wrl1Lc+fOpXHjxlG/fv2obdu2VLduXbp16xarMnl4fnUkEgmNHTtW9vutVKkSTZs2jR4+fEhCobBU+vTDKQqut5zZY9++fUsdOnSQ/VMGDBhAfn5+JflOC+TNmzc0dOhQWdlVXV1dmj59OmtlUF++fEnVq1cnANS7d2+5jJlZM5rVq1fTly9fSiwrOTmZlixZQtra2nnOZtq1a6eUwu95ER4ezkpBJB4eZbNr1y6ytrYmBwcHcnd3LxVTU05+e0VBJNXg165doxo1asjqS8+ePZvi4+MV/FoL5+vXrzRnzhxZNTYVFRUaOHAgvXjxosRtJycn04QJE2jZsmVERPTt2zfasmWLrH521taoUSNat25diZXUt2/faOTIkfmm0La2tqaRI0fS3r17ycvLi0QiUYk/YxaZmZnk6elJx44do9mzZ1PHjh2pfPnypKWlRZ8/f2ZNDg+PMhCLxeTj4/PDVUTkFUU2BAIB7dixQ1aS1NjYmHbu3Mlp/vjExETasmULValSRTawtmrVii5dulTiATU0NDTXvqCgINq4cSM1adJEbjBv0qQJbdy4kYKDgxWW9+rVK7K3t5e1OW7cOLKzs5PNnrI2fX196tixIy1ZsoRu376tkEJ+9+4ddejQIU9TGwCqVasWrV27lk6cOEGurq7k7+/PWmnZ/IiKiqKjR4/SggULfrgfOg9PSeAVRR7ExcXRrFmzKMtDqmbNmnT9+nVOf/xCoZDOnDlDdnZ2ck/iu3btopSUFE5kBgQE0Lp166hRo0Zyg2zz5s1py5Yt9O3bt2K3KZFI6NKlS1StWjW6d+8eEUmL6jx69IjWrl1LvXr1IhMTEzl5DMNQnTp1aMKECXT48GH69OlTkb7r9PR0OnLkSK7+F7SVL1+eGjduTH379qVp06bRhg0b6PTp0/T48WMKCgoq1kOBRCKhDx8+0Nq1a6lly5ay9a7jx48X+3tjg7weDHh42IBXFAXg5+dHAwYMIABkYGCglEpzEomEHj16RL1795YNPHv37uVc7pcvX2j16tVUv3592aCqSDXALDIzMykwMDDPYxKJhD5//kzHjh2jSZMmUd26dXNVX/P09CyyLIlEQk+fPqXBgwfLZhgTJkygY8eO0apVq2jSpEnUvXt3qlu3rlyFtby2adOmFSrr7t27NH36dLK0tCyygipsU7TCnUgkosuXL1OHDh3y7LtAIKDg4GB68uQJnTlzhjZs2EDTpk2jvn37kp2dHVWoUIHOnTunkOzikpCQQC4uLrR06VLq3LkzDRgwQClyswgMDKT169dT48aNlb5+FRgYSDNnzix2BUU2ePr0KQ0aNKhEM+qCFAWfPfY7jx8/xtevXzFs2DCOepU3nz9/xp49e7B69WrOq3Flx9fXF05OTrCwsMDIkSOVIjMxMREvX77Es2fP4OHhgUuXLimUGiE0NBT79+/Hx48fcfHixTzPSUlJQUhICL59+ya3hYSEYMiQIRg9enShMrLXyM7IyJA7bm1tDU3N4uVF0tbWxqtXr4p8flxcHA4dOoQ9e/bI6i+MGzcOBgYGcp8pIiJCLjdRTgwNDbFlyxaMGTOmWP3N4uXLl7CysoKJiYncfiLCly9f8OzZMzx//hzPnj3Dx48fkX1MsbS0hL+/P6dBot++fYOTkxOcnJzw8uVL2f4jR44U+n9mg1evXmHz5s24cOECxGIxrKys8OXLF87Sn2cnMDAQjo6OcHJyAiBNaDh06FCF2uJrZvP8kmRmZiplnSA1NZWuXbtGEyZMIDMzMwLYq0GdF+/evaMJEybk62mWfdPX1ycbGxvq3LkzjRs3jpYtW0aHDh0iFxcX8vb2pqSkJIX7ERwcTMOGDSNDQ0MSiUSUkpJCDx8+pNWrV1PPnj3J2Ng4l3nR1taWJk6cSEeOHCFfX1/O/j8hISG0bds2atGihVwfGjRoQGvXruXUq5FIuiB95coVudimypUr08aNGykhIYFT2UTSmdu8efNIQ0ODAFD16tXpypUrJfq+wZueeHjYQSwW0+vXr2nv3r2cDIKurq7Uv39/sra2zmWqA0CdO3emmzdv0vv37zkbkJKSkmjhwoWkpaVFAKhs2bLUqFGjXA4LBgYG1LlzZ1q6dCm5uLiw2p+goCCKjIyU2xcWFkY7d+6k1q1by/Wjbt26tGrVKvL19WVNfn6kpqbS3r17Ze7p+O5ZeOrUKU4dYrIQCoW0Z88e2RqgkZERbdu2jZXgXl5R8PD8hCQnJ9Pz589p3759NHXqVGrVqhWZmZnluy5UUkQiER08eJAqVKiQ5+ylRo0aNGrUKNq/fz+9f/+eVRfo7GSlrImNjaWIiAjas2cPtW3bVk5x2tjY0PLly8nb25uTPuQkIiKCFi9eLDeL6tmzJz18+FAps1qJREI3btyg2rVrEwBSU1OjmTNnsrqmyisKHp5fBIlEQmlpaay3e/fuXapXr16+Jq6VK1eyLjMnCQkJNGzYMAKkMUd//PEHqaioyCmqxYsX04cPHziRLxKJcg36Hz9+pHHjxslyu2lqatKECROUpqCIpKbITp06yb6Hfv36cRI/VJCi4PMR8/D8RDAMw3pGUSJC7dq1cfLkScTFxcltsbGxiIuLw4cPHxASEoLKlSuzKjuLJ0+eYPjw4QgODgYgzZp8//59WFtbY9CgQRg8eDDq1q3L2aJ4amoqRo0ahfPnz4OI8PDhQ2zevBk3b94EAJiYmMDR0RFTp05F+fLlOelDTiIiIrBkyRIcOnQIEokEDRs2xJYtW9CuXTulyJcjPw3yM2/8jIKH5+dAIBDQwoUL5WYOyDaDUMbCcHh4ONnZ2ZGOjg6dOHGCGjRoIOtDzZo1af/+/ZzM4vLj/+3df4gcZx3H8feHYkjOlKQxamOsWqEUf6AYtNQazGnENAex+oehIrbBgjSQYP84w4kghZqIgv4hqKC1SSNq/rFqKC1aPaXEpE3PpLn87I+EGs/ERGtzVYK/6tc/5rk43u3OztzN7Ma7zwuWnZ2Z59nPPfPcPjuzs7MXL16Mbdu2xcKFCwOya8bt2LGj8ct84NNjzexyMzY2xsaNGxkdHWXBggWXfnMlP93f38/mzZsby3Ds2DEGBgYu7clMWLVqFYODgwwMDDR6muv4+DiLFi0Csjftu3btYmhoiNOnT9PX18eWLVsYHBzsyqnzPj3WzGyS4eHhS7/zMnFbvHhx7Nu3ryvPf+LEiVi3bl1EROzZs+fSJcclxYYNG2JsbKwrOSZQsEfR/DdCzMwuMzt37mTNmjWMj4//z/wLFy6wdevWKV+wrNupU6dYvXo1Bw8eZP369axcuZL9+/fT39/PyMgI27dvZ/ny5Y1mqGJWHnqS9Efgtx1XbG0p8Kca49TFuapxrmqcq5rZmOv1EfHKVgtm5UAxE5JGot1xuh5yrmqcqxrnqmau5fKhJzMzK+SBwszMCnmgmOpbvQ7QhnNV41zVOFc1cyqXP6MwM7NC3qMwM7NCHijMzKzQnBwoJH1U0lFJ/5bU9lQySTdLekrSs5KGcvOXSHpE0jPp/qqacnWsV9L1kp7M3V6UdFdadrek3+eWDXQrV1rvOUmH03OPVC3fRC5J10j6paTjaZt/OrestvZq11dyyyXpa2n5qKQVZcvORIlcH095RiXtlfT23LKW27NLufoljee2zefLlu1Cts/kch2R9JKkJWlZI20m6T5J5yUdabO82f7V7ivbs/kGvAm4HvgV8M4261wBnATeCMwDDgFvTsu+DAyl6SHgSzXlqlRvyvgHsi/KANwNDDbQXqVyAc8BS2f6d9WZC1gGrEjTVwJP57ZjLe1V1Fdy6wwADwMCbgQeL1u24Vw3AVel6bUTuYq2Z5dy9QMPTqds09kmrb8OGO5Cm70XWAEcabO80f41J/coIuJ4RDzVYbUbgGcj4lRE/APYBdySlt0C3J+m7wc+XFO0qvWuBk5GxHS/hV7WTP/enrVXRJyNiANp+i/AcaDuayMU9ZV81p2ReQxYLGlZybKN5YqIvRHxQnr4GNDMdcQr5mqobBP1fwz4QY3P31JEPAr8uWCVRvvXnBwoSloO/C73eIz/vsC8OiLOQvZCBNR1gfqq9d7K1E66Ke163lfXIZ4KuQL4maTfSPrUNMo3lQsASW8A3gE8nptdR3sV9ZVO65QpO11V676D7F3phHbbs1u53i3pkKSHJb2lYtmmsyGpD7gZ+GFudlNt1kmj/WvW/nCRpJ8DV7dY9LmI+EmZKlrMm/G5xEW5KtYzD/gQ8Nnc7G8C95DlvAf4CvDJLuZ6T0SckfQq4BFJJ9I7oWmrsb0Wkv1D3xURL6bZ026vydW3mDe5r7Rbp5F+1uE5p64ovY9soFiZm1379qyQ6wDZIdW/ps+OfgxcV7Js09kmrAN+HRH5d/pNtVknjfavWTtQRMQHZljFGHBN7vFrgTNp+pykZRFxNu3ena8jl6Qq9a4FDkTEuVzdl6YlfRt4sJu5IuJMuj8v6Udku72P0uP2kvQyskHiexHxQK7uabfXJEV9pdM680qUna4yuZD0NuBeYG1EPD8xv2B7Np4rN5gTEQ9J+oakpWXKNp0tZ8oefYNt1kmj/cuHntp7ArhO0rXp3futwO60bDdwe5q+HSizh1JGlXqnHBtNL5YTPgK0PEOiiVySXi7pyolp4IO55+9Ze0kS8B3geER8ddKyutqrqK/ks96Wzk65ERhPh8vKlJ2ujnVLeh3wAPCJiHg6N79oe3Yj19Vp2yHpBrLXqufLlG06W8q0CFhFrs813GadNNu/6v50/v/hRvaiMAb8HTgH/DTNfw3wUG69AbKzZE6SHbKamP8K4BfAM+l+SU25WtbbIlcf2T/NoknlvwscBkZTZ1jWrVxkZ1UcSrejl0t7kR1KidQmT6bbQN3t1aqvAHcCd6ZpAV9Pyw+TO9uuXT+rqY065boXeCHXNiOdtmeXcm1Kz3uI7EP2m7rRXmWypccbgF2TyjXWZmRvCs8C/yR77bqjm/3Ll/AwM7NCPvRkZmaFPFCYmVkhDxRmZlbIA4WZmRXyQGFmZoU8UJiZWSEPFGZmVsgDhVnDJL0rXXhwfvr27lFJb+11LrOy/IU7sy6Q9AVgPrAAGIuIL/Y4kllpHijMuiBdZ+cJ4G9kl6N4qceRzErzoSez7lgCLCT7lb35Pc5iVon3KMy6QNJusl8Xu5bs4oObehzJrLRZ+3sUZpcLSbcB/4qI70u6Atgr6f0RMdzrbGZleI/CzMwK+TMKMzMr5IHCzMwKeaAwM7NCHijMzKyQBwozMyvkgcLMzAp5oDAzs0L/AWz7E/6lkc/+AAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "x, y = np.meshgrid(np.arange(-1.,1.1, 0.2), np.arange(-1.1,1.1, 0.2))\n", "U = y\n", "V = -x\n", "\n", "# plotting the slope field\n", "fig, ax = plt.subplots()\n", "q = ax.quiver(x, y, U, V, angles='xy')\n", "ax.quiverkey(q, X=0.8, Y=1.02, U=2,\n", " label='dy/dx', labelpos='E')\n", "\n", "# plotting the simulated trajectories\n", "ax.plot(data1[:, 0], data1[:, 1])\n", "ax.plot(data2[:, 0], data2[:, 1])\n", "\n", "ax.set_xlabel('x')\n", "ax.set_ylabel('y')\n", "\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.7.9" } }, "nbformat": 4, "nbformat_minor": 4 }